Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4UrbanMscModel.hh

Version: [ ReleaseNotes ] [ 1.0 ] [ 1.1 ] [ 2.0 ] [ 3.0 ] [ 3.1 ] [ 3.2 ] [ 4.0 ] [ 4.0.p1 ] [ 4.0.p2 ] [ 4.1 ] [ 4.1.p1 ] [ 5.0 ] [ 5.0.p1 ] [ 5.1 ] [ 5.1.p1 ] [ 5.2 ] [ 5.2.p1 ] [ 5.2.p2 ] [ 6.0 ] [ 6.0.p1 ] [ 6.1 ] [ 6.2 ] [ 6.2.p1 ] [ 6.2.p2 ] [ 7.0 ] [ 7.0.p1 ] [ 7.1 ] [ 7.1.p1 ] [ 8.0 ] [ 8.0.p1 ] [ 8.1 ] [ 8.1.p1 ] [ 8.1.p2 ] [ 8.2 ] [ 8.2.p1 ] [ 8.3 ] [ 8.3.p1 ] [ 8.3.p2 ] [ 9.0 ] [ 9.0.p1 ] [ 9.0.p2 ] [ 9.1 ] [ 9.1.p1 ] [ 9.1.p2 ] [ 9.1.p3 ] [ 9.2 ] [ 9.2.p1 ] [ 9.2.p2 ] [ 9.2.p3 ] [ 9.2.p4 ] [ 9.3 ] [ 9.3.p1 ] [ 9.3.p2 ] [ 9.4 ] [ 9.4.p1 ] [ 9.4.p2 ] [ 9.4.p3 ] [ 9.4.p4 ] [ 9.5 ] [ 9.5.p1 ] [ 9.5.p2 ] [ 9.6 ] [ 9.6.p1 ] [ 9.6.p2 ] [ 9.6.p3 ] [ 9.6.p4 ] [ 10.0 ] [ 10.0.p1 ] [ 10.0.p2 ] [ 10.0.p3 ] [ 10.0.p4 ] [ 10.1 ] [ 10.1.p1 ] [ 10.1.p2 ] [ 10.1.p3 ] [ 10.2 ] [ 10.2.p1 ] [ 10.2.p2 ] [ 10.2.p3 ] [ 10.3 ] [ 10.3.p1 ] [ 10.3.p2 ] [ 10.3.p3 ] [ 10.4 ] [ 10.4.p1 ] [ 10.4.p2 ] [ 10.4.p3 ] [ 10.5 ] [ 10.5.p1 ] [ 10.6 ] [ 10.6.p1 ] [ 10.6.p2 ] [ 10.6.p3 ] [ 10.7 ] [ 10.7.p1 ] [ 10.7.p2 ] [ 10.7.p3 ] [ 10.7.p4 ] [ 11.0 ] [ 11.0.p1 ] [ 11.0.p2 ] [ 11.0.p3, ] [ 11.0.p4 ] [ 11.1 ] [ 11.1.1 ] [ 11.1.2 ] [ 11.1.3 ] [ 11.2 ] [ 11.2.1 ] [ 11.2.2 ] [ 11.3.0 ]

Diff markup

