82 std::vector<G4RDVEMDataSet*>*
set =
new std::vector<G4RDVEMDataSet*>;
86 G4int nOfBins = energyVector.size();
92 for (
size_t m=0;
m<numOfCouples;
m++) {
101 G4cout <<
"eIonisation CS for " <<
m <<
"th material "
103 <<
" eEl= " << nElements <<
G4endl;
111 for (
G4int i=0; i<nElements; i++) {
117 G4double density = nAtomsPerVolume[i];
122 energies->push_back(e);
129 value += cross * p * density;
131 if(verbose>0 &&
m == 0 && e>=1. && e<=0.) {
132 G4cout <<
"G4eIonCrossSH: e(MeV)= " << e/
MeV
134 <<
" cross= " << cross
136 <<
" value= " << value
137 <<
" tcut(MeV)= " << tcut/
MeV
138 <<
" rho= " << density
145 cs->push_back(value);
151 set->push_back(setForMat);