G4Decay Class Reference

#include <G4Decay.hh>

Inheritance diagram for G4Decay:

Public Member Functions
	G4Decay (const G4String &processName="Decay")
virtual	~G4Decay ()
virtual G4VParticleChange *	PostStepDoIt (const G4Track &aTrack, const G4Step &aStep)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &aTrack, const G4Step &aStep)
virtual void	BuildPhysicsTable (const G4ParticleDefinition &)
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual void	StartTracking (G4Track *)
virtual void	EndTracking ()
void	SetExtDecayer (G4VExtDecayer *)
const G4VExtDecayer *	GetExtDecayer () const
G4double	GetRemainderLifeTime () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
Public Member Functions inherited from G4VRestDiscreteProcess
	G4VRestDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)
	G4VRestDiscreteProcess (G4VRestDiscreteProcess &)
virtual	~G4VRestDiscreteProcess ()
virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &, G4double, G4double, G4double &, G4GPILSelection *)
virtual G4VParticleChange *	AlongStepDoIt (const G4Track &, const G4Step &)
Public Member Functions inherited from G4VProcess
	G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
	G4VProcess (const G4VProcess &right)
virtual	~G4VProcess ()
G4int	operator== (const G4VProcess &right) const
G4int	operator!= (const G4VProcess &right) const
G4double	GetCurrentInteractionLength () const
void	SetPILfactor (G4double value)
G4double	GetPILfactor () const
G4double	AlongStepGPIL (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
G4double	AtRestGPIL (const G4Track &track, G4ForceCondition *condition)
G4double	PostStepGPIL (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual void	PreparePhysicsTable (const G4ParticleDefinition &)
virtual G4bool	StorePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
virtual G4bool	RetrievePhysicsTable (const G4ParticleDefinition *, const G4String &, G4bool)
const G4String &	GetPhysicsTableFileName (const G4ParticleDefinition *, const G4String &directory, const G4String &tableName, G4bool ascii=false)
const G4String &	GetProcessName () const
G4ProcessType	GetProcessType () const
void	SetProcessType (G4ProcessType)
G4int	GetProcessSubType () const
void	SetProcessSubType (G4int)
virtual void	SetProcessManager (const G4ProcessManager *)
virtual const G4ProcessManager *	GetProcessManager ()
virtual void	ResetNumberOfInteractionLengthLeft ()
G4double	GetNumberOfInteractionLengthLeft () const
G4double	GetTotalNumberOfInteractionLengthTraversed () const
G4bool	isAtRestDoItIsEnabled () const
G4bool	isAlongStepDoItIsEnabled () const
G4bool	isPostStepDoItIsEnabled () const
virtual void	DumpInfo () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	SetMasterProcess (G4VProcess *masterP)
const G4VProcess *	GetMasterProcess () const
virtual void	BuildWorkerPhysicsTable (const G4ParticleDefinition &part)
virtual void	PrepareWorkerPhysicsTable (const G4ParticleDefinition &)

Protected Member Functions
virtual G4VParticleChange *	DecayIt (const G4Track &aTrack, const G4Step &aStep)
virtual void	DaughterPolarization (const G4Track &aTrack, G4DecayProducts *products)
virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
virtual G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *condition)
Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
void	ClearNumberOfInteractionLengthLeft ()

Protected Attributes
G4int	verboseLevel
const G4double	HighestValue
G4double	fRemainderLifeTime
G4ParticleChangeForDecay	fParticleChangeForDecay
G4VExtDecayer *	pExtDecayer
Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager
G4VParticleChange *	pParticleChange
G4ParticleChange	aParticleChange
G4double	theNumberOfInteractionLengthLeft
G4double	currentInteractionLength
G4double	theInitialNumberOfInteractionLength
G4String	theProcessName
G4String	thePhysicsTableFileName
G4ProcessType	theProcessType
G4int	theProcessSubType
G4double	thePILfactor
G4bool	enableAtRestDoIt
G4bool	enableAlongStepDoIt
G4bool	enablePostStepDoIt
G4int	verboseLevel

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)

Detailed Description

