G4MuonicAtomDecay Class Reference

#include <G4MuonicAtomDecay.hh>

Inheritance diagram for G4MuonicAtomDecay:

Public Member Functions
	G4MuonicAtomDecay (G4HadronicInteraction *hiptr=nullptr, const G4String &processName="MuonicAtomDecay")
virtual	~G4MuonicAtomDecay ()
virtual G4bool	IsApplicable (const G4ParticleDefinition &)
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel () const
virtual void	ProcessDescription (std::ostream &outFile) const
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)
virtual G4double	AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
virtual G4VParticleChange *	AtRestDoIt (const G4Track &theTrack, const G4Step &theStep)
virtual G4VParticleChange *	PostStepDoIt (const G4Track &theTrack, const G4Step &theStep)
virtual G4double	GetMeanFreePath (const G4Track &aTrack, G4double previousStepSize, G4ForceCondition *condition)
virtual G4double	GetMeanLifeTime (const G4Track &aTrack, G4ForceCondition *condition)
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 const G4VProcess *	GetCreatorProcess () const
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 &)

Additional Inherited Members
Static Public Member Functions inherited from G4VProcess
static const G4String &	GetProcessTypeName (G4ProcessType)
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 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

Detailed Description

Definition at line 50 of file G4MuonicAtomDecay.hh.

Constructor & Destructor Documentation

G4MuonicAtomDecay()

G4MuonicAtomDecay::G4MuonicAtomDecay ( G4HadronicInteraction *  hiptr = nullptr, const G4String &  processName = "MuonicAtomDecay"  )

explicit

Definition at line 67 of file G4MuonicAtomDecay.cc.

  : G4VRestDiscreteProcess(name, fDecay),
    fMuMass(G4MuonMinus::MuonMinus()->GetPDGMass()),
    cmptr(hiptr),
    verboseLevel(0)
{
  // This is not a hadronic process; assume it is a kind of decay
  enableAtRestDoIt = true;
  enablePostStepDoIt = true; // it is a streach; fixme
  theProcessSubType = 221; // (see G4DecayProcessType.hh) fixme
  if (!cmptr) {
    // cmptr = new G4CascadeInterface(); // Bertini - Pointer owned by InteractionRegistry
    cmptr = new G4MuMinusCapturePrecompound(); // Precompound
  }
}

~G4MuonicAtomDecay()

G4MuonicAtomDecay::~G4MuonicAtomDecay ( )

virtual

Definition at line 86 of file G4MuonicAtomDecay.cc.

{}

Member Function Documentation

AtRestDoIt()

virtual G4VParticleChange * G4MuonicAtomDecay::AtRestDoIt ( const G4Track &  theTrack, const G4Step &  theStep  )

inlinevirtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 86 of file G4MuonicAtomDecay.hh.

  {return DecayIt(theTrack, theStep);}

AtRestGetPhysicalInteractionLength()

G4double G4MuonicAtomDecay::AtRestGetPhysicalInteractionLength ( const G4Track &  aTrack, G4ForceCondition *  condition  )

virtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 113 of file G4MuonicAtomDecay.cc.

{
  *condition = NotForced;
   // check if this is the beginning of tracking
  if (theNumberOfInteractionLengthLeft < 0.) {
    ResetNumberOfInteractionLengthLeft();
  }
  return theNumberOfInteractionLengthLeft*GetMeanLifeTime(aTrack, condition);
}

GetMeanFreePath()

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

virtual

Implements G4VRestDiscreteProcess.

Definition at line 585 of file G4MuonicAtomDecay.cc.

{
  // based on G4Decay::GetMeanFreePath check; fixme
 
  // 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;
    const G4double HighestValue = 20.0; //
    if ( rKineticEnergy > HighestValue) {
      // gamma >>  1
      pathlength = ( rKineticEnergy + 1.0)* aCtau;
    } else if ( rKineticEnergy < DBL_MIN ) {
      // too slow particle
#ifdef G4VERBOSE
      if (GetVerboseLevel()>1) {
        G4cout << "G4MuonicAtomDecay::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;
}

GetMeanLifeTime()

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

virtual

Implements G4VRestDiscreteProcess.

Definition at line 135 of file G4MuonicAtomDecay.cc.

{
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   G4ParticleDefinition* aParticleDef = aParticle->GetDefinition();
   G4MuonicAtom* muatom = static_cast<G4MuonicAtom*>(aParticleDef);
   G4double meanlife = muatom->GetPDGLifeTime();
#ifdef G4VERBOSE
   if (GetVerboseLevel()>1) {
     G4cout << "mean life time: "<< meanlife/ns << "[ns]" << G4endl;
   }
#endif
   return  meanlife;
}

Referenced by AtRestGetPhysicalInteractionLength().

GetVerboseLevel()

G4int G4MuonicAtomDecay::GetVerboseLevel ( ) const

inline

Definition at line 74 of file G4MuonicAtomDecay.hh.

{return verboseLevel;}

Referenced by GetMeanFreePath(), and GetMeanLifeTime().

IsApplicable()

G4bool G4MuonicAtomDecay::IsApplicable ( const G4ParticleDefinition &  a )

virtual

Reimplemented from G4VProcess.

Definition at line 91 of file G4MuonicAtomDecay.cc.

{
  return ( a.GetParticleType() == "MuonicAtom" );
}

PostStepDoIt()

virtual G4VParticleChange * G4MuonicAtomDecay::PostStepDoIt ( const G4Track &  theTrack, const G4Step &  theStep  )

inlinevirtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 90 of file G4MuonicAtomDecay.hh.

  {return DecayIt(theTrack, theStep);}

PostStepGetPhysicalInteractionLength()

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

virtual

Reimplemented from G4VRestDiscreteProcess.

Definition at line 126 of file G4MuonicAtomDecay.cc.

{
  *condition = NotForced;
  return DBL_MAX; // this will need to be changed in future; fixme
}

ProcessDescription()

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

virtual

Reimplemented from G4VProcess.

Definition at line 443 of file G4MuonicAtomDecay.cc.

{
  outFile << "MuonicAtom process where Muon decays in orbit or is captured by the nucleus." <<G4endl;
}

SetVerboseLevel()

void G4MuonicAtomDecay::SetVerboseLevel ( G4int  value )

inline

Definition at line 71 of file G4MuonicAtomDecay.hh.

{verboseLevel = value;}

---

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

• G4MuonicAtomDecay.hh
• G4MuonicAtomDecay.cc