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 7.0.p1)


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