MinEkineCuts Class Reference

#include <MinEkineCuts.hh>

Inheritance diagram for MinEkineCuts:

Public Member Functions
	MinEkineCuts (const G4String &processName="MinEkineCuts")
virtual	~MinEkineCuts ()
virtual G4double	PostStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4ForceCondition *condition)

---

Detailed Description

Definition at line 45 of file MinEkineCuts.hh.

---

Constructor & Destructor Documentation

MinEkineCuts::MinEkineCuts

(

const G4String &

processName = "MinEkineCuts"

)

Definition at line 43 of file MinEkineCuts.cc.

References fUserDefined, G4cout, G4endl, G4VProcess::GetProcessName(), G4VProcess::SetProcessType(), and G4VProcess::verboseLevel.

  : SpecialCuts(aName)
{
   if (verboseLevel>1) {
     G4cout << GetProcessName() << " is created "<< G4endl;
   }
   SetProcessType(fUserDefined);
}

MinEkineCuts::~MinEkineCuts

(

)

[virtual]

Definition at line 52 of file MinEkineCuts.cc.

{}

---

Member Function Documentation

G4double MinEkineCuts::PostStepGetPhysicalInteractionLength	(	const G4Track &	track,
		G4double	previousStepSize,
		G4ForceCondition *	condition
	)			[virtual]

Reimplemented from SpecialCuts.

Definition at line 59 of file MinEkineCuts.cc.

References DBL_MAX, G4Track::GetDefinition(), G4Track::GetDynamicParticle(), G4DynamicParticle::GetKineticEnergy(), G4VPhysicalVolume::GetLogicalVolume(), G4Track::GetMaterialCutsCouple(), G4ParticleDefinition::GetPDGCharge(), G4EnergyLossTables::GetRange(), G4LogicalVolume::GetUserLimits(), G4UserLimits::GetUserMinEkine(), G4Track::GetVolume(), and NotForced.

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

   G4double     proposedStep = DBL_MAX;
   // get the pointer to UserLimits
   G4UserLimits* pUserLimits = aTrack.GetVolume()->GetLogicalVolume()->GetUserLimits();
   const G4DynamicParticle* aParticle = aTrack.GetDynamicParticle();
   G4ParticleDefinition* aParticleDef = aTrack.GetDefinition();
  
   if (pUserLimits && aParticleDef->GetPDGCharge() != 0.0) {
     //min kinetic energy
     G4double temp = DBL_MAX;
     G4double    eKine     = aParticle->GetKineticEnergy();
     const G4MaterialCutsCouple* couple = aTrack.GetMaterialCutsCouple();
     G4double eMin = pUserLimits->GetUserMinEkine(aTrack);

     G4double    rangeNow = DBL_MAX;

     rangeNow = G4EnergyLossTables::GetRange(aParticleDef,eKine,couple);

     if (eKine < eMin ) {
       proposedStep = 0.;
     } else {
       // charged particles only
       G4double rangeMin = G4EnergyLossTables::GetRange(aParticleDef,eMin,couple);
       temp = rangeNow - rangeMin;
       if (proposedStep > temp) proposedStep = temp;        
     }
   }
   return proposedStep;
}

---

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

• MinEkineCuts.hh
• MinEkineCuts.cc

---