G4PolarizedGammaConversionXS Class Reference

#include <G4PolarizedGammaConversionXS.hh>

Inheritance diagram for G4PolarizedGammaConversionXS:

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

Protected Member Functions
void	SetXmax (G4double xmax)
void	SetXmin (G4double xmin)
void	SetYmin (G4double ymin)

Protected Attributes
G4double	fA
G4double	fCoul
G4double	fXmax
G4double	fXmin
G4double	fYmin
G4double	fZ

Private Attributes
G4StokesVector	fFinalElectronPolarization
G4StokesVector	fFinalPositronPolarization

Static Private Attributes
static G4double	SCRN [2][19]

Detailed Description

Definition at line 38 of file G4PolarizedGammaConversionXS.hh.

Constructor & Destructor Documentation

G4PolarizedGammaConversionXS() [1/2]

G4PolarizedGammaConversionXS::G4PolarizedGammaConversionXS

(

)

Definition at line 43 of file G4PolarizedGammaConversionXS.cc.

{
  fFinalElectronPolarization = G4StokesVector::ZERO;
  fFinalPositronPolarization = G4StokesVector::ZERO;
}

References fFinalElectronPolarization, fFinalPositronPolarization, and G4StokesVector::ZERO.

~G4PolarizedGammaConversionXS()

G4PolarizedGammaConversionXS::~G4PolarizedGammaConversionXS ( )

override

Definition at line 49 of file G4PolarizedGammaConversionXS.cc.

{}

G4PolarizedGammaConversionXS() [2/2]

G4PolarizedGammaConversionXS::G4PolarizedGammaConversionXS ( const G4PolarizedGammaConversionXS &  )

delete

Member Function Documentation

GetPol2()

G4StokesVector G4PolarizedGammaConversionXS::GetPol2 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 176 of file G4PolarizedGammaConversionXS.cc.

{
  // electron/positron
  return fFinalElectronPolarization;
}

References fFinalElectronPolarization.

Referenced by G4PolarizedGammaConversionModel::SampleSecondaries().

GetPol3()

G4StokesVector G4PolarizedGammaConversionXS::GetPol3 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 182 of file G4PolarizedGammaConversionXS.cc.

{
  // photon
  return fFinalPositronPolarization;
}

References fFinalPositronPolarization.

Referenced by G4PolarizedGammaConversionModel::SampleSecondaries().

GetXmax()

G4double G4VPolarizedXS::GetXmax ( G4double  y )

virtualinherited

Reimplemented in G4PolarizedAnnihilationXS.

Definition at line 84 of file G4VPolarizedXS.cc.

{ return fXmax; }

References G4VPolarizedXS::fXmax.

GetXmin()

G4double G4VPolarizedXS::GetXmin ( G4double  y )

virtualinherited

Reimplemented in G4PolarizedAnnihilationXS.

Definition at line 81 of file G4VPolarizedXS.cc.

{ return fXmin; }

References G4VPolarizedXS::fXmin.

GetYmin()

G4double G4VPolarizedXS::GetYmin ( )

inlineinherited

Definition at line 65 of file G4VPolarizedXS.hh.

{ return fYmin; }

References G4VPolarizedXS::fYmin.

Initialize()

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

overridevirtual

Implements G4VPolarizedXS.

Definition at line 51 of file G4PolarizedGammaConversionXS.cc.

{
  G4double aLept1E = aGammaE - aLept0E;
 
  G4double Stokes_P3 = beamPol.z();
 
  G4double Lept0E  = aLept0E / CLHEP::electron_mass_c2 + 1.;
  G4double Lept0E2 = Lept0E * Lept0E;
  G4double GammaE  = aGammaE / CLHEP::electron_mass_c2;
  G4double Lept1E  = aLept1E / CLHEP::electron_mass_c2 - 1.;
  G4double Lept1E2 = Lept1E * Lept1E;
 
  // *******  Gamma Transvers Momentum
  G4double TMom = std::sqrt(Lept0E2 - 1.) * sintheta;
  G4double u    = TMom;
  G4double u2   = u * u;
  G4double Xsi  = 1. / (1. + u2);
  G4double Xsi2 = Xsi * Xsi;
 
  G4double delta =
    12. * std::pow(fZ, 1. / 3.) * Lept0E * Lept1E * Xsi / (121. * GammaE);
  G4double GG = 0.;
 
  if(delta < 0.5)
  {
    GG = std::log(2. * Lept0E * Lept1E / GammaE) - 2. - fCoul;
  }
  else if(delta < 120.)
  {
    for(G4int j = 1; j < 19; ++j)
    {
      if(SCRN[0][j] >= delta)
      {
        GG = std::log(2. * Lept0E * Lept1E / GammaE) - 2. - fCoul -
             (SCRN[1][j - 1] + (delta - SCRN[0][j - 1]) *
                                 (SCRN[1][j] - SCRN[1][j - 1]) /
                                 (SCRN[0][j] - SCRN[0][j - 1]));
        break;
      }
    }
  }
  else
  {
    G4double alpha_sc = (111. * std::pow(fZ, -1. / 3.)) / Xsi;
    GG                = std::log(alpha_sc) - 2. - fCoul;
  }
 
  if(GG < -1.)
    GG = -1.;
 
  G4double I_Lepton = (Lept0E2 + Lept1E2) * (3 + 2 * GG) +
                      2. * Lept0E * Lept1E * (1. + 4. * u2 * Xsi2 * GG);
 
  G4double L_Lepton1 = GammaE *
                       ((Lept0E - Lept1E) * (3. + 2. * GG) +
                        2 * Lept1E * (1. + 4. * u2 * Xsi2 * GG)) /
                       I_Lepton;
 
  G4double T_Lepton1 =
    4. * GammaE * Lept1E * Xsi * u * (1. - 2. * Xsi) * GG / I_Lepton;
 
  G4double Stokes_S1 = (Stokes_P3 * T_Lepton1);
  G4double Stokes_S2 = 0.;
  G4double Stokes_S3 = (Stokes_P3 * L_Lepton1);
 
  fFinalElectronPolarization.setX(Stokes_S1);
  fFinalElectronPolarization.setY(Stokes_S2);
  fFinalElectronPolarization.setZ(Stokes_S3);
 
  if(fFinalElectronPolarization.mag2() > 1.)
  {
    G4ExceptionDescription ed;
    ed << "\t" << fFinalElectronPolarization << "\t GG\t" << GG << "\t delta\t"
       << delta << "\n";
    G4Exception("G4PolarizedGammaConversionXS::Initialize", "pol022",
                JustWarning, ed);
    fFinalElectronPolarization.setX(0.);
    fFinalElectronPolarization.setY(0.);
    fFinalElectronPolarization.setZ(Stokes_S3);
    if(Stokes_S3 > 1.)
      fFinalElectronPolarization.setZ(1.);
  }
 
  G4double L_Lepton2 = GammaE *
                       ((Lept1E - Lept0E) * (3. + 2. * GG) +
                        2 * Lept0E * (1. + 4. * u2 * Xsi2 * GG)) /
                       I_Lepton;
 
  G4double T_Lepton2 =
    4. * GammaE * Lept0E * Xsi * u * (1. - 2. * Xsi) * GG / I_Lepton;
 
  G4double Stokes_SS1 = (Stokes_P3 * T_Lepton2);
  G4double Stokes_SS2 = 0.;
  G4double Stokes_SS3 = (Stokes_P3 * L_Lepton2);
 
  fFinalPositronPolarization.SetPhoton();
 
  fFinalPositronPolarization.setX(Stokes_SS1);
  fFinalPositronPolarization.setY(Stokes_SS2);
  fFinalPositronPolarization.setZ(Stokes_SS3);
 
  if(fFinalPositronPolarization.mag2() > 1.)
  {
    G4ExceptionDescription ed;
    ed << "\t" << fFinalPositronPolarization << "\t GG\t" << GG << "\t delta\t"
       << delta << "\n";
    G4Exception("G4PolarizedGammaConversionXS::Initialize", "pol023",
                JustWarning, ed);
  }
}

References CLHEP::electron_mass_c2, G4VPolarizedXS::fCoul, fFinalElectronPolarization, fFinalPositronPolarization, G4VPolarizedXS::fZ, G4Exception(), JustWarning, CLHEP::Hep3Vector::mag2(), SCRN, G4StokesVector::SetPhoton(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), and CLHEP::Hep3Vector::z().

Referenced by G4PolarizedGammaConversionModel::SampleSecondaries().

operator=()

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

delete

SetMaterial()

void G4VPolarizedXS::SetMaterial ( G4double  A, G4double  Z, G4double  coul  )

inlineinherited

Definition at line 71 of file G4VPolarizedXS.hh.

  {
    fA    = A;
    fZ    = Z;
    fCoul = coul;
  }

References A, G4VPolarizedXS::fA, G4VPolarizedXS::fCoul, G4VPolarizedXS::fZ, and Z.

Referenced by G4PolarizedBremsstrahlungModel::SampleSecondaries(), G4PolarizedGammaConversionModel::SampleSecondaries(), and G4PolarizedPhotoElectricModel::SampleSecondaries().

SetXmax()

void G4VPolarizedXS::SetXmax ( G4double  xmax )

inlineprotectedinherited

Definition at line 84 of file G4VPolarizedXS.hh.

{ fXmax = xmax; }

References G4VPolarizedXS::fXmax.

Referenced by G4PolarizedIonisationMollerXS::G4PolarizedIonisationMollerXS(), and G4PolarizedIonisationBhabhaXS::Initialize().

SetXmin()

void G4VPolarizedXS::SetXmin ( G4double  xmin )

inlineprotectedinherited

Definition at line 83 of file G4VPolarizedXS.hh.

{ fXmin = xmin; }

References G4VPolarizedXS::fXmin.

SetYmin()

void G4VPolarizedXS::SetYmin ( G4double  ymin )

inlineprotectedinherited

Definition at line 85 of file G4VPolarizedXS.hh.

{ fYmin = ymin; }

References G4VPolarizedXS::fYmin.

Referenced by G4PolarizedComptonXS::G4PolarizedComptonXS().

TotalXSection()

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

virtualinherited

Reimplemented in G4PolarizedAnnihilationXS, G4PolarizedComptonXS, G4PolarizedIonisationBhabhaXS, and G4PolarizedIonisationMollerXS.

Definition at line 86 of file G4VPolarizedXS.cc.

{
  G4ExceptionDescription ed;
  ed << "WARNING virtual function G4VPolarizedXS::TotalXSection() "
        "called.\n";
  G4Exception("G4VPolarizedXS::TotalXSection", "pol032", FatalException, ed);
  return 0.;
}

References FatalException, and G4Exception().

Referenced by G4PolarizedIonisationModel::ComputeCrossSectionPerElectron().

XSection()

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

overridevirtual

Implements G4VPolarizedXS.

Definition at line 164 of file G4PolarizedGammaConversionXS.cc.

{
  G4ExceptionDescription ed;
  ed << "ERROR dummy routine G4PolarizedGammaConversionXS::XSection "
        "called \n";
  G4Exception("G4PolarizedGammaConversionXS::Initialize", "pol024",
              FatalException, ed);
  return 0.;
}

References FatalException, and G4Exception().

Field Documentation

fA

G4double G4VPolarizedXS::fA

protectedinherited

Definition at line 90 of file G4VPolarizedXS.hh.

Referenced by G4VPolarizedXS::SetMaterial().

fCoul

G4double G4VPolarizedXS::fCoul

protectedinherited

Definition at line 91 of file G4VPolarizedXS.hh.

Referenced by G4PolarizedBremsstrahlungXS::Initialize(), Initialize(), and G4VPolarizedXS::SetMaterial().

fFinalElectronPolarization

G4StokesVector G4PolarizedGammaConversionXS::fFinalElectronPolarization

private

Definition at line 62 of file G4PolarizedGammaConversionXS.hh.

Referenced by G4PolarizedGammaConversionXS(), GetPol2(), and Initialize().

fFinalPositronPolarization

G4StokesVector G4PolarizedGammaConversionXS::fFinalPositronPolarization

private

Definition at line 63 of file G4PolarizedGammaConversionXS.hh.

Referenced by G4PolarizedGammaConversionXS(), GetPol3(), and Initialize().

fXmax

G4double G4VPolarizedXS::fXmax

protectedinherited

Definition at line 88 of file G4VPolarizedXS.hh.

Referenced by G4VPolarizedXS::GetXmax(), and G4VPolarizedXS::SetXmax().

fXmin

G4double G4VPolarizedXS::fXmin

protectedinherited

Definition at line 88 of file G4VPolarizedXS.hh.

Referenced by G4VPolarizedXS::GetXmin(), and G4VPolarizedXS::SetXmin().

fYmin

G4double G4VPolarizedXS::fYmin

protectedinherited

Definition at line 88 of file G4VPolarizedXS.hh.

Referenced by G4VPolarizedXS::GetYmin(), and G4VPolarizedXS::SetYmin().

fZ

G4double G4VPolarizedXS::fZ

protectedinherited

Definition at line 90 of file G4VPolarizedXS.hh.

Referenced by G4PolarizedBremsstrahlungXS::Initialize(), Initialize(), G4VPolarizedXS::SetMaterial(), and G4PolarizedComptonXS::TotalXSection().

SCRN

G4double G4PolarizedGammaConversionXS::SCRN

staticprivate

Initial value:

= {
  { 0.5, 1.0, 2.0, 4.0, 8.0, 15.0, 20.0, 25.0, 30.0, 35.0, 40.0, 45.0, 50.0,
    60.0, 70.0, 80.0, 90.0, 100.0, 120.0 },
  { 0.0145, 0.0490, 0.1400, 0.3312, 0.6758, 1.126, 1.367, 1.564, 1.731, 1.875,
    2.001, 2.114, 2.216, 2.393, 2.545, 2.676, 2.793, 2.897, 3.078 }
}

Definition at line 60 of file G4PolarizedGammaConversionXS.hh.

Referenced by Initialize().

---

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

• source/processes/electromagnetic/polarisation/include/G4PolarizedGammaConversionXS.hh
• source/processes/electromagnetic/polarisation/src/G4PolarizedGammaConversionXS.cc