62 G4cout <<
"Born excitation model is constructed " <<
G4endl;
91 G4cout <<
"Calling G4DNABornExcitationModel2::Initialise()" <<
G4endl;
96 G4Exception(
"G4DNABornExcitationModel2::Initialise",
"em0001",
97 FatalException,
"Model already initialized for another particle type.");
102 std::ostringstream fullFileName;
103 char *path = getenv(
"G4LEDATA");
107 G4Exception(
"G4DNABornExcitationModel2::Initialise",
"G4LEDATA-CHECK",
111 fullFileName << path;
115 fullFileName <<
"/dna/bornExcitation-e.dat";
121 fullFileName <<
"/dna/bornExcitation-p.dat";
133 for(
size_t level = 0; level<
fTableData->size(); ++level)
136 (*fTableData)(level)->SetSpline(
true);
139 size_t finalBin_i = 2000;
147 for(
size_t energy_i = 0; energy_i < finalBin_i; ++energy_i)
152 for(
size_t level = 0; level<
fTableData->size(); ++level)
154 finalXS += (*fTableData)(level)->
Value(energy);
168 G4cout <<
"Born excitation model is initialized " <<
G4endl
195 G4cout <<
"Calling CrossSectionPerVolume() of G4DNABornExcitationModel2"
222 G4cout <<
"__________________________________" <<
G4endl;
223 G4cout <<
"G4DNABornExcitationModel2 - XS INFO START" <<
G4endl;
225 G4cout <<
"Cross section per water molecule (cm^2)=" << sigma/
cm/
cm <<
G4endl;
226 G4cout <<
"Cross section per water molecule (cm^-1)=" << sigma*waterDensity/(1./
cm) << G4endl;
227 G4cout <<
"G4DNABornExcitationModel2 - XS INFO END" <<
G4endl;
230 return sigma*waterDensity;
244 G4cout <<
"Calling SampleSecondaries() of G4DNABornExcitationModel2"
252 G4double newEnergy = k - excitationEnergy;
279 G4Exception(
"G4DNABornExcitationModel2::GetPartialCrossSection",
282 "Model initialized for another particle type.");
306 partialXS = (*fTableData)(i)->
Value(k);
307 if (partialXS > value)