66 G4double re2 = classic_electr_radius * classic_electr_radius;
69 G4double gmo2 = (gamma - 1.)*(gamma - 1.);
71 G4double pref = gamma2*re2/(gmo2*(gamma + 1.0));
80 if (flag==0) polarized=
false;
84 phi0+= ((1. - 2.*gamma)/gamma2)*(1./e + 1./(1.-e));
85 phi0+= 1./(e*e) + 1./((1. - e)*(1. - e));
90 if (flag<=1) usephi=0.;
94 G4double xx = (gamma - f*e*gmo*(3 + gamma))/(4*f*e*gamma2);
95 G4double yy = (-1 + f*e*gmo2 + 2*gamma)/(4*f*e*gamma2);
96 G4double zz = (-(e*gmo*(3 + gamma)) + e2*gmo*(3 + gamma) +
97 gamma*(-1 + 2*gamma))/(4*f*e*gamma2);
99 phi0 += xx*pol0.x()*pol1.x() + yy*pol0.y()*pol1.y() + zz*pol0.z()*pol1.z();
104 G4double xz = -((-1 + 2*e)*gmo)/(2*sqrttwo*gamma2*
105 std::sqrt(-((f*e)/gpo)));
108 G4double zx = -((-1 + 2*e)*gmo)/(2*sqrttwo*gamma2*
109 std::sqrt(-((f*e)/gpo)));
111 phi0+=yx*pol0.y()*pol1.x() + xy*pol0.x()*pol1.y();
112 phi0+=zx*pol0.z()*pol1.x() + xz*pol0.x()*pol1.z();
113 phi0+=zy*pol0.z()*pol1.y() + yz*pol0.y()*pol1.z();
127 G4double xxP1K1 = (std::sqrt(gpo/(1 + e2*gmo + gamma - 2*e*gamma))*
128 (gamma - e*gpo))/(4*e2*gamma);
130 G4double xzP1K1 = (-1 + 2*e*gamma)/(2*sqrttwo*f*gamma*
131 std::sqrt(e*e2*(1 + e + gamma - e*gamma)));
133 G4double yyP1K1 = (-gamma2 + e*(-1 + gamma*(2 + gamma)))/(4*f*e2*gamma2);
135 G4double zxP1K1 = (1 + 2*e2*gmo - 2*e*gamma)/(2*sqrttwo*f*e*gamma*
136 std::sqrt(e*(1 + e + gamma - e*gamma)));
138 G4double zzP1K1 = (-gamma + e*(1 - 2*e*gmo + gamma))/(4*f*e2*gamma*
139 std::sqrt(1 - (2*e)/(f*gpo)));
140 phi2[0] += xxP1K1*pol0.x() + xyP1K1*pol0.y() + xzP1K1*pol0.z();
141 phi2[1] += yxP1K1*pol0.x() + yyP1K1*pol0.y() + yzP1K1*pol0.z();
142 phi2[2] += zxP1K1*pol0.x() + zyP1K1*pol0.y() + zzP1K1*pol0.z();
146 G4double xxP1K2 = ((1 + e*(-3 + gamma))*std::sqrt(gpo/(1 + e2*gmo + gamma -
147 2*e*gamma)))/(4*f*e*gamma);
149 G4double xzP1K2 = (-2 + 2*e + gamma)/(2*sqrttwo*f2*gamma*
150 std::sqrt(e*(1 + e + gamma - e*gamma)));
152 G4double yyP1K2 = (1 - 2*gamma + e*(-1 + gamma*(2 + gamma)))/(4*f2*e*gamma2);
154 G4double zxP1K2 = (2*e*(1 + e*gmo - 2*gamma) + gamma)/(2*sqrttwo*f2*gamma*
155 std::sqrt(e*(1 + e + gamma - e*gamma)));
157 G4double zzP1K2 = (1 - 2*gamma + e*(-1 - 2*e*gmo + 3*gamma))/
158 (4*f2*e*gamma*std::sqrt(1 - (2*e)/(f*gpo)));
159 phi2[0] += xxP1K2*pol1.x() + xyP1K2*pol1.y() + xzP1K2*pol1.z();
160 phi2[1] += yxP1K2*pol1.x() + yyP1K2*pol1.y() + yzP1K2*pol1.z();
161 phi2[2] += zxP1K2*pol1.x() + zyP1K2*pol1.y() + zzP1K2*pol1.z();
171 G4double xxP2K1 = (-1 + e + e*gamma)/(4*f2*gamma*
172 std::sqrt((e*(2 + e*gmo))/gpo));
174 G4double xzP2K1 = -((1 + 2*f*gamma)*std::sqrt(f/(-2 + e - e*gamma)))/
175 (2*sqrttwo*f2*e*gamma);
177 G4double yyP2K1 = (1 - 2*gamma + e*(-1 + gamma*(2 + gamma)))/(4*f2*e*gamma2);
179 G4double zxP2K1 = (1 + 2*e*(-2 + e + gamma - e*gamma))/(2*sqrttwo*f*e*
180 std::sqrt(-(f*(2 + e*gmo)))*gamma);
182 G4double zzP2K1 = (std::sqrt((e*gpo)/(2 + e*gmo))*
183 (-3 + e*(5 + 2*e*gmo - 3*gamma) + 2*gamma))/(4*f2*e*gamma);
185 phi3[0] += xxP2K1*pol0.x() + xyP2K1*pol0.y() + xzP2K1*pol0.z();
186 phi3[1] += yxP2K1*pol0.x() + yyP2K1*pol0.y() + yzP2K1*pol0.z();
187 phi3[2] += zxP2K1*pol0.x() + zyP2K1*pol0.y() + zzP2K1*pol0.z();
192 G4double xxP2K2 = (-2 - e*(-3 + gamma) + gamma)/
193 (4*f*e*gamma* std::sqrt((e*(2 + e*gmo))/gpo));
195 G4double xzP2K2 = ((-2*e + gamma)*std::sqrt(f/(-2 + e - e*gamma)))/
196 (2*sqrttwo*f*e2*gamma);
198 G4double yyP2K2 = (-gamma2 + e*(-1 + gamma*(2 + gamma)))/(4*f*e2*gamma2);
200 G4double zxP2K2 = (gamma + 2*e*(-1 + e - e*gamma))/
201 (2*sqrttwo*e2* std::sqrt(-(f*(2 + e*gmo)))*gamma);
203 G4double zzP2K2 = (std::sqrt((e*gpo)/(2 + e*gmo))*
204 (-2 + e*(3 + 2*e*gmo - gamma) + gamma))/(4*f*e2*gamma);
205 phi3[0] += xxP2K2*pol1.x() + xyP2K2*pol1.y() + xzP2K2*pol1.z();
206 phi3[1] += yxP2K2*pol1.x() + yyP2K2*pol1.y() + yzP2K2*pol1.z();
207 phi3[2] += zxP2K2*pol1.x() + zyP2K2*pol1.y() + zzP2K2*pol1.z();
236 if (xmax != 1./2.)
G4cout<<
" warning xmax expected to be 1/2 but is "<<xmax<<
G4endl;
239 G4double re2 = classic_electr_radius * classic_electr_radius;
241 G4double gmo2 = (gamma - 1.)*(gamma - 1.);
242 G4double logMEM = std::log(1./x - 1.);
246 sigma0 += (gmo2/gamma2)*(0.5 - x);
247 sigma0 += ((1. - 2.*gamma)/gamma2)*logMEM;
248 sigma0 += 1./x - 1./(1. -
x);
251 sigma2 += ((gamma2 + 2.*gamma - 3.)/gamma2)*(0.5 - x);
252 sigma2 += (1./gamma - 2.)*logMEM;
255 sigma3 += (2.*(1. - gamma)/gamma2)*(0.5 -
x);
256 sigma3 += (1. - 3.*gamma)/(2.*gamma2)*logMEM;
258 xs+=pref*(sigma0 + sigma2*pol0.z()*pol1.z() + sigma3*(pol0.x()*pol1.x()+pol0.y()*pol1.y()));
G4double XSection(const G4StokesVector &pol2, const G4StokesVector &pol3)
CLHEP::Hep3Vector G4ThreeVector
void Initialize(G4double x, G4double y, G4double phi, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0)
G4PolarizedMollerCrossSection()
G4GLOB_DLL std::ostream G4cout
static const double twopi
void SetXmax(G4double xmax)
const G4double x[NPOINTSGL]
G4double TotalXSection(G4double xmin, G4double xmax, G4double y, const G4StokesVector &pol0, const G4StokesVector &pol1)
virtual ~G4PolarizedMollerCrossSection()