Geant4 Cross Reference |
1 // 1 // 2 // ******************************************* 2 // ******************************************************************** 3 // * License and Disclaimer << 3 // * DISCLAIMER * 4 // * 4 // * * 5 // * The Geant4 software is copyright of th << 5 // * The following disclaimer summarizes all the specific disclaimers * 6 // * the Geant4 Collaboration. It is provided << 6 // * of contributors to this software. The specific disclaimers,which * 7 // * conditions of the Geant4 Software License << 7 // * govern, are listed with their locations in: * 8 // * LICENSE and available at http://cern.ch/ << 8 // * http://cern.ch/geant4/license * 9 // * include a list of copyright holders. << 10 // * 9 // * * 11 // * Neither the authors of this software syst 10 // * Neither the authors of this software system, nor their employing * 12 // * institutes,nor the agencies providing fin 11 // * institutes,nor the agencies providing financial support for this * 13 // * work make any representation or warran 12 // * work make any representation or warranty, express or implied, * 14 // * regarding this software system or assum 13 // * regarding this software system or assume any liability for its * 15 // * use. Please see the license in the file << 14 // * use. * 16 // * for the full disclaimer and the limitatio << 17 // * 15 // * * 18 // * This code implementation is the result << 16 // * This code implementation is the intellectual property of the * 19 // * technical work of the GEANT4 collaboratio << 17 // * GEANT4 collaboration. * 20 // * By using, copying, modifying or distri << 18 // * By copying, distributing or modifying the Program (or any work * 21 // * any work based on the software) you ag << 19 // * based on the Program) you indicate your acceptance of this * 22 // * use in resulting scientific publicati << 20 // * statement, and all its terms. * 23 // * acceptance of all terms of the Geant4 Sof << 24 // ******************************************* 21 // ******************************************************************** 25 // 22 // 26 // 23 // 27 // ------------------------------------------- 24 // ------------------------------------------------------------------- 28 // 25 // 29 // GEANT4 Class header file 26 // GEANT4 Class header file 30 // 27 // >> 28 // 31 // File name: G4EmModelManager 29 // File name: G4EmModelManager 32 // 30 // 33 // Author: Vladimir Ivanchenko 31 // Author: Vladimir Ivanchenko 34 // 32 // 35 // Creation date: 07.05.2002 33 // Creation date: 07.05.2002 36 // 34 // 37 // Modifications: 35 // Modifications: 38 // 36 // 39 // 03-12-02 V.Ivanchenko fix a bug in model se 37 // 03-12-02 V.Ivanchenko fix a bug in model selection 40 // 20-01-03 Migrade to cut per region (V.Ivanc 38 // 20-01-03 Migrade to cut per region (V.Ivanchenko) 41 // 27-01-03 Make models region aware (V.Ivanch 39 // 27-01-03 Make models region aware (V.Ivanchenko) 42 // 13-02-03 The set of models is defined for r 40 // 13-02-03 The set of models is defined for region (V.Ivanchenko) 43 // 26-03-03 Add GetDEDXDispersion (V.Ivanchenk 41 // 26-03-03 Add GetDEDXDispersion (V.Ivanchenko) 44 // 13-04-03 Add startFromNull (V.Ivanchenko) 42 // 13-04-03 Add startFromNull (V.Ivanchenko) 45 // 13-05-03 Add calculation of precise range ( << 46 // 21-07-03 Add UpdateEmModel method (V.Ivanch << 47 // 03-11-03 Substitute STL vector for G4Region << 48 // 11-04-05 Remove access to fluctuation model << 49 // 10-01-06 PreciseRange -> CSDARange (V.Ivant << 50 // 20-01-06 Introduce G4EmTableType and reduci << 51 // 13-05-06 Add GetModel by index method (VI) << 52 // 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 // 43 // 58 // Class Description: 44 // Class Description: 59 // 45 // 60 // It is the unified energy loss process it ca 46 // It is the unified energy loss process it calculates the continuous 61 // energy loss for charged particles using a s 47 // energy loss for charged particles using a set of Energy Loss 62 // models valid for different energy regions. 48 // models valid for different energy regions. There are a possibility 63 // to create and access to dE/dx and range tab 49 // to create and access to dE/dx and range tables, or to calculate 64 // that information on fly. 50 // that information on fly. 65 51 66 // ------------------------------------------- 52 // ------------------------------------------------------------------- 67 // 53 // 68 54 >> 55 69 #ifndef G4EmModelManager_h 56 #ifndef G4EmModelManager_h 70 #define G4EmModelManager_h 1 57 #define G4EmModelManager_h 1 71 58 72 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 73 60 74 #include "globals.hh" 61 #include "globals.hh" 75 #include "G4DataVector.hh" 62 #include "G4DataVector.hh" 76 #include "G4EmTableType.hh" << 77 #include "G4EmProcessSubType.hh" << 78 #include "G4Region.hh" << 79 << 80 #include "G4VEmModel.hh" << 81 #include "G4VEmFluctuationModel.hh" << 82 #include "G4DynamicParticle.hh" << 83 #include <iostream> << 84 63 85 class G4RegionModels 64 class G4RegionModels 86 { 65 { 87 66 88 friend class G4EmModelManager; 67 friend class G4EmModelManager; 89 68 >> 69 public: >> 70 90 private: 71 private: 91 72 92 G4RegionModels(G4int nMod, std::vector<G4int << 73 G4RegionModels(G4int nMod, G4std::vector<G4int>& list, G4DataVector& lowE); 93 G4DataVector& lowE, const G4R << 94 74 95 ~G4RegionModels(); 75 ~G4RegionModels(); 96 76 97 inline G4int SelectIndex(G4double e) const { << 77 G4int SelectIndex(G4double e) const { 98 G4int idx = 0; 78 G4int idx = 0; 99 if (nModelsForRegion>1) { 79 if (nModelsForRegion>1) { 100 idx = nModelsForRegion; 80 idx = nModelsForRegion; 101 // Loop checking, 03-Aug-2015, Vladimir << 81 do {idx--;} while (idx && e <= lowKineticEnergy[idx]); 102 do {--idx;} while (idx > 0 && e <= lowKi << 103 } 82 } 104 return theListOfModelIndexes[idx]; 83 return theListOfModelIndexes[idx]; 105 }; 84 }; 106 85 107 inline G4int ModelIndex(G4int n) const { << 86 G4int ModelIndex(G4int n) const { 108 return theListOfModelIndexes[n]; 87 return theListOfModelIndexes[n]; 109 }; 88 }; 110 89 111 inline G4int NumberOfModels() const { << 90 G4int NumberOfModels() const { 112 return nModelsForRegion; 91 return nModelsForRegion; 113 }; 92 }; 114 93 115 inline G4double LowEdgeEnergy(G4int n) const << 116 return lowKineticEnergy[n]; << 117 }; << 118 << 119 inline const G4Region* Region() const { << 120 return theRegion; << 121 }; << 122 << 123 G4RegionModels(G4RegionModels &) = delete; << 124 G4RegionModels & operator=(const G4RegionMod << 125 << 126 const G4Region* theRegion; << 127 G4int nModelsForRegion; 94 G4int nModelsForRegion; 128 G4int* theListOfModelIndexes; 95 G4int* theListOfModelIndexes; 129 G4double* lowKineticEnergy; 96 G4double* lowKineticEnergy; 130 97 131 }; 98 }; 132 99 133 //....oooOO0OOooo........oooOO0OOooo........oo << 100 #include "G4VEmModel.hh" >> 101 #include "G4VEmFluctuationModel.hh" >> 102 #include "G4DynamicParticle.hh" 134 103 135 class G4Region; 104 class G4Region; 136 class G4ParticleDefinition; 105 class G4ParticleDefinition; >> 106 class G4DataVector; 137 class G4PhysicsVector; 107 class G4PhysicsVector; 138 class G4MaterialCutsCouple; 108 class G4MaterialCutsCouple; 139 109 140 //....oooOO0OOooo........oooOO0OOooo........oo 110 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 141 111 142 class G4EmModelManager 112 class G4EmModelManager 143 { 113 { 144 public: 114 public: 145 115 146 G4EmModelManager(); 116 G4EmModelManager(); 147 117 148 ~G4EmModelManager(); << 118 ~G4EmModelManager(); 149 119 150 void Clear(); 120 void Clear(); 151 121 152 const G4DataVector* Initialise(const G4Parti << 122 const G4DataVector* Initialise(const G4ParticleDefinition*, 153 const G4Parti << 123 const G4ParticleDefinition*, 154 G4int verb); << 124 G4double, >> 125 G4int); 155 126 156 void FillDEDXVector(G4PhysicsVector*, const << 127 const G4DataVector* Cuts() const; 157 G4EmTableType t = fRestr << 158 128 159 void FillLambdaVector(G4PhysicsVector*, cons << 129 const G4DataVector* SubCutoff() const; 160 G4bool startFromNull = << 161 G4EmTableType t = fRes << 162 130 163 void AddEmModel(G4int, G4VEmModel*, G4VEmFlu << 131 void FillDEDXVector(G4PhysicsVector*, const G4MaterialCutsCouple*); 164 const G4Region* r); << 165 132 166 // Get model pointer from the model list << 133 void FillLambdaVector(G4PhysicsVector*, const G4MaterialCutsCouple*, 167 G4VEmModel* GetModel(G4int idx, G4bool ver = << 134 G4bool startFromNull = true); 168 135 169 // Get model pointer from the model list for << 136 void FillSubLambdaVector(G4PhysicsVector*, const G4MaterialCutsCouple*, 170 // no check on material cuts couple index << 137 G4bool startFromNull = true); 171 G4VEmModel* GetRegionModel(G4int idx, std::s << 172 138 173 // total number of models for material cut c << 139 G4VEmModel* SelectModel(G4double& energy, size_t& index); 174 // no check on material cuts couple index << 175 G4int NumberOfRegionModels(std::size_t index << 176 140 177 // Automatic documentation << 141 G4double SampleFluctuations(const G4Material* material, 178 void DumpModelList(std::ostream& out, G4int << 142 const G4DynamicParticle* dp, >> 143 G4double& tmax, >> 144 G4double& length, >> 145 G4double& eloss, >> 146 G4double& preStepKinEnergy, >> 147 size_t& index); >> 148 >> 149 G4double GetDEDXDispersion( const G4Material* material, >> 150 const G4DynamicParticle* dp, >> 151 G4double& tmax, >> 152 G4double& length, >> 153 size_t& index); 179 154 180 // Select model for given material cuts coup << 155 void AddEmModel(G4int, G4VEmModel*, G4VEmFluctuationModel*, const G4Region*); 181 inline G4VEmModel* SelectModel(G4double ener << 182 156 183 // Access to cuts << 157 private: 184 inline const G4DataVector* Cuts() const; << 185 158 186 // Set flag of fluorescence << 159 // hide assignment operator 187 inline void SetFluoFlag(G4bool val); << 188 160 189 // total number of models << 161 G4EmModelManager(G4EmModelManager &); 190 inline G4int NumberOfModels() const; << 162 G4EmModelManager & operator=(const G4EmModelManager &right); 191 163 192 // hide assignment operator << 164 // ===================================================================== 193 G4EmModelManager(G4EmModelManager &) = delet << 194 G4EmModelManager & operator=(const G4EmModel << 195 165 196 private: 166 private: 197 167 198 const G4ParticleDefinition* particle = nullp << 168 G4DataVector theCuts; 199 const G4DataVector* theCuts = nullpt << 169 G4DataVector theSubCuts; 200 G4DataVector* theCutsNew = nul << 170 201 << 171 G4std::vector<G4VEmModel*> models; 202 // may be changed in run time << 172 G4std::vector<G4VEmFluctuationModel*> flucModels; 203 G4RegionModels* currRegionModel << 173 G4std::vector<const G4Region*> regions; 204 G4VEmModel* currModel = null << 174 G4std::vector<G4int> orderOfModels; 205 << 175 G4DataVector upperEkin; 206 G4int nEmModels = 0; << 176 207 G4int nRegions = 0; << 177 G4int nEmModels; 208 << 178 G4int nRegions; 209 G4int verboseLevel = 0 << 179 G4int nCouples; 210 G4bool severalModels = << 180 211 G4bool fluoFlag = false << 181 G4std::vector<G4int> idxOfRegionModels; 212 << 182 G4std::vector<G4RegionModels*> setOfRegionModels; 213 std::vector<G4VEmModel*> models; << 183 214 std::vector<G4VEmFluctuationModel*> flucMod << 184 G4double minSubRange; 215 std::vector<const G4Region*> regions << 185 216 std::vector<G4int> orderOf << 186 const G4ParticleDefinition* particle; 217 std::vector<G4int> isUsed; << 187 const G4ParticleDefinition* secondaryParticle; >> 188 >> 189 G4int verboseLevel; >> 190 >> 191 // cash >> 192 G4int currentIdx; 218 193 219 std::vector<G4int> idxOfRegionMod << 220 std::vector<G4RegionModels*> setOfRegionMod << 221 }; 194 }; 222 195 223 //....oooOO0OOooo........oooOO0OOooo........oo 196 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 224 //....oooOO0OOooo........oooOO0OOooo........oo 197 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 225 198 226 inline << 199 inline G4VEmModel* G4EmModelManager::SelectModel(G4double& kinEnergy, size_t& index) 227 G4VEmModel* G4EmModelManager::SelectModel(G4do << 228 { 200 { 229 if(severalModels) { << 201 currentIdx = (setOfRegionModels[idxOfRegionModels[index]])->SelectIndex(kinEnergy); 230 if(nRegions > 1) { << 202 return models[currentIdx]; 231 currRegionModel = setOfRegionModels[idxO << 232 } << 233 currModel = models[currRegionModel->Select << 234 } << 235 return currModel; << 236 } 203 } 237 204 >> 205 238 //....oooOO0OOooo........oooOO0OOooo........oo 206 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 239 207 240 inline const G4DataVector* G4EmModelManager::C << 208 inline G4double G4EmModelManager::SampleFluctuations( >> 209 const G4Material* material, >> 210 const G4DynamicParticle* dp, >> 211 G4double& tmax, >> 212 G4double& length, >> 213 G4double& eloss, >> 214 G4double& kinEnergy, >> 215 size_t& index) >> 216 241 { 217 { 242 return theCuts; << 218 currentIdx = (setOfRegionModels[idxOfRegionModels[index]])->SelectIndex(kinEnergy); >> 219 G4VEmFluctuationModel* fm = flucModels[currentIdx]; >> 220 G4double x = eloss; >> 221 if(fm) x = fm->SampleFluctuations(material,dp,tmax,length,eloss); >> 222 return x; 243 } 223 } 244 224 245 //....oooOO0OOooo........oooOO0OOooo........oo 225 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 246 226 247 inline void G4EmModelManager::SetFluoFlag(G4bo << 227 inline G4double G4EmModelManager::GetDEDXDispersion( >> 228 const G4Material* material, >> 229 const G4DynamicParticle* dp, >> 230 G4double& tmax, >> 231 G4double& length, >> 232 size_t& index) >> 233 { >> 234 currentIdx = (setOfRegionModels[idxOfRegionModels[index]]) >> 235 ->SelectIndex(dp->GetKineticEnergy()); >> 236 G4VEmFluctuationModel* fm = flucModels[currentIdx]; >> 237 G4double x = 0.0; >> 238 if(fm) x = fm->Dispersion(material,dp,tmax,length); >> 239 return x; >> 240 } >> 241 >> 242 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... >> 243 >> 244 inline const G4DataVector* G4EmModelManager::Cuts() const 248 { 245 { 249 fluoFlag = val; << 246 return &theCuts; 250 } 247 } 251 248 252 //....oooOO0OOooo........oooOO0OOooo........oo 249 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 253 250 254 inline G4int G4EmModelManager::NumberOfModels( << 251 inline const G4DataVector* G4EmModelManager::SubCutoff() const 255 { 252 { 256 return nEmModels; << 253 return &theSubCuts; 257 } 254 } 258 255 259 //....oooOO0OOooo........oooOO0OOooo........oo 256 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 260 257 261 #endif 258 #endif 262 259 263 260