4 // See the class description in the header file.
6 #include "TG4SpecialCuts.h"
7 #include "TG4CutVector.h"
10 #include <G4UserLimits.hh>
12 #include <G4EnergyLossTables.hh>
14 TG4SpecialCuts::TG4SpecialCuts(TG3ParticleWSP particle, TG4CutVector* cutVector,
15 const G4String& processName)
16 : G4UserSpecialCuts(processName),
22 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForGamma;
23 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForGamma;
25 case kElectron: case kEplus:
26 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForElectron;
27 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForElectron;
30 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForHadron;
31 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForHadron;
34 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForNeutralHadron;
35 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForNeutralHadron;
38 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForMuon;
39 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForMuon;
42 fPtrMinEkineInCutVector = &TG4CutVector::GetMinEkineForOther;
43 fPtrMinEkineInLimits = &TG4Limits::GetMinEkineForOther;
45 case kNofParticlesWSP:
46 TG4Globals::Exception("TG4SpecialCuts: Wrong particle specification.");
51 TG4SpecialCuts::TG4SpecialCuts() {
55 TG4SpecialCuts::TG4SpecialCuts(const TG4SpecialCuts& right) {
57 TG4Globals::Exception(
58 "TG4SpecialCuts is protected from copying.");
61 TG4SpecialCuts::~TG4SpecialCuts() {
67 TG4SpecialCuts& TG4SpecialCuts::operator=(const TG4SpecialCuts& right)
69 // check assignement to self
70 if (this == &right) return *this;
72 TG4Globals::Exception(
73 "TG4SpecialCuts is protected from assigning.");
80 G4double TG4SpecialCuts::PostStepGetPhysicalInteractionLength(
81 const G4Track& track, G4double previousStepSize,
82 G4ForceCondition* condition)
84 // Returns the Step-size (actual length) which is allowed
89 *condition = NotForced;
91 G4double proposedStep = DBL_MAX;
93 = track.GetVolume()->GetLogicalVolume()->GetUserLimits();
97 = (limits->GetUserMaxTrackLength(track) - track.GetTrackLength());
98 if (proposedStep < 0.) return 0.;
102 = (track.GetDynamicParticle()->GetTotalMomentum()) /
103 (track.GetTotalEnergy());
105 = (limits->GetUserMaxTime(track) - track.GetGlobalTime());
106 G4double temp = beta*c_light*time;
107 if (temp < 0.) return 0.;
108 if (proposedStep > temp) proposedStep = temp;
110 // min remaining range
111 G4ParticleDefinition* particle = track.GetDefinition();
112 G4double kinEnergy = track.GetKineticEnergy();
113 G4Material* material = track.GetMaterial();
115 = G4EnergyLossTables::GetRange(particle, kinEnergy, material);
116 temp = (rangeNow - limits->GetUserMinRange(track));
117 if (temp < 0.) return 0.;
118 if (proposedStep > temp) proposedStep = temp;
120 // min kinetic energy (from limits)
121 // the kin energy cut can be applied only in case
122 // G4EnergyLossTables are defined for the particle
123 if (G4EnergyLossTables::GetDEDXTable(particle)) {
124 TG4Limits* tg4Limits = dynamic_cast<TG4Limits*>(limits);
126 G4String text = "TG4SpecialCuts::PostStepGetPhysicalInteractionLength:\n";
127 text = text + " Unknown limits type.";
128 TG4Globals::Exception(text);
131 = (tg4Limits->*fPtrMinEkineInLimits)(track);
133 = G4EnergyLossTables::GetRange(particle, minEkine, material);
134 temp = rangeNow - minR;
135 if (temp < 0.) return 0.;
136 if (proposedStep > temp) proposedStep = temp;
140 else if (fCutVector) {
141 // min kinetic energy (from cut vector)
142 // the kin energy cut can be applied only in case
143 // G4EnergyLossTables are defined for the particle
144 G4ParticleDefinition* particle = track.GetDefinition();
145 if (G4EnergyLossTables::GetDEDXTable(particle)) {
146 G4double kinEnergy = track.GetKineticEnergy();
147 G4Material* material = track.GetMaterial();
149 = G4EnergyLossTables::GetRange(particle, kinEnergy, material);
151 = (fCutVector->*fPtrMinEkineInCutVector)(track);
153 = G4EnergyLossTables::GetRange(particle, minEkine, material);
154 G4double temp = rangeNow - minR;
155 if (temp < 0.) return 0.;
156 if (proposedStep > temp) proposedStep = temp;
162 G4VParticleChange* TG4SpecialCuts::PostStepDoIt(const G4Track& track,
165 // Kills the current particle, if requested by G4UserLimits.
168 aParticleChange.Initialize(track);
169 aParticleChange.SetEnergyChange(0.) ;
170 aParticleChange.SetLocalEnergyDeposit(track.GetKineticEnergy()) ;
171 aParticleChange.SetStatusChange(fStopAndKill);
172 return &aParticleChange;