Geant4 Cross Reference

Cross-Referencing   Geant4
Geant4/processes/electromagnetic/utils/include/G4EnergyLossTables.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/utils/include/G4EnergyLossTables.hh (Version 11.3.0) and /processes/electromagnetic/utils/include/G4EnergyLossTables.hh (Version 11.0.p3,)


** Warning: Cannot open xref database.

  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 // -------------------------------------------    
 30                                                   
 31 #ifndef included_G4EnergyLossTables               
 32 #define included_G4EnergyLossTables               
 33                                                   
 34 #include <map>                                    
 35 #include "globals.hh"                             
 36                                                   
 37 #include "G4PhysicsTable.hh"                      
 38 #include "G4ParticleDefinition.hh"                
 39 #include "G4Material.hh"                          
 40 #include "G4ios.hh"                               
 41                                                   
 42 //--------------------------------------------    
 43 // A utility class, containing the energy loss    
 44 // for each particle                              
 45 //                                                
 46 // Energy loss processes have to register thei    
 47 // class. The responsibility of creating and d    
 48 // remains with the energy loss classes.          
 49 // -------------------------------------------    
 50 //                                                
 51 // P. Urban, 06/04/1998                           
 52 // L. Urban, 27/05/1988 , modifications + new     
 53 // L.Urban , 13/10/98 , revision                  
 54 // L.Urban,  26/10/98 , revision, Interpolate     
 55 // L.Urban , 08/02/99,  cache mechanism           
 56 // L.Urban , 12/04/99 , bug fixed                 
 57 // don't use the helper class.                    
 58 // It can't be hidden for Rogue Wave uses it.     
 59 // 10.11.99: moved from RWT hash dictionary to    
 60 // 26.10.01: all static functions movev from .    
 61 // 15.01.03 Add interfaces required for "cut p    
 62 // 12.03.03 Add warnings to obsolete interface    
 63 // 12.04.03 move exception to new method (V.Iv    
 64 //                                                
 65 // -------------------------------------------    
 66                                                   
 67 //....oooOO0OOooo........oooOO0OOooo........oo    
 68                                                   
 69 class G4EnergyLossTablesHelper {                  
 70                                                   
 71 friend class G4EnergyLossTables;                  
 72   // the only instances are within the class G    
 73                                                   
 74 public:                                           
 75   G4EnergyLossTablesHelper();                     
 76                                                   
 77 private:                                          
 78   G4EnergyLossTablesHelper(const G4PhysicsTabl    
 79          const G4PhysicsTable* aRangeTable,       
 80                            const G4PhysicsTabl    
 81                            const G4PhysicsTabl    
 82                            const G4PhysicsTabl    
 83          G4double aLowestKineticEnergy,           
 84          G4double aHighestKineticEnergy,          
 85          G4double aMassRatio,                     
 86                            G4int aNumberOfBins    
 87   // data to be stored in the dictionary          
 88   const G4PhysicsTable* theDEDXTable;             
 89   const G4PhysicsTable* theRangeTable;            
 90   const G4PhysicsTable* theInverseRangeTable;     
 91   const G4PhysicsTable* theLabTimeTable;          
 92   const G4PhysicsTable* theProperTimeTable;       
 93   G4double theLowestKineticEnergy;                
 94   G4double theHighestKineticEnergy;               
 95   G4double theMassRatio;                          
 96   G4int theNumberOfBins;                          
 97 };                                                
 98                                                   
 99 //....oooOO0OOooo........oooOO0OOooo........oo    
100                                                   
101 class G4MaterialCutsCouple;                       
102                                                   
103 class G4EnergyLossTables {                        
104                                                   
105 public:                                           
106                                                   
107   // get the table for a given particle           
108   // (0 if the table was not found)               
109   static const G4PhysicsTable* GetDEDXTable(      
110     const G4ParticleDefinition* p);               
111   static const G4PhysicsTable* GetRangeTable(     
112     const G4ParticleDefinition* p);               
113   static const G4PhysicsTable* GetInverseRange    
114     const G4ParticleDefinition* p);               
115   static const G4PhysicsTable* GetLabTimeTable    
116     const G4ParticleDefinition* p);               
117   static const G4PhysicsTable* GetProperTimeTa    
118     const G4ParticleDefinition* p);               
119                                                   
120   // get the DEDX or the range for a given par    
121   static G4double GetDEDX(                        
122     const G4ParticleDefinition *aParticle,        
123     G4double KineticEnergy,                       
124     const G4Material *aMaterial);                 
125   static G4double GetRange(                       
126     const G4ParticleDefinition *aParticle,        
127     G4double KineticEnergy,                       
128     const G4Material *aMaterial);                 
129   static G4double GetLabTime(                     
130     const G4ParticleDefinition *aParticle,        
131     G4double KineticEnergy,                       
132     const G4Material *aMaterial);                 
133   static G4double GetDeltaLabTime(                
134     const G4ParticleDefinition *aParticle,        
135     G4double KineticEnergyStart,                  
136     G4double KineticEnergyEnd,                    
137     const G4Material *aMaterial);                 
138   static G4double GetProperTime(                  
139     const G4ParticleDefinition *aParticle,        
140     G4double KineticEnergy,                       
141     const G4Material *aMaterial);                 
142   static G4double GetDeltaProperTime(             
143     const G4ParticleDefinition *aParticle,        
144     G4double KineticEnergyStart,                  
145     G4double KineticEnergyEnd,                    
146     const G4Material *aMaterial);                 
147                                                   
148   static G4double GetPreciseDEDX(                 
149     const G4ParticleDefinition *aParticle,        
150     G4double KineticEnergy,                       
151     const G4Material *aMaterial);                 
152   static G4double GetPreciseRangeFromEnergy(      
153     const G4ParticleDefinition *aParticle,        
154     G4double KineticEnergy,                       
155     const G4Material *aMaterial);                 
156   static G4double GetPreciseEnergyFromRange(      
157     const G4ParticleDefinition *aParticle,        
158     G4double range,                               
159     const G4Material *aMaterial);                 
160                                                   
161   // get the DEDX or the range for a given par    
162   static G4double GetDEDX(                        
163     const G4ParticleDefinition *aParticle,        
164     G4double KineticEnergy,                       
165     const G4MaterialCutsCouple *couple,           
166     G4bool check = true);                         
167   static G4double GetRange(                       
168     const G4ParticleDefinition *aParticle,        
169     G4double KineticEnergy,                       
170     const G4MaterialCutsCouple *couple,           
171     G4bool check = true);                         
172                                                   
173   static G4double GetPreciseDEDX(                 
174     const G4ParticleDefinition *aParticle,        
175     G4double KineticEnergy,                       
176     const G4MaterialCutsCouple *couple);          
177   static G4double GetPreciseRangeFromEnergy(      
178     const G4ParticleDefinition *aParticle,        
179     G4double KineticEnergy,                       
180     const G4MaterialCutsCouple *couple);          
181   static G4double GetPreciseEnergyFromRange(      
182     const G4ParticleDefinition *aParticle,        
183     G4double range,                               
184     const G4MaterialCutsCouple *couple,           
185     G4bool check = true);                         
186                                                   
187   // to be called only by energy loss processe    
188   static void Register(                           
189     const G4ParticleDefinition* p,                
190     const G4PhysicsTable* tDEDX,                  
191     const G4PhysicsTable* tRange,                 
192     const G4PhysicsTable* tInverseRange,          
193     const G4PhysicsTable* tLabTime,               
194     const G4PhysicsTable* tProperTime,            
195     G4double lowestKineticEnergy,                 
196     G4double highestKineticEnergy,                
197     G4double massRatio,                           
198     G4int NumberOfBins);                          
199                                                   
200 public:                                           
201   typedef const G4ParticleDefinition* K;          
202                                                   
203 private:                                          
204                                                   
205   static void CPRWarning();                       
206   static void ParticleHaveNoLoss(const G4Parti    
207                                  const G4Strin    
208                                                   
209   /*                                              
210   typedef std::map<K,G4EnergyLossTablesHelper,    
211   static G4ThreadLocal helper_map *dict;          
212                                                   
213   static G4EnergyLossTablesHelper GetTables(co    
214                                                   
215   static G4ThreadLocal G4EnergyLossTablesHelpe    
216   static G4ThreadLocal G4EnergyLossTablesHelpe    
217   static G4ThreadLocal G4ParticleDefinition* l    
218   static G4ThreadLocal G4double QQPositron ;      
219   static G4ThreadLocal G4double Chargesquare ;    
220   static G4ThreadLocal G4int oldIndex ;           
221   static G4ThreadLocal G4double rmin,rmax,Thig    
222   static G4ThreadLocal G4int  let_counter;        
223   static G4ThreadLocal G4int  let_max_num_warn    
224   static G4ThreadLocal G4bool first_loss;         
225   */                                              
226   using helper_map = std::map<K, G4EnergyLossT    
227   static helper_map *dict;                        
228                                                   
229   static G4EnergyLossTablesHelper GetTables(co    
230                                                   
231   static G4EnergyLossTablesHelper *t;             
232   static G4EnergyLossTablesHelper *null_loss;     
233   static G4ParticleDefinition* lastParticle ;     
234   static G4double QQPositron ;                    
235   static G4double Chargesquare ;                  
236   static G4int oldIndex ;                         
237   static G4double rmin,rmax,Thigh ;               
238   static G4int  let_counter;                      
239   static G4int  let_max_num_warnings;             
240   static G4bool first_loss;                       
241                                                   
242 };                                                
243                                                   
244 //....oooOO0OOooo........oooOO0OOooo........oo    
245                                                   
246 #endif                                            
247