G4PolarizedAnnihilationXS Class Reference

#include <G4PolarizedAnnihilationXS.hh>

Inheritance diagram for G4PolarizedAnnihilationXS:

Public Member Functions
	G4PolarizedAnnihilationXS ()
virtual	~G4PolarizedAnnihilationXS () override
virtual void	Initialize (G4double eps, G4double gamma, G4double phi, const G4StokesVector &p0, const G4StokesVector &p1, G4int flag=0) override
G4double	DiceEpsilon ()
virtual G4double	XSection (const G4StokesVector &pol2, const G4StokesVector &pol3) override
virtual G4double	TotalXSection (G4double xmin, G4double xmax, G4double y, const G4StokesVector &pol0, const G4StokesVector &pol1) override
virtual G4StokesVector	GetPol2 () override
virtual G4StokesVector	GetPol3 () override
virtual G4double	GetXmin (G4double y) override
virtual G4double	GetXmax (G4double y) override
G4double	getVar (G4int)
void	getCoeff ()
G4PolarizedAnnihilationXS &	operator= (const G4PolarizedAnnihilationXS &right)=delete
	G4PolarizedAnnihilationXS (const G4PolarizedAnnihilationXS &)=delete
Public Member Functions inherited from G4VPolarizedXS
	G4VPolarizedXS ()
virtual	~G4VPolarizedXS ()
G4double	GetYmin ()
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 47 of file G4PolarizedAnnihilationXS.hh.

Constructor & Destructor Documentation

G4PolarizedAnnihilationXS() [1/2]

G4PolarizedAnnihilationXS::G4PolarizedAnnihilationXS

(

)

Definition at line 39 of file G4PolarizedAnnihilationXS.cc.

  : polxx(0.)
  , polyy(0.)
  , polzz(0.)
  , polxz(0.)
  , polzx(0.)
  , polxy(0.)
  , polyx(0.)
  , polyz(0.)
  , polzy(0.)
  , fPhi0(0.)
{
  fPhi2  = G4ThreeVector(0., 0., 0.);
  fPhi3  = G4ThreeVector(0., 0., 0.);
  fDice  = 0.;
  fPolXS = 0.;
  fUnpXS = 0.;
  ISPxx = ISPyy = ISPzz = ISPnd = 0.;
}

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

~G4PolarizedAnnihilationXS()

G4PolarizedAnnihilationXS::~G4PolarizedAnnihilationXS ( )

overridevirtual

Definition at line 60 of file G4PolarizedAnnihilationXS.cc.

{}

G4PolarizedAnnihilationXS() [2/2]

G4PolarizedAnnihilationXS::G4PolarizedAnnihilationXS ( const G4PolarizedAnnihilationXS & )

delete

Member Function Documentation

DiceEpsilon()

G4double G4PolarizedAnnihilationXS::DiceEpsilon

(

)

Definition at line 294 of file G4PolarizedAnnihilationXS.cc.

{ return fDice; }

getCoeff()

void G4PolarizedAnnihilationXS::getCoeff

(

)

GetPol2()

G4StokesVector G4PolarizedAnnihilationXS::GetPol2 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 251 of file G4PolarizedAnnihilationXS.cc.

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

GetPol3()

G4StokesVector G4PolarizedAnnihilationXS::GetPol3 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 258 of file G4PolarizedAnnihilationXS.cc.

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

getVar()

G4double G4PolarizedAnnihilationXS::getVar

(

G4int

choice

)

Definition at line 297 of file G4PolarizedAnnihilationXS.cc.

{
  if(choice == -1)
    return fPolXS / fUnpXS;
  if(choice == 0)
    return fUnpXS;
  if(choice == 1)
    return ISPxx;
  if(choice == 2)
    return ISPyy;
  if(choice == 3)
    return ISPzz;
  if(choice == 4)
    return ISPnd;
  return 0;
}

GetXmax()

G4double G4PolarizedAnnihilationXS::GetXmax ( G4double y )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 288 of file G4PolarizedAnnihilationXS.cc.

{
  return 0.5 * (1. + std::sqrt((y - 1.) / (y + 1.)));
}

GetXmin()

G4double G4PolarizedAnnihilationXS::GetXmin ( G4double y )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 282 of file G4PolarizedAnnihilationXS.cc.

{
  return 0.5 * (1. - std::sqrt((y - 1.) / (y + 1.)));
}

