Geant4 Cross Reference |
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: G4EmModelManager.hh,v 1.25 2008/10/13 14:56:56 vnivanch Exp $ >> 27 // GEANT4 tag $Name: geant4-09-02 $ 26 // 28 // 27 // ------------------------------------------- 29 // ------------------------------------------------------------------- 28 // 30 // 29 // GEANT4 Class header file 31 // GEANT4 Class header file 30 // 32 // 31 // File name: G4EmModelManager 33 // File name: G4EmModelManager 32 // 34 // 33 // Author: Vladimir Ivanchenko 35 // Author: Vladimir Ivanchenko 34 // 36 // 35 // Creation date: 07.05.2002 37 // Creation date: 07.05.2002 36 // 38 // 37 // Modifications: 39 // Modifications: 38 // 40 // 39 // 03-12-02 V.Ivanchenko fix a bug in model se 41 // 03-12-02 V.Ivanchenko fix a bug in model selection 40 // 20-01-03 Migrade to cut per region (V.Ivanc 42 // 20-01-03 Migrade to cut per region (V.Ivanchenko) 41 // 27-01-03 Make models region aware (V.Ivanch 43 // 27-01-03 Make models region aware (V.Ivanchenko) 42 // 13-02-03 The set of models is defined for r 44 // 13-02-03 The set of models is defined for region (V.Ivanchenko) 43 // 26-03-03 Add GetDEDXDispersion (V.Ivanchenk 45 // 26-03-03 Add GetDEDXDispersion (V.Ivanchenko) 44 // 13-04-03 Add startFromNull (V.Ivanchenko) 46 // 13-04-03 Add startFromNull (V.Ivanchenko) 45 // 13-05-03 Add calculation of precise range ( 47 // 13-05-03 Add calculation of precise range (V.Ivanchenko) 46 // 21-07-03 Add UpdateEmModel method (V.Ivanch 48 // 21-07-03 Add UpdateEmModel method (V.Ivanchenko) 47 // 03-11-03 Substitute STL vector for G4Region 49 // 03-11-03 Substitute STL vector for G4RegionModels (V.Ivanchenko) 48 // 11-04-05 Remove access to fluctuation model 50 // 11-04-05 Remove access to fluctuation models (V.Ivanchenko) 49 // 10-01-06 PreciseRange -> CSDARange (V.Ivant 51 // 10-01-06 PreciseRange -> CSDARange (V.Ivantchenko) 50 // 20-01-06 Introduce G4EmTableType and reduci 52 // 20-01-06 Introduce G4EmTableType and reducing number of methods (VI) 51 // 13-05-06 Add GetModel by index method (VI) 53 // 13-05-06 Add GetModel by index method (VI) 52 // 15-03-07 Add maxCutInRange (V.Ivanchenko) 54 // 15-03-07 Add maxCutInRange (V.Ivanchenko) 53 // 08-04-08 Simplify Select method for only on << 54 // 03-08-09 Removed unused members and simplif << 55 // model is used (VI) << 56 // 14-07-11 Use pointer to the vector of cuts << 57 // 55 // 58 // Class Description: 56 // Class Description: 59 // 57 // 60 // It is the unified energy loss process it ca 58 // It is the unified energy loss process it calculates the continuous 61 // energy loss for charged particles using a s 59 // energy loss for charged particles using a set of Energy Loss 62 // models valid for different energy regions. 60 // models valid for different energy regions. There are a possibility 63 // to create and access to dE/dx and range tab 61 // to create and access to dE/dx and range tables, or to calculate 64 // that information on fly. 62 // that information on fly. 65 63 66 // ------------------------------------------- 64 // ------------------------------------------------------------------- 67 // 65 // 68 66 >> 67 69 #ifndef G4EmModelManager_h 68 #ifndef G4EmModelManager_h 70 #define G4EmModelManager_h 1 69 #define G4EmModelManager_h 1 71 70 72 //....oooOO0OOooo........oooOO0OOooo........oo 71 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 73 72 74 #include "globals.hh" 73 #include "globals.hh" 75 #include "G4DataVector.hh" 74 #include "G4DataVector.hh" 76 #include "G4EmTableType.hh" 75 #include "G4EmTableType.hh" 77 #include "G4EmProcessSubType.hh" 76 #include "G4EmProcessSubType.hh" 78 #include "G4Region.hh" 77 #include "G4Region.hh" 79 78 80 #include "G4VEmModel.hh" << 81 #include "G4VEmFluctuationModel.hh" << 82 #include "G4DynamicParticle.hh" << 83 #include <iostream> << 84 << 85 class G4RegionModels 79 class G4RegionModels 86 { 80 { 87 81 88 friend class G4EmModelManager; 82 friend class G4EmModelManager; 89 83 >> 84 public: >> 85 90 private: 86 private: 91 87 92 G4RegionModels(G4int nMod, std::vector<G4int 88 G4RegionModels(G4int nMod, std::vector<G4int>& indx, 93 G4DataVector& lowE, const G4R << 89 G4DataVector& lowE, const G4Region* reg); 94 90 95 ~G4RegionModels(); 91 ~G4RegionModels(); 96 92 97 inline G4int SelectIndex(G4double e) const { << 93 G4int SelectIndex(G4double e) const { 98 G4int idx = 0; 94 G4int idx = 0; 99 if (nModelsForRegion>1) { 95 if (nModelsForRegion>1) { 100 idx = nModelsForRegion; 96 idx = nModelsForRegion; 101 // Loop checking, 03-Aug-2015, Vladimir << 97 do {idx--;} while (idx && e <= lowKineticEnergy[idx]); 102 do {--idx;} while (idx > 0 && e <= lowKi << 103 } 98 } 104 return theListOfModelIndexes[idx]; 99 return theListOfModelIndexes[idx]; 105 }; 100 }; 106 101 107 inline G4int ModelIndex(G4int n) const { << 102 G4int ModelIndex(G4int n) const { 108 return theListOfModelIndexes[n]; 103 return theListOfModelIndexes[n]; 109 }; 104 }; 110 105 111 inline G4int NumberOfModels() const { << 106 G4int NumberOfModels() const { 112 return nModelsForRegion; 107 return nModelsForRegion; 113 }; 108 }; 114 109 115 inline G4double LowEdgeEnergy(G4int n) const << 110 G4double LowEdgeEnergy(G4int n) const { 116 return lowKineticEnergy[n]; 111 return lowKineticEnergy[n]; 117 }; 112 }; 118 113 119 inline const G4Region* Region() const { << 114 const G4Region* Region() const { 120 return theRegion; 115 return theRegion; 121 }; 116 }; 122 117 123 G4RegionModels(G4RegionModels &) = delete; << 124 G4RegionModels & operator=(const G4RegionMod << 125 << 126 const G4Region* theRegion; 118 const G4Region* theRegion; 127 G4int nModelsForRegion; 119 G4int nModelsForRegion; 128 G4int* theListOfModelIndexes; 120 G4int* theListOfModelIndexes; 129 G4double* lowKineticEnergy; 121 G4double* lowKineticEnergy; 130 122 131 }; 123 }; 132 124 133 //....oooOO0OOooo........oooOO0OOooo........oo << 125 #include "G4VEmModel.hh" >> 126 #include "G4VEmFluctuationModel.hh" >> 127 #include "G4DynamicParticle.hh" 134 128 135 class G4Region; 129 class G4Region; 136 class G4ParticleDefinition; 130 class G4ParticleDefinition; >> 131 class G4DataVector; 137 class G4PhysicsVector; 132 class G4PhysicsVector; 138 class G4MaterialCutsCouple; 133 class G4MaterialCutsCouple; 139 134 140 //....oooOO0OOooo........oooOO0OOooo........oo 135 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 141 136 142 class G4EmModelManager 137 class G4EmModelManager 143 { 138 { 144 public: 139 public: 145 140 146 G4EmModelManager(); 141 G4EmModelManager(); 147 142 148 ~G4EmModelManager(); 143 ~G4EmModelManager(); 149 144 150 void Clear(); 145 void Clear(); 151 146 152 const G4DataVector* Initialise(const G4Parti << 147 const G4DataVector* Initialise(const G4ParticleDefinition*, 153 const G4Parti << 148 const G4ParticleDefinition*, 154 G4int verb); << 149 G4double, >> 150 G4int); >> 151 >> 152 const G4DataVector* Cuts() const; >> 153 >> 154 const G4DataVector* SubCutoff() const; 155 155 156 void FillDEDXVector(G4PhysicsVector*, const 156 void FillDEDXVector(G4PhysicsVector*, const G4MaterialCutsCouple*, 157 G4EmTableType t = fRestr << 157 G4EmTableType t = fRestricted); 158 158 159 void FillLambdaVector(G4PhysicsVector*, cons 159 void FillLambdaVector(G4PhysicsVector*, const G4MaterialCutsCouple*, 160 G4bool startFromNull = << 160 G4bool startFromNull = true, G4EmTableType t = fRestricted); 161 G4EmTableType t = fRes << 162 161 163 void AddEmModel(G4int, G4VEmModel*, G4VEmFlu << 162 G4VEmModel* SelectModel(G4double& energy, size_t& index); 164 const G4Region* r); << 165 163 166 // Get model pointer from the model list << 164 G4VEmModel* GetModel(G4int, G4bool ver = false); 167 G4VEmModel* GetModel(G4int idx, G4bool ver = << 168 165 169 // Get model pointer from the model list for << 166 G4int NumberOfModels() const; 170 // no check on material cuts couple index << 171 G4VEmModel* GetRegionModel(G4int idx, std::s << 172 167 173 // total number of models for material cut c << 168 void AddEmModel(G4int, G4VEmModel*, G4VEmFluctuationModel*, const G4Region*); 174 // no check on material cuts couple index << 169 175 G4int NumberOfRegionModels(std::size_t index << 170 void UpdateEmModel(const G4String&, G4double, G4double); 176 171 177 // Automatic documentation << 172 void DumpModelList(G4int verb); 178 void DumpModelList(std::ostream& out, G4int << 179 173 180 // Select model for given material cuts coup << 174 private: 181 inline G4VEmModel* SelectModel(G4double ener << 182 175 183 // Access to cuts << 176 // hide assignment operator 184 inline const G4DataVector* Cuts() const; << 185 177 186 // Set flag of fluorescence << 178 G4EmModelManager(G4EmModelManager &); 187 inline void SetFluoFlag(G4bool val); << 179 G4EmModelManager & operator=(const G4EmModelManager &right); 188 180 189 // total number of models << 181 // ===================================================================== 190 inline G4int NumberOfModels() const; << 191 << 192 // hide assignment operator << 193 G4EmModelManager(G4EmModelManager &) = delet << 194 G4EmModelManager & operator=(const G4EmModel << 195 182 196 private: 183 private: 197 184 198 const G4ParticleDefinition* particle = nullp << 185 G4DataVector theCuts; 199 const G4DataVector* theCuts = nullpt << 186 G4DataVector theSubCuts; 200 G4DataVector* theCutsNew = nul << 187 201 << 188 std::vector<G4VEmModel*> models; 202 // may be changed in run time << 189 std::vector<G4VEmFluctuationModel*> flucModels; 203 G4RegionModels* currRegionModel << 190 std::vector<const G4Region*> regions; 204 G4VEmModel* currModel = null << 191 std::vector<G4int> orderOfModels; 205 << 192 206 G4int nEmModels = 0; << 193 G4int nEmModels; 207 G4int nRegions = 0; << 194 G4int nRegions; 208 << 195 G4int nCouples; 209 G4int verboseLevel = 0 << 196 210 G4bool severalModels = << 197 G4int* idxOfRegionModels; 211 G4bool fluoFlag = false << 198 G4RegionModels** setOfRegionModels; 212 << 199 213 std::vector<G4VEmModel*> models; << 200 G4double minSubRange; 214 std::vector<G4VEmFluctuationModel*> flucMod << 201 G4double maxCutInRange; 215 std::vector<const G4Region*> regions << 202 G4double maxSubCutInRange; 216 std::vector<G4int> orderOf << 203 217 std::vector<G4int> isUsed; << 204 const G4ParticleDefinition* particle; >> 205 const G4ParticleDefinition* secondaryParticle; >> 206 const G4ParticleDefinition* theGamma; >> 207 const G4ParticleDefinition* thePositron; >> 208 >> 209 G4int verboseLevel; >> 210 >> 211 // cash >> 212 G4int currentIdx; 218 213 219 std::vector<G4int> idxOfRegionMod << 220 std::vector<G4RegionModels*> setOfRegionMod << 221 }; 214 }; 222 215 223 //....oooOO0OOooo........oooOO0OOooo........oo 216 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 224 //....oooOO0OOooo........oooOO0OOooo........oo 217 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 225 218 226 inline << 219 inline G4VEmModel* G4EmModelManager::SelectModel(G4double& kinEnergy, 227 G4VEmModel* G4EmModelManager::SelectModel(G4do << 220 size_t& index) 228 { 221 { 229 if(severalModels) { << 222 currentIdx = 230 if(nRegions > 1) { << 223 (setOfRegionModels[idxOfRegionModels[index]])->SelectIndex(kinEnergy); 231 currRegionModel = setOfRegionModels[idxO << 224 return models[currentIdx]; 232 } << 233 currModel = models[currRegionModel->Select << 234 } << 235 return currModel; << 236 } 225 } 237 226 238 //....oooOO0OOooo........oooOO0OOooo........oo 227 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 239 228 240 inline const G4DataVector* G4EmModelManager::C 229 inline const G4DataVector* G4EmModelManager::Cuts() const 241 { 230 { 242 return theCuts; << 231 return &theCuts; 243 } 232 } 244 233 245 //....oooOO0OOooo........oooOO0OOooo........oo 234 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 246 235 247 inline void G4EmModelManager::SetFluoFlag(G4bo << 236 inline const G4DataVector* G4EmModelManager::SubCutoff() const 248 { 237 { 249 fluoFlag = val; << 238 return &theSubCuts; 250 } 239 } 251 240 252 //....oooOO0OOooo........oooOO0OOooo........oo 241 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 253 242 254 inline G4int G4EmModelManager::NumberOfModels( 243 inline G4int G4EmModelManager::NumberOfModels() const 255 { 244 { 256 return nEmModels; 245 return nEmModels; 257 } 246 } 258 247 259 //....oooOO0OOooo........oooOO0OOooo........oo 248 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 260 249 261 #endif 250 #endif 262 251 263 252