38 G4cout <<
"Regular transparent X-ray TR radiator EM process is called"
52 out <<
"Process describing radiator of X-ray transition radiation.\n";
59 G4double result, sum = 0., tmp, cof1, cof2, cofMin, aMa, bMb, sigma;
92 result = (k - cof1) * (k - cof1) * (k + cof2) * (k + cof2);
97 0.5 * std::sin(tmp) * std::sin(tmp) * std::abs(k - cofMin) / result;
101 sum += std::sin(tmp) * std::sin(tmp) * std::abs(k - cofMin) / result;
104 result = 4. * (cof1 + cof2) * (cof1 + cof2) * sum /
energy;
105 result *= (1. - std::exp(-
fPlateNumber * sigma)) / (1. - std::exp(-sigma));
128 G4complex Ha(Qa * std::cos(aZa), -Qa * std::sin(aZa));
129 G4complex Hb(Qb * std::cos(bZb), -Qb * std::sin(bZb));
133 1.0 / ((1. -
Q) * (1. -
Q) +
134 4. *
Q * std::sin(0.5 * (aZa + bZb)) * std::sin(0.5 * (aZa + bZb)));
137 G4complex F2 = (1.0 - Ha) * (1.0 - Ha) * Hb * (1.0 - Hs) * (1.0 - Hs) *
140 return 2.0 * std::real(R);
G4double D(G4double temp)
static constexpr double pi
std::complex< G4double > G4complex
G4GLOB_DLL std::ostream G4cout
G4double GetPlateLinearPhotoAbs(G4double)
G4double GetGasFormationZone(G4double, G4double, G4double)
G4complex OneInterfaceXTRdEdx(G4double energy, G4double gamma, G4double varAngle)
G4double GetGasCompton(G4double)
G4double GetPlateFormationZone(G4double, G4double, G4double)
G4double GetGasLinearPhotoAbs(G4double)
G4double GetPlateCompton(G4double)
G4XTRTransparentRegRadModel(G4LogicalVolume *anEnvelope, G4Material *, G4Material *, G4double, G4double, G4int, const G4String &processName="XTRTransparentRegRadModel")
~G4XTRTransparentRegRadModel()
G4double GetStackFactor(G4double energy, G4double gamma, G4double varAngle) override
void ProcessDescription(std::ostream &) const override
G4double SpectralXTRdEdx(G4double energy) override
G4double energy(const ThreeVector &p, const G4double m)