Initialize()

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

overridevirtual

Implements G4VPolarizedXS.

Definition at line 77 of file G4PolarizedAnnihilationXS.cc.

{
  G4double diffXSFactor = re2 / (gam - 1.);
  DefineCoefficients(pol0, pol1);
 
  // prepare dicing
  fDice = 0.;
  G4double symmXS =
    0.125 * ((-1. / sqr(gam + 1.)) / sqr(eps) +
             ((sqr(gam) + 4. * gam - 1.) / sqr(gam + 1.)) / eps - 1.);
 
  G4ThreeVector epsVector(1. / sqr(eps), 1. / eps, 1.);
  G4ThreeVector oneEpsVector(1. / sqr(1. - eps), 1. / (1. - eps), 1.);
  G4ThreeVector sumEpsVector(epsVector + oneEpsVector);
  G4ThreeVector difEpsVector(epsVector - oneEpsVector);
  G4ThreeVector calcVector(0., 0., 0.);
 
  // temporary variables
  G4double helpVar2 = 0., helpVar1 = 0.;
 
  // unpolarised contribution
  helpVar1   = (gam * gam + 4. * gam + 1.) / sqr(gam + 1.);
  helpVar2   = -1. / sqr(gam + 1.);
  calcVector = G4ThreeVector(helpVar2, helpVar1, -1.);
  fUnpXS     = 0.125 * calcVector * sumEpsVector;
 
  // initial particles polarised contribution
  helpVar2   = 1. / sqr(gam + 1.);
  helpVar1   = -(gam * gam + 4. * gam + 1.) / sqr(gam + 1.);
  calcVector = G4ThreeVector(helpVar2, helpVar1, 0.5 * (gam + 3.));
  ISPxx      = 0.25 * (calcVector * sumEpsVector) / (gam - 1.);
 
  helpVar1   = 1. / sqr(gam + 1.);
  calcVector = G4ThreeVector(-helpVar1, 2. * gam * helpVar1, -1.);
  ISPyy      = 0.125 * calcVector * sumEpsVector;
 
  helpVar1   = 1. / (gam - 1.);
  helpVar2   = 1. / sqr(gam + 1.);
  calcVector = G4ThreeVector(
    -(gam * gam + 1.) * helpVar2,
    (gam * gam * (gam + 1.) + 7. * gam + 3.) * helpVar2, -(gam + 3.));
  ISPzz = 0.125 * helpVar1 * (calcVector * sumEpsVector);
 
  helpVar1   = std::sqrt(std::fabs(eps * (1. - eps) * 2. * (gam + 1.) - 1.));
  calcVector = G4ThreeVector(-1. / (gam * gam - 1.), 2. / (gam - 1.), 0.);
  ISPnd      = 0.125 * (calcVector * difEpsVector) * helpVar1;
 
  fPolXS = ISPxx * polxx + ISPyy * polyy + ISPzz * polzz;
  fPolXS += ISPnd * (polzx + polxz);
  fPhi0 = fUnpXS + fPolXS;
  fDice = symmXS;
 
  if(polzz != 0.)
  {
    fDice *= (1. + (polzz * ISPzz / fUnpXS));
    if(fDice < 0.)
      fDice = 0.;
  }
  // prepare final state coefficients
  if(flag == 2)
  {
    // circular polarisation
    G4double circ1 = 0., circ2 = 0., circ3 = 0.;
    helpVar1 = 8. * sqr(1. - eps) * sqr(eps) * (gam - 1.) * sqr(gam + 1.) /
               std::sqrt(gam * gam - 1.);
    helpVar2 = sqr(gam + 1.) * sqr(eps) * (-2. * eps + 3.) -
               (gam * gam + 3. * gam + 2.) * eps;
    circ1 = helpVar2 + gam;
    circ1 /= helpVar1;
    circ2 = helpVar2 + 1.;
    circ2 /= helpVar1;
    helpVar1 = std::sqrt(std::fabs(eps * (1. - eps) * 2. * (gam + 1.) - 1.));
    helpVar1 /= std::sqrt(gam * gam - 1.);
    calcVector = G4ThreeVector(1., -2. * gam, 0.);
    circ3      = 0.125 * (calcVector * sumEpsVector) / (gam + 1.);
    circ3 *= helpVar1;
 
    fPhi2.setZ(circ2 * pol1.z() + circ1 * pol0.z() +
               circ3 * (pol1.x() + pol0.x()));
    fPhi3.setZ(-circ1 * pol1.z() - circ2 * pol0.z() -
               circ3 * (pol1.x() + pol0.x()));
 
    // common to both linear polarisation
    calcVector          = G4ThreeVector(-1., 2. * gam, 0.);
    G4double linearZero = 0.125 * (calcVector * sumEpsVector) / sqr(gam + 1.);
 
    //        Linear Polarisation #1
    helpVar1 = std::sqrt(std::fabs(2. * (gam + 1.) * (1. - eps) * eps - 1.)) /
               ((gam + 1.) * eps * (1. - eps));
    helpVar2 = helpVar1 * helpVar1;
 
    // photon 1
    G4double diagContrib = 0.125 * helpVar2 * (polxx + polyy - polzz);
    G4double nonDiagContrib =
      0.125 * helpVar1 * (-polxz / (1. - eps) + polzx / eps);
 
    fPhi2.setX(linearZero + diagContrib + nonDiagContrib);
 
    // photon 2
    nonDiagContrib = 0.125 * helpVar1 * (polxz / eps - polzx / (1. - eps));
 
    fPhi3.setX(linearZero + diagContrib + nonDiagContrib);
 
    //        Linear Polarisation #2
    helpVar1 =
      std::sqrt(gam * gam - 1.) * (2. * (gam + 1.) * eps * (1. - eps) - 1.);
    helpVar1 /= 8. * sqr(1. - eps) * sqr(eps) * sqr(gam + 1.) * (gam - 1.);
    helpVar2 = std::sqrt((gam * gam - 1.) *
                         std::fabs(2. * (gam + 1.) * eps * (1. - eps) - 1.));
    helpVar2 /= 8. * sqr(1. - eps) * sqr(eps) * sqr(gam + 1.) * (gam - 1.);
 
    G4double contrib21 = (-polxy + polyx) * helpVar1;
    G4double contrib32 =
      -(eps * (gam + 1.) - 1.) * polyz + (eps * (gam + 1.) - gam) * polzy;
 
    contrib32 *= helpVar2;
    fPhi2.setY(contrib21 + contrib32);
 
    contrib32 =
      -(eps * (gam + 1.) - gam) * polyz + (eps * (gam + 1.) - 1.) * polzy;
    contrib32 *= helpVar2;
    fPhi3.setY(contrib21 + contrib32);
  }
  fPhi0 *= diffXSFactor;
  fPhi2 *= diffXSFactor;
  fPhi3 *= diffXSFactor;
}

operator=()

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

delete

TotalXSection()

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

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 220 of file G4PolarizedAnnihilationXS.cc.

{
  G4double totalXSFactor = pi * re2 / (gam + 1.);  // atomic number ignored
  DefineCoefficients(pol0, pol1);
 
  G4double xs = 0.;
 
  G4double gam2     = gam * gam;
  G4double sqrtgam1 = std::sqrt(gam2 - 1.);
  G4double logMEM   = std::log(gam + sqrtgam1);
  G4double unpME    = (gam * (gam + 4.) + 1.) * logMEM;
  unpME += -(gam + 3.) * sqrtgam1;
  unpME /= 4. * (gam2 - 1.);
  G4double longPart = (3 + gam * (gam * (gam + 1.) + 7.)) * logMEM;
  longPart += -(5. + gam * (3 * gam + 4.)) * sqrtgam1;
  longPart /= 4. * sqr(gam - 1.) * (gam + 1.);
  G4double tranPart = -(5 * gam + 1.) * logMEM;
  tranPart += (gam + 5.) * sqrtgam1;
  tranPart /= 4. * sqr(gam - 1.) * (gam + 1.);
 
  xs += unpME;
  xs += polzz * longPart;
  xs += (polxx + polyy) * tranPart;
 
  return xs * totalXSFactor;
}

XSection()

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

overridevirtual

Implements G4VPolarizedXS.

Definition at line 211 of file G4PolarizedAnnihilationXS.cc.

{
  G4double xs = fPhi0 + pol2 * fPhi2 + pol3 * fPhi3;
  return xs;
}

---

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

• G4PolarizedAnnihilationXS.hh
• G4PolarizedAnnihilationXS.cc