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: G4VEmModel.cc,v 1.37 2010-10-14 16:27:35 vnivanch Exp $ >> 27 // GEANT4 tag $Name: geant4-09-04-patch-01 $ >> 28 // 26 // ------------------------------------------- 29 // ------------------------------------------------------------------- 27 // 30 // 28 // GEANT4 Class file 31 // GEANT4 Class file 29 // 32 // 30 // 33 // 31 // File name: G4VEmModel 34 // File name: G4VEmModel 32 // 35 // 33 // Author: Vladimir Ivanchenko 36 // Author: Vladimir Ivanchenko 34 // 37 // 35 // Creation date: 25.07.2005 38 // Creation date: 25.07.2005 36 // 39 // 37 // Modifications: 40 // Modifications: 38 // 25.10.2005 Set default highLimit=100.TeV (V 41 // 25.10.2005 Set default highLimit=100.TeV (V.Ivanchenko) 39 // 06.02.2006 add method ComputeMeanFreePath() 42 // 06.02.2006 add method ComputeMeanFreePath() (mma) 40 // 16.02.2009 Move implementations of virtual 43 // 16.02.2009 Move implementations of virtual methods to source (VI) 41 // 44 // 42 // 45 // 43 // Class Description: 46 // Class Description: 44 // 47 // 45 // Abstract interface to energy loss models 48 // Abstract interface to energy loss models 46 49 47 // ------------------------------------------- 50 // ------------------------------------------------------------------- 48 // 51 // 49 52 50 #include "G4VEmModel.hh" 53 #include "G4VEmModel.hh" 51 #include "G4ElementData.hh" << 52 #include "G4LossTableManager.hh" 54 #include "G4LossTableManager.hh" 53 #include "G4LossTableBuilder.hh" << 54 #include "G4ProductionCutsTable.hh" 55 #include "G4ProductionCutsTable.hh" 55 #include "G4ParticleChangeForLoss.hh" 56 #include "G4ParticleChangeForLoss.hh" 56 #include "G4ParticleChangeForGamma.hh" 57 #include "G4ParticleChangeForGamma.hh" 57 #include "G4EmParameters.hh" << 58 #include "G4SystemOfUnits.hh" << 59 #include "G4EmUtility.hh" << 60 #include "G4Log.hh" << 61 #include "Randomize.hh" << 62 #include <iostream> << 63 58 64 //....oooOO0OOooo........oooOO0OOooo........oo 59 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 65 //....oooOO0OOooo........oooOO0OOooo........oo 60 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 66 61 67 G4VEmModel::G4VEmModel(const G4String& nam): 62 G4VEmModel::G4VEmModel(const G4String& nam): 68 inveplus(1.0/CLHEP::eplus), << 63 flucModel(0),anglModel(0), name(nam), lowLimit(0.1*keV), highLimit(100.0*TeV), 69 lowLimit(0.1*CLHEP::keV), << 64 eMinActive(0.0),eMaxActive(DBL_MAX), 70 highLimit(100.0*CLHEP::TeV), << 65 polarAngleLimit(0.0),secondaryThreshold(DBL_MAX),theLPMflag(false), 71 polarAngleLimit(CLHEP::pi), << 66 pParticleChange(0), 72 name(nam) << 67 currentCouple(0),currentElement(0), >> 68 nsec(5),flagDeexcitation(false) 73 { 69 { 74 xsec.resize(nsec); 70 xsec.resize(nsec); 75 fEmManager = G4LossTableManager::Instance(); << 71 nSelectors = 0; 76 fEmManager->Register(this); << 72 G4LossTableManager::Instance()->Register(this); 77 isMaster = fEmManager->IsMaster(); << 78 << 79 G4LossTableBuilder* bld = fEmManager->GetTab << 80 theDensityFactor = bld->GetDensityFactors(); << 81 theDensityIdx = bld->GetCoupleIndexes(); << 82 } 73 } 83 74 84 //....oooOO0OOooo........oooOO0OOooo........oo 75 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 85 76 86 G4VEmModel::~G4VEmModel() 77 G4VEmModel::~G4VEmModel() 87 { 78 { 88 if(localElmSelectors) { << 79 G4LossTableManager::Instance()->DeRegister(this); 89 for(G4int i=0; i<nSelectors; ++i) { << 80 G4int n = elmSelectors.size(); 90 delete (*elmSelectors)[i]; << 81 if(n > 0) { >> 82 for(G4int i=0; i<n; ++i) { >> 83 delete elmSelectors[i]; 91 } 84 } 92 delete elmSelectors; << 93 } 85 } 94 delete anglModel; 86 delete anglModel; 95 << 96 if(localTable && xSectionTable != nullptr) { << 97 xSectionTable->clearAndDestroy(); << 98 delete xSectionTable; << 99 xSectionTable = nullptr; << 100 } << 101 fEmManager->DeRegister(this); << 102 } 87 } 103 88 104 //....oooOO0OOooo........oooOO0OOooo........oo 89 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 105 90 106 G4ParticleChangeForLoss* G4VEmModel::GetPartic 91 G4ParticleChangeForLoss* G4VEmModel::GetParticleChangeForLoss() 107 { 92 { 108 G4ParticleChangeForLoss* p = nullptr; << 93 G4ParticleChangeForLoss* p = 0; 109 if (pParticleChange != nullptr) { << 94 if (pParticleChange) { 110 p = static_cast<G4ParticleChangeForLoss*>( 95 p = static_cast<G4ParticleChangeForLoss*>(pParticleChange); 111 } else { 96 } else { 112 p = new G4ParticleChangeForLoss(); 97 p = new G4ParticleChangeForLoss(); 113 pParticleChange = p; 98 pParticleChange = p; 114 } 99 } 115 if(fTripletModel != nullptr) { fTripletModel << 116 return p; 100 return p; 117 } 101 } 118 102 119 //....oooOO0OOooo........oooOO0OOooo........oo 103 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 120 104 121 G4ParticleChangeForGamma* G4VEmModel::GetParti 105 G4ParticleChangeForGamma* G4VEmModel::GetParticleChangeForGamma() 122 { 106 { 123 G4ParticleChangeForGamma* p = nullptr; << 107 G4ParticleChangeForGamma* p = 0; 124 if (pParticleChange != nullptr) { << 108 if (pParticleChange) { 125 p = static_cast<G4ParticleChangeForGamma*> 109 p = static_cast<G4ParticleChangeForGamma*>(pParticleChange); 126 } else { 110 } else { 127 p = new G4ParticleChangeForGamma(); 111 p = new G4ParticleChangeForGamma(); 128 pParticleChange = p; 112 pParticleChange = p; 129 } 113 } 130 if(fTripletModel != nullptr) { fTripletModel << 131 return p; 114 return p; 132 } 115 } 133 116 134 //....oooOO0OOooo........oooOO0OOooo........oo 117 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 135 118 136 void G4VEmModel::InitialiseElementSelectors(co << 119 void G4VEmModel::InitialiseElementSelectors(const G4ParticleDefinition* p, 137 co << 120 const G4DataVector& cuts) 138 { << 139 if(highLimit <= lowLimit) { return; } << 140 G4EmUtility::InitialiseElementSelectors(this << 141 localElmSelectors = true; << 142 } << 143 << 144 //....oooOO0OOooo........oooOO0OOooo........oo << 145 << 146 void G4VEmModel::InitialiseLocal(const G4Parti << 147 {} << 148 << 149 //....oooOO0OOooo........oooOO0OOooo........oo << 150 << 151 void G4VEmModel::InitialiseForMaterial(const G << 152 const G << 153 { 121 { 154 if(material != nullptr) { << 122 // initialise before run 155 G4int n = (G4int)material->GetNumberOfElem << 123 G4LossTableManager* man = G4LossTableManager::Instance(); 156 for(G4int i=0; i<n; ++i) { << 124 G4bool spline = man->SplineFlag(); 157 G4int Z = material->GetElement(i)->GetZa << 125 158 InitialiseForElement(part, Z); << 126 // two times less bins because probability functon is normalized >> 127 // so correspondingly is more smooth >> 128 G4int nbins = G4int(man->GetNumberOfBinsPerDecade() >> 129 * std::log10(highLimit/lowLimit) / 6.0); >> 130 if(nbins < 5) { nbins = 5; } >> 131 >> 132 G4ProductionCutsTable* theCoupleTable= >> 133 G4ProductionCutsTable::GetProductionCutsTable(); >> 134 G4int numOfCouples = theCoupleTable->GetTableSize(); >> 135 >> 136 // prepare vector >> 137 if(numOfCouples > nSelectors) { >> 138 elmSelectors.reserve(numOfCouples); >> 139 for(G4int i=nSelectors; i<numOfCouples; ++i) { elmSelectors.push_back(0); } >> 140 nSelectors = numOfCouples; >> 141 } >> 142 >> 143 // initialise vector >> 144 for(G4int i=0; i<numOfCouples; ++i) { >> 145 currentCouple = theCoupleTable->GetMaterialCutsCouple(i); >> 146 const G4Material* material = currentCouple->GetMaterial(); >> 147 G4int idx = currentCouple->GetIndex(); >> 148 >> 149 // selector already exist check if should be deleted >> 150 G4bool create = true; >> 151 if(elmSelectors[i]) { >> 152 if(material == elmSelectors[i]->GetMaterial()) { create = false; } >> 153 else { delete elmSelectors[i]; } 159 } 154 } 160 } << 155 if(create) { >> 156 elmSelectors[i] = new G4EmElementSelector(this,material,nbins, >> 157 lowLimit,highLimit,spline); >> 158 } >> 159 elmSelectors[i]->Initialise(p, cuts[idx]); >> 160 //elmSelectors[i]->Dump(p); >> 161 } 161 } 162 } 162 163 163 //....oooOO0OOooo........oooOO0OOooo........oo 164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 164 165 165 void G4VEmModel::InitialiseForElement(const G4 << 166 {} << 167 << 168 //....oooOO0OOooo........oooOO0OOooo........oo << 169 << 170 G4double G4VEmModel::ComputeDEDXPerVolume(cons 166 G4double G4VEmModel::ComputeDEDXPerVolume(const G4Material*, 171 cons << 167 const G4ParticleDefinition*, 172 G4do << 168 G4double,G4double) 173 { 169 { 174 return 0.0; 170 return 0.0; 175 } 171 } 176 172 177 //....oooOO0OOooo........oooOO0OOooo........oo 173 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 178 174 179 G4double G4VEmModel::CrossSectionPerVolume(con << 175 G4double G4VEmModel::CrossSectionPerVolume(const G4Material* material, 180 con << 176 const G4ParticleDefinition* p, 181 G4d << 177 G4double ekin, 182 G4d << 178 G4double emin, 183 G4d << 179 G4double emax) 184 { << 180 { 185 SetupForMaterial(p, mat, ekin); << 181 SetupForMaterial(p, material, ekin); 186 const G4double* theAtomNumDensityVector = ma << 182 G4double cross = 0.0; 187 G4int nelm = (G4int)mat->GetNumberOfElements << 183 const G4ElementVector* theElementVector = material->GetElementVector(); >> 184 const G4double* theAtomNumDensityVector = material->GetVecNbOfAtomsPerVolume(); >> 185 G4int nelm = material->GetNumberOfElements(); 188 if(nelm > nsec) { 186 if(nelm > nsec) { 189 xsec.resize(nelm); 187 xsec.resize(nelm); 190 nsec = nelm; 188 nsec = nelm; 191 } 189 } 192 G4double cross = 0.0; << 190 for (G4int i=0; i<nelm; i++) { 193 for (G4int i=0; i<nelm; ++i) { << 194 cross += theAtomNumDensityVector[i]* 191 cross += theAtomNumDensityVector[i]* 195 ComputeCrossSectionPerAtom(p,mat->GetEle << 192 ComputeCrossSectionPerAtom(p,(*theElementVector)[i],ekin,emin,emax); 196 xsec[i] = cross; 193 xsec[i] = cross; 197 } 194 } 198 return cross; 195 return cross; 199 } 196 } 200 197 201 //....oooOO0OOooo........oooOO0OOooo........oo 198 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 202 199 203 G4double G4VEmModel::GetPartialCrossSection(co << 200 const G4Element* G4VEmModel::SelectRandomAtom(const G4Material* material, 204 co << 201 const G4ParticleDefinition* pd, 205 G4 << 202 G4double kinEnergy, 206 { << 203 G4double tcut, 207 return 0.0; << 204 G4double tmax) 208 } << 205 { 209 << 206 const G4ElementVector* theElementVector = material->GetElementVector(); 210 //....oooOO0OOooo........oooOO0OOooo........oo << 207 G4int n = material->GetNumberOfElements() - 1; 211 << 208 currentElement = (*theElementVector)[n]; 212 void G4VEmModel::StartTracking(G4Track*) << 209 if (n > 0) { 213 {} << 210 G4double x = G4UniformRand()* 214 << 211 G4VEmModel::CrossSectionPerVolume(material,pd,kinEnergy,tcut,tmax); 215 //....oooOO0OOooo........oooOO0OOooo........oo << 216 << 217 const G4Element* G4VEmModel::SelectRandomAtom( << 218 << 219 << 220 << 221 << 222 { << 223 G4int n = (G4int)mat->GetNumberOfElements(); << 224 fCurrentElement = mat->GetElement(0); << 225 if (n > 1) { << 226 const G4double x = G4UniformRand()* << 227 G4VEmModel::CrossSectionPerVolume(mat,pd << 228 for(G4int i=0; i<n; ++i) { 212 for(G4int i=0; i<n; ++i) { 229 if (x <= xsec[i]) { 213 if (x <= xsec[i]) { 230 fCurrentElement = mat->GetElement(i); << 214 currentElement = (*theElementVector)[i]; 231 break; << 215 break; 232 } 216 } 233 } 217 } 234 } 218 } 235 return fCurrentElement; << 219 return currentElement; 236 } << 237 << 238 //....oooOO0OOooo........oooOO0OOooo........oo << 239 << 240 const G4Element* G4VEmModel::GetCurrentElement << 241 { << 242 const G4Element* elm = fCurrentElement; << 243 if(nullptr == elm && nullptr != mat) { << 244 elm = G4EmUtility::SampleRandomElement(mat << 245 } << 246 return elm; << 247 } << 248 << 249 //....oooOO0OOooo........oooOO0OOooo........oo << 250 << 251 G4int G4VEmModel::SelectRandomAtomNumber(const << 252 { << 253 const G4Element* elm = GetCurrentElement(mat << 254 return (nullptr == elm) ? 0 : elm->GetZasInt << 255 } << 256 << 257 //....oooOO0OOooo........oooOO0OOooo........oo << 258 << 259 const G4Isotope* G4VEmModel::GetCurrentIsotope << 260 { << 261 const G4Isotope* iso = nullptr; << 262 const G4Element* el = elm; << 263 if(nullptr == el && nullptr != fCurrentCoupl << 264 el = GetCurrentElement(fCurrentCouple->Get << 265 } << 266 if(nullptr != el) { << 267 iso = G4EmUtility::SampleRandomIsotope(el) << 268 } << 269 return iso; << 270 } << 271 << 272 //....oooOO0OOooo........oooOO0OOooo........oo << 273 << 274 G4int G4VEmModel::SelectIsotopeNumber(const G4 << 275 { << 276 auto iso = GetCurrentIsotope(elm); << 277 return (nullptr != iso) ? iso->GetN() : 0; << 278 } 220 } 279 221 280 //....oooOO0OOooo........oooOO0OOooo........oo 222 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 281 223 282 G4double G4VEmModel::ComputeCrossSectionPerAto 224 G4double G4VEmModel::ComputeCrossSectionPerAtom(const G4ParticleDefinition*, 283 << 225 G4double, G4double, G4double, 284 << 226 G4double, G4double) 285 { << 286 return 0.0; << 287 } << 288 << 289 //....oooOO0OOooo........oooOO0OOooo........oo << 290 << 291 G4double << 292 G4VEmModel::ComputeCrossSectionPerShell(const << 293 G4int, << 294 G4doub << 295 { 227 { 296 return 0.0; 228 return 0.0; 297 } 229 } 298 230 299 //....oooOO0OOooo........oooOO0OOooo........oo 231 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 300 232 301 void G4VEmModel::DefineForRegion(const G4Regio 233 void G4VEmModel::DefineForRegion(const G4Region*) 302 {} 234 {} 303 235 304 //....oooOO0OOooo........oooOO0OOooo........oo 236 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 305 237 306 void G4VEmModel::FillNumberOfSecondaries(G4int << 238 G4double G4VEmModel::MinEnergyCut(const G4ParticleDefinition*, 307 G4int << 239 const G4MaterialCutsCouple*) 308 { 240 { 309 numberOfTriplets = 0; << 241 return 0.0; 310 numberOfRecoil = 0; << 311 } 242 } 312 243 313 //....oooOO0OOooo........oooOO0OOooo........oo 244 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 314 245 315 G4double G4VEmModel::ChargeSquareRatio(const G 246 G4double G4VEmModel::ChargeSquareRatio(const G4Track& track) 316 { 247 { 317 return GetChargeSquareRatio(track.GetParticl 248 return GetChargeSquareRatio(track.GetParticleDefinition(), 318 track.GetMateria << 249 track.GetMaterial(), track.GetKineticEnergy()); 319 } 250 } 320 251 321 //....oooOO0OOooo........oooOO0OOooo........oo 252 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 322 253 323 G4double G4VEmModel::GetChargeSquareRatio(cons 254 G4double G4VEmModel::GetChargeSquareRatio(const G4ParticleDefinition* p, 324 cons << 255 const G4Material*, G4double) 325 { 256 { 326 const G4double q = p->GetPDGCharge()*inveplu << 257 G4double q = p->GetPDGCharge()/CLHEP::eplus; 327 return q*q; 258 return q*q; 328 } 259 } 329 260 330 //....oooOO0OOooo........oooOO0OOooo........oo 261 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 331 262 332 G4double G4VEmModel::GetParticleCharge(const G 263 G4double G4VEmModel::GetParticleCharge(const G4ParticleDefinition* p, 333 const G << 264 const G4Material*, G4double) 334 { 265 { 335 return p->GetPDGCharge(); 266 return p->GetPDGCharge(); 336 } 267 } 337 268 338 //....oooOO0OOooo........oooOO0OOooo........oo 269 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 339 270 340 void G4VEmModel::CorrectionsAlongStep(const G4 271 void G4VEmModel::CorrectionsAlongStep(const G4MaterialCutsCouple*, 341 const G4 << 272 const G4DynamicParticle*, 342 const G4 << 273 G4double&,G4double&,G4double) 343 {} 274 {} 344 275 345 //....oooOO0OOooo........oooOO0OOooo........oo 276 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 346 277 347 G4double G4VEmModel::Value(const G4MaterialCut << 278 void G4VEmModel::SampleDeexcitationAlongStep(const G4Material*, 348 const G4ParticleDef << 279 const G4Track&, 349 { << 280 G4double& ) 350 SetCurrentCouple(couple); << 281 {} 351 return pFactor*e*e*CrossSectionPerVolume(pBa << 352 } << 353 << 354 //....oooOO0OOooo........oooOO0OOooo........oo << 355 << 356 G4double G4VEmModel::MinPrimaryEnergy(const G4 << 357 const G4 << 358 G4double << 359 { << 360 return 0.0; << 361 } << 362 << 363 //....oooOO0OOooo........oooOO0OOooo........oo << 364 << 365 G4double G4VEmModel::MinEnergyCut(const G4Part << 366 const G4Mate << 367 { << 368 return 0.0; << 369 } << 370 282 371 //....oooOO0OOooo........oooOO0OOooo........oo 283 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 372 284 373 G4double G4VEmModel::MaxSecondaryEnergy(const 285 G4double G4VEmModel::MaxSecondaryEnergy(const G4ParticleDefinition*, 374 G4doub << 286 G4double kineticEnergy) 375 { 287 { 376 return kineticEnergy; 288 return kineticEnergy; 377 } 289 } 378 290 379 //....oooOO0OOooo........oooOO0OOooo........oo 291 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 380 292 381 void G4VEmModel::SetupForMaterial(const G4Part << 293 void G4VEmModel::SetupForMaterial(const G4ParticleDefinition*, 382 const G4Mate << 294 const G4Material*, G4double) 383 { << 295 {} 384 GetChargeSquareRatio(p, mat, ekin); << 385 } << 386 296 387 //....oooOO0OOooo........oooOO0OOooo........oo 297 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 388 298 389 void 299 void 390 G4VEmModel::SetParticleChange(G4VParticleChang 300 G4VEmModel::SetParticleChange(G4VParticleChange* p, G4VEmFluctuationModel* f) 391 { 301 { 392 if(p != nullptr && pParticleChange != p) { p << 302 if(p && pParticleChange != p) { pParticleChange = p; } 393 if(flucModel != f) { flucModel = f; } << 303 flucModel = f; 394 } << 395 << 396 //....oooOO0OOooo........oooOO0OOooo........oo << 397 << 398 void G4VEmModel::SetCrossSectionTable(G4Physic << 399 { << 400 xSectionTable = p; << 401 localTable = isLocal; << 402 } << 403 << 404 //....oooOO0OOooo........oooOO0OOooo........oo << 405 << 406 void G4VEmModel::SetLPMFlag(G4bool) << 407 { << 408 if (G4EmParameters::Instance()->Verbose() > << 409 G4ExceptionDescription ed; << 410 ed << "The obsolete method SetLPMFlag(..) << 411 << " is called. Please, use G4EmParamet << 412 << " instead"; << 413 G4Exception("G4VEmModel::SetLPMFlag", "em0 << 414 } << 415 } << 416 << 417 //....oooOO0OOooo........oooOO0OOooo........oo << 418 << 419 void G4VEmModel::SetMasterThread(G4bool) << 420 {} << 421 << 422 //....oooOO0OOooo........oooOO0OOooo........oo << 423 << 424 void G4VEmModel::ModelDescription(std::ostrea << 425 { << 426 outFile << "The description for this model h << 427 } 304 } 428 305 429 //....oooOO0OOooo........oooOO0OOooo........oo 306 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 430 307