Definition at line 63 of file G4Decay.hh.

Constructor & Destructor Documentation

G4Decay::G4Decay

(

const G4String &

processName = "Decay"

)

Definition at line 63 of file G4Decay.cc.

References DECAY, fParticleChangeForDecay, G4cout, G4endl, GetVerboseLevel(), G4VProcess::pParticleChange, and G4VProcess::SetProcessSubType().

                               :G4VRestDiscreteProcess(processName, fDecay),
                verboseLevel(1),
                                HighestValue(20.0),
                fRemainderLifeTime(-1.0),
                                pExtDecayer(0)
{
  // set Process Sub Type
  SetProcessSubType(static_cast<int>(DECAY));

#ifdef G4VERBOSE
  if (GetVerboseLevel()>1) {
    G4cout << "G4Decay  constructor " << "  Name:" << processName << G4endl;
  }
#endif

  pParticleChange = &fParticleChangeForDecay;
}

G4Decay::~G4Decay ( )

virtual

Definition at line 82 of file G4Decay.cc.

References pExtDecayer.

{
  if (pExtDecayer) {
    delete pExtDecayer;
  }
}

Member Function Documentation

G4VParticleChange * G4Decay::AtRestDoIt ( const G4Track &  aTrack, const G4Step &  aStep  )

inlinevirtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 191 of file G4Decay.hh.

References DecayIt().

{
  return DecayIt(aTrack, aStep);
}

G4double G4Decay::AtRestGetPhysicalInteractionLength ( const G4Track &  track, G4ForceCondition *  condition  )

virtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 438 of file G4Decay.cc.

References DBL_MIN, fRemainderLifeTime, G4Track::GetDynamicParticle(), GetMeanLifeTime(), G4DynamicParticle::GetPreAssignedDecayProperTime(), G4Track::GetProperTime(), NotForced, and G4VProcess::theNumberOfInteractionLengthLeft.

{
     // condition is set to "Not Forced"
  *condition = NotForced;

  G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime();
  if (pTime >= 0.) {
    fRemainderLifeTime = pTime - track.GetProperTime();
    if (fRemainderLifeTime <= 0.0) fRemainderLifeTime = DBL_MIN;
  } else {
    fRemainderLifeTime = 
      theNumberOfInteractionLengthLeft * GetMeanLifeTime(track, condition);
  }
  return fRemainderLifeTime;
}

void G4Decay::BuildPhysicsTable ( const G4ParticleDefinition &  )

virtual

Reimplemented from G4VProcess.

Definition at line 176 of file G4Decay.cc.

{
  return;
}

void G4Decay::DaughterPolarization ( const G4Track &  aTrack, G4DecayProducts *  products  )

protectedvirtual

Reimplemented in G4PionDecayMakeSpin.

Definition at line 348 of file G4Decay.cc.

Referenced by DecayIt().

{
}

G4VParticleChange * G4Decay::DecayIt ( const G4Track &  aTrack, const G4Step &  aStep  )

protectedvirtual

Reimplemented in G4DecayWithSpin.

Definition at line 181 of file G4Decay.cc.

References G4VParticleChange::AddSecondary(), G4DecayProducts::Boost(), G4VProcess::ClearNumberOfInteractionLengthLeft(), python.hepunit::cm, DaughterPolarization(), G4VDecayChannel::DecayIt(), G4DecayProducts::DumpInfo(), G4DecayProducts::entries(), FatalException, fParticleChangeForDecay, fRemainderLifeTime, fStopAndKill, fStopButAlive, G4cout, G4endl, G4Exception(), G4ParticleDefinition::GetDecayTable(), G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4Track::GetGlobalTime(), G4DynamicParticle::GetKineticEnergy(), G4Track::GetLocalTime(), G4DynamicParticle::GetMass(), G4DynamicParticle::GetMomentumDirection(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGStable(), G4Track::GetPosition(), G4DynamicParticle::GetPreAssignedDecayProducts(), G4DynamicParticle::GetTotalEnergy(), G4Track::GetTouchableHandle(), G4Track::GetTrackStatus(), GetVerboseLevel(), G4VDecayChannel::GetVerboseLevel(), G4VExtDecayer::ImportDecayProducts(), G4ParticleChangeForDecay::Initialize(), G4DecayProducts::IsChecked(), JustWarning, python.hepunit::MeV, ns, pExtDecayer, G4DecayProducts::PopProducts(), G4VParticleChange::ProposeLocalEnergyDeposit(), G4ParticleChangeForDecay::ProposeLocalTime(), G4VParticleChange::ProposeTrackStatus(), G4DecayTable::SelectADecayChannel(), G4Track::SetGoodForTrackingFlag(), G4VParticleChange::SetNumberOfSecondaries(), G4Track::SetTouchableHandle(), G4VDecayChannel::SetVerboseLevel(), test::x, and z.

Referenced by AtRestDoIt(), G4DecayWithSpin::DecayIt(), and PostStepDoIt().

{
  // The DecayIt() method returns by pointer a particle-change object.
  // Units are expressed in GEANT4 internal units.

  //   Initialize ParticleChange
  //     all members of G4VParticleChange are set to equal to 
  //     corresponding member in G4Track
  fParticleChangeForDecay.Initialize(aTrack);

  // get particle 
  const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
  const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();

  // check if  the particle is stable
  if (aParticleDef->GetPDGStable()) return &fParticleChangeForDecay ;
 

  //check if thePreAssignedDecayProducts exists
  const G4DecayProducts* o_products = (aParticle->GetPreAssignedDecayProducts());
  G4bool isPreAssigned = (o_products != 0);   
  G4DecayProducts* products = 0;

  // decay table
  G4DecayTable   *decaytable = aParticleDef->GetDecayTable();
 
  // check if external decayer exists
  G4bool isExtDecayer = (decaytable == 0) && (pExtDecayer !=0);

  // Error due to NO Decay Table 
  if ( (decaytable == 0) && !isExtDecayer &&!isPreAssigned ){
    if (GetVerboseLevel()>0) {
      G4cout <<  "G4Decay::DoIt  : decay table not defined  for ";
      G4cout << aParticle->GetDefinition()->GetParticleName()<< G4endl;
    }
    G4Exception( "G4Decay::DecayIt ",
                 "DECAY101",JustWarning, 
                 "Decay table is not defined");

    fParticleChangeForDecay.SetNumberOfSecondaries(0);
    // Kill the parent particle
    fParticleChangeForDecay.ProposeTrackStatus( fStopAndKill ) ;
    fParticleChangeForDecay.ProposeLocalEnergyDeposit(0.0); 
    
    ClearNumberOfInteractionLengthLeft();
    return &fParticleChangeForDecay ;
  }

  if (isPreAssigned) {
    // copy decay products 
    products = new G4DecayProducts(*o_products); 
  } else if ( isExtDecayer ) {
    // decay according to external decayer
    products = pExtDecayer->ImportDecayProducts(aTrack);
  } else {
    // decay acoording to decay table
    // choose a decay channel
    G4VDecayChannel *decaychannel = decaytable->SelectADecayChannel();
    if (decaychannel == 0 ){
      // decay channel not found
      G4Exception("G4Decay::DoIt", "DECAY003", FatalException,
                  " can not determine decay channel ");
    } else {
      // execute DecayIt() 
#ifdef G4VERBOSE
      G4int temp = decaychannel->GetVerboseLevel();
      if (GetVerboseLevel()>1) {
    G4cout << "G4Decay::DoIt  : selected decay channel  addr:" << decaychannel <<G4endl;
    decaychannel->SetVerboseLevel(GetVerboseLevel());
      }
#endif
      products = decaychannel->DecayIt(aParticle->GetMass());
#ifdef G4VERBOSE
      if (GetVerboseLevel()>1) {
    decaychannel->SetVerboseLevel(temp);
      }
#endif
#ifdef G4VERBOSE
      if (GetVerboseLevel()>2) {
    if (! products->IsChecked() ) products->DumpInfo();
      }
#endif
    }
  }
  
  // get parent particle information ...................................
  G4double   ParentEnergy  = aParticle->GetTotalEnergy();
  G4double   ParentMass    = aParticle->GetMass();
  if (ParentEnergy < ParentMass) {
    if (GetVerboseLevel()>0) {
      G4cout << "G4Decay::DoIt  : Total Energy is less than its mass" << G4endl;
      G4cout << " Particle: " << aParticle->GetDefinition()->GetParticleName();
      G4cout << " Energy:"    << ParentEnergy/MeV << "[MeV]";
      G4cout << " Mass:"    << ParentMass/MeV << "[MeV]";
      G4cout << G4endl;
    }
    G4Exception( "G4Decay::DecayIt ",
                 "DECAY102",JustWarning, 
                 "Total Energy is less than its mass");
    ParentEnergy = ParentMass;
  }

  G4ThreeVector ParentDirection(aParticle->GetMomentumDirection());

  //boost all decay products to laboratory frame
  G4double energyDeposit = 0.0;
  G4double finalGlobalTime = aTrack.GetGlobalTime();
  G4double finalLocalTime = aTrack.GetLocalTime();
  if (aTrack.GetTrackStatus() == fStopButAlive ){
    // AtRest case
    finalGlobalTime += fRemainderLifeTime;
    finalLocalTime += fRemainderLifeTime;
    energyDeposit += aParticle->GetKineticEnergy();
    if (isPreAssigned) products->Boost( ParentEnergy, ParentDirection);
  } else {
    // PostStep case
    if (!isExtDecayer) products->Boost( ParentEnergy, ParentDirection);
  }
   // set polarization for daughter particles
   DaughterPolarization(aTrack, products);


  //add products in fParticleChangeForDecay
  G4int numberOfSecondaries = products->entries();
  fParticleChangeForDecay.SetNumberOfSecondaries(numberOfSecondaries);
#ifdef G4VERBOSE
  if (GetVerboseLevel()>1) {
    G4cout << "G4Decay::DoIt  : Decay vertex :";
    G4cout << " Time: " << finalGlobalTime/ns << "[ns]";
    G4cout << " X:" << (aTrack.GetPosition()).x() /cm << "[cm]";
    G4cout << " Y:" << (aTrack.GetPosition()).y() /cm << "[cm]";
    G4cout << " Z:" << (aTrack.GetPosition()).z() /cm << "[cm]";
    G4cout << G4endl;
    G4cout << "G4Decay::DoIt  : decay products in Lab. Frame" << G4endl;
    products->DumpInfo();
  }
#endif
  G4int index;
  G4ThreeVector currentPosition;
  const G4TouchableHandle thand = aTrack.GetTouchableHandle();
  for (index=0; index < numberOfSecondaries; index++)
  {
     // get current position of the track
     currentPosition = aTrack.GetPosition();
     // create a new track object
     G4Track* secondary = new G4Track( products->PopProducts(),
                      finalGlobalTime ,
                      currentPosition );
     // switch on good for tracking flag
     secondary->SetGoodForTrackingFlag();
     secondary->SetTouchableHandle(thand);
     // add the secondary track in the List
     fParticleChangeForDecay.AddSecondary(secondary);
  }
  delete products;

  // Kill the parent particle
  fParticleChangeForDecay.ProposeTrackStatus( fStopAndKill ) ;
  fParticleChangeForDecay.ProposeLocalEnergyDeposit(energyDeposit); 
  fParticleChangeForDecay.ProposeLocalTime( finalLocalTime );

  // Clear NumberOfInteractionLengthLeft
  ClearNumberOfInteractionLengthLeft();

  return &fParticleChangeForDecay ;
} 

void G4Decay::EndTracking ( )

virtual

Reimplemented from G4VProcess.

Definition at line 362 of file G4Decay.cc.

References G4VProcess::ClearNumberOfInteractionLengthLeft(), and G4VProcess::currentInteractionLength.

{
  // Clear NumberOfInteractionLengthLeft
  ClearNumberOfInteractionLengthLeft();

  currentInteractionLength = -1.0;
}

const G4VExtDecayer * G4Decay::GetExtDecayer ( ) const

inline

Definition at line 201 of file G4Decay.hh.

References pExtDecayer.

{
  return pExtDecayer;
}

G4double G4Decay::GetMeanFreePath ( const G4Track &  aTrack, G4double  previousStepSize, G4ForceCondition *  condition  )

protectedvirtual

Implements G4VRestDiscreteProcess.

Definition at line 130 of file G4Decay.cc.

References python.hepunit::c_light, DBL_MAX, DBL_MIN, G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4DynamicParticle::GetKineticEnergy(), G4DynamicParticle::GetMass(), G4ParticleDefinition::GetParticleName(), G4ParticleDefinition::GetPDGLifeTime(), G4ParticleDefinition::GetPDGStable(), G4DynamicParticle::GetTotalMomentum(), GetVerboseLevel(), python.hepunit::GeV, and HighestValue.

Referenced by PostStepGetPhysicalInteractionLength().

{
   // get particle 
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
   G4double aMass = aParticle->GetMass();
   G4double aLife = aParticleDef->GetPDGLifeTime();


    // returns the mean free path in GEANT4 internal units
   G4double pathlength;
   G4double aCtau = c_light * aLife;

   // check if the particle is stable?
   if (aParticleDef->GetPDGStable()) {
     pathlength = DBL_MAX;

   //check if the particle has very short life time ?
   } else if (aCtau < DBL_MIN) { 
     pathlength =  DBL_MIN;
 
   } else {
    //calculate the mean free path 
    // by using normalized kinetic energy (= Ekin/mass)
     G4double   rKineticEnergy = aParticle->GetKineticEnergy()/aMass; 
     if ( rKineticEnergy > HighestValue) {
       // gamma >>  1
       pathlength = ( rKineticEnergy + 1.0)* aCtau;
     } else if ( rKineticEnergy < DBL_MIN ) {
       // too slow particle
#ifdef G4VERBOSE
       if (GetVerboseLevel()>1) {
     G4cout << "G4Decay::GetMeanFreePath()   !!particle stops!!";
         G4cout << aParticleDef->GetParticleName() << G4endl;
     G4cout << "KineticEnergy:" << aParticle->GetKineticEnergy()/GeV <<"[GeV]";
       }
#endif
       pathlength = DBL_MIN;
     } else {
       // beta <1 
       pathlength = (aParticle->GetTotalMomentum())/aMass*aCtau ;
     }
   }
  return  pathlength;
}

G4double G4Decay::GetMeanLifeTime ( const G4Track &  aTrack, G4ForceCondition *  condition  )

protectedvirtual

Implements G4VRestDiscreteProcess.

Definition at line 101 of file G4Decay.cc.

References G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4ParticleDefinition::GetPDGLifeTime(), G4ParticleDefinition::GetPDGStable(), GetVerboseLevel(), and ns.

Referenced by AtRestGetPhysicalInteractionLength().

{
   // returns the mean free path in GEANT4 internal units
   G4double meanlife;

   // get particle 
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   const G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
   G4double aLife = aParticleDef->GetPDGLifeTime();

   // check if the particle is stable?
   if (aParticleDef->GetPDGStable()) {
    //1000000 times the life time of the universe
     meanlife = 1e24 * s;
    
   } else {
     meanlife = aLife;
   }

#ifdef G4VERBOSE
   if (GetVerboseLevel()>1) {
     G4cout << "mean life time: "<< meanlife/ns << "[ns]" << G4endl;
   }
#endif

   return  meanlife;
}

G4double G4Decay::GetRemainderLifeTime ( ) const

inline

Definition at line 207 of file G4Decay.hh.

References fRemainderLifeTime.

{
  return fRemainderLifeTime;
}

G4int G4Decay::GetVerboseLevel ( ) const

inline

Definition at line 188 of file G4Decay.hh.

References verboseLevel.

Referenced by DecayIt(), G4Decay(), GetMeanFreePath(), and GetMeanLifeTime().

{ return verboseLevel; }

G4bool G4Decay::IsApplicable ( const G4ParticleDefinition &  aParticleType )

virtual

Reimplemented from G4VProcess.

Definition at line 89 of file G4Decay.cc.

References G4ParticleDefinition::GetPDGLifeTime(), G4ParticleDefinition::GetPDGMass(), and python.hepunit::MeV.

Referenced by XrayFluoPhysicsList::AddDecay(), PhysicsList::AddDecay(), OpNovicePhysicsList::ConstructDecay(), eRositaPhysicsList::ConstructGeneral(), RE06PhysicsList::ConstructGeneral(), Par01PhysicsList::ConstructGeneral(), ExGflashPhysicsList::ConstructGeneral(), B03PhysicsList::ConstructGeneral(), ExP02PhysicsList::ConstructGeneral(), ExN02PhysicsList::ConstructGeneral(), ExP01PhysicsList::ConstructGeneral(), F01PhysicsList::ConstructGeneral(), F02PhysicsList::ConstructGeneral(), F03PhysicsList::ConstructGeneral(), UltraPhysicsList::ConstructGeneral(), DMXPhysicsList::ConstructGeneral(), XrayTelPhysicsList::ConstructGeneral(), G4AdjointPhysicsList::ConstructGeneral(), LXeGeneralPhysics::ConstructProcess(), GammaRayTelGeneralPhysics::ConstructProcess(), and G4DecayPhysics::ConstructProcess().

{
   // check if the particle is stable?
   if (aParticleType.GetPDGLifeTime() <0.0) {
     return false;
   } else if (aParticleType.GetPDGMass() <= 0.0*MeV) {
     return false;
   } else {
     return true; 
   }
}

G4VParticleChange * G4Decay::PostStepDoIt ( const G4Track &  aTrack, const G4Step &  aStep  )

virtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 468 of file G4Decay.cc.

References DecayIt(), fParticleChangeForDecay, fStopAndKill, fStopButAlive, G4Track::GetTrackStatus(), and G4ParticleChangeForDecay::Initialize().

{
  if ( (aTrack.GetTrackStatus() == fStopButAlive ) ||
       (aTrack.GetTrackStatus() == fStopAndKill )   ){
    fParticleChangeForDecay.Initialize(aTrack);
    return &fParticleChangeForDecay;
  } else {
    return DecayIt(aTrack, aStep);
  }
}

G4double G4Decay::PostStepGetPhysicalInteractionLength ( const G4Track &  track, G4double  previousStepSize, G4ForceCondition *  condition  )

virtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 371 of file G4Decay.cc.

References python.hepunit::c_light, python.hepunit::cm, G4VProcess::currentInteractionLength, DBL_MAX, DBL_MIN, G4DynamicParticle::DumpInfo(), fRemainderLifeTime, G4cout, G4endl, G4DynamicParticle::GetDefinition(), G4Track::GetDynamicParticle(), G4DynamicParticle::GetMass(), G4Track::GetMaterial(), GetMeanFreePath(), G4Material::GetName(), G4ParticleDefinition::GetPDGLifeTime(), G4DynamicParticle::GetPreAssignedDecayProperTime(), G4Track::GetProperTime(), G4DynamicParticle::GetTotalMomentum(), NotForced, python.hepunit::perMillion, G4VProcess::SubtractNumberOfInteractionLengthLeft(), G4VProcess::theNumberOfInteractionLengthLeft, and verboseLevel.

{
   // condition is set to "Not Forced"
  *condition = NotForced;

   // pre-assigned Decay time
  G4double pTime = track.GetDynamicParticle()->GetPreAssignedDecayProperTime();
  G4double aLife = track.GetDynamicParticle()->GetDefinition()->GetPDGLifeTime();

  if (pTime < 0.) {
    // normal case 
    if ( previousStepSize > 0.0){
      // subtract NumberOfInteractionLengthLeft 
      SubtractNumberOfInteractionLengthLeft(previousStepSize);
      if(theNumberOfInteractionLengthLeft<0.){
    theNumberOfInteractionLengthLeft=perMillion;
      }
      fRemainderLifeTime = theNumberOfInteractionLengthLeft*aLife;
    }
    // get mean free path
    currentInteractionLength = GetMeanFreePath(track, previousStepSize, condition);
    
#ifdef G4VERBOSE
    if ((currentInteractionLength <=0.0) || (verboseLevel>2)){
      G4cout << "G4Decay::PostStepGetPhysicalInteractionLength " << G4endl;
      track.GetDynamicParticle()->DumpInfo();
      G4cout << " in Material  " << track.GetMaterial()->GetName() <<G4endl;
      G4cout << "MeanFreePath = " << currentInteractionLength/cm << "[cm]" <<G4endl;
    }
#endif

    G4double value;
    if (currentInteractionLength <DBL_MAX) {
      value = theNumberOfInteractionLengthLeft * currentInteractionLength;
      //fRemainderLifeTime = theNumberOfInteractionLengthLeft*aLife;
    } else {
      value = DBL_MAX;
    }

    return value;

  } else {
    //pre-assigned Decay time case
    // reminder proper time
    fRemainderLifeTime = pTime - track.GetProperTime();
    if (fRemainderLifeTime <= 0.0) fRemainderLifeTime = DBL_MIN;
    
    G4double  rvalue=0.0; 
    // use pre-assigned Decay time to determine PIL
    if (aLife>0.0) {
      // ordinary particle
      rvalue = (fRemainderLifeTime/aLife)*GetMeanFreePath(track, previousStepSize, condition);
    } else {
     // shortlived particle
      rvalue = c_light * fRemainderLifeTime;
     // by using normalized kinetic energy (= Ekin/mass)
     G4double   aMass =  track.GetDynamicParticle()->GetMass();
     rvalue *= track.GetDynamicParticle()->GetTotalMomentum()/aMass;
    }
    return rvalue;
  }
}

void G4Decay::SetExtDecayer

(

G4VExtDecayer *

val

)

Definition at line 458 of file G4Decay.cc.

References DECAY_External, pExtDecayer, and G4VProcess::SetProcessSubType().

Referenced by P6DExtDecayerPhysics::ConstructProcess().

{
  pExtDecayer = val;

  // set Process Sub Type
  if ( pExtDecayer !=0 ) {
    SetProcessSubType(static_cast<int>(DECAY_External));
  }
}

void G4Decay::SetVerboseLevel ( G4int  value )

inline

Definition at line 185 of file G4Decay.hh.

References verboseLevel.

{ verboseLevel = value; }

void G4Decay::StartTracking ( G4Track *  )

virtual

Reimplemented from G4VProcess.

Definition at line 354 of file G4Decay.cc.

References G4VProcess::currentInteractionLength, fRemainderLifeTime, and G4VProcess::ResetNumberOfInteractionLengthLeft().

{
  currentInteractionLength = -1.0;
  ResetNumberOfInteractionLengthLeft();
 
  fRemainderLifeTime = -1.0;
}

Field Documentation

G4ParticleChangeForDecay G4Decay::fParticleChangeForDecay

protected

Definition at line 178 of file G4Decay.hh.

Referenced by DecayIt(), G4Decay(), and PostStepDoIt().

G4double G4Decay::fRemainderLifeTime

protected

Definition at line 175 of file G4Decay.hh.

Referenced by AtRestGetPhysicalInteractionLength(), G4DecayWithSpin::DecayIt(), DecayIt(), GetRemainderLifeTime(), PostStepGetPhysicalInteractionLength(), and StartTracking().

const G4double G4Decay::HighestValue

protected

Definition at line 172 of file G4Decay.hh.

Referenced by GetMeanFreePath().

G4VExtDecayer* G4Decay::pExtDecayer

protected

Definition at line 181 of file G4Decay.hh.

Referenced by DecayIt(), GetExtDecayer(), SetExtDecayer(), and ~G4Decay().

G4int G4Decay::verboseLevel

protected

Definition at line 164 of file G4Decay.hh.

Referenced by GetVerboseLevel(), PostStepGetPhysicalInteractionLength(), and SetVerboseLevel().

---

The documentation for this class was generated from the following files:

• G4Decay.hh
• G4Decay.cc