Changeset 1192 for trunk/source/processes/electromagnetic/lowenergy/src/G4DNAMillerGreenExcitationModel.cc
- Timestamp:
- Nov 19, 2009, 2:53:25 PM (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/source/processes/electromagnetic/lowenergy/src/G4DNAMillerGreenExcitationModel.cc
r1058 r1192 24 24 // ******************************************************************** 25 25 // 26 // $Id: G4DNAMillerGreenExcitationModel.cc,v 1. 3 2009/02/16 11:00:11sincerti Exp $27 // GEANT4 tag $Name: geant4-09-03-beta-cand-01$26 // $Id: G4DNAMillerGreenExcitationModel.cc,v 1.6 2009/08/13 11:32:47 sincerti Exp $ 27 // GEANT4 tag $Name: emlowen-V09-02-64 $ 28 28 // 29 29 … … 49 49 // 4 = entering in methods 50 50 51 G4cout << "Miller & Green excitation model is constructed " << G4endl; 52 51 if( verboseLevel>0 ) 52 { 53 G4cout << "Miller & Green excitation model is constructed " << G4endl; 54 } 53 55 } 54 56 … … 186 188 187 189 // 188 189 G4cout << "Miller & Green excitation model is initialized " << G4endl 190 << "Energy range: " 191 << LowEnergyLimit() / eV << " eV - " 192 << HighEnergyLimit() / keV << " keV for " 193 << particle->GetParticleName() 194 << G4endl; 195 190 if( verboseLevel>0 ) 191 { 192 G4cout << "Miller & Green excitation model is initialized " << G4endl 193 << "Energy range: " 194 << LowEnergyLimit() / eV << " eV - " 195 << HighEnergyLimit() / keV << " keV for " 196 << particle->GetParticleName() 197 << G4endl; 198 } 199 196 200 if(!isInitialised) 197 201 { … … 224 228 const G4Material* material = couple->GetMaterial(); 225 229 226 size_t j = material->GetNumberOfElements(); 227 while (j>0) 230 if (material->GetName() == "G4_WATER") 228 231 { 229 j--; 230 const G4Element* element(material->GetElement(j)); 231 if (element->GetZ() == 8.) 232 { 233 G4double density = material->GetAtomicNumDensityVector()[j]; 234 if (density > 0.) 235 { 236 flagMaterialIsWater = true; 237 densityWater = density; 238 239 if (verboseLevel > 3) 240 G4cout << "Water material is found with density(cm^-3)=" << density/(cm*cm*cm) << G4endl; 241 } 242 } 232 G4double density = material->GetAtomicNumDensityVector()[1]; 233 flagMaterialIsWater = true; 234 densityWater = density; 235 236 if (verboseLevel > 3) 237 G4cout << "****** Water material is found with density(cm^-3)=" << density/(cm*cm*cm) << G4endl; 243 238 } 244 239 245 240 } 246 } // if(numOfCouples>0) 241 242 } // if(numOfCouples>0) 247 243 248 244 } // if (theCoupleTable)
Note: See TracChangeset
for help on using the changeset viewer.