G4UserSpecialCuts Class Reference

#include <G4UserSpecialCuts.hh>

Inheritance diagram for G4UserSpecialCuts:

Public Member Functions
	G4UserSpecialCuts (const G4String &processName="UserSpecialCut")
virtual	~G4UserSpecialCuts ()
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 &)

---

Detailed Description

Definition at line 47 of file G4UserSpecialCuts.hh.

---

Constructor & Destructor Documentation

G4UserSpecialCuts::G4UserSpecialCuts

(

const G4String &

processName = "UserSpecialCut"

)

Definition at line 53 of file G4UserSpecialCuts.cc.

References G4cout, G4endl, G4VProcess::GetProcessName(), G4LossTableManager::Instance(), G4VProcess::SetProcessSubType(), USER_SPECIAL_CUTS, and G4VProcess::verboseLevel.

  : G4VProcess(aName, fGeneral  )
{
  // set Process Sub Type
  SetProcessSubType(static_cast<int>(USER_SPECIAL_CUTS));

  if (verboseLevel>0)
    {
      G4cout << GetProcessName() << " is created "<< G4endl;
    }
  theLossTableManager = G4LossTableManager::Instance();
}

G4UserSpecialCuts::~G4UserSpecialCuts

(

)

[virtual]

Definition at line 68 of file G4UserSpecialCuts.cc.

{
}

---

Member Function Documentation

virtual G4VParticleChange* G4UserSpecialCuts::AlongStepDoIt	(	const G4Track &	,
		const G4Step &
	)			[inline, virtual]

Implements G4VProcess.

Definition at line 90 of file G4UserSpecialCuts.hh.

                              {return 0;};

virtual G4double G4UserSpecialCuts::AlongStepGetPhysicalInteractionLength	(	const G4Track &	,
		G4double	,
		G4double	,
		G4double &	,
		G4GPILSelection *
	)			[inline, virtual]

Implements G4VProcess.

Definition at line 81 of file G4UserSpecialCuts.hh.

                             { return -1.0; };

virtual G4VParticleChange* G4UserSpecialCuts::AtRestDoIt	(	const G4Track &	,
		const G4Step &
	)			[inline, virtual]

Implements G4VProcess.

Definition at line 75 of file G4UserSpecialCuts.hh.

                             {return 0;};

virtual G4double G4UserSpecialCuts::AtRestGetPhysicalInteractionLength	(	const G4Track &	,
		G4ForceCondition *
	)			[inline, virtual]

Implements G4VProcess.

Definition at line 69 of file G4UserSpecialCuts.hh.

                             { return -1.0; };

G4VParticleChange * G4UserSpecialCuts::PostStepDoIt	(	const G4Track &	,
		const G4Step &
	)			[virtual]

Implements G4VProcess.

Definition at line 139 of file G4UserSpecialCuts.cc.

References G4VProcess::aParticleChange, fStopAndKill, G4Track::GetKineticEnergy(), G4ParticleChange::Initialize(), G4ParticleChange::ProposeEnergy(), G4VParticleChange::ProposeLocalEnergyDeposit(), and G4VParticleChange::ProposeTrackStatus().

{
   aParticleChange.Initialize(aTrack);
   aParticleChange.ProposeEnergy(0.) ;
   aParticleChange.ProposeLocalEnergyDeposit(aTrack.GetKineticEnergy()) ;
   aParticleChange.ProposeTrackStatus(fStopAndKill);
   return &aParticleChange;
}

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

Implements G4VProcess.

Definition at line 82 of file G4UserSpecialCuts.cc.

References DBL_MAX, DBL_MIN, G4Track::GetKineticEnergy(), G4VPhysicalVolume::GetLogicalVolume(), G4ParticleDefinition::GetPDGCharge(), G4ParticleDefinition::GetPDGMass(), G4LossTableManager::GetRange(), G4LogicalVolume::GetUserLimits(), G4UserLimits::GetUserMaxTime(), G4UserLimits::GetUserMaxTrackLength(), G4UserLimits::GetUserMinRange(), G4Track::GetVolume(), and NotForced.

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

   G4double ProposedStep = DBL_MAX;
   G4UserLimits* pUserLimits =
                 aTrack.GetVolume()->GetLogicalVolume()->GetUserLimits();
   if (pUserLimits)
   {
     // check max kinetic energy first
     //
     G4double Ekine = aTrack.GetKineticEnergy();
     if(Ekine <= pUserLimits->GetUserMinEkine(aTrack)) { return 0.; }

     // max track length
     //
     ProposedStep = (pUserLimits->GetUserMaxTrackLength(aTrack)
                   - aTrack.GetTrackLength());
     if (ProposedStep < 0.) { return 0.; }

     // max time limit
     //
     G4double tlimit = pUserLimits->GetUserMaxTime(aTrack);
     if(tlimit < DBL_MAX) {
       G4double beta  = (aTrack.GetDynamicParticle()->GetTotalMomentum())
         /(aTrack.GetTotalEnergy());
       G4double dTime = (tlimit - aTrack.GetGlobalTime());
       G4double temp  = beta*c_light*dTime;
       if (temp < 0.) { return 0.; }
       if (ProposedStep > temp) { ProposedStep = temp; }
     }
                 
     // min remaining range 
     // (only for charged particle except for chargedGeantino)
     //
     G4double Rmin = pUserLimits->GetUserMinRange(aTrack);
     if (Rmin > DBL_MIN) {
       G4ParticleDefinition* Particle = aTrack.GetDefinition();
       if ( (Particle->GetPDGCharge() != 0.) && (Particle->GetPDGMass() > 0.0))
       {
         const G4MaterialCutsCouple* couple = aTrack.GetMaterialCutsCouple();
         G4double RangeNow = theLossTableManager->GetRange(Particle,Ekine,couple);
         G4double temp = RangeNow - Rmin;
         if (temp < 0.) { return 0.; }
         if (ProposedStep > temp) { ProposedStep = temp; }
       }         
     }
   }
   return ProposedStep;
}

---

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

• G4UserSpecialCuts.hh
• G4UserSpecialCuts.cc

---