Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/lowenergy/src/G4hICRU49Nuclear.cc

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/lowenergy/src/G4hICRU49Nuclear.cc (Version 11.3.0) and /processes/electromagnetic/lowenergy/src/G4hICRU49Nuclear.cc (Version 1.1)


  1 //                                                  1 
  2 // *******************************************    
  3 // * License and Disclaimer                       
  4 // *                                              
  5 // * The  Geant4 software  is  copyright of th    
  6 // * the Geant4 Collaboration.  It is provided    
  7 // * conditions of the Geant4 Software License    
  8 // * LICENSE and available at  http://cern.ch/    
  9 // * include a list of copyright holders.         
 10 // *                                              
 11 // * Neither the authors of this software syst    
 12 // * institutes,nor the agencies providing fin    
 13 // * work  make  any representation or  warran    
 14 // * regarding  this  software system or assum    
 15 // * use.  Please see the license in the file     
 16 // * for the full disclaimer and the limitatio    
 17 // *                                              
 18 // * This  code  implementation is the result     
 19 // * technical work of the GEANT4 collaboratio    
 20 // * By using,  copying,  modifying or  distri    
 21 // * any work based  on the software)  you  ag    
 22 // * use  in  resulting  scientific  publicati    
 23 // * acceptance of all terms of the Geant4 Sof    
 24 // *******************************************    
 25 //                                                
 26 //                                                
 27 // -------------------------------------------    
 28 //                                                
 29 // GEANT4 Class file                              
 30 //                                                
 31 //                                                
 32 // File name:     G4hICRU49Nuclear                
 33 //                                                
 34 // Author:        V.Ivanchenko (Vladimir.Ivanc    
 35 //                                                
 36 // Creation date: 20 July 2000                    
 37 //                                                
 38 // Modifications:                                 
 39 // 20/07/2000  V.Ivanchenko First implementati    
 40 // 10/05/2001  V.Ivanchenko Clean up againist     
 41 //                                                
 42 // Class Description:                             
 43 //                                                
 44 // Nuclear stopping power parametrised accordi    
 45 // ICRU Report N49, 1993. Moliere model.          
 46 // G.Moliere "Theorie der Streuung schneller g    
 47 // Einzelstreuungam abbgeschirmten Coulomb-Fel    
 48 // 133 (1947).                                    
 49 //                                                
 50 // Class Description: End                         
 51 //                                                
 52 // -------------------------------------------    
 53 //                                                
 54 //....oooOO0OOooo........oooOO0OOooo........oo    
 55                                                   
 56 #include "G4hICRU49Nuclear.hh"                    
 57                                                   
 58 #include "globals.hh"                             
 59 #include "G4SystemOfUnits.hh"                     
 60 #include "G4UnitsTable.hh"                        
 61 #include "Randomize.hh"                           
 62                                                   
 63 //....oooOO0OOooo........oooOO0OOooo........oo    
 64                                                   
 65 G4hICRU49Nuclear::G4hICRU49Nuclear():G4VhNucle    
 66 {;}                                               
 67                                                   
 68 //....oooOO0OOooo........oooOO0OOooo........oo    
 69                                                   
 70 G4hICRU49Nuclear::~G4hICRU49Nuclear()             
 71 {;}                                               
 72                                                   
 73 //....oooOO0OOooo........oooOO0OOooo........oo    
 74                                                   
 75 G4double G4hICRU49Nuclear::NuclearStoppingPowe    
 76                                 G4double z1, G    
 77                                 G4double m1, G    
 78 {                                                 
 79   G4double energy = kineticEnergy/keV ;  // en    
 80   G4double ionloss = 0.0;                         
 81                                                   
 82   G4double rm = (m1 + m2Local) * ( std::pow(z1    
 83                                                   
 84   G4double er = 32.536 * m2Local * energy / (     
 85                                                   
 86   static const G4double a[104][2] = {             
 87   { 1.0E+8, 5.831E-8},                            
 88   { 8.0E+7, 7.288E-8},                            
 89   { 6.0E+7, 9.719E-8},                            
 90   { 5.0E+7, 1.166E-7},                            
 91   { 4.0E+7, 1.457E-7},                            
 92   { 3.0E+7, 1.942E-7},                            
 93   { 2.0E+7, 2.916E-7},                            
 94   { 1.5E+7, 3.887E-7},                            
 95                                                   
 96   { 1.0E+7, 5.833E-7},                            
 97   { 8.0E+6, 7.287E-7},                            
 98   { 6.0E+6, 9.712E-7},                            
 99   { 5.0E+6, 1.166E-6},                            
100   { 4.0E+6, 1.457E-6},                            
101   { 3.0E+6, 1.941E-6},                            
102   { 2.0E+6, 2.911E-6},                            
103   { 1.5E+6, 3.878E-6},                            
104                                                   
105   { 1.0E+6, 5.810E-6},                            
106   { 8.0E+5, 7.262E-6},                            
107   { 6.0E+5, 9.663E-6},                            
108   { 5.0E+5, 1.157E-5},                            
109   { 4.0E+5, 1.442E-5},                            
110   { 3.0E+5, 1.913E-5},                            
111   { 2.0E+5, 2.845E-5},                            
112   { 1.5E+5, 3.762E-5},                            
113                                                   
114   { 1.0E+5, 5.554E-5},                            
115   { 8.0E+4, 6.866E-5},                            
116   { 6.0E+4, 9.020E-5},                            
117   { 5.0E+4, 1.070E-4},                            
118   { 4.0E+4, 1.319E-4},                            
119   { 3.0E+4, 1.722E-4},                            
120   { 2.0E+4, 2.499E-4},                            
121   { 1.5E+4, 3.248E-4},                            
122                                                   
123   { 1.0E+4, 4.688E-4},                            
124   { 8.0E+3, 5.729E-4},                            
125   { 6.0E+3, 7.411E-4},                            
126   { 5.0E+3, 8.718E-4},                            
127   { 4.0E+3, 1.063E-3},                            
128   { 3.0E+3, 1.370E-3},                            
129   { 2.0E+3, 1.955E-3},                            
130   { 1.5E+3, 2.511E-3},                            
131                                                   
132   { 1.0E+3, 3.563E-3},                            
133   { 8.0E+2, 4.314E-3},                            
134   { 6.0E+2, 5.511E-3},                            
135   { 5.0E+2, 6.430E-3},                            
136   { 4.0E+2, 7.756E-3},                            
137   { 3.0E+2, 9.855E-3},                            
138   { 2.0E+2, 1.375E-2},                            
139   { 1.5E+2, 1.736E-2},                            
140                                                   
141   { 1.0E+2, 2.395E-2},                            
142   { 8.0E+1, 2.850E-2},                            
143   { 6.0E+1, 3.552E-2},                            
144   { 5.0E+1, 4.073E-2},                            
145   { 4.0E+1, 4.802E-2},                            
146   { 3.0E+1, 5.904E-2},                            
147   { 1.5E+1, 9.426E-2},                            
148                                                   
149   { 1.0E+1, 1.210E-1},                            
150   { 8.0E+0, 1.377E-1},                            
151   { 6.0E+0, 1.611E-1},                            
152   { 5.0E+0, 1.768E-1},                            
153   { 4.0E+0, 1.968E-1},                            
154   { 3.0E+0, 2.235E-1},                            
155   { 2.0E+0, 2.613E-1},                            
156   { 1.5E+0, 2.871E-1},                            
157                                                   
158   { 1.0E+0, 3.199E-1},                            
159   { 8.0E-1, 3.354E-1},                            
160   { 6.0E-1, 3.523E-1},                            
161   { 5.0E-1, 3.609E-1},                            
162   { 4.0E-1, 3.693E-1},                            
163   { 3.0E-1, 3.766E-1},                            
164   { 2.0E-1, 3.803E-1},                            
165   { 1.5E-1, 3.788E-1},                            
166                                                   
167   { 1.0E-1, 3.711E-1},                            
168   { 8.0E-2, 3.644E-1},                            
169   { 6.0E-2, 3.530E-1},                            
170   { 5.0E-2, 3.444E-1},                            
171   { 4.0E-2, 3.323E-1},                            
172   { 3.0E-2, 3.144E-1},                            
173   { 2.0E-2, 2.854E-1},                            
174   { 1.5E-2, 2.629E-1},                            
175                                                   
176   { 1.0E-2, 2.298E-1},                            
177   { 8.0E-3, 2.115E-1},                            
178   { 6.0E-3, 1.883E-1},                            
179   { 5.0E-3, 1.741E-1},                            
180   { 4.0E-3, 1.574E-1},                            
181   { 3.0E-3, 1.372E-1},                            
182   { 2.0E-3, 1.116E-1},                            
183   { 1.5E-3, 9.559E-2},                            
184                                                   
185   { 1.0E-3, 7.601E-2},                            
186   { 8.0E-4, 6.668E-2},                            
187   { 6.0E-4, 5.605E-2},                            
188   { 5.0E-4, 5.008E-2},                            
189   { 4.0E-4, 4.352E-2},                            
190   { 3.0E-4, 3.617E-2},                            
191   { 2.0E-4, 2.768E-2},                            
192   { 1.5E-4, 2.279E-2},                            
193                                                   
194   { 1.0E-4, 1.723E-2},                            
195   { 8.0E-5, 1.473E-2},                            
196   { 6.0E-5, 1.200E-2},                            
197   { 5.0E-5, 1.052E-2},                            
198   { 4.0E-5, 8.950E-3},                            
199   { 3.0E-5, 7.246E-3},                            
200   { 2.0E-5, 5.358E-3},                            
201   { 1.5E-5, 4.313E-3},                            
202   { 0.0, 3.166E-3}                                
203   };                                              
204                                                   
205   for (G4int i=1; i<104; i++)                     
206     {                                             
207       if (er > a[i][0]) {                         
208   ionloss =                                       
209          (a[i][1]-a[i-1][1])*(er-a[i-1][0])/(a    
210   break;                                          
211       }                                           
212     }                                             
213                                                   
214   // Stragling                                    
215   if(lossFlucFlag) {                              
216     G4double sig = 4.0 * m1 * m2Local / ((m1 +    
217                   (4.0 + 0.197*std::pow(er,-1.    
218                                                   
219     ionloss *= G4RandGauss::shoot(1.0,sig) ;      
220   }                                               
221                                                   
222   ionloss *= 8.462 * z1 * z2 * m1 / rm ; // Re    
223                                                   
224   if ( ionloss < 0.0) ionloss = 0.0 ;             
225                                                   
226   return ionloss;                                 
227 }                                                 
228                                                   
229                                                   
230                                                   
231                                                   
232                                                   
233