Geant4 Cross Reference |
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 #ifndef G4EmSaturation_h 28 #define G4EmSaturation_h 1 29 30 // ------------------------------------------- 31 // 32 // GEANT4 Class header file 33 // 34 // 35 // File name: G4EmSaturation 36 // 37 // Author: Vladimir Ivanchenko 38 // 39 // Creation date: 18.02.2008 40 // 41 // Modifications: 42 // 43 // 44 // Class Description: 45 // Compution on saturation effect, which red 46 // deposition at the step. Default implement 47 // account Birks effect. Birks coefficients 48 // from G4 database on materials are provide 49 // 50 // This class assumed to be G4ThreadLocal, b 51 // cache value for material 52 // 53 // ------------------------------------------- 54 55 //....oooOO0OOooo........oooOO0OOooo........oo 56 //....oooOO0OOooo........oooOO0OOooo........oo 57 58 #include "globals.hh" 59 #include "G4Step.hh" 60 #include "G4ParticleDefinition.hh" 61 #include <vector> 62 63 //....oooOO0OOooo........oooOO0OOooo........oo 64 65 class G4NistManager; 66 class G4MaterialCutsCouple; 67 class G4Material; 68 69 class G4EmSaturation 70 { 71 public: 72 73 explicit G4EmSaturation(G4int verb); 74 virtual ~G4EmSaturation(); 75 76 // this method may be overwritten in the der 77 // which implements alternative algorithm of 78 virtual G4double VisibleEnergyDeposition(con 79 const G4MaterialCutsCouple*, 80 G4double length, 81 G4double edepTotal, 82 G4double edepNIEL = 0.0) const; 83 84 // activate default model 85 void InitialiseG4Saturation(); 86 87 // find and Birks coefficient 88 G4double FindG4BirksCoefficient(const G4Mate 89 90 // dump coeffitients used in run time 91 void DumpBirksCoefficients(); 92 93 // dump G4 list 94 void DumpG4BirksCoefficients(); 95 96 // this method should not be overwritten 97 inline G4double VisibleEnergyDepositionAtASt 98 99 inline void SetVerbose(G4int); 100 101 // hide assignment operator 102 G4EmSaturation & operator=(const G4EmSaturat 103 G4EmSaturation(const G4EmSaturation&) = dele 104 105 private: 106 107 void InitialiseBirksCoefficient(const G4Mate 108 109 void InitialiseG4materials(); 110 111 const G4ParticleDefinition* electron; 112 const G4ParticleDefinition* proton; 113 G4NistManager* nist; 114 115 G4int verbose; 116 G4int nG4Birks; 117 G4int nWarnings; 118 119 static size_t nMaterials; 120 121 // list of materials used in run time 122 static std::vector<G4double> massFactors; 123 static std::vector<G4double> effCharges; 124 125 // list of G4 materials 126 static std::vector<G4double> g4MatData; 127 static std::vector<G4String> g4MatNames; 128 }; 129 130 //....oooOO0OOooo........oooOO0OOooo........oo 131 132 inline void G4EmSaturation::SetVerbose(G4int v 133 { 134 verbose = val; 135 } 136 137 //....oooOO0OOooo........oooOO0OOooo........oo 138 139 inline G4double G4EmSaturation::VisibleEnergyD 140 const G4Step* step) const 141 { 142 return VisibleEnergyDeposition(step->GetTrac 143 step->GetTrac 144 step->GetStepLength(), 145 step->GetTota 146 step->GetNonI 147 } 148 149 //....oooOO0OOooo........oooOO0OOooo........oo 150 151 #endif 152 153