52 for (
G4int k =0; k<size; k++)
71 std::vector<G4int> numberOfIonisations;
75 for (
G4int i = 0; i<numberOfElements; i++)
81 G4int ionisations = 0;
82 if(averageNumberOfIonisations > 0.0) {
86 numberOfIonisations.push_back(ionisations);
89 return numberOfIonisations;
101 G4double averageEnergy = energy - eLoss/2.;
103 size_t indexInMaterialTable = couple->
GetIndex();
std::vector< G4int > GenerateNumberOfIonisations(const G4MaterialCutsCouple *couple, G4double incidentEnergy, G4double eLoss) const
virtual G4double FindValue(G4double x, G4int componentId=0) const =0
G4long G4Poisson(G4double mean)
G4double AverageNOfIonisations(const G4MaterialCutsCouple *couple, G4int index, G4double energy, G4double eLoss) const
std::vector< G4VEMDataSet * > xsis
G4double energy(const ThreeVector &p, const G4double m)
size_t GetNumberOfElements() const
void AddXsiTable(G4VEMDataSet *p)
const G4Material * GetMaterial() const