82 if(t0 >= tm)
return 0.0;
91 for (
size_t i=0; i<=length; i++) {
92 p.push_back(theBRparam->
Parameter(i, Z, e));
96 G4double y = IntSpectrum(z0, 1.0, p);
101 <<
"; tMax(MeV)= " << tmax/
MeV
104 <<
"; xp[0]= " << xp[0]
129 if(t0 >= tm)
return 0.0;
139 for (
size_t i=0; i<=length; i++) {
140 p.push_back(theBRparam->
Parameter(i, Z, e));
143 G4double x = AverageValue(t0, tm, p);
144 G4double y = IntSpectrum(z0, 1.0, p);
150 x += p[0]*(t0 - zmin - c*(std::atan(t0/c) - std::atan(zmin/c)));
156 <<
"; tMax(MeV)= " << tmax/
MeV
157 <<
"; e(MeV)= " << e/
MeV
183 if(t0 >= tm)
return 0.0;
190 for (
size_t i=0; i<=length; i++) {
191 p.push_back(theBRparam->
Parameter(i, Z, e));
193 G4double amaj =
std::max(p[length], 1. - (p[1] - p[0])*xp[0]/(xp[1] - xp[0]) );
202 fun = Function(tgam, p);
205 G4cout <<
"WARNING in G4eBremsstrahlungSpectrum::SampleEnergy:"
206 <<
" Majoranta " << amaj
230 G4double k = (p[1] - p[0])/(xp[1] - xp[0]);
231 sum += (1. - k*xp[0])*std::log(x2/x1) + k*(x2 - x1);
234 for (
size_t i=0; i<length-1; i++) {
240 sum += z2 - z1 + std::log(x2/x1)*(z1*x2 - z2*x1)/(x2 - x1);
243 if(sum < 0.0) sum = 0.0;
258 G4double k = (p[1] - p[0])/(xp[1] - xp[0]);
259 sum += (z2 - z1)*(1. - k*xp[0]);
262 sum += 0.5*k*(z2 - z1);
265 for (
size_t i=0; i<length-1; i++) {
271 sum += 0.5*(z2 - z1)*(x2 + x1) + z1*x2 - z2*x1;
274 if(sum < 0.0) sum = 0.0;
284 f = p[0] + (p[1] - p[0])*(x - xp[0])/(xp[1] - xp[0]);
288 for (
size_t i=0; i<length-1; i++) {
291 f = p[i] + (p[i+1] - p[i])*(x - xp[i])/(xp[i+1] - xp[i]);
311 return kineticEnergy;
~G4eBremsstrahlungSpectrum()
G4double Excitation(G4int Z, G4double kineticEnergy) const
G4eBremsstrahlungSpectrum(const G4DataVector &bins, const G4String &name)
G4GLOB_DLL std::ostream G4cout
static constexpr double eV
G4double Parameter(G4int parameterIndex, G4int Z, G4double energy) const
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
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
T min(const T t1, const T t2)
brief Return the smallest of the two arguments
static constexpr double MeV
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