59 G4cout <<
"G4LivermoreRayleighModel is constructed " <<
G4endl;
85 G4cout <<
"Calling Initialise() of G4LivermoreRayleighModel." <<
G4endl
98 char* path = std::getenv(
"G4LEDATA");
103 for(
G4int i=0; i<numOfCouples; ++i)
111 for (
G4int j=0; j<nelm; ++j)
113 G4int Z = (*theElementVector)[j]->GetZasInt();
141 G4cout <<
"Calling ReadData() of G4LivermoreRayleighModel"
147 const char* datadir = path;
151 datadir = std::getenv(
"G4LEDATA");
154 G4Exception(
"G4LivermoreRayleighModelModel::ReadData()",
"em0006",
156 "Environment variable G4LEDATA not defined");
168 std::ostringstream ostCS;
169 ostCS << datadir <<
"/livermore/rayl/re-cs-" << Z <<
".dat";
170 std::ifstream finCS(ostCS.str().c_str());
172 if( !finCS .is_open() )
175 ed <<
"G4LivermoreRayleighModel data file <" << ostCS.str().c_str()
176 <<
"> is not opened!" <<
G4endl;
178 ed,
"G4LEDATA version should be G4EMLOW6.27 or later.");
184 G4cout <<
"File " << ostCS.str()
185 <<
" is opened by G4LivermoreRayleighModel" <<
G4endl;
201 G4cout <<
"G4LivermoreRayleighModel::ComputeCrossSectionPerAtom()"
211 if(intZ < 1 || intZ >
maxZ) {
return xs; }
220 if(!pv) {
return xs; }
227 }
else if(e >= pv->
Energy(0)) {
233 G4cout <<
"****** DEBUG: tcs value for Z=" << Z <<
" at energy (MeV)="
235 G4cout <<
" cs (Geant4 internal unit)=" << xs <<
G4endl;
236 G4cout <<
" -> first E*E*cs value in CS data file (iu) =" << (*pv)[0]
238 G4cout <<
" -> last E*E*cs value in CS data file (iu) =" << (*pv)[
n]
240 G4cout <<
"*********************************************************"
249 std::vector<G4DynamicParticle*>*,
255 G4cout <<
"Calling SampleSecondaries() of G4LivermoreRayleighModel"