Geant4  10.02.p03
Public Member Functions | Private Member Functions | Private Attributes | Static Private Attributes | List of all members
G4PolarizedPairProductionCrossSection Class Reference

#include <G4PolarizedPairProductionCrossSection.hh>

Inheritance diagram for G4PolarizedPairProductionCrossSection:
Collaboration diagram for G4PolarizedPairProductionCrossSection:

Public Member Functions

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

Private Member Functions

void InitializeMe ()
 

Private Attributes

G4StokesVector theFinalElectronPolarization
 
G4StokesVector theFinalPositronPolarization
 

Static Private Attributes

static G4bool scrnInitialized =false
 
static G4double SCRN [3][20]
 

Additional Inherited Members

- Protected Member Functions inherited from G4VPolarizedCrossSection
void SetXmin (G4double xmin)
 
void SetXmax (G4double xmax)
 
void SetYmin (G4double ymin)
 
- Protected Attributes inherited from G4VPolarizedCrossSection
G4double fXmin
 
G4double fXmax
 
G4double fYmin
 
G4double theA
 
G4double theZ
 
G4double fCoul
 

Detailed Description

Definition at line 53 of file G4PolarizedPairProductionCrossSection.hh.

Constructor & Destructor Documentation

◆ G4PolarizedPairProductionCrossSection()

G4PolarizedPairProductionCrossSection::G4PolarizedPairProductionCrossSection ( )

Definition at line 74 of file G4PolarizedPairProductionCrossSection.cc.

75 {
76  InitializeMe();
77 }
Here is the call graph for this function:

Member Function Documentation

◆ GetPol2()

G4StokesVector G4PolarizedPairProductionCrossSection::GetPol2 ( )
virtual

Reimplemented from G4VPolarizedCrossSection.

Definition at line 195 of file G4PolarizedPairProductionCrossSection.cc.

196 {
197  // electron/positron
198  return theFinalElectronPolarization;
199 }

◆ GetPol3()

G4StokesVector G4PolarizedPairProductionCrossSection::GetPol3 ( )
virtual

Reimplemented from G4VPolarizedCrossSection.

Definition at line 200 of file G4PolarizedPairProductionCrossSection.cc.

201 {
202  // photon
203  return theFinalPositronPolarization;;
204 }

◆ Initialize()

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

Reimplemented from G4VPolarizedCrossSection.

Definition at line 80 of file G4PolarizedPairProductionCrossSection.cc.

