69 G4double gmo2 = (gamma - 1.)*(gamma - 1.);
71 G4double pref = gamma2*re2/(gmo2*(gamma + 1.0));
78 G4bool polarized=(!pol0.IsZero())||(!pol1.IsZero());
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.;
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();
126 if (!pol0.IsZero()) {
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();
145 if (!pol1.IsZero()) {
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();
168 if (!pol0.IsZero()) {
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();
190 if (!pol1.IsZero()) {
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();
CLHEP::Hep3Vector G4ThreeVector
int classic_electr_radius
Public Member Functions inherited from G4VPolarizedCrossSection