BgsGentleGenocide.cc

Go to the documentation of this file.

//-----------------------------------------------------------------------------
// File and Version Information:
// $Id: BgsGentleGenocide.cc,v 1.1 2007/10/12 07:30:20 caogf Exp $
//
// Description:
//  Class BgsGentleGenocide
//
//      A process that prevents the associated particle from being tracked
//      if it has traveled more than the specified number of steps and
//      falls below the specified kinetic energy.
//
//      To fix an apparent "loophole" in hadronic ionization for tritium
//
// Environment:
//  Software developed for the BaBar Detector at the SLAC B-Factory.
//
// Author List:
//  David Williams
//
// Copyright Information:
//  Copyright (C) 2001         SLAC
//
//-----------------------------------------------------------------------------
 
#include "PhySim/BgsGentleGenocide.hh"
#include "PhySim/BaBar.hh"
 
//
// Constructor
//
BgsGentleGenocide::BgsGentleGenocide( G4double theMaxKineticEnergy, G4int theMinNumberSteps,
                                      const G4String& name, G4ProcessType type )
    : BgsGenocide( name, type )
    , maxKineticEnergy( theMaxKineticEnergy )
    , minNumberSteps( theMinNumberSteps ) {
  ;
}
 
//
// Destructor
//
BgsGentleGenocide::~BgsGentleGenocide() { ; }
 
//
// Decide what to do at each step
//
G4double BgsGentleGenocide::PostStepGetPhysicalInteractionLength(
    const G4Track& track, G4double previousStepSize, G4ForceCondition* condition ) {
  *condition = NotForced;
 
  //
  // Do nothing if we haven't reached the minimum number of steps
  //
  if ( track.GetCurrentStepNumber() < minNumberSteps ) return DBL_MAX;
 
  //
  // Do nothing if we are above the threshold kinetic energy
  //
  const G4DynamicParticle* particle = track.GetDynamicParticle();
 
  if ( particle->GetKineticEnergy() > maxKineticEnergy ) return DBL_MAX;
 
  //
  // Die!
  //
  //    const G4ThreeVector &pos = track.GetPosition();
 
  const G4VProcess* creator            = track.GetCreatorProcess();
  G4String          creatorProcessName = "(none?)";
  if ( creator ) creatorProcessName = creator->GetProcessName();
 
  /*
  G4cout << "BgsGentleGenocide: rogue "
         << particle->GetDefinition()->GetParticleName()
         << " terminated, created by "
         << creatorProcessName << G4endl;
 
  G4cout << " KE = "
         << particle->GetKineticEnergy()/keV << " keV, "
         << " global coords (cm): ("
         << pos.x()/cm << ","
         << pos.y()/cm << ","
         << pos.z()/cm << ")" << G4endl;
  */
 
  *condition = NotForced;
  return 0;
}