52 : gammaPol2(false), electronPol3(false)
61 phi0 = 0.; polXS = 0.; unpXS = 0.;
64 polxx = polyy = polzz = polxz = polzx = polyz = polzy = polxy = polyx = 0.;
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.;
91 DefineCoefficients(pol0,pol1);
92 diffXSFactor = re2/(4.*
X);
95 unpXS = (eps2 + 1. - eps*sinT2)/(2.*eps);
97 polXS = -sinT2*pol0.
x() + (1. - eps)*sinT*polzx + ((eps2 - 1.)/eps)*cosT*polzz;
100 phi0 = unpXS + polXS;
104 G4double PHI21 = -sinT2 + 0.5*(cos2T + 3.)*pol0.
x() - ((1. - eps)/eps)*sinT*polzx;
106 G4double PHI22 = cosT*pol0.
y() + ((1. - eps)/(2.*eps))*sinT*polzy;
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();
109 PHI23 += (eps - 1.)*(-sinT2*polxz + sinT*polyy - 0.5*sin2T*polxx);
114 G4double PHI32 = -sinT2*polxy + ((1. - eps)/eps)*sinT*polyz + 0.5*(cos2T + 3.)*pol1.
y();
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;
128 PHI31add += (1. - eps)*sinT2*polyy;
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);
148 phi3 =
G4ThreeVector(PHI31 + PHI31add, PHI32, PHI33 + PHI33add);
151 unpXS *= diffXSFactor;
152 polXS *= diffXSFactor;
153 phi0 *= diffXSFactor;
154 phi2 *= diffXSFactor;
155 phi3 *= diffXSFactor;
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);
218 return 1./phi0 * phi2;
227 return 1./phi0 * phi3;
233 void G4PolarizedComptonCrossSection::DefineCoefficients(
const G4StokesVector & pol0,
236 polxx=pol0.
x()*pol1.
x();
237 polyy=pol0.
y()*pol1.
y();
238 polzz=pol0.
z()*pol1.
z();
240 polxz=pol0.
x()*pol1.
z();
241 polzx=pol0.
z()*pol1.
x();
243 polyz=pol0.
y()*pol1.
z();
244 polzy=pol0.
z()*pol1.
y();
246 polxy=pol0.
x()*pol1.
y();
247 polyx=pol0.
y()*pol1.
x();
G4double XSection(const G4StokesVector &pol2, const G4StokesVector &pol3)
CLHEP::Hep3Vector G4ThreeVector
void SetYmin(G4double ymin)
G4double TotalXSection(G4double xmin, G4double xmax, G4double y, const G4StokesVector &pol0, const G4StokesVector &pol1)
virtual void Initialize(G4double eps, G4double X, G4double phi, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0)
virtual ~G4PolarizedComptonCrossSection()
G4PolarizedComptonCrossSection()
static const G4StokesVector ZERO