86 {
87  G4double aLept1E = aGammaE - aLept0E;
88 
89  G4double Stokes_P3 = beamPol.z() ;
90  // **************************************************************************
91 
92  G4double m0_c2 = electron_mass_c2;
93  G4double Lept0E = aLept0E/m0_c2+1., Lept0E2 = Lept0E * Lept0E ;
94  G4double GammaE = aGammaE/m0_c2;
95  G4double Lept1E = aLept1E/m0_c2-1., Lept1E2 = Lept1E * Lept1E ;
96 
97 
98  // const G4Element* theSelectedElement = theModel->SelectedAtom();
99 
100  // ******* Gamma Transvers Momentum
101 
102  G4double TMom = std::sqrt(Lept0E2 -1.)* sintheta, u = TMom , u2 =u * u ;
103  G4double Xsi = 1./(1.+u2) , Xsi2 = Xsi * Xsi ;
104 
105  // G4double theZ = theSelectedElement->GetZ();
106 
107  // G4double fCoul = theSelectedElement->GetfCoulomb();
108  G4double delta = 12. * std::pow(theZ, 1./3.) * Lept0E * Lept1E * Xsi / (121. * GammaE);
109  G4double GG=0.;
110 
111  if(delta < 0.5) {
112  GG = std::log(2.* Lept0E * Lept1E / GammaE) - 2. - fCoul;
113  }
114  else if ( delta < 120) {
115  for (G4int j=2; j<=19; j++) {
116  if(SCRN[1][j] >= delta) {
117  GG =std::log(2 * Lept0E * Lept1E / GammaE) - 2 - fCoul
118  -(SCRN[2][j-1]+(delta-SCRN[1][j-1])*(SCRN[2][j]-SCRN[2][j-1])/(SCRN[1][j]-SCRN[1][j-1]));
119  break;
120  }
121  }
122  }
123  else {
124  G4double alpha_sc = (111 * std::pow(theZ, -1./3.)) / Xsi;
125  GG = std::log(alpha_sc)- 2 - fCoul;
126  }
127 
128  if(GG<-1) GG=-1; // *KL* do we need this ?!
129 
130 
131  G4double I_Lepton = (Lept0E2 + Lept1E2)*(3+2*GG) + 2 * Lept0E * Lept1E * (1 + 4 * u2 * Xsi2 * GG);
132 
133  // G4double D_Lepton1 = -8 * Lept0E * Lept1E * u2 * Xsi2 * GG / I_Lepton;
134 
135  G4double L_Lepton1 = GammaE * ((Lept0E - Lept1E) * (3. + 2. * GG)+2 * Lept1E * (1 + 4 * u2 * Xsi2 * GG))/I_Lepton;
136 
137  G4double T_Lepton1 = 4 * GammaE * Lept1E * Xsi * u * (1. - 2. * Xsi) * GG / I_Lepton ;
138 
139 
140  G4double Stokes_S1 = (Stokes_P3 * T_Lepton1) ;
141  G4double Stokes_S2 = 0;
142  G4double Stokes_S3 = (Stokes_P3 * L_Lepton1) ;
143 
144 
145  theFinalElectronPolarization.setX(Stokes_S1);
146  theFinalElectronPolarization.setY(Stokes_S2);
147  theFinalElectronPolarization.setZ(Stokes_S3);
148 
149  if(theFinalElectronPolarization.mag2()>1) {
150  G4cout<<" WARNING in pol-conv theFinalElectronPolarization \n";
151  G4cout
152  <<"\t"<<theFinalElectronPolarization
153  <<"\t GG\t"<<GG
154  <<"\t delta\t"<<delta
155  <<G4endl;
156  theFinalElectronPolarization.setX(0);
157  theFinalElectronPolarization.setY(0);
158  theFinalElectronPolarization.setZ(Stokes_S3);
159  if(Stokes_S3>1) theFinalElectronPolarization.setZ(1);
160  }
161 
162 
163  G4double L_Lepton2 = GammaE * ((Lept1E - Lept0E) * (3. + 2. * GG)+2 * Lept0E * (1 + 4 * u2 * Xsi2 * GG))/I_Lepton;
164 
165  G4double T_Lepton2 = 4 * GammaE * Lept0E * Xsi * u * (1. - 2. * Xsi) * GG / I_Lepton ;
166 
167  G4double Stokes_SS1 = (Stokes_P3 * T_Lepton2) ;
168  G4double Stokes_SS2 = 0;
169  G4double Stokes_SS3 = (Stokes_P3 * L_Lepton2) ;
170 
171  theFinalPositronPolarization.SetPhoton();
172 
173  theFinalPositronPolarization.setX(Stokes_SS1);
174  theFinalPositronPolarization.setY(Stokes_SS2);
175  theFinalPositronPolarization.setZ(Stokes_SS3);
176 
177  if(theFinalPositronPolarization.mag2()>1) {
178  G4cout<<" WARNING in pol-conv theFinalPositronPolarization \n";
179  G4cout
180  <<"\t"<<theFinalPositronPolarization
181  <<"\t GG\t"<<GG
182  <<"\t delta\t"<<delta
183  <<G4endl;
184  }
185 }
G4int
int G4int
Definition: G4Types.hh:78
CLHEP::Hep3Vector::setY
void setY(double)
CLHEP::Hep3Vector::setZ
void setZ(double)
CLHEP::Hep3Vector::setX
void setX(double)
CLHEP::Hep3Vector::mag2
double mag2() const
G4cout
G4GLOB_DLL std::ostream G4cout
python.hepunit.electron_mass_c2
float electron_mass_c2
Definition: hepunit.py:274
G4endl
#define G4endl
Definition: G4ios.hh:61
G4double
double G4double
Definition: G4Types.hh:76
Here is the call graph for this function:

◆ InitializeMe()

void G4PolarizedPairProductionCrossSection::InitializeMe ( )
private

Definition at line 47 of file G4PolarizedPairProductionCrossSection.cc.

48 {
49  if (!scrnInitialized) {
50  SCRN [1][1]= 0.5 ; SCRN [2][1] = 0.0145;
51  SCRN [1][2]= 1.0 ; SCRN [2][2] = 0.0490;
52  SCRN [1][3]= 2.0 ; SCRN [2][3] = 0.1400;
53  SCRN [1][4]= 4.0 ; SCRN [2][4] = 0.3312;
54  SCRN [1][5]= 8.0 ; SCRN [2][5] = 0.6758;
55  SCRN [1][6]= 15.0 ; SCRN [2][6] = 1.126;
56  SCRN [1][7]= 20.0 ; SCRN [2][7] = 1.367;
57  SCRN [1][8]= 25.0 ; SCRN [2][8] = 1.564;
58  SCRN [1][9]= 30.0 ; SCRN [2][9] = 1.731;
59  SCRN [1][10]= 35.0 ; SCRN [2][10]= 1.875;
60  SCRN [1][11]= 40.0 ; SCRN [2][11]= 2.001;
61  SCRN [1][12]= 45.0 ; SCRN [2][12]= 2.114;
62  SCRN [1][13]= 50.0 ; SCRN [2][13]= 2.216;
63  SCRN [1][14]= 60.0 ; SCRN [2][14]= 2.393;
64  SCRN [1][15]= 70.0 ; SCRN [2][15]= 2.545;
65  SCRN [1][16]= 80.0 ; SCRN [2][16]= 2.676;
66  SCRN [1][17]= 90.0 ; SCRN [2][17]= 2.793;
67  SCRN [1][18]= 100.0 ; SCRN [2][18]= 2.897;
68  SCRN [1][19]= 120.0 ; SCRN [2][19]= 3.078;
69 
70  scrnInitialized=true;
71  }
72 }
Here is the caller graph for this function:

◆ XSection()

G4double G4PolarizedPairProductionCrossSection::XSection ( const G4StokesVector &  pol2,
const G4StokesVector &  pol3 
)
virtual

Implements G4VPolarizedCrossSection.

Definition at line 187 of file G4PolarizedPairProductionCrossSection.cc.

189 {
190  G4cout<<"ERROR dummy routine G4PolarizedPairProductionCrossSection::XSection called \n";
191  return 0;
192 }
G4cout
G4GLOB_DLL std::ostream G4cout

Member Data Documentation

◆ SCRN

G4double G4PolarizedPairProductionCrossSection::SCRN
staticprivate

Definition at line 75 of file G4PolarizedPairProductionCrossSection.hh.

◆ scrnInitialized

G4bool G4PolarizedPairProductionCrossSection::scrnInitialized =false
staticprivate

Definition at line 74 of file G4PolarizedPairProductionCrossSection.hh.

◆ theFinalElectronPolarization

G4StokesVector G4PolarizedPairProductionCrossSection::theFinalElectronPolarization
private

Definition at line 69 of file G4PolarizedPairProductionCrossSection.hh.

◆ theFinalPositronPolarization

G4StokesVector G4PolarizedPairProductionCrossSection::theFinalPositronPolarization
private

Definition at line 70 of file G4PolarizedPairProductionCrossSection.hh.

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