Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/standard/src/G4CoulombScattering.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/standard/src/G4CoulombScattering.cc (Version 11.3.0) and /processes/electromagnetic/standard/src/G4CoulombScattering.cc (Version 4.0)


  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:     G4CoulombScattering             
 33 //                                                
 34 // Author:        Vladimir Ivanchenko             
 35 //                                                
 36 // Creation date: 22.08.2004                      
 37 //                                                
 38 // Modifications:                                 
 39 // 01.08.06 V.Ivanchenko add choice between G4    
 40 //          G4CoulombScatteringModel              
 41 //                                                
 42                                                   
 43 //                                                
 44 // -------------------------------------------    
 45 //                                                
 46 //....oooOO0OOooo........oooOO0OOooo........oo    
 47 //....oooOO0OOooo........oooOO0OOooo........oo    
 48                                                   
 49 #include "G4CoulombScattering.hh"                 
 50 #include "G4SystemOfUnits.hh"                     
 51 #include "G4eCoulombScatteringModel.hh"           
 52 #include "G4IonCoulombScatteringModel.hh"         
 53 #include "G4Proton.hh"                            
 54 #include "G4EmParameters.hh"                      
 55                                                   
 56 //....oooOO0OOooo........oooOO0OOooo........oo    
 57                                                   
 58 G4CoulombScattering::G4CoulombScattering(const    
 59   : G4VEmProcess(nam),                            
 60     q2Max(CLHEP::TeV*CLHEP::TeV),                 
 61     isCombined(comb)                              
 62 {                                                 
 63   SetBuildTableFlag(true);                        
 64   SetStartFromNullFlag(false);                    
 65   SetSplineFlag(false);                           
 66   SetCrossSectionType(fEmOnePeak);                
 67   SetSecondaryParticle(G4Proton::Proton());       
 68   SetProcessSubType(fCoulombScattering);          
 69 }                                                 
 70                                                   
 71 //....oooOO0OOooo........oooOO0OOooo........oo    
 72                                                   
 73 G4CoulombScattering::G4CoulombScattering(const    
 74   : G4CoulombScattering(nam, true)                
 75 {}                                                
 76                                                   
 77 //....oooOO0OOooo........oooOO0OOooo........oo    
 78                                                   
 79 G4CoulombScattering::G4CoulombScattering(G4boo    
 80   : G4CoulombScattering("CoulombScat", comb)      
 81 {}                                                
 82                                                   
 83 //....oooOO0OOooo........oooOO0OOooo........oo    
 84                                                   
 85 G4CoulombScattering::~G4CoulombScattering() =     
 86                                                   
 87 //....oooOO0OOooo........oooOO0OOooo........oo    
 88                                                   
 89 G4bool G4CoulombScattering::IsApplicable(const    
 90 {                                                 
 91   return (p.GetPDGCharge() != 0.0);               
 92 }                                                 
 93                                                   
 94 //....oooOO0OOooo........oooOO0OOooo........oo    
 95                                                   
 96 void G4CoulombScattering::InitialiseProcess(co    
 97 {                                                 
 98   // second initialisation not allowed for the    
 99   // this means that polar angle limit change     
100   // after first initialisation                   
101   if(isInitialised) { return; }                   
102                                                   
103   G4EmParameters* param = G4EmParameters::Inst    
104   G4double a = param->FactorForAngleLimit()*CL    
105   q2Max = 0.5*a*a;                                
106   G4double theta = param->MscThetaLimit();        
107                                                   
108   // restricted or non-restricted cross sectio    
109   if(isCombined) {                                
110     if(theta == CLHEP::pi) {                      
111       // for restriced single scattering chang    
112       SetCrossSectionType(fEmIncreasing);         
113       SetStartFromNullFlag(true);                 
114     }                                             
115   } else {                                        
116     SetSplineFlag(true);                          
117     SetCrossSectionType(fEmDecreasing);           
118   }                                               
119   isInitialised = true;                           
120   G4double mass = p->GetPDGMass();                
121   G4String name = p->GetParticleName();           
122                                                   
123   G4bool ion = false;                             
124   if (mass > CLHEP::GeV || p->GetParticleType(    
125     SetBuildTableFlag(false);                     
126     ion = true;                                   
127     if(name != "GenericIon") { SetVerboseLevel    
128   } else {                                        
129     if(name != "e-" && name != "e+" &&            
130        name != "mu+" && name != "mu-" && name     
131        name != "kaon+" && name != "proton" ) {    
132   }                                               
133   /*                                              
134   G4cout << "### G4CoulombScattering::Initiali    
135      << p->GetParticleName()                      
136    << " Emin(MeV)= " << MinKinEnergy()/MeV        
137    << " Emax(TeV)= " << MaxKinEnergy()/TeV        
138    << " nbins= " << LambdaBinning()               
139    << " theta= " << theta                         
140    << " mass(MeV)= " << mass                      
141    << " isCombined=" << isCombined                
142    << " ion=" << ion                              
143    << G4endl;                                     
144   */                                              
145   if(nullptr == EmModel(0)) {                     
146     if(ion) { SetEmModel(new G4IonCoulombScatt    
147     else { SetEmModel(new G4eCoulombScattering    
148   }                                               
149   G4VEmModel* model = EmModel(0);                 
150   G4double emin = std::max(param->MinKinEnergy    
151   G4double emax = std::min(param->MaxKinEnergy    
152   model->SetPolarAngleLimit(theta);               
153   model->SetLowEnergyLimit(emin);                 
154   model->SetHighEnergyLimit(emax);                
155   AddEmModel(1, model);                           
156 }                                                 
157                                                   
158 //....oooOO0OOooo........oooOO0OOooo........oo    
159                                                   
160 G4double G4CoulombScattering::MinPrimaryEnergy    
161                  const G4Material* mat)           
162 {                                                 
163   // Pure Coulomb scattering                      
164   G4double emin = 0.0;                            
165                                                   
166   // Coulomb scattering combined with multiple    
167   G4double theta = G4EmParameters::Instance()-    
168                                                   
169   if(0.0 < theta) {                               
170     G4double p2 = q2Max*mat->GetIonisation()->    
171     G4double mass = part->GetPDGMass();           
172     emin = p2/(std::sqrt(p2 + mass*mass) + mas    
173   }                                               
174                                                   
175   return emin;                                    
176 }                                                 
177                                                   
178 //....oooOO0OOooo........oooOO0OOooo........oo    
179                                                   
180 void G4CoulombScattering::StreamProcessInfo(st    
181 {                                                 
182   G4double tetmin = G4EmParameters::Instance()    
183   outFile << "      ";                            
184   if(tetmin > 179.) { outFile << "ThetaMin(p)"    
185   else              { outFile << tetmin; }        
186   outFile << " < Theta(degree) < 180";            
187                                                   
188   if(q2Max < DBL_MAX) {                           
189     outFile << ", pLimit(GeV^1)= " << std::sqr    
190   }                                               
191   outFile << G4endl;                              
192 }                                                 
193                                                   
194 //....oooOO0OOooo........oooOO0OOooo........oo    
195                                                   
196 void G4CoulombScattering::ProcessDescription(s    
197 {                                                 
198   out <<                                          
199   "  Coulomb scattering. Simulation of elastic    
200   "    events individually. May be used in com    
201   "    scattering, where Coulomb scattering is    
202   "    collisions and multiple scattering for     
203   G4VEmProcess::ProcessDescription(out);          
204 }                                                 
205                                                   
206 //....oooOO0OOooo........oooOO0OOooo........oo    
207