89 std::vector<G4VEMDataSet*>*
93 std::vector<G4VEMDataSet*>*
set =
new std::vector<G4VEMDataSet*>;
101 G4int nOfBins = energyVector.size();
107 for (
size_t mLocal=0; mLocal<numOfCouples; mLocal++) {
115 G4double tcut = (*energyCuts)[mLocal];
120 for (
G4int i=0; i<nElements; i++) {
130 G4double density = nAtomsPerVolume[i];
135 energies->push_back(e);
137 log_energies->push_back(std::log10(e));
145 value *= elemCs*density;
147 cs->push_back(value);
149 if (value==0.) value=1e-300;
150 log_cs->push_back(std::log10(value));
160 set->push_back(setForMat);
173 if(energy > cutEnergy)