Differences between /processes/electromagnetic/standard/include/G4UrbanMscModel.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4UrbanMscModel.hh (Version 11.2)


  1 //                                                  1 //
  2 // *******************************************      2 // ********************************************************************
  3 // * License and Disclaimer                         3 // * License and Disclaimer                                           *
  4 // *                                                4 // *                                                                  *
  5 // * The  Geant4 software  is  copyright of th      5 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
  6 // * the Geant4 Collaboration.  It is provided      6 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
  7 // * conditions of the Geant4 Software License      7 // * conditions of the Geant4 Software License,  included in the file *
  8 // * LICENSE and available at  http://cern.ch/      8 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
  9 // * include a list of copyright holders.           9 // * include a list of copyright holders.                             *
 10 // *                                               10 // *                                                                  *
 11 // * Neither the authors of this software syst     11 // * Neither the authors of this software system, nor their employing *
 12 // * institutes,nor the agencies providing fin     12 // * institutes,nor the agencies providing financial support for this *
 13 // * work  make  any representation or  warran     13 // * work  make  any representation or  warranty, express or implied, *
 14 // * regarding  this  software system or assum     14 // * regarding  this  software system or assume any liability for its *
 15 // * use.  Please see the license in the file      15 // * use.  Please see the license in the file  LICENSE  and URL above *
 16 // * for the full disclaimer and the limitatio     16 // * for the full disclaimer and the limitation of liability.         *
 17 // *                                               17 // *                                                                  *
 18 // * This  code  implementation is the result      18 // * This  code  implementation is the result of  the  scientific and *
 19 // * technical work of the GEANT4 collaboratio     19 // * technical work of the GEANT4 collaboration.                      *
 20 // * By using,  copying,  modifying or  distri     20 // * By using,  copying,  modifying or  distributing the software (or *
 21 // * any work based  on the software)  you  ag     21 // * any work based  on the software)  you  agree  to acknowledge its *
 22 // * use  in  resulting  scientific  publicati     22 // * use  in  resulting  scientific  publications,  and indicate your *
 23 // * acceptance of all terms of the Geant4 Sof     23 // * acceptance of all terms of the Geant4 Software license.          *
 24 // *******************************************     24 // ********************************************************************
 25 //                                                 25 //
 26 //                                                 26 //
 27 // -------------------------------------------     27 // -------------------------------------------------------------------
 28 //                                                 28 //
 29 //                                                 29 //
 30 // GEANT4 Class header file                        30 // GEANT4 Class header file
 31 //                                                 31 //
 32 //                                                 32 //
 33 // File name:     G4UrbanMscModel                  33 // File name:     G4UrbanMscModel
 34 //                                                 34 //
 35 // Author:        Laszlo Urban                     35 // Author:        Laszlo Urban
 36 //                                                 36 //
 37 // Creation date: 19.02.2013                       37 // Creation date: 19.02.2013
 38 //                                                 38 //
 39 // Created from G4UrbanMscModel96                  39 // Created from G4UrbanMscModel96 
 40 //                                                 40 //
 41 // New parametrization for theta0                  41 // New parametrization for theta0
 42 // Correction for very small step length           42 // Correction for very small step length
 43 //                                                 43 //
 44 // Class Description:                              44 // Class Description:
 45 //                                                 45 //
 46 // Implementation of the model of multiple sca     46 // Implementation of the model of multiple scattering based on
 47 // H.W.Lewis Phys Rev 78 (1950) 526 and L.Urba     47 // H.W.Lewis Phys Rev 78 (1950) 526 and L.Urban model
 48                                                    48 
 49 // -------------------------------------------     49 // -------------------------------------------------------------------
 50 //                                                 50 //
 51                                                    51 
 52 #ifndef G4UrbanMscModel_h                          52 #ifndef G4UrbanMscModel_h
 53 #define G4UrbanMscModel_h 1                        53 #define G4UrbanMscModel_h 1
 54                                                    54 
 55 //....oooOO0OOooo........oooOO0OOooo........oo     55 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 56                                                    56 
 57 #include <CLHEP/Units/SystemOfUnits.h>             57 #include <CLHEP/Units/SystemOfUnits.h>
 58                                                    58 
 59 #include "G4VMscModel.hh"                          59 #include "G4VMscModel.hh"
 60 #include "G4MscStepLimitType.hh"                   60 #include "G4MscStepLimitType.hh"
 61 #include "G4Log.hh"                                61 #include "G4Log.hh"
 62 #include "G4Exp.hh"                                62 #include "G4Exp.hh"
 63                                                    63 
 64 class G4ParticleChangeForMSC;                      64 class G4ParticleChangeForMSC;
 65 class G4SafetyHelper;                              65 class G4SafetyHelper;
 66                                                    66 
 67 //....oooOO0OOooo........oooOO0OOooo........oo     67 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 68                                                    68 
 69 class G4UrbanMscModel : public G4VMscModel         69 class G4UrbanMscModel : public G4VMscModel
 70 {                                                  70 {
 71                                                    71 
 72 public:                                            72 public:
 73                                                    73 
 74   explicit G4UrbanMscModel(const G4String& nam     74   explicit G4UrbanMscModel(const G4String& nam = "UrbanMsc");
 75                                                    75 
 76   ~G4UrbanMscModel() override;                     76   ~G4UrbanMscModel() override;
 77                                                    77 
 78   void Initialise(const G4ParticleDefinition*,     78   void Initialise(const G4ParticleDefinition*, 
 79       const G4DataVector&) override;               79       const G4DataVector&) override;
 80                                                    80 
 81   void StartTracking(G4Track*) override;           81   void StartTracking(G4Track*) override;
 82                                                    82 
 83   G4double                                         83   G4double 
 84   ComputeCrossSectionPerAtom(const G4ParticleD     84   ComputeCrossSectionPerAtom(const G4ParticleDefinition* particle,
 85            G4double KineticEnergy,                 85            G4double KineticEnergy,
 86            G4double AtomicNumber,                  86            G4double AtomicNumber,
 87            G4double AtomicWeight=0.,               87            G4double AtomicWeight=0., 
 88            G4double cut =0.,                       88            G4double cut =0.,
 89            G4double emax=DBL_MAX) override;        89            G4double emax=DBL_MAX) override;
 90                                                    90 
 91   G4ThreeVector& SampleScattering(const G4Thre     91   G4ThreeVector& SampleScattering(const G4ThreeVector&, 
 92           G4double safety) override;               92           G4double safety) override;
 93                                                    93 
 94   G4double ComputeTruePathLengthLimit(const G4     94   G4double ComputeTruePathLengthLimit(const G4Track& track,
 95                     G4double& currentMinimalSt     95                     G4double& currentMinimalStep) override;
 96                                                    96 
 97   G4double ComputeGeomPathLength(G4double true     97   G4double ComputeGeomPathLength(G4double truePathLength) override;
 98                                                    98 
 99   G4double ComputeTrueStepLength(G4double geom     99   G4double ComputeTrueStepLength(G4double geomStepLength) override;
100                                                   100 
101   G4double ComputeTheta0(G4double truePathLeng    101   G4double ComputeTheta0(G4double truePathLength, G4double KineticEnergy);
102                                                   102 
103   inline void SetDisplacementAlgorithm96(const    103   inline void SetDisplacementAlgorithm96(const G4bool val);
104                                                   104 
105   inline void SetPositronCorrection(const G4bo    105   inline void SetPositronCorrection(const G4bool val);
106                                                   106 
107   //  hide assignment operator                    107   //  hide assignment operator
108   G4UrbanMscModel & operator=(const  G4UrbanMs    108   G4UrbanMscModel & operator=(const  G4UrbanMscModel &right) = delete;
109   G4UrbanMscModel(const  G4UrbanMscModel&) = d    109   G4UrbanMscModel(const  G4UrbanMscModel&) = delete;
110                                                   110 
111 private:                                          111 private:
112                                                   112 
113   G4double SampleCosineTheta(G4double trueStep    113   G4double SampleCosineTheta(G4double trueStepLength, G4double KineticEnergy);
114                                                   114 
115   void SampleDisplacement(G4double sinTheta, G    115   void SampleDisplacement(G4double sinTheta, G4double phi);
116                                                   116 
117   void SampleDisplacementNew(G4double sinTheta    117   void SampleDisplacementNew(G4double sinTheta, G4double phi);
118                                                   118 
119   void InitialiseModelCache();                    119   void InitialiseModelCache();
120                                                   120 
121   inline void SetParticle(const G4ParticleDefi    121   inline void SetParticle(const G4ParticleDefinition*);
122                                                   122 
123   inline G4double Randomizetlimit();              123   inline G4double Randomizetlimit();
124                                                   124   
125   inline G4double SimpleScattering();             125   inline G4double SimpleScattering();
126                                                   126 
127   inline G4double ComputeStepmin();               127   inline G4double ComputeStepmin();
128                                                   128 
129   inline G4double ComputeTlimitmin();             129   inline G4double ComputeTlimitmin();
130                                                   130 
131   CLHEP::HepRandomEngine* rndmEngineMod;          131   CLHEP::HepRandomEngine* rndmEngineMod;
132                                                   132 
133   const G4ParticleDefinition* particle = nullp    133   const G4ParticleDefinition* particle = nullptr;
134   const G4ParticleDefinition* positron;           134   const G4ParticleDefinition* positron;
135   G4ParticleChangeForMSC* fParticleChange = nu    135   G4ParticleChangeForMSC* fParticleChange = nullptr;
136   const G4MaterialCutsCouple* couple = nullptr    136   const G4MaterialCutsCouple* couple = nullptr;
137                                                   137 
138   G4double mass;                                  138   G4double mass;
139   G4double charge,chargeSquare;                   139   G4double charge,chargeSquare;
140   G4double masslimite,fr;                         140   G4double masslimite,fr;
141                                                   141 
142   G4double taubig;                                142   G4double taubig;
143   G4double tausmall;                              143   G4double tausmall;
144   G4double taulim;                                144   G4double taulim;
145   G4double currentTau;                            145   G4double currentTau;
146   G4double tlimit;                                146   G4double tlimit;
147   G4double tlimitmin;                             147   G4double tlimitmin;
148   G4double tlimitminfix,tlimitminfix2;            148   G4double tlimitminfix,tlimitminfix2;
149   G4double tgeom;                                 149   G4double tgeom;
150                                                   150 
151   G4double geombig;                               151   G4double geombig;
152   G4double geommin;                               152   G4double geommin;
153   G4double geomlimit;                             153   G4double geomlimit;
154   G4double skindepth;                             154   G4double skindepth;
155   G4double smallstep;                             155   G4double smallstep;
156                                                   156 
157   G4double presafety;                             157   G4double presafety;
158                                                   158 
159   G4double lambda0;                               159   G4double lambda0;
160   G4double lambdaeff;                             160   G4double lambdaeff;
161   G4double tPathLength;                           161   G4double tPathLength;
162   G4double zPathLength;                           162   G4double zPathLength;
163   G4double par1,par2,par3;                        163   G4double par1,par2,par3;
164                                                   164 
165   G4double stepmin;                               165   G4double stepmin;
166                                                   166 
167   G4double currentKinEnergy;                      167   G4double currentKinEnergy;
168   G4double currentLogKinEnergy;                   168   G4double currentLogKinEnergy;
169   G4double currentRange;                          169   G4double currentRange; 
170   G4double rangeinit;                             170   G4double rangeinit;
171   G4double currentRadLength;                      171   G4double currentRadLength;
172                                                   172 
173   G4double drr,finalr;                            173   G4double drr,finalr;
174                                                   174 
175   G4double tlow;                                  175   G4double tlow;
176   G4double invmev;                                176   G4double invmev;
177   G4double xmeanth = 0.0;                         177   G4double xmeanth = 0.0;
178   G4double x2meanth = 1./3.;                      178   G4double x2meanth = 1./3.;
179   G4double rndmarray[2];                          179   G4double rndmarray[2];
180                                                   180 
181   struct mscData {                                181   struct mscData {
182     G4double Z23, sqrtZ, factmin;                 182     G4double Z23, sqrtZ, factmin;
183     G4double coeffth1, coeffth2;                  183     G4double coeffth1, coeffth2;
184     G4double coeffc1, coeffc2, coeffc3, coeffc    184     G4double coeffc1, coeffc2, coeffc3, coeffc4;
185     G4double stepmina, stepminb;                  185     G4double stepmina, stepminb;
186     G4double doverra, doverrb;                    186     G4double doverra, doverrb;
187     G4double posa, posb, posc, posd, pose;        187     G4double posa, posb, posc, posd, pose;
188   };                                              188   };
189   static std::vector<mscData*> msc;               189   static std::vector<mscData*> msc;
190                                                   190 
191   // index of G4MaterialCutsCouple                191   // index of G4MaterialCutsCouple
192   G4int idx = 0;                                  192   G4int idx = 0;
193                                                   193 
194   G4bool firstStep = true;                        194   G4bool firstStep = true;
195   G4bool insideskin = false;                      195   G4bool insideskin = false;
196                                                   196 
197   G4bool latDisplasmentbackup = false;            197   G4bool latDisplasmentbackup = false;
198   G4bool dispAlg96 = true;                        198   G4bool dispAlg96 = true;
199   G4bool fPosiCorrection = true;                  199   G4bool fPosiCorrection = true;
200   G4bool isFirstInstance = false;                 200   G4bool isFirstInstance = false;
201 };                                                201 };
202                                                   202 
203 //....oooOO0OOooo........oooOO0OOooo........oo    203 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
204 //....oooOO0OOooo........oooOO0OOooo........oo    204 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
205                                                   205 
206 inline                                            206 inline
207 void G4UrbanMscModel::SetParticle(const G4Part    207 void G4UrbanMscModel::SetParticle(const G4ParticleDefinition* p)
208 {                                                 208 {
209   if (p != particle) {                            209   if (p != particle) {
210     particle = p;                                 210     particle = p;
211     mass = p->GetPDGMass();                       211     mass = p->GetPDGMass();
212     charge = p->GetPDGCharge()/CLHEP::eplus;      212     charge = p->GetPDGCharge()/CLHEP::eplus;
213     chargeSquare = charge*charge;                 213     chargeSquare = charge*charge;
214   }                                               214   }
215 }                                                 215 }
216                                                   216 
217 //....oooOO0OOooo........oooOO0OOooo........oo    217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
218                                                   218 
219 inline G4double G4UrbanMscModel::Randomizetlim    219 inline G4double G4UrbanMscModel::Randomizetlimit()
220 {                                                 220 {
221   G4double res = tlimitmin;                       221   G4double res = tlimitmin;
222   if(tlimit > tlimitmin)                          222   if(tlimit > tlimitmin)
223   {                                               223   {
224     res = G4RandGauss::shoot(rndmEngineMod,tli    224     res = G4RandGauss::shoot(rndmEngineMod,tlimit,0.1*(tlimit-tlimitmin));
225     res = std::max(res, tlimitmin);               225     res = std::max(res, tlimitmin);
226   }                                               226   }
227   return res;                                     227   return res;
228 }                                                 228 }
229                                                   229 
230 //....oooOO0OOooo........oooOO0OOooo........oo    230 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
231                                                   231 
232 inline G4double G4UrbanMscModel::SimpleScatter    232 inline G4double G4UrbanMscModel::SimpleScattering()
233 {                                                 233 {
234   // 'large angle scattering'                     234   // 'large angle scattering'
235   // 2 model functions with correct xmean and     235   // 2 model functions with correct xmean and x2mean
236   const G4double a = (2.*xmeanth+9.*x2meanth-3    236   const G4double a = (2.*xmeanth+9.*x2meanth-3.)/(2.*xmeanth-3.*x2meanth+1.);
237   const G4double prob = (a+2.)*xmeanth/a;         237   const G4double prob = (a+2.)*xmeanth/a;
238                                                   238 
239   // sampling                                     239   // sampling
240   rndmEngineMod->flatArray(2, rndmarray);         240   rndmEngineMod->flatArray(2, rndmarray);
241   return (rndmarray[1] < prob) ?                  241   return (rndmarray[1] < prob) ? 
242     -1.+2.*G4Exp(G4Log(rndmarray[0])/(a+1.)) :    242     -1.+2.*G4Exp(G4Log(rndmarray[0])/(a+1.)) : -1.+2.*rndmarray[0];
243 }                                                 243 }
244                                                   244 
245 //....oooOO0OOooo........oooOO0OOooo........oo    245 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
246                                                   246 
247 inline G4double G4UrbanMscModel::ComputeStepmi    247 inline G4double G4UrbanMscModel::ComputeStepmin()
248 {                                                 248 {
249   // define stepmin using estimation of the ra    249   // define stepmin using estimation of the ratio 
250   // of lambda_elastic/lambda_transport           250   // of lambda_elastic/lambda_transport
251   const G4double rat = currentKinEnergy*invmev    251   const G4double rat = currentKinEnergy*invmev;
252   return lambda0*msc[idx]->factmin/               252   return lambda0*msc[idx]->factmin/
253     (0.002 + rat*(msc[idx]->stepmina + msc[idx    253     (0.002 + rat*(msc[idx]->stepmina + msc[idx]->stepminb*rat));
254 }                                                 254 }
255                                                   255 
256 //....oooOO0OOooo........oooOO0OOooo........oo    256 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
257                                                   257 
258 inline G4double G4UrbanMscModel::ComputeTlimit    258 inline G4double G4UrbanMscModel::ComputeTlimitmin()
259 {                                                 259 {
260   G4double x = (particle == positron) ?           260   G4double x = (particle == positron) ? 
261     0.7*msc[idx]->sqrtZ*stepmin : 0.87*msc[idx    261     0.7*msc[idx]->sqrtZ*stepmin : 0.87*msc[idx]->Z23*stepmin;
262   if(currentKinEnergy < tlow) { x *= 0.5*(1.+c    262   if(currentKinEnergy < tlow) { x *= 0.5*(1.+currentKinEnergy/tlow); }
263   return std::max(x, tlimitminfix);               263   return std::max(x, tlimitminfix);
264 }                                                 264 }
265                                                   265 
266 //....oooOO0OOooo........oooOO0OOooo........oo    266 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
267                                                   267 
268 inline void G4UrbanMscModel::SetDisplacementAl    268 inline void G4UrbanMscModel::SetDisplacementAlgorithm96(const G4bool val)
269 {                                                 269 {
270   dispAlg96 = val;                                270   dispAlg96 = val;
271 }                                                 271 }
272                                                   272 
273 //....oooOO0OOooo........oooOO0OOooo........oo    273 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
274                                                   274 
275 inline void G4UrbanMscModel::SetPositronCorrec    275 inline void G4UrbanMscModel::SetPositronCorrection(const G4bool val)
276 {                                                 276 {
277   fPosiCorrection = val;                          277   fPosiCorrection = val;
278 }                                                 278 }
279                                                   279 
280 //....oooOO0OOooo........oooOO0OOooo........oo    280 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
281                                                   281 
282 #endif                                            282 #endif
283                                                   283 
284                                                   284