G4PolarizedComptonXS Class Reference

#include <G4PolarizedComptonXS.hh>

Inheritance diagram for G4PolarizedComptonXS:

Public Member Functions
	G4PolarizedComptonXS ()
	~G4PolarizedComptonXS () override
void	Initialize (G4double eps, G4double X, G4double phi, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0) override
G4double	XSection (const G4StokesVector &pol2, const G4StokesVector &pol3) override
G4double	TotalXSection (G4double xmin, G4double xmax, G4double y, const G4StokesVector &pol0, const G4StokesVector &pol1) override
G4StokesVector	GetPol2 () override
G4StokesVector	GetPol3 () override
G4PolarizedComptonXS &	operator= (const G4PolarizedComptonXS &right)=delete
	G4PolarizedComptonXS (const G4PolarizedComptonXS &)=delete
Public Member Functions inherited from G4VPolarizedXS
	G4VPolarizedXS ()
virtual	~G4VPolarizedXS ()
G4double	GetYmin ()
virtual G4double	GetXmin (G4double y)
virtual G4double	GetXmax (G4double y)
void	SetMaterial (G4double A, G4double Z, G4double coul)
G4VPolarizedXS &	operator= (const G4VPolarizedXS &right)=delete
	G4VPolarizedXS (const G4VPolarizedXS &)=delete

Additional Inherited Members
Protected Member Functions inherited from G4VPolarizedXS
void	SetXmin (G4double xmin)
void	SetXmax (G4double xmax)
void	SetYmin (G4double ymin)
Protected Attributes inherited from G4VPolarizedXS
G4double	fXmin
G4double	fXmax
G4double	fYmin
G4double	fA
G4double	fZ
G4double	fCoul

Detailed Description

Definition at line 45 of file G4PolarizedComptonXS.hh.

Constructor & Destructor Documentation

G4PolarizedComptonXS() [1/2]

G4PolarizedComptonXS::G4PolarizedComptonXS

(

)

Definition at line 41 of file G4PolarizedComptonXS.cc.

{
  SetYmin(0.);
 
  fPhi0  = 0.;
  fPolXS = 0.;
  fUnpXS = 0.;
  fPhi2  = G4ThreeVector(0., 0., 0.);
  fPhi3  = G4ThreeVector(0., 0., 0.);
  polxx = polyy = polzz = polxz = polzx = polyz = polzy = polxy = polyx = 0.;
}

Referenced by G4PolarizedComptonXS(), and operator=().

~G4PolarizedComptonXS()

G4PolarizedComptonXS::~G4PolarizedComptonXS ( )

override

Definition at line 54 of file G4PolarizedComptonXS.cc.

{}

G4PolarizedComptonXS() [2/2]

G4PolarizedComptonXS::G4PolarizedComptonXS ( const G4PolarizedComptonXS & )

delete

Member Function Documentation

GetPol2()

G4StokesVector G4PolarizedComptonXS::GetPol2 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 214 of file G4PolarizedComptonXS.cc.

{
  // Note, mean polarization can not contain correlation effects.
  return G4StokesVector(1. / fPhi0 * fPhi2);
}

GetPol3()

G4StokesVector G4PolarizedComptonXS::GetPol3 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 221 of file G4PolarizedComptonXS.cc.

{
  // Note, mean polarization can not contain correlation effects.
  return G4StokesVector(1. / fPhi0 * fPhi3);
}

Initialize()

void G4PolarizedComptonXS::Initialize ( G4double eps, G4double X, G4double phi, const G4StokesVector & p0, const G4StokesVector & p1, G4int flag = 0 )

overridevirtual

Implements G4VPolarizedXS.

Definition at line 57 of file G4PolarizedComptonXS.cc.

{
  G4double cosT = 1. - (1. / eps - 1.) / X;
  if(cosT > 1. + 1.e-8)
    cosT = 1.;
  else if(cosT < -1. - 1.e-8)
    cosT = -1.;
  G4double cosT2 = cosT * cosT;
  G4double cosT3 = cosT2 * cosT;
  G4double sinT2 = 1. - cosT2;
  if(sinT2 > 1. + 1.e-8)
    sinT2 = 1.;
  else if(sinT2 < 0.)
    sinT2 = 0.;
  G4double sinT  = std::sqrt(sinT2);
  G4double cos2T = 2. * cosT2 - 1.;
  G4double sin2T = 2. * sinT * cosT;
  G4double eps2  = sqr(eps);
  DefineCoefficients(pol0, pol1);
  G4double diffXSFactor = re2 / (4. * X);
 
  // unpolarized Cross Section
  fUnpXS = (eps2 + 1. - eps * sinT2) / (2. * eps);
  // initial polarization dependence
  fPolXS = -sinT2 * pol0.x() + (1. - eps) * sinT * polzx +
           ((eps2 - 1.) / eps) * cosT * polzz;
  fPolXS *= 0.5;
 
  fPhi0 = fUnpXS + fPolXS;
 
  if(flag == 2)
  {
    // polarization of outgoing photon
    G4double phi21 = -sinT2 + 0.5 * (cos2T + 3.) * pol0.x() -
                     ((1. - eps) / eps) * sinT * polzx;
    phi21 *= 0.5;
    G4double phi22 = cosT * pol0.y() + ((1. - eps) / (2. * eps)) * sinT * polzy;
    G4double phi23 = ((eps2 + 1.) / eps) * cosT * pol0.z() -
                     ((1. - eps) / eps) * (eps * cosT2 + 1.) * pol1.z();
    phi23 += 0.5 * (1. - eps) * sin2T * pol1.x();
    phi23 += (eps - 1.) * (-sinT2 * polxz + sinT * polyy - 0.5 * sin2T * polxx);
    phi23 *= 0.5;
 
    fPhi2 = G4ThreeVector(phi21, phi22, phi23);
 
    // polarization of outgoing electron
    G4double phi32 = -sinT2 * polxy + ((1. - eps) / eps) * sinT * polyz +
                     0.5 * (cos2T + 3.) * pol1.y();
    phi32 *= 0.5;
 
    G4double phi31    = 0.;
    G4double phi31add = 0.;
    G4double phi33    = 0.;
    G4double phi33add = 0.;
 
    if((1. - eps) > 1.e-12)
    {
      G4double helpVar = std::sqrt(eps2 - 2. * cosT * eps + 1.);
 
      phi31 = (1. - eps) * (1. + cosT) * sinT * pol0.z();
      phi31 +=
        (-eps * cosT3 + eps * cosT2 + (eps - 2.) * cosT + eps) * pol1.x();
      phi31 += -(eps * cosT2 - eps * cosT + cosT + 1.) * sinT * pol1.z();
      phi31 /= 2. * helpVar;
 
      phi31add = -eps * sqr(1. - cosT) * (1. + cosT) * polxx;
      phi31add += (1. - eps) * sinT2 * polyy;
      phi31add += -(-eps2 + cosT * (cosT * eps - eps + 1.) * eps + eps - 1.) *
                  sinT * polxz / eps;
      phi31add /= 2. * helpVar;
 
      phi33 = ((1. - eps) / eps) *
              (-eps * cosT2 + eps * (eps + 1.) * cosT - 1.) * pol0.z();
      phi33 += -(eps * cosT2 + (1. - eps) * eps * cosT + 1.) * sinT * pol1.x();
      phi33 +=
        -(-eps2 * cosT3 + eps * (eps2 - eps + 1.) * cosT2 - cosT + eps2) *
        pol1.z() / eps;
      phi33 /= -2. * helpVar;
 
      phi33add = (eps * (eps - cosT - 1.) * cosT + 1.) * sinT * polxx;
      phi33add += -(-eps2 + cosT * eps + eps - 1.) * sinT2 * polxz;
      phi33add += (eps - 1.) * (cosT - eps) * sinT * polyy;
      phi33add /= -2. * helpVar;
    }
    else
    {
      phi31 = -pol1.z() -
              (X - 1.) * std::sqrt(1. - eps) * pol1.x() / std::sqrt(2. * X);
      phi31add = -(-X * X * pol1.z() - 2. * X * (2. * pol0.z() - pol1.z()) -
                   (4. * pol0.x() + 5.) * pol1.z()) *
                 (1. - eps) / (4. * X);
 
      phi33 = pol1.x() -
              (X - 1.) * std::sqrt(1. - eps) * pol1.z() / std::sqrt(2. * X);
      phi33add = -(X * X - 2. * X + 4. * pol0.x() + 5.) * (1. - eps) *
                 pol1.x() / (4. * X);
    }
    fPhi3 = G4ThreeVector(phi31 + phi31add, phi32, phi33 + phi33add);
  }
  fUnpXS *= diffXSFactor;
  fPolXS *= diffXSFactor;
  fPhi0 *= diffXSFactor;
  fPhi2 *= diffXSFactor;
  fPhi3 *= diffXSFactor;
}

operator=()

G4PolarizedComptonXS & G4PolarizedComptonXS::operator= ( const G4PolarizedComptonXS & right )

delete

TotalXSection()

G4double G4PolarizedComptonXS::TotalXSection ( G4double xmin, G4double xmax, G4double y, const G4StokesVector & pol0, const G4StokesVector & pol1 )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 193 of file G4PolarizedComptonXS.cc.

{
  G4double k1 = 1. + 2. * k0;
 
  G4double unit = fZ * CLHEP::pi * CLHEP::classic_electr_radius *
                  CLHEP::classic_electr_radius;
 
  G4double pre = unit / (sqr(k0) * sqr(1. + 2. * k0));
 
  G4double xs_0 = ((k0 - 2.) * k0 - 2.) * sqr(k1) * std::log(k1) +
                  2. * k0 * (k0 * (k0 + 1.) * (k0 + 8.) + 2.);
  G4double xs_pol = (k0 + 1.) * sqr(k1) * std::log(k1) -
                    2. * k0 * (5. * sqr(k0) + 4. * k0 + 1.);
 
  return pre * (xs_0 / k0 + pol0.p3() * pol1.z() * xs_pol);
}

XSection()

G4double G4PolarizedComptonXS::XSection ( const G4StokesVector & pol2, const G4StokesVector & pol3 )

overridevirtual

Implements G4VPolarizedXS.

Definition at line 166 of file G4PolarizedComptonXS.cc.

{
  G4bool gammaPol2    = !(pol2 == G4StokesVector::ZERO);
  G4bool electronPol3 = !(pol3 == G4StokesVector::ZERO);
 
  G4double phi = 0.;
  // polarization independent part
  phi += fPhi0;
 
  if(gammaPol2)
  {
    // part depending on the polarization of the final photon
    phi += fPhi2 * pol2;
  }
 
  if(electronPol3)
  {
    // part depending on the polarization of the final electron
    phi += fPhi3 * pol3;
  }
 
  // return cross section.
  return phi;
}

---

The documentation for this class was generated from the following files:

• G4PolarizedComptonXS.hh
• G4PolarizedComptonXS.cc