Geant4  10.03.p01
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | Protected Member Functions | List of all members
G4ContinuousGainOfEnergy Class Reference

#include <G4ContinuousGainOfEnergy.hh>

Inheritance diagram for G4ContinuousGainOfEnergy:
Inheritance graph
[legend]
Collaboration diagram for G4ContinuousGainOfEnergy:
Collaboration graph
[legend]

Public Member Functions

 G4ContinuousGainOfEnergy (const G4String &name="EnergyGain", G4ProcessType type=fElectromagnetic)
 
virtual ~G4ContinuousGainOfEnergy ()
 
void PreparePhysicsTable (const G4ParticleDefinition &)
 
void BuildPhysicsTable (const G4ParticleDefinition &)
 
G4VParticleChange * AlongStepDoIt (const G4Track &, const G4Step &)
 
void SetLossFluctuations (G4bool val)
 
void SetIsIntegral (G4bool val)
 
void SetDirectEnergyLossProcess (G4VEnergyLossProcess *aProcess)
 
void SetDirectParticle (G4ParticleDefinition *p)
 
- Public Member Functions inherited from G4VContinuousProcess
 G4VContinuousProcess (const G4String &, G4ProcessType aType=fNotDefined)
 
 G4VContinuousProcess (G4VContinuousProcess &)
 
virtual ~G4VContinuousProcess ()
 
virtual G4double AlongStepGetPhysicalInteractionLength (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &proposedSafety, G4GPILSelection *selection)
 
virtual G4double PostStepGetPhysicalInteractionLength (const G4Track &, G4double, G4ForceCondition *)
 
virtual G4double AtRestGetPhysicalInteractionLength (const G4Track &, G4ForceCondition *)
 
virtual G4VParticleChange * AtRestDoIt (const G4Track &, const G4Step &)
 
virtual G4VParticleChange * PostStepDoIt (const G4Track &, const G4Step &)
 
- Public Member Functions inherited from G4VProcess
 G4VProcess (const G4String &aName="NoName", G4ProcessType aType=fNotDefined)
 
 G4VProcess (const G4VProcess &right)
 
virtual ~G4VProcess ()
 
G4int operator== (const G4VProcess &right) const
 
G4int operator!= (const G4VProcess &right) const
 
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 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
 
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 G4double GetContinuousStepLimit (const G4Track &track, G4double previousStepSize, G4double currentMinimumStep, G4double &currentSafety)
 
- Protected Member Functions inherited from G4VContinuousProcess
void SetGPILSelection (G4GPILSelection selection)
 
G4GPILSelection GetGPILSelection () const
 
- Protected Member Functions inherited from G4VProcess
void SubtractNumberOfInteractionLengthLeft (G4double previousStepSize)
 
void ClearNumberOfInteractionLengthLeft ()
 

Additional Inherited Members

- Static Public Member Functions inherited from G4VProcess
static const G4String & GetProcessTypeName (G4ProcessType)
 
- Protected Attributes inherited from G4VProcess
const G4ProcessManager * aProcessManager
 
G4VParticleChange * pParticleChange
 
G4ParticleChange aParticleChange
 
G4double theNumberOfInteractionLengthLeft
 
G4double currentInteractionLength
 
G4double theInitialNumberOfInteractionLength
 
G4String theProcessName
 
G4String thePhysicsTableFileName
 
G4ProcessType theProcessType
 
G4int theProcessSubType
 
G4double thePILfactor
 
G4bool enableAtRestDoIt
 
G4bool enableAlongStepDoIt
 
G4bool enablePostStepDoIt
 
G4int verboseLevel
 

Detailed Description

Definition at line 69 of file G4ContinuousGainOfEnergy.hh.

Constructor & Destructor Documentation

G4ContinuousGainOfEnergy::G4ContinuousGainOfEnergy ( const G4String &  name = "EnergyGain",
G4ProcessType  type = fElectromagnetic 
)

Definition at line 45 of file G4ContinuousGainOfEnergy.cc.

46  : G4VContinuousProcess(name, type)
47 {
48 
49 
50  linLossLimit=0.05;
51  lossFluctuationArePossible =true;
52  lossFluctuationFlag=true;
53  is_integral = false;
54 
55  //Will be properly set in SetDirectParticle()
56  IsIon=false;
57  massRatio =1.;
58  chargeSqRatio=1.;
59  preStepChargeSqRatio=1.;
60 
61  //Some initialization
62  currentCoupleIndex=9999999;
63  currentCutInRange=0.;
64  currentMaterialIndex=9999999;
65  currentTcut=0.;
66  preStepKinEnergy=0.;
67  preStepRange=0.;
68  preStepScaledKinEnergy=0.;
69 
70  currentCouple=0;
71 }
G4ContinuousGainOfEnergy::~G4ContinuousGainOfEnergy ( )
virtual

Definition at line 75 of file G4ContinuousGainOfEnergy.cc.

76 {
77 
78 }

Member Function Documentation

G4VParticleChange * G4ContinuousGainOfEnergy::AlongStepDoIt ( const G4Track &  track,
const G4Step &  step 
)
virtual

Reimplemented from G4VContinuousProcess.

Definition at line 115 of file G4ContinuousGainOfEnergy.cc.

117 {
118 
119  //Caution in this method the step length should be the true step length
120  // A problem is that this is compute by the multiple scattering that does not know the energy at the end of the adjoint step. This energy is used during the
121  //Forward sim. Nothing we can really do against that at this time. This is inherent to the MS method
122  //
123 
124 
125 
126  aParticleChange.Initialize(track);
127 
128  // Get the actual (true) Step length
129  //----------------------------------
130  G4double length = step.GetStepLength();
131  G4double degain = 0.0;
132 
133 
134 
135  // Compute this for weight change after continuous energy loss
136  //-------------------------------------------------------------
137  G4double DEDX_before = theDirectEnergyLossProcess->GetDEDX(preStepKinEnergy, currentCouple);
138 
139 
140 
141  // For the fluctuation we generate a new dynamic particle with energy =preEnergy+egain
142  // and then compute the fluctuation given in the direct case.
143  //-----------------------------------------------------------------------
144  G4DynamicParticle* dynParticle = new G4DynamicParticle();
145  *dynParticle = *(track.GetDynamicParticle());
146  dynParticle->SetDefinition(theDirectPartDef);
147  G4double Tkin = dynParticle->GetKineticEnergy();
148 
149 
150  size_t n=1;
151  if (is_integral ) n=10;
152  n=1;
153  G4double dlength= length/n;
154  for (size_t i=0;i<n;i++) {
155  if (Tkin != preStepKinEnergy && IsIon) {
156  chargeSqRatio = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,Tkin);
157  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,chargeSqRatio);
158 
159  }
160 
161  G4double r = theDirectEnergyLossProcess->GetRange(Tkin, currentCouple);
162  if( dlength <= linLossLimit * r ) {
163  degain = DEDX_before*dlength;
164  }
165  else {
166  G4double x = r + dlength;
167  //degain = theDirectEnergyLossProcess->GetKineticEnergy(x,currentCouple) - theDirectEnergyLossProcess->GetKineticEnergy(r,currentCouple);
168  G4double E = theDirectEnergyLossProcess->GetKineticEnergy(x,currentCouple);
169  if (IsIon){
170  chargeSqRatio = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,E);
171  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,chargeSqRatio);
172  G4double x1= theDirectEnergyLossProcess->GetRange(E, currentCouple);
173 
174  // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
175  G4int ii=0;
176  const G4int iimax = 100;
177  while (std::abs(x-x1)>0.01*x) {
178  E = theDirectEnergyLossProcess->GetKineticEnergy(x,currentCouple);
179  chargeSqRatio = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,E);
180  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,chargeSqRatio);
181  x1= theDirectEnergyLossProcess->GetRange(E, currentCouple);
182  ++ii;
183  if(ii >= iimax) { break; }
184  }
185  }
186 
187  degain=E-Tkin;
188 
189 
190 
191  }
192  //G4cout<<degain<<G4endl;
193  G4double tmax = currentModel->MaxSecondaryKinEnergy(dynParticle);
194  tmax = std::min(tmax,currentTcut);
195 
196 
197  dynParticle->SetKineticEnergy(Tkin+degain);
198 
199  // Corrections, which cannot be tabulated for ions
200  //----------------------------------------
201  G4double esecdep=0;//not used in most models
202  currentModel->CorrectionsAlongStep(currentCouple, dynParticle, degain,esecdep, dlength);
203 
204  // Sample fluctuations
205  //-------------------
206 
207 
208  G4double deltaE =0.;
209  if (lossFluctuationFlag ) {
210  deltaE = currentModel->GetModelOfFluctuations()->
211  SampleFluctuations(currentCouple,dynParticle,tmax,dlength,degain)-degain;
212  }
213 
214  G4double egain=degain+deltaE;
215  if (egain <=0) egain=degain;
216  Tkin+=egain;
217  dynParticle->SetKineticEnergy(Tkin);
218  }
219 
220 
221 
222 
223 
224  delete dynParticle;
225 
226  if (IsIon){
227  chargeSqRatio = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,Tkin);
228  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,chargeSqRatio);
229 
230  }
231 
232  G4double DEDX_after = theDirectEnergyLossProcess->GetDEDX(Tkin, currentCouple);
233 
234 
235  G4double weight_correction=DEDX_after/DEDX_before;
236 
237 
238  aParticleChange.ProposeEnergy(Tkin);
239 
240 
241  //Caution!!!
242  // It is important to select the weight of the post_step_point
243  // as the current weight and not the weight of the track, as t
244  // the weight of the track is changed after having applied all
245  // the along_step_do_it.
246 
247  // G4double new_weight=weight_correction*track.GetWeight(); //old
248  G4double new_weight=weight_correction*step.GetPostStepPoint()->GetWeight();
249  aParticleChange.SetParentWeightByProcess(false);
250  aParticleChange.ProposeParentWeight(new_weight);
251 
252 
253  return &aParticleChange;
254 
255 }
G4VEnergyLossProcess::GetKineticEnergy
G4double GetKineticEnergy(G4double &range, const G4MaterialCutsCouple *)
Definition: G4VEnergyLossProcess.hh:837
G4VEmModel::MaxSecondaryKinEnergy
G4double MaxSecondaryKinEnergy(const G4DynamicParticle *dynParticle)
Definition: G4VEmModel.hh:481
G4StepPoint::GetWeight
G4double GetWeight() const
G4VEmModel::CorrectionsAlongStep
virtual void CorrectionsAlongStep(const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
Definition: G4VEmModel.cc:370
G4VEnergyLossProcess::SetDynamicMassCharge
void SetDynamicMassCharge(G4double massratio, G4double charge2ratio)
Definition: G4VEnergyLossProcess.hh:652
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4Step::GetStepLength
G4double GetStepLength() const
G4Track::GetDynamicParticle
const G4DynamicParticle * GetDynamicParticle() const
G4VEnergyLossProcess::GetDEDX
G4double GetDEDX(G4double &kineticEnergy, const G4MaterialCutsCouple *)
Definition: G4VEnergyLossProcess.hh:771
G4VParticleChange::ProposeParentWeight
void ProposeParentWeight(G4double finalWeight)
G4VEmModel::GetModelOfFluctuations
G4VEmFluctuationModel * GetModelOfFluctuations()
Definition: G4VEmModel.hh:609
G4int
int G4int
Definition: G4Types.hh:78
G4VParticleChange::SetParentWeightByProcess
void SetParentWeightByProcess(G4bool)
Definition: G4VParticleChange.cc:522
G4DynamicParticle::SetKineticEnergy
void SetKineticEnergy(G4double aEnergy)
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:353
G4ParticleChange::Initialize
virtual void Initialize(const G4Track &)
Definition: G4ParticleChange.cc:228
G4Step::GetPostStepPoint
G4StepPoint * GetPostStepPoint() const
G4ParticleChange::ProposeEnergy
void ProposeEnergy(G4double finalEnergy)
G4VProcess::aParticleChange
G4ParticleChange aParticleChange
Definition: G4VProcess.hh:289
G4INCL::Math::min
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
Definition: G4INCLGlobals.hh:117
G4double
double G4double
Definition: G4Types.hh:76
G4DynamicParticle::SetDefinition
void SetDefinition(const G4ParticleDefinition *aParticleDefinition)
Definition: G4DynamicParticle.cc:305
G4VEnergyLossProcess::GetRange
G4double GetRange(G4double &kineticEnergy, const G4MaterialCutsCouple *)
Definition: G4VEnergyLossProcess.hh:791

Here is the call graph for this function:

void G4ContinuousGainOfEnergy::BuildPhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 92 of file G4ContinuousGainOfEnergy.cc.

93 {//theDirectEnergyLossProcess->BuildPhysicsTable(part);
94 ;
95 }
G4double G4ContinuousGainOfEnergy::GetContinuousStepLimit ( const G4Track &  track,
G4double  previousStepSize,
G4double  currentMinimumStep,
G4double &  currentSafety 
)
protectedvirtual

Implements G4VContinuousProcess.

Definition at line 268 of file G4ContinuousGainOfEnergy.cc.

270 {
271  G4double x = DBL_MAX;
272  x=.1*mm;
273 
274 
275  DefineMaterial(track.GetMaterialCutsCouple());
276 
277  preStepKinEnergy = track.GetKineticEnergy();
278  preStepScaledKinEnergy = track.GetKineticEnergy()*massRatio;
279  currentModel = theDirectEnergyLossProcess->SelectModelForMaterial(preStepScaledKinEnergy,currentCoupleIndex);
280  G4double emax_model=currentModel->HighEnergyLimit();
281  if (IsIon) {
282  chargeSqRatio = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,preStepKinEnergy);
283  preStepChargeSqRatio = chargeSqRatio;
284  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,preStepChargeSqRatio);
285  }
286 
287 
288  G4double maxE =1.1*preStepKinEnergy;
289  /*if (preStepKinEnergy< 0.05*MeV) maxE =2.*preStepKinEnergy;
290  else if (preStepKinEnergy< 0.1*MeV) maxE =1.5*preStepKinEnergy;
291  else if (preStepKinEnergy< 0.5*MeV) maxE =1.25*preStepKinEnergy;*/
292 
293  if (preStepKinEnergy < currentTcut) maxE = std::min(currentTcut,maxE);
294 
295  maxE=std::min(emax_model*1.001,maxE);
296 
297  preStepRange = theDirectEnergyLossProcess->GetRange(preStepKinEnergy, currentCouple);
298 
299  if (IsIon) {
300  G4double chargeSqRatioAtEmax = currentModel->GetChargeSquareRatio(theDirectPartDef,currentMaterial,maxE);
301  theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,chargeSqRatioAtEmax);
302  }
303 
304  G4double r1 = theDirectEnergyLossProcess->GetRange(maxE, currentCouple);
305 
306  if (IsIon) theDirectEnergyLossProcess->SetDynamicMassCharge(massRatio,preStepChargeSqRatio);
307 
308 
309 
310  x=r1-preStepRange;
311  x=std::max(r1-preStepRange,0.001*mm);
312 
313  return x;
314 
315 
316 }
mm
static constexpr double mm
Definition: G4SIunits.hh:115
G4VEnergyLossProcess::SetDynamicMassCharge
void SetDynamicMassCharge(G4double massratio, G4double charge2ratio)
Definition: G4VEnergyLossProcess.hh:652
G4VEmModel::HighEnergyLimit
G4double HighEnergyLimit() const
Definition: G4VEmModel.hh:633
G4Track::GetMaterialCutsCouple
const G4MaterialCutsCouple * GetMaterialCutsCouple() const
G4Track::GetKineticEnergy
G4double GetKineticEnergy() const
G4VEmModel::GetChargeSquareRatio
virtual G4double GetChargeSquareRatio(const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
Definition: G4VEmModel.cc:353
G4INCL::Math::max
T max(const T t1, const T t2)
brief Return the largest of the two arguments
Definition: G4INCLGlobals.hh:112
G4INCL::Math::min
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
Definition: G4INCLGlobals.hh:117
G4double
double G4double
Definition: G4Types.hh:76
DBL_MAX
#define DBL_MAX
Definition: templates.hh:83
G4VEnergyLossProcess::SelectModelForMaterial
G4VEmModel * SelectModelForMaterial(G4double kinEnergy, size_t &idx) const
Definition: G4VEnergyLossProcess.hh:627
G4VEnergyLossProcess::GetRange
G4double GetRange(G4double &kineticEnergy, const G4MaterialCutsCouple *)
Definition: G4VEnergyLossProcess.hh:791

Here is the call graph for this function:

void G4ContinuousGainOfEnergy::PreparePhysicsTable ( const G4ParticleDefinition &  )
virtual

Reimplemented from G4VProcess.

Definition at line 82 of file G4ContinuousGainOfEnergy.cc.

84 {//theDirectEnergyLossProcess->PreparePhysicsTable(part);
85 
86 ;
87 }
void G4ContinuousGainOfEnergy::SetDirectEnergyLossProcess ( G4VEnergyLossProcess *  aProcess)
inline

Definition at line 112 of file G4ContinuousGainOfEnergy.hh.

112 {theDirectEnergyLossProcess=aProcess;};
void G4ContinuousGainOfEnergy::SetDirectParticle ( G4ParticleDefinition *  p)

Definition at line 99 of file G4ContinuousGainOfEnergy.cc.

100 {theDirectPartDef=p;
101  if (theDirectPartDef->GetParticleType()== "nucleus") {
102  IsIon=true;
103  massRatio = proton_mass_c2/theDirectPartDef->GetPDGMass();
104  G4double q=theDirectPartDef->GetPDGCharge();
105  chargeSqRatio=q*q;
106 
107 
108  }
109 
110 }
CLHEP::proton_mass_c2
static constexpr double proton_mass_c2
Definition: PhysicalConstants.h:80
p
const char * p
Definition: xmltok.h:285
G4ParticleDefinition::GetParticleType
const G4String & GetParticleType() const
Definition: G4ParticleDefinition.hh:141
G4ParticleDefinition::GetPDGMass
G4double GetPDGMass() const
Definition: G4ParticleDefinition.hh:122
G4double
double G4double
Definition: G4Types.hh:76
G4ParticleDefinition::GetPDGCharge
G4double GetPDGCharge() const
Definition: G4ParticleDefinition.hh:124

Here is the call graph for this function:

void G4ContinuousGainOfEnergy::SetIsIntegral ( G4bool  val)
inline

Definition at line 110 of file G4ContinuousGainOfEnergy.hh.

110 {is_integral= val;}
void G4ContinuousGainOfEnergy::SetLossFluctuations ( G4bool  val)

Definition at line 258 of file G4ContinuousGainOfEnergy.cc.

259 {
260  if(val && !lossFluctuationArePossible) return;
261  lossFluctuationFlag = val;
262 }
The documentation for this class was generated from the following files: