76 theParam(spec),verbose(0)
93 std::vector<G4VEMDataSet*>*
set =
new std::vector<G4VEMDataSet*>;
101 G4int nOfBins = energyVector.size();
107 for (
size_t mLocal=0; mLocal<numOfCouples; mLocal++) {
117 G4cout <<
"eIonisation CS for " << mLocal <<
"th material "
119 <<
" eEl= " << nElements <<
G4endl;
122 G4double tcut = (*energyCuts)[mLocal];
127 for (
G4int i=0; i<nElements; i++) {
138 G4double density = nAtomsPerVolume[i];
143 energies->push_back(e);
144 log_energies->push_back(std::log10(e));
152 value += cross * p * density;
154 if(
verbose>0 && mLocal == 0 && e>=1. && e<=0.)
156 G4cout <<
"G4eIonCrossSH: e(MeV)= " << e/
MeV
158 <<
" cross= " << cross
160 <<
" value= " << value
161 <<
" tcut(MeV)= " << tcut/
MeV
162 <<
" rho= " << density
168 if (value == 0.) value = 1e-300;
169 log_value = std::log10(value);
171 cs->push_back(value);
172 log_cs->push_back(log_value);
182 set->push_back(setForMat);
194 if(energy > cutEnergy)