Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/include/G4WentzelOKandVIxSection.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/G4WentzelOKandVIxSection.hh (Version 11.3.0) and /processes/electromagnetic/standard/include/G4WentzelOKandVIxSection.hh (Version 9.4)


  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 // $Id: G4WentzelOKandVIxSection.hh,v 1.6 2010/06/25 09:41:38 gunter Exp $
                                                   >>  27 // GEANT4 tag $Name: geant4-09-04-beta-01 $
 26 //                                                 28 //
 27 // -------------------------------------------     29 // -------------------------------------------------------------------
 28 //                                                 30 //
 29 //                                                 31 //
 30 // GEANT4 Class header file                        32 // GEANT4 Class header file
 31 //                                                 33 //
 32 //                                                 34 //
 33 // File name:     G4WentzelOKandVIxSection         35 // File name:     G4WentzelOKandVIxSection
 34 //                                                 36 //
 35 // Authors:       V.Ivanchenko and O.Kadri         37 // Authors:       V.Ivanchenko and O.Kadri 
 36 //                                                 38 //
 37 // Creation date: 21.05.2010                       39 // Creation date: 21.05.2010 
 38 //                                                 40 //
 39 // Modifications:                                  41 // Modifications:
 40 //                                                 42 //
 41 //                                                 43 //
 42 // Class Description:                              44 // Class Description:
 43 //                                                 45 //
 44 // Implementation of the computation of total      46 // Implementation of the computation of total and transport cross sections,
 45 // sample scattering angle for the single scat     47 // sample scattering angle for the single scattering case.
 46 // to be used by single and multiple scatterin     48 // to be used by single and multiple scattering models. References:
 47 // 1) G.Wentzel, Z. Phys. 40 (1927) 590.           49 // 1) G.Wentzel, Z. Phys. 40 (1927) 590.
 48 // 2) J.M. Fernandez-Varea et al., NIM B73 (19     50 // 2) J.M. Fernandez-Varea et al., NIM B73 (1993) 447.
 49 //                                                 51 //
 50 // -------------------------------------------     52 // -------------------------------------------------------------------
 51 //                                                 53 //
 52                                                    54 
 53 #ifndef G4WentzelOKandVIxSection_h                 55 #ifndef G4WentzelOKandVIxSection_h
 54 #define G4WentzelOKandVIxSection_h 1               56 #define G4WentzelOKandVIxSection_h 1
 55                                                    57 
 56 //....oooOO0OOooo........oooOO0OOooo........oo     58 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 57                                                    59 
 58 #include "globals.hh"                              60 #include "globals.hh"
 59 #include "G4Material.hh"                           61 #include "G4Material.hh"
 60 #include "G4Element.hh"                            62 #include "G4Element.hh"
 61 #include "G4ElementVector.hh"                      63 #include "G4ElementVector.hh"
 62 #include "G4NistManager.hh"                        64 #include "G4NistManager.hh"
 63 #include "G4NuclearFormfactorType.hh"          << 
 64 #include "G4ThreeVector.hh"                        65 #include "G4ThreeVector.hh"
 65 #include "G4Pow.hh"                                66 #include "G4Pow.hh"
 66                                                    67 
 67 class G4ParticleDefinition;                        68 class G4ParticleDefinition;
 68 class G4ScreeningMottCrossSection;             << 
 69                                                    69 
 70 //....oooOO0OOooo........oooOO0OOooo........oo     70 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
 71                                                    71 
 72 class G4WentzelOKandVIxSection                     72 class G4WentzelOKandVIxSection 
 73 {                                                  73 {
 74                                                    74 
 75 public:                                            75 public:
 76                                                    76 
 77   explicit G4WentzelOKandVIxSection(G4bool com <<  77   G4WentzelOKandVIxSection();
 78                                                    78 
 79   virtual ~G4WentzelOKandVIxSection();             79   virtual ~G4WentzelOKandVIxSection();
 80                                                    80 
 81   void Initialise(const G4ParticleDefinition*,     81   void Initialise(const G4ParticleDefinition*, G4double CosThetaLim);
 82                                                    82 
 83   void SetupParticle(const G4ParticleDefinitio     83   void SetupParticle(const G4ParticleDefinition*);
 84                                                    84 
 85   // return cos(ThetaMax) for msc and cos(thet     85   // return cos(ThetaMax) for msc and cos(thetaMin) for single scattering
 86   // cut = DBL_MAX means no scattering off ele     86   // cut = DBL_MAX means no scattering off electrons 
 87   virtual G4double SetupKinematic(G4double kin <<  87   G4double SetupTarget(G4int Z, G4double cut = DBL_MAX);
 88   G4double SetupTarget(G4int Z, G4double cut); << 
 89                                                    88 
 90   G4double ComputeTransportCrossSectionPerAtom     89   G4double ComputeTransportCrossSectionPerAtom(G4double CosThetaMax);
 91                                                    90  
 92   G4ThreeVector& SampleSingleScattering(G4doub <<  91   G4ThreeVector SampleSingleScattering(G4double CosThetaMin,
 93           G4double CosThetaMax,                <<  92                G4double CosThetaMax,
 94           G4double elecRatio);                 <<  93                G4double elecRatio = 0.0);
 95                                                << 
 96   G4double ComputeSecondTransportMoment(G4doub << 
 97                                                    94 
 98   inline G4double ComputeNuclearCrossSection(G     95   inline G4double ComputeNuclearCrossSection(G4double CosThetaMin,
 99                G4double CosThetaMax);              96                G4double CosThetaMax);
100                                                    97  
101   inline G4double ComputeElectronCrossSection(     98   inline G4double ComputeElectronCrossSection(G4double CosThetaMin,
102                 G4double CosThetaMax);             99                 G4double CosThetaMax);
103                                                << 100  
                                                   >> 101   inline G4double SetupKinematic(G4double kinEnergy, const G4Material* mat);
                                                   >> 102   
104   inline void SetTargetMass(G4double value);      103   inline void SetTargetMass(G4double value);
105                                                   104 
                                                   >> 105   inline void SetRelativisticMass(G4double value);
                                                   >> 106 
106   inline G4double GetMomentumSquare() const;      107   inline G4double GetMomentumSquare() const;
107                                                   108 
108   inline G4double GetCosThetaNuc() const;         109   inline G4double GetCosThetaNuc() const;
109                                                   110 
110   inline G4double GetCosThetaElec() const;        111   inline G4double GetCosThetaElec() const;
111                                                   112 
112   //  hide assignment operator                 << 113 private:
113   G4WentzelOKandVIxSection & operator=         << 
114   (const G4WentzelOKandVIxSection &right) = de << 
115   G4WentzelOKandVIxSection(const  G4WentzelOKa << 
116                                                << 
117 protected:                                     << 
118                                                   114 
119   void ComputeMaxElectronScattering(G4double c    115   void ComputeMaxElectronScattering(G4double cut);
120                                                   116 
121   void InitialiseA();                          << 117   //  hide assignment operator
122                                                << 118   G4WentzelOKandVIxSection & operator=(const  G4WentzelOKandVIxSection &right);
123   inline G4double FlatFormfactor(G4double x);  << 119   G4WentzelOKandVIxSection(const  G4WentzelOKandVIxSection&);
124                                                   120 
125   const G4ParticleDefinition* theProton;          121   const G4ParticleDefinition* theProton;
126   const G4ParticleDefinition* theElectron;        122   const G4ParticleDefinition* theElectron;
127   const G4ParticleDefinition* thePositron;        123   const G4ParticleDefinition* thePositron;
128   const G4ParticleDefinition* particle = nullp << 124   const G4Material* currentMaterial;
129   const G4Material* currentMaterial = nullptr; << 
130                                                   125 
131   G4NistManager* fNistManager;                 << 126   G4NistManager*  fNistManager;
132   G4Pow*         fG4pow;                       << 127   G4Pow*          fG4pow;
133                                                   128 
134   G4ScreeningMottCrossSection* fMottXSection = << 129   G4double numlimit;
135                                                   130 
136   G4ThreeVector temp;                          << 131   G4double elecXSRatio;
                                                   >> 132 
                                                   >> 133   // integer parameters
                                                   >> 134   G4int    nwarnings;
                                                   >> 135   G4int    nwarnlimit;
137                                                   136 
138   // single scattering parameters                 137   // single scattering parameters
139   G4double coeff;                                 138   G4double coeff;
140   G4double cosTetMaxElec = 1.0;                << 139   G4double cosTetMaxElec;
141   G4double cosTetMaxNuc = 1.0;                 << 140   G4double cosTetMaxNuc;
142                                                << 141   G4double cosThetaMax;
143   // for the combined mode it is cos(thetaMax) << 142   G4double alpha2;
144   // for single scattering it is cos(thetaMin) << 143 
145   G4double cosThetaMax = 1.0;                  << 144   // projectile
146                                                << 145   const G4ParticleDefinition* particle;
147   G4double chargeSquare = 0.0;                 << 146 
148   G4double charge3 = 0.0;                      << 147   G4double chargeSquare;
149   G4double spin = 0.0;                         << 148   G4double charge3;
150   G4double mass = 0.0;                         << 149   G4double spin;
151   G4double tkin = 0.0;                         << 150   G4double mass;
152   G4double mom2 = 0.0;                         << 151   G4double tkin;
153   G4double momCM2 = 0.0;                       << 152   G4double mom2;
154   G4double invbeta2 = 1.0;                     << 153   G4double invbeta2;
155   G4double kinFactor = 1.0;                    << 154   G4double kinFactor;
156   G4double etag = DBL_MAX;                     << 155   G4double etag;
157   G4double ecut = DBL_MAX;                     << 156   G4double ecut;
                                                   >> 157   G4double lowEnergyLimit;
158                                                   158 
159   // target                                       159   // target
                                                   >> 160   G4int    targetZ;
160   G4double targetMass;                            161   G4double targetMass;
161   G4double screenZ = 0.0;                      << 162   G4double screenZ;
162   G4double formfactA = 0.0;                    << 163   G4double formfactA;
163   G4double factorA2 = 0.0;                     << 164   G4double factorA2;
164   G4double factB = 0.0;                        << 165   G4double factB;
165   G4double factD = 0.0;                        << 166   G4double factD;
166   G4double fMottFactor = 1.0;                  << 
167   G4double gam0pcmp = 1.0;                     << 
168   G4double pcmp2 = 1.0;                        << 
169                                                << 
170   // integer parameters                        << 
171   G4int targetZ = 0;                           << 
172   G4int nwarnings = 0;                         << 
173                                                   167 
174   G4NuclearFormfactorType fNucFormfactor = fEx << 
175                                                << 
176   G4bool isCombined;                           << 
177                                                << 
178   static G4double ScreenRSquareElec[100];      << 
179   static G4double ScreenRSquare[100];             168   static G4double ScreenRSquare[100];
180   static G4double FormFactor[100];                169   static G4double FormFactor[100];
181 };                                                170 };
182                                                   171 
                                                   >> 172 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo....
                                                   >> 173 
                                                   >> 174 inline G4double 
                                                   >> 175 G4WentzelOKandVIxSection::SetupKinematic(G4double ekin, const G4Material* mat)
                                                   >> 176 {
                                                   >> 177   if(ekin != tkin || mat != currentMaterial) { 
                                                   >> 178     currentMaterial = mat;
                                                   >> 179     tkin  = ekin;
                                                   >> 180     mom2  = tkin*(tkin + 2.0*mass);
                                                   >> 181     invbeta2 = 1.0 +  mass*mass/mom2;
                                                   >> 182     factB = spin/invbeta2;
                                                   >> 183     cosTetMaxNuc = cosThetaMax;
                                                   >> 184     if(std::fabs(cosThetaMax) < 1.0) {
                                                   >> 185       cosTetMaxNuc = 
                                                   >> 186   std::max(cosThetaMax,1.-factorA2*mat->GetIonisation()->GetInvA23()/mom2);
                                                   >> 187     }
                                                   >> 188   } 
                                                   >> 189   return cosTetMaxNuc;
                                                   >> 190 }
                                                   >> 191 
183 //....oooOO0OOooo........oooOO0OOooo........oo    192 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
184                                                   193 
185 inline void G4WentzelOKandVIxSection::SetTarge    194 inline void G4WentzelOKandVIxSection::SetTargetMass(G4double value)
186 {                                                 195 {
187   targetMass = value;                             196   targetMass = value;
188   factD = std::sqrt(mom2)/value;                  197   factD = std::sqrt(mom2)/value;
189 }                                                 198 }
190                                                   199 
191 //....oooOO0OOooo........oooOO0OOooo........oo    200 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
192                                                   201 
                                                   >> 202 inline void G4WentzelOKandVIxSection::SetRelativisticMass(G4double value)
                                                   >> 203 {
                                                   >> 204   mass = value;
                                                   >> 205 }
                                                   >> 206 
                                                   >> 207 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
                                                   >> 208 
193 inline G4double G4WentzelOKandVIxSection::GetM    209 inline G4double G4WentzelOKandVIxSection::GetMomentumSquare() const
194 {                                                 210 {
195   return mom2;                                    211   return mom2;
196 }                                                 212 }
197                                                   213 
198 //....oooOO0OOooo........oooOO0OOooo........oo    214 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
199                                                   215 
200 inline G4double G4WentzelOKandVIxSection::GetC    216 inline G4double G4WentzelOKandVIxSection::GetCosThetaNuc() const
201 {                                                 217 {
202   return cosTetMaxNuc;                            218   return cosTetMaxNuc;
203 }                                                 219 }
204                                                   220 
205 //....oooOO0OOooo........oooOO0OOooo........oo    221 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
206                                                   222 
207 inline G4double G4WentzelOKandVIxSection::GetC    223 inline G4double G4WentzelOKandVIxSection::GetCosThetaElec() const
208 {                                                 224 {
209   return cosTetMaxElec;                           225   return cosTetMaxElec;
210 }                                                 226 }
211                                                   227 
212 //....oooOO0OOooo........oooOO0OOooo........oo    228 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
213                                                   229 
214 inline G4double                                   230 inline G4double 
215 G4WentzelOKandVIxSection::ComputeNuclearCrossS    231 G4WentzelOKandVIxSection::ComputeNuclearCrossSection(G4double cosTMin,
216                  G4double cosTMax)                232                  G4double cosTMax)
217 {                                                 233 {
218   return targetZ*kinFactor*fMottFactor*(cosTMi << 234   G4double xsec = 0.0;
219     ((1.0 - cosTMin + screenZ)*(1.0 - cosTMax  << 235   if(cosTMax < cosTMin) {
                                                   >> 236     xsec = targetZ*kinFactor*(cosTMin - cosTMax)/
                                                   >> 237       ((1.0 - cosTMin + screenZ)*(1.0 - cosTMax + screenZ));
                                                   >> 238   }
                                                   >> 239   return xsec;
220 }                                                 240 }
221                                                   241 
222 //....oooOO0OOooo........oooOO0OOooo........oo    242 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
223                                                   243 
224 inline G4double                                   244 inline G4double 
225 G4WentzelOKandVIxSection::ComputeElectronCross    245 G4WentzelOKandVIxSection::ComputeElectronCrossSection(G4double cosTMin,
226                   G4double cosTMax)               246                   G4double cosTMax)
227 {                                                 247 {
                                                   >> 248   G4double xsec = 0.0;
228   G4double cost1 = std::max(cosTMin,cosTetMaxE    249   G4double cost1 = std::max(cosTMin,cosTetMaxElec);
229   G4double cost2 = std::max(cosTMax,cosTetMaxE    250   G4double cost2 = std::max(cosTMax,cosTetMaxElec);
230   return (cost1 <= cost2) ? 0.0 : kinFactor*fM << 251   if(cost1 > cost2) {
231     ((1.0 - cost1 + screenZ)*(1.0 - cost2 + sc << 252     xsec = kinFactor*(cost1 - cost2)/((1.0 - cost1 + screenZ)*(1.0 - cost2 + screenZ));
                                                   >> 253   }
                                                   >> 254   return xsec;
232 }                                                 255 }
233                                                << 
234 inline G4double G4WentzelOKandVIxSection::Flat << 
235 {                                              << 
236   return 3.0*(std::sin(x) - x*std::cos(x))/(x* << 
237 }                                              << 
238                                                << 
239 //....oooOO0OOooo........oooOO0OOooo........oo    256 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
240                                                   257 
241 #endif                                            258 #endif
242                                                   259 
243                                                   260