G4PolarizedBremsstrahlungXS Class Reference

#include <G4PolarizedBremsstrahlungXS.hh>

Inheritance diagram for G4PolarizedBremsstrahlungXS:

Public Member Functions
	G4PolarizedBremsstrahlungXS ()
	G4PolarizedBremsstrahlungXS (const G4PolarizedBremsstrahlungXS &)=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
G4PolarizedBremsstrahlungXS &	operator= (const G4PolarizedBremsstrahlungXS &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
	~G4PolarizedBremsstrahlungXS () 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	fFinalGammaPolarization
G4StokesVector	fFinalLeptonPolarization

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

Detailed Description

Definition at line 44 of file G4PolarizedBremsstrahlungXS.hh.

Constructor & Destructor Documentation

G4PolarizedBremsstrahlungXS() [1/2]

G4PolarizedBremsstrahlungXS::G4PolarizedBremsstrahlungXS

(

)

Definition at line 45 of file G4PolarizedBremsstrahlungXS.cc.

{
  fFinalLeptonPolarization = G4StokesVector::ZERO;
  fFinalGammaPolarization  = G4StokesVector::ZERO;
}

References fFinalGammaPolarization, fFinalLeptonPolarization, and G4StokesVector::ZERO.

~G4PolarizedBremsstrahlungXS()

G4PolarizedBremsstrahlungXS::~G4PolarizedBremsstrahlungXS ( )

override

Definition at line 51 of file G4PolarizedBremsstrahlungXS.cc.

{}

G4PolarizedBremsstrahlungXS() [2/2]

G4PolarizedBremsstrahlungXS::G4PolarizedBremsstrahlungXS ( const G4PolarizedBremsstrahlungXS &  )

delete

Member Function Documentation

GetPol2()

G4StokesVector G4PolarizedBremsstrahlungXS::GetPol2 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 188 of file G4PolarizedBremsstrahlungXS.cc.

{
  // electron/positron
  return fFinalLeptonPolarization;
}

References fFinalLeptonPolarization.

GetPol3()

G4StokesVector G4PolarizedBremsstrahlungXS::GetPol3 ( )

overridevirtual

Reimplemented from G4VPolarizedXS.

Definition at line 194 of file G4PolarizedBremsstrahlungXS.cc.

{
  // photon
  return fFinalGammaPolarization;
}

References fFinalGammaPolarization.

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 G4PolarizedBremsstrahlungXS::Initialize ( G4double  eps, G4double  X, G4double  phi, const G4StokesVector &  p0, const G4StokesVector &  p1, G4int  flag = 0  )

overridevirtual

Implements G4VPolarizedXS.

Definition at line 53 of file G4PolarizedBremsstrahlungXS.cc.

{
  G4double aLept1E = aLept0E - aGammaE;
 
  G4double Stokes_S1 = beamPol.x();
  G4double Stokes_S2 = beamPol.y();
  G4double Stokes_S3 = beamPol.z();
 
  G4double Lept0E  = aLept0E / electron_mass_c2 + 1.;
  G4double Lept0E2 = Lept0E * Lept0E;
  G4double GammaE  = aGammaE / electron_mass_c2;
  G4double GammaE2 = GammaE * GammaE;
  G4double Lept1E  = aLept1E / electron_mass_c2 + 1.;
  G4double Lept1E2 = Lept1E * Lept1E;
 
  // *******  Gamma Transverse 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_Lept = (Lept0E2 + Lept1E2) * (3. + 2. * GG) -
                    2 * Lept0E * Lept1E * (1. + 4. * u2 * Xsi2 * GG);
  G4double F_Lept =
    Lept1E * 4. * GammaE * u * Xsi * (1. - 2 * Xsi) * GG / I_Lept;
  G4double E_Lept =
    Lept0E * 4. * GammaE * u * Xsi * (2. * Xsi - 1.) * GG / I_Lept;
  G4double M_Lept =
    4. * Lept0E * Lept1E * (1. + GG - 2. * Xsi2 * u2 * GG) / I_Lept;
  G4double P_Lept =
    GammaE2 * (1. + 8. * GG * (Xsi - 0.5) * (Xsi - 0.5)) / I_Lept;
 
  G4double Stokes_SS1 = M_Lept * Stokes_S1 + E_Lept * Stokes_S3;
  G4double Stokes_SS2 = M_Lept * Stokes_S2;
  G4double Stokes_SS3 = (M_Lept + P_Lept) * Stokes_S3 + F_Lept * Stokes_S1;
 
  fFinalLeptonPolarization.setX(Stokes_SS1);
  fFinalLeptonPolarization.setY(Stokes_SS2);
  fFinalLeptonPolarization.setZ(Stokes_SS3);
 
  if(fFinalLeptonPolarization.mag2() > 1.)
  {
    G4ExceptionDescription ed;
    ed << " WARNING in pol-brem fFinalLeptonPolarization \n";
    ed << "\t" << fFinalLeptonPolarization << "\t GG\t" << GG << "\t delta\t"
       << delta;
    G4Exception("G4PolarizedBremsstrahlungXS::Initialize", "pol014",
                JustWarning, ed);
    fFinalLeptonPolarization.setX(0);
    fFinalLeptonPolarization.setY(0);
    fFinalLeptonPolarization.setZ(Stokes_SS3);
    if(Stokes_SS3 > 1)
      fFinalLeptonPolarization.setZ(1);
  }
 
  G4double I_Gamma = (Lept0E2 + Lept1E2) * (3. + 2. * GG) -
                     2. * Lept0E * Lept1E * (1. + 4. * u2 * Xsi2 * GG);
  G4double D_Gamma = 8. * Lept0E * Lept1E * u2 * Xsi2 * GG / I_Gamma;
  G4double L_Gamma = GammaE *
                     ((Lept0E + Lept1E) * (3. + 2. * GG) -
                      2. * Lept1E * (1. + 4. * u2 * Xsi2 * GG)) /
                     I_Gamma;
  G4double T_Gamma =
    4. * GammaE * Lept1E * Xsi * u * (2. * Xsi - 1.) * GG / I_Gamma;
 
  G4double Stokes_P1 = D_Gamma;
  G4double Stokes_P2 = 0.;
  G4double Stokes_P3 = (Stokes_S3 * L_Gamma + Stokes_S1 * T_Gamma);
 
  fFinalGammaPolarization.SetPhoton();
 
  fFinalGammaPolarization.setX(Stokes_P1);
  fFinalGammaPolarization.setY(Stokes_P2);
  fFinalGammaPolarization.setZ(Stokes_P3);
 
  if(fFinalGammaPolarization.mag2() > 1.)
  {
    G4ExceptionDescription ed;
    ed << " WARNING in pol-brem fFinalGammaPolarization \n";
    ed << "\t" << fFinalGammaPolarization << "\t GG\t" << GG << "\t delta\t"
       << delta;
    G4Exception("G4PolarizedBremsstrahlungXS::Initialize", "pol015",
                JustWarning, ed);
  }
}

References source.hepunit::electron_mass_c2, G4VPolarizedXS::fCoul, fFinalGammaPolarization, fFinalLeptonPolarization, G4VPolarizedXS::fZ, G4Exception(), JustWarning, CLHEP::Hep3Vector::mag2(), SCRN, G4StokesVector::SetPhoton(), CLHEP::Hep3Vector::setX(), CLHEP::Hep3Vector::setY(), CLHEP::Hep3Vector::setZ(), CLHEP::Hep3Vector::x(), CLHEP::Hep3Vector::y(), and CLHEP::Hep3Vector::z().

operator=()

G4PolarizedBremsstrahlungXS & G4PolarizedBremsstrahlungXS::operator= ( const G4PolarizedBremsstrahlungXS &  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 G4PolarizedBremsstrahlungXS::XSection ( const G4StokesVector &  pol2, const G4StokesVector &  pol3  )

overridevirtual

Implements G4VPolarizedXS.

Definition at line 175 of file G4PolarizedBremsstrahlungXS.cc.

{
  G4ExceptionDescription ed;
  ed << "ERROR dummy routine G4PolarizedBremsstrahlungXS::XSection "
        "called.\n";
  G4Exception("G4PolarizedBremsstrahlungXS::XSection", "pol016", 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 Initialize(), G4PolarizedGammaConversionXS::Initialize(), and G4VPolarizedXS::SetMaterial().

fFinalGammaPolarization

G4StokesVector G4PolarizedBremsstrahlungXS::fFinalGammaPolarization

private

Definition at line 68 of file G4PolarizedBremsstrahlungXS.hh.

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

fFinalLeptonPolarization

G4StokesVector G4PolarizedBremsstrahlungXS::fFinalLeptonPolarization

private

Definition at line 67 of file G4PolarizedBremsstrahlungXS.hh.

Referenced by G4PolarizedBremsstrahlungXS(), GetPol2(), 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 Initialize(), G4PolarizedGammaConversionXS::Initialize(), G4VPolarizedXS::SetMaterial(), and G4PolarizedComptonXS::TotalXSection().

SCRN

G4double G4PolarizedBremsstrahlungXS::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 65 of file G4PolarizedBremsstrahlungXS.hh.

Referenced by Initialize().

---

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

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