G4PolarizedPEEffectCrossSection Class Reference

#include <G4PolarizedPEEffectCrossSection.hh>

Inheritance diagram for G4PolarizedPEEffectCrossSection:

Collaboration diagram for G4PolarizedPEEffectCrossSection:

Public Member Functions
	G4PolarizedPEEffectCrossSection ()
virtual	~G4PolarizedPEEffectCrossSection ()
virtual void	Initialize (G4double aGammaE, G4double aLept0E, G4double sintheta, const G4StokesVector &beamPol, const G4StokesVector &, G4int flag=0)
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 Attributes
G4StokesVector	theFinalElectronPolarization

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 50 of file G4PolarizedPEEffectCrossSection.hh.

Constructor & Destructor Documentation

G4PolarizedPEEffectCrossSection()

G4PolarizedPEEffectCrossSection::G4PolarizedPEEffectCrossSection

(

)

Definition at line 48 of file G4PolarizedPEEffectCrossSection.cc.

  {
   cout<<"G4PolarizedPEEffectCrossSection() init\n";

}

~G4PolarizedPEEffectCrossSection()

G4PolarizedPEEffectCrossSection::~G4PolarizedPEEffectCrossSection ( )

virtual

Definition at line 55 of file G4PolarizedPEEffectCrossSection.cc.

{}

Member Function Documentation

GetPol2()

G4StokesVector G4PolarizedPEEffectCrossSection::GetPol2 ( )

virtual

Reimplemented from G4VPolarizedCrossSection.

Definition at line 150 of file G4PolarizedPEEffectCrossSection.cc.

{
  return theFinalElectronPolarization;
}

Here is the caller graph for this function:

GetPol3()

G4StokesVector G4PolarizedPEEffectCrossSection::GetPol3 ( )

virtual

Reimplemented from G4VPolarizedCrossSection.

Definition at line 157 of file G4PolarizedPEEffectCrossSection.cc.

{
  return G4StokesVector();
}

Initialize()

void G4PolarizedPEEffectCrossSection::Initialize ( G4double  aGammaE, G4double  aLept0E, G4double  sintheta, const G4StokesVector &  beamPol, const G4StokesVector &  , G4int  flag = 0  )

virtual

Reimplemented from G4VPolarizedCrossSection.

Definition at line 59 of file G4PolarizedPEEffectCrossSection.cc.

{
  //   cout<<"G4PolarizedPEEffectCrossSection::Initialize()\n";

// G4StokesVector PolarizedPhotoElectricEffect::Transfer_G4StokesVector(
//     G4double aGammaE,                 // Incoming Primary Gamma Energy.
//     G4ThreeVector aGammaDir,          // Incoming Primary Gamma Direction.
//     G4StokesVector beamPol,          // Incoming Primary Gamma  polarization.
//     G4double aLept0E,                 // The Lepton  e- of interest Total energy.
//     G4ThreeVector aParticl_01_Dir,    // The Lepton  e- of interest direction.
//     G4double cos_aTetha_Angle            // The lepton of interest Scattering angle.
//     )

// ***********************************************************
// ************  added by Karim   Polarization transfer to e- in PhotoelectricEffect.
// ************ 
// ***********************************************************
    G4double Gfactor   = aLept0E/electron_mass_c2+1;
    G4double Gfactor_2 = Gfactor * Gfactor;

    G4double BETA    = sqrt(1. - 1./(Gfactor_2));

    G4double Stokes_P3  = beamPol.z()   ;                       
    
    G4double m0_c2  = electron_mass_c2; 
    G4double Lept0E = aLept0E/m0_c2+1,   Lept0E2 = Lept0E * Lept0E ;
    G4double GammaE = aGammaE/m0_c2;


//     G4double cosTheta = cos_aTetha_Angle;
//     G4double sinTheta = sqrt(1- cos_aTetha_Angle * cos_aTetha_Angle);
    G4double cosTheta = std::sqrt(1. - sinTheta*sinTheta);

    G4double D_Lepton0 = (1./GammaE) * ((2./(GammaE*Lept0E*(1-BETA*cosTheta)))-1.);   

    G4double I_Lepton0 = 1.0+D_Lepton0;

    G4double A_Lepton0 = (Lept0E/(Lept0E+1))*(2.0/(GammaE*Lept0E) 
                          + BETA*cosTheta 
                          +(2.0/((GammaE*Lept0E2)*(1.0-BETA*cosTheta)))) / I_Lepton0 ;
    
    G4double B_Lepton0 = (Lept0E/(Lept0E+1.0)) * BETA * sinTheta * (2.0/(GammaE*Lept0E*(1-BETA*cosTheta))-1.0)/I_Lepton0;   

    G4double Stokes_S1 = (Stokes_P3 * B_Lepton0) ;
    G4double Stokes_S2 = 0.;
    G4double Stokes_S3 = (Stokes_P3 * A_Lepton0) ; 

    
    theFinalElectronPolarization.setX(Stokes_S1);
    theFinalElectronPolarization.setY(Stokes_S2);
    theFinalElectronPolarization.setZ(Stokes_S3);

    if((theFinalElectronPolarization.x()*theFinalElectronPolarization.x()
        + theFinalElectronPolarization.y()* theFinalElectronPolarization.y()
    + theFinalElectronPolarization.z()* theFinalElectronPolarization.z())>1)
    
    {
    cout<<"Warning: PhotoelectricEffect Problem in pol-transfer photon to lepton:Px2 + Py2 + Pz2 > 1"<<endl;
    cout<<"Polarization transfer forced to be total and similar as incoming Photo"<<endl;
    // *KL* Surprising if it arrives (never seen it up to now)
    theFinalElectronPolarization = beamPol; // suplement de securite
//  cout<<"PhotoEffect okay :"
//            <<"\t"<<(aLept0E-m0_c2)/aGammaE
//        <<"\t"<<aGammaE
//        <<"\t"<<aLept0E
//        <<"\t"<<cos_aTetha_Angle
//        <<"\t"<<beamPol
//        <<"\t"<<theFinalElectronPolarization
//        <<"\t"<<A_Lepton0
//        <<endl;
    }
}

Here is the call graph for this function:

Here is the caller graph for this function:

XSection()

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

virtual

Implements G4VPolarizedCrossSection.

Definition at line 141 of file G4PolarizedPEEffectCrossSection.cc.

{
  cout<<"ERROR dummy routine G4PolarizedPEEffectCrossSection::XSection() called\n";
  return 0;
} 

Member Data Documentation

theFinalElectronPolarization

G4StokesVector G4PolarizedPEEffectCrossSection::theFinalElectronPolarization

private

Definition at line 68 of file G4PolarizedPEEffectCrossSection.hh.

---

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

• G4PolarizedPEEffectCrossSection.hh
• G4PolarizedPEEffectCrossSection.cc