#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) override

virtual G4VParticleChange *	AtRestDoIt (const G4Track &aTrack, const G4Step &aStep) override

virtual void	BuildPhysicsTable (const G4ParticleDefinition &) override

virtual G4bool	IsApplicable (const G4ParticleDefinition &) override

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition) override

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition) override

virtual void	StartTracking (G4Track *) override

virtual void	EndTracking () override

void	SetExtDecayer (G4VExtDecayer *)

const G4VExtDecayer *	GetExtDecayer () const

G4double	GetRemainderLifeTime () const

virtual void	ProcessDescription (std::ostream &outFile) const override

Public Member Functions inherited from G4VRestDiscreteProcess
	G4VRestDiscreteProcess (const G4String &, G4ProcessType aType=fNotDefined)

	G4VRestDiscreteProcess (G4VRestDiscreteProcess &)

virtual	~G4VRestDiscreteProcess ()

G4VRestDiscreteProcess &	operator= (const G4VRestDiscreteProcess &)=delete

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

virtual G4VParticleChange *	PostStepDoIt (const G4Track &, const G4Step &)

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)

virtual G4VParticleChange *	AtRestDoIt (const G4Track &, const G4Step &)

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 ()

G4VProcess &	operator= (const G4VProcess &)=delete

G4bool	operator== (const G4VProcess &right) const

G4bool	operator!= (const G4VProcess &right) const

virtual G4VParticleChange *	PostStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AlongStepDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4VParticleChange *	AtRestDoIt (const G4Track &track, const G4Step &stepData)=0

virtual G4double	AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)=0

virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &track, G4ForceCondition *condition)=0

virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)=0

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 G4bool	IsApplicable (const G4ParticleDefinition &)

virtual void	BuildPhysicsTable (const G4ParticleDefinition &)

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	StartTracking (G4Track *)

virtual void	EndTracking ()

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

virtual void	ProcessDescription (std::ostream &outfile) 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) override

virtual G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *condition) override

virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)=0

virtual G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *condition)=0

Protected Member Functions inherited from G4VProcess
void	SubtractNumberOfInteractionLengthLeft (G4double prevStepSize)

void	ClearNumberOfInteractionLengthLeft ()

Protected Attributes
G4int	verboseLevel

const G4double	HighestValue

G4double	fRemainderLifeTime

G4ParticleChangeForDecay	fParticleChangeForDecay

G4VExtDecayer *	pExtDecayer

Protected Attributes inherited from G4VProcess
const G4ProcessManager *	aProcessManager = nullptr

G4VParticleChange *	pParticleChange = nullptr

G4ParticleChange	aParticleChange

G4double	theNumberOfInteractionLengthLeft = -1.0

G4double	currentInteractionLength = -1.0

G4double	theInitialNumberOfInteractionLength = -1.0

G4String	theProcessName

G4String	thePhysicsTableFileName

G4ProcessType	theProcessType = fNotDefined

G4int	theProcessSubType = -1

G4double	thePILfactor = 1.0

G4int	verboseLevel = 0

G4bool	enableAtRestDoIt = true

G4bool	enableAlongStepDoIt = true

G4bool	enablePostStepDoIt = true

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

Detailed Description

Definition at line 62 of file G4Decay.hh.

Constructor & Destructor Documentation

◆ G4Decay()

G4Decay::G4Decay ( const G4String & processName = "Decay" )

Definition at line 62 of file G4Decay.cc.

                               :G4VRestDiscreteProcess(processName, fDecay),
                verboseLevel(1),
                                HighestValue(20.0),
                fRemainderLifeTime(-1.0),
                                pExtDecayer(nullptr)
{
  // 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::~G4Decay ( )

virtual

Definition at line 81 of file G4Decay.cc.

{
  if (pExtDecayer != nullptr) {
    delete pExtDecayer;
  }
}

Member Function Documentation

◆ AtRestDoIt()

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

inlineoverridevirtual

Reimplemented from G4VRestDiscreteProcess.

Reimplemented in G4DecayWithSpin.

Definition at line 183 of file G4Decay.hh.

{
  return DecayIt(aTrack, aStep);
}

◆ AtRestGetPhysicalInteractionLength()

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

overridevirtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 477 of file G4Decay.cc.

{
     // 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;
}

◆ BuildPhysicsTable()

void G4Decay::BuildPhysicsTable ( const G4ParticleDefinition & )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 175 of file G4Decay.cc.

{
  return;
}

◆ DaughterPolarization()

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

protectedvirtual

Reimplemented in G4PionDecayMakeSpin.

Definition at line 386 of file G4Decay.cc.

{
  // empty implementation
}

Referenced by DecayIt().

◆ DecayIt()

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

protectedvirtual

Definition at line 180 of file G4Decay.cc.

{
  // 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 != nullptr);   
  G4DecayProducts* products = nullptr;
 
  // decay table
  G4DecayTable   *decaytable = aParticleDef->GetDecayTable();
 
  // check if external decayer exists
  G4bool isExtDecayer = (decaytable == nullptr) && (pExtDecayer != nullptr);
 
  // Error due to NO Decay Table 
  if ( (decaytable == nullptr) && !isExtDecayer && !isPreAssigned ){
    if (GetVerboseLevel()>0) {
      G4cout <<  "G4Decay::DoIt  : decay table not defined  for ";
      G4cout << aParticle->GetDefinition()->GetParticleName()<< G4endl;
    }
    G4ExceptionDescription ed;
    ed << "For " << aParticle->GetDefinition()->GetParticleName()
       << " decay probability exist but decay table is not defined "
       << "- the particle will be killed;\n"
       << "    isExtDecayer: " << isExtDecayer
       << "; isPreAssigned: " << isPreAssigned;    
    G4Exception( "G4Decay::DecayIt ",
                 "DECAY101",JustWarning, ed);
 
    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 according to decay table.
    // Keep trying to choose a candidate decay channel if the dynamic mass
    // of the decaying particle is below the sum of the PDG masses of the
    // candidate daughter particles.
    // This is needed because the decay table used in Geant4 is based on
    // the assumption of nominal PDG masses, but a wide resonance can have
    // a dynamic masses well below its nominal PDG masses, and therefore
    // some of its decay channels can be below the kinematical threshold. 
    // Note that, for simplicity, we ignore here the possibility that
    // one or more of the candidate daughter particles can be, in turn,
    // wide resonance. However, if this is the case, and the channel is
    // accepted, then the masses of the resonance daughter particles will
    // be sampled by taking into account their widths.
    G4VDecayChannel* decaychannel = nullptr;
    G4double massParent = aParticle->GetMass();
    decaychannel = decaytable->SelectADecayChannel(massParent);
    if ( decaychannel == nullptr) {
      // decay channel not found
       G4ExceptionDescription ed;
      ed << "Can not determine decay channel for " 
     << aParticleDef->GetParticleName() << G4endl 
     << "  mass of dynamic particle: " 
     << massParent/GeV << " (GEV)" << G4endl
     << "  dacay table has " << decaytable->entries() 
     << " entries" << G4endl;
      G4double checkedmass=massParent;
      if (massParent < 0.) {
    checkedmass=aParticleDef->GetPDGMass();
    ed << "Using PDG mass ("<<checkedmass/GeV 
       << "(GeV)) in IsOKWithParentMass" << G4endl; 
      }
      for (G4int ic =0;ic <decaytable->entries();++ic) {
    G4VDecayChannel * dc= decaytable->GetDecayChannel(ic);
    ed << ic << ": BR " << dc->GetBR() << ", IsOK? " 
       << dc->IsOKWithParentMass(checkedmass)
       << ", --> "; 
    G4int ndaughters=dc->GetNumberOfDaughters();
    for (G4int id=0;id<ndaughters;++id) {
      if (id>0) ed << " + ";   // seperator, except for first
      ed << dc->GetDaughterName(id);
    }
    ed << G4endl;
      }
      G4Exception("G4Decay::DoIt", "DECAY003", FatalException,ed);
    } 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) {
    G4ExceptionDescription ed;
    ed << "Total Energy is less than its mass - increased the energy"
       << "\n Particle: " << aParticle->GetDefinition()->GetParticleName()
       << "\n Energy:"    << ParentEnergy/MeV << "[MeV]"
       << "\n Mass:"      << ParentMass/MeV << "[MeV]";
    G4Exception( "G4Decay::DecayIt ",
                 "DECAY102",JustWarning, ed);
    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 ;
} 

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

◆ EndTracking()

void G4Decay::EndTracking ( )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 401 of file G4Decay.cc.

{
  // Clear NumberOfInteractionLengthLeft
  ClearNumberOfInteractionLengthLeft();
 
  currentInteractionLength = -1.0;
}

◆ GetExtDecayer()

const G4VExtDecayer * G4Decay::GetExtDecayer ( ) const

inline

Definition at line 193 of file G4Decay.hh.

{
  return pExtDecayer;
}

◆ GetMeanFreePath()

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

overrideprotectedvirtual

Implements G4VRestDiscreteProcess.

Definition at line 129 of file G4Decay.cc.

{
   // 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;
}

Referenced by PostStepGetPhysicalInteractionLength().

◆ GetMeanLifeTime()

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

overrideprotectedvirtual

Implements G4VRestDiscreteProcess.

Definition at line 100 of file G4Decay.cc.

{
   // 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;
}

Referenced by AtRestGetPhysicalInteractionLength().

◆ GetRemainderLifeTime()

G4double G4Decay::GetRemainderLifeTime ( ) const

inline

Definition at line 199 of file G4Decay.hh.

{
  return fRemainderLifeTime;
}

◆ IsApplicable()

G4bool G4Decay::IsApplicable ( const G4ParticleDefinition & aParticleType )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 88 of file G4Decay.cc.

{
   // 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; 
   }
}

Referenced by LBE::ConstructGeneral(), and G4DecayPhysics::ConstructProcess().

◆ PostStepDoIt()

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

overridevirtual

Reimplemented from G4VRestDiscreteProcess.

Reimplemented in G4DecayWithSpin.

Definition at line 507 of file G4Decay.cc.

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

◆ PostStepGetPhysicalInteractionLength()

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

overridevirtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 410 of file G4Decay.cc.

{
   // 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 = 0.0;
    
    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;
  }
}

◆ ProcessDescription()

void G4Decay::ProcessDescription ( std::ostream & outFile ) const

overridevirtual

Reimplemented from G4VProcess.

Reimplemented in G4DecayWithSpin, and G4PionDecayMakeSpin.

Definition at line 521 of file G4Decay.cc.

{
  outFile << GetProcessName() << ": Decay of particles. \n"
      << "kinematics of daughters are dertermined by DecayChannels " 
          << " or by PreAssignedDecayProducts\n";
}

◆ SetExtDecayer()

void G4Decay::SetExtDecayer ( G4VExtDecayer * val )

Definition at line 497 of file G4Decay.cc.

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

◆ StartTracking()

void G4Decay::StartTracking ( G4Track * )

overridevirtual

Reimplemented from G4VProcess.

Definition at line 393 of file G4Decay.cc.

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

Member Data Documentation

◆ fParticleChangeForDecay

G4ParticleChangeForDecay G4Decay::fParticleChangeForDecay

protected

Definition at line 176 of file G4Decay.hh.

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

◆ fRemainderLifeTime

G4double G4Decay::fRemainderLifeTime

protected

Definition at line 173 of file G4Decay.hh.

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

◆ HighestValue

const G4double G4Decay::HighestValue

protected

Definition at line 170 of file G4Decay.hh.

Referenced by GetMeanFreePath().

◆ pExtDecayer

G4VExtDecayer* G4Decay::pExtDecayer

protected

Definition at line 179 of file G4Decay.hh.

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

◆ verboseLevel

G4int G4Decay::verboseLevel

protected

Definition at line 162 of file G4Decay.hh.

Referenced by PostStepGetPhysicalInteractionLength().

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

Public Member Functions

Protected Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ G4Decay()

◆ ~G4Decay()

Member Function Documentation

◆ AtRestDoIt()

◆ AtRestGetPhysicalInteractionLength()

◆ BuildPhysicsTable()

◆ DaughterPolarization()

◆ DecayIt()

◆ EndTracking()

◆ GetExtDecayer()

◆ GetMeanFreePath()

◆ GetMeanLifeTime()

◆ GetRemainderLifeTime()

◆ IsApplicable()

◆ PostStepDoIt()

◆ PostStepGetPhysicalInteractionLength()

◆ ProcessDescription()

◆ SetExtDecayer()

◆ StartTracking()

Member Data Documentation

◆ fParticleChangeForDecay

◆ fRemainderLifeTime

◆ HighestValue

◆ pExtDecayer

◆ verboseLevel