81 if(t0 >= tm)
return 0.0;
90 for (
size_t i=0; i<=
length; i++) {
100 <<
"; tMax(MeV)= " << tmax/
MeV
103 <<
"; xp[0]= " <<
xp[0]
128 if(t0 >= tm)
return 0.0;
138 for (
size_t i=0; i<=
length; i++) {
149 x += p[0]*(t0 - zmin - c*(std::atan(t0/c) - std::atan(zmin/c)));
155 <<
"; tMax(MeV)= " << tmax/
MeV
156 <<
"; e(MeV)= " << e/
MeV
182 if(t0 >= tm)
return 0.0;
189 for (
size_t i=0; i<=
length; i++) {
204 G4cout <<
"WARNING in G4eBremsstrahlungSpectrum::SampleEnergy:"
205 <<
" Majoranta " << amaj
230 sum += (1. - k*
xp[0])*std::log(x2/x1) + k*(x2 - x1);
233 for (
size_t i=0; i<
length-1; i++) {
239 sum += z2 - z1 + std::log(x2/x1)*(z1*x2 - z2*x1)/(x2 - x1);
242 if(sum < 0.0) sum = 0.0;
258 sum += (z2 - z1)*(1. - k*
xp[0]);
261 sum += 0.5*k*(z2 - z1);
264 for (
size_t i=0; i<
length-1; i++) {
270 sum += 0.5*(z2 - z1)*(x2 + x1) + z1*x2 - z2*x1;
273 if(sum < 0.0) sum = 0.0;
283 f = p[0] + (p[1] - p[0])*(x -
xp[0])/(
xp[1] -
xp[0]);
287 for (
size_t i=0; i<
length-1; i++) {
290 f = p[i] + (p[i+1] - p[i])*(x -
xp[i])/(
xp[i+1] -
xp[i]);
310 return kineticEnergy;
G4BremsstrahlungParameters * theBRparam
~G4eBremsstrahlungSpectrum()
G4double Function(G4double x, const G4DataVector &p) const
G4double Excitation(G4int Z, G4double kineticEnergy) const
G4eBremsstrahlungSpectrum(const G4DataVector &bins, const G4String &name)
G4GLOB_DLL std::ostream G4cout
G4double Parameter(G4int parameterIndex, G4int Z, G4double energy) const
G4double AverageValue(G4double xMin, G4double xMax, const G4DataVector &p) const
G4double SampleEnergy(G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
G4double AverageEnergy(G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
T max(const T t1, const T t2)
brief Return the largest of the two arguments
const G4double x[NPOINTSGL]
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
G4double IntSpectrum(G4double xMin, G4double xMax, const G4DataVector &p) const
G4double Probability(G4int Z, G4double tMin, G4double tMax, G4double kineticEnergy, G4int shell=0, const G4ParticleDefinition *pd=0) const
G4double MaxEnergyOfSecondaries(G4double kineticEnergy, G4int Z=0, const G4ParticleDefinition *pd=0) const
G4double ParameterC(G4int index) const