52 : gammaPol2(false), electronPol3(false)
80 if(cosT > 1.+1.
e-8) cosT = 1.;
81 if(cosT < -1.-1.
e-8) cosT = -1.;
85 if(sinT2 > 1. + 1.
e-8) sinT2 = 1.;
86 if(sinT2 < 0.) sinT2 = 0.;
95 unpXS = (eps2 + 1. - eps*sinT2)/(2.*eps);
104 G4double PHI21 = -sinT2 + 0.5*(cos2T + 3.)*pol0.
x() - ((1. -
eps)/eps)*sinT*
polzx;
107 G4double PHI23 = ((eps2 + 1.)/
eps)*cosT*pol0.
z() - ((1. -
eps)/eps)*(eps*cosT2 + 1.)*pol1.
z();
108 PHI23 += 0.5*(1. -
eps)*sin2T*pol1.
x();
117 G4double PHI31 = 0., PHI31add = 0., PHI33 = 0., PHI33add = 0.;
119 if ((1. - eps) > 1.
e-12){
120 G4double helpVar = std::sqrt(eps2 - 2.*cosT*eps + 1.);
122 PHI31 = (1. -
eps)*(1. + cosT)*sinT*pol0.
z();
123 PHI31 += (-eps*cosT3 + eps*cosT2 + (eps - 2.)*cosT +
eps)*pol1.
x();
124 PHI31 += -(eps*cosT2 - eps*cosT + cosT + 1.)*sinT*pol1.
z();
127 PHI31add = -eps*
sqr(1. - cosT)*(1. + cosT)*
polxx;
129 PHI31add += -(-eps2 + cosT*(cosT*eps - eps + 1.)*eps + eps - 1.)*sinT*
polxz/
eps;
130 PHI31add /= 2.*helpVar;
132 PHI33 = ((1. -
eps)/eps)*(-eps*cosT2 + eps*(eps + 1.)*cosT - 1.)*pol0.
z();
133 PHI33 += -(eps*cosT2 + (1. -
eps)*eps*cosT + 1.)*sinT*pol1.
x();
134 PHI33 += -(-eps2*cosT3 + eps*(eps2 - eps + 1.)*cosT2 - cosT + eps2)*pol1.
z()/
eps;
135 PHI33 /= -2.*helpVar;
137 PHI33add = (eps*(eps - cosT - 1.)*cosT + 1.)*sinT*
polxx;
138 PHI33add += -(-eps2 + cosT*eps + eps - 1.)*sinT2*
polxz;
139 PHI33add += (eps - 1.)*(cosT -
eps)*sinT*
polyy;
140 PHI33add /= -2.*helpVar;
142 PHI31 = -pol1.
z() - (X - 1.)*std::sqrt(1. - eps)*pol1.
x()/std::sqrt(2.*X);
143 PHI31add = -(-X*X*pol1.
z() - 2.*X*(2.*pol0.
z() - pol1.
z()) - (4.*pol0.
x() + 5.)*pol1.
z())*(1. - eps)/(4.*
X);
145 PHI33 = pol1.
x() - (X - 1.)*std::sqrt(1. - eps)*pol1.
z()/std::sqrt(2.*X);
146 PHI33add = -(X*X - 2.*X + 4.*pol0.
x() + 5.)*(1. -
eps)*pol1.
x()/(4.*
X);
204 G4double xs_0 = ((k0 - 2.)*k0 -2.)*
sqr(k1)*std::log(k1) + 2.*k0*(k0*(k0 + 1.)*(k0 + 8.) + 2.);
205 G4double xs_pol = (k0 + 1.)*
sqr(k1)*std::log(k1) - 2.*k0*(5.*
sqr(k0) + 4.*k0 + 1.);
207 return pre*(xs_0/k0 + pol0.
p3()*pol1.
z()*xs_pol);
CLHEP::Hep3Vector G4ThreeVector
virtual void Initialize(G4double eps, G4double X, G4double phi, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0) override
static constexpr double hbarc
G4StokesVector GetPol2() override
static constexpr double classic_electr_radius
G4double XSection(const G4StokesVector &pol2, const G4StokesVector &pol3) override
G4PolarizedComptonCrossSection()
G4StokesVector GetPol3() override
virtual ~G4PolarizedComptonCrossSection()
static const G4StokesVector ZERO
G4double TotalXSection(G4double xmin, G4double xmax, G4double y, const G4StokesVector &pol0, const G4StokesVector &pol1) override
static constexpr double pi
static const G4double eps
void SetYmin(G4double ymin)
void DefineCoefficients(const G4StokesVector &pol0, const G4StokesVector &pol1)