G4PolarizedPEEffectCrossSection.cc

Go to the documentation of this file.
00001 //
00002 // ********************************************************************
00003 // * License and Disclaimer                                           *
00004 // *                                                                  *
00005 // * The  Geant4 software  is  copyright of the Copyright Holders  of *
00006 // * the Geant4 Collaboration.  It is provided  under  the terms  and *
00007 // * conditions of the Geant4 Software License,  included in the file *
00008 // * LICENSE and available at  http://cern.ch/geant4/license .  These *
00009 // * include a list of copyright holders.                             *
00010 // *                                                                  *
00011 // * Neither the authors of this software system, nor their employing *
00012 // * institutes,nor the agencies providing financial support for this *
00013 // * work  make  any representation or  warranty, express or implied, *
00014 // * regarding  this  software system or assume any liability for its *
00015 // * use.  Please see the license in the file  LICENSE  and URL above *
00016 // * for the full disclaimer and the limitation of liability.         *
00017 // *                                                                  *
00018 // * This  code  implementation is the result of  the  scientific and *
00019 // * technical work of the GEANT4 collaboration.                      *
00020 // * By using,  copying,  modifying or  distributing the software (or *
00021 // * any work based  on the software)  you  agree  to acknowledge its *
00022 // * use  in  resulting  scientific  publications,  and indicate your *
00023 // * acceptance of all terms of the Geant4 Software license.          *
00024 // ********************************************************************
00025 //
00026 // $Id: G4PolarizedPEEffectCrossSection.cc 69847 2013-05-16 09:36:18Z gcosmo $
00027 //
00028 // GEANT4 Class file
00029 //
00030 //
00031 // File name:     G4PolarizedPEEffectCrossSection
00032 //
00033 // Author:        Karim Laihem 
00034 //
00035 // Creation date: 15.03.2007
00036 //
00037 // Modifications:
00038 //   19-03-07 Modified to fit in g4.8.2 framework (A.Schaelicke)
00039 //
00040 // Class Description:
00041 //
00042 #include "G4PolarizedPEEffectCrossSection.hh"
00043 #include "G4PhysicalConstants.hh"
00044 
00045 using namespace std;
00046 
00047 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00048 G4PolarizedPEEffectCrossSection::G4PolarizedPEEffectCrossSection()
00049   {
00050    cout<<"G4PolarizedPEEffectCrossSection() init\n";
00051 
00052 }
00053 
00054 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00055 G4PolarizedPEEffectCrossSection::~G4PolarizedPEEffectCrossSection()
00056 {}
00057 
00058 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00059 void G4PolarizedPEEffectCrossSection::Initialize(G4double aGammaE, 
00060                                                  G4double aLept0E, 
00061                                                  G4double sinTheta,
00062                                                  const G4StokesVector & beamPol,
00063                                                  const G4StokesVector & /*p1*/,
00064                                                  G4int /*flag*/)
00065 {
00066   //   cout<<"G4PolarizedPEEffectCrossSection::Initialize()\n";
00067 
00068 // G4StokesVector PolarizedPhotoElectricEffect::Transfer_G4StokesVector(
00069 //     G4double aGammaE,                 // Incoming Primary Gamma Energy.
00070 //     G4ThreeVector aGammaDir,          // Incoming Primary Gamma Direction.
00071 //     G4StokesVector beamPol,          // Incoming Primary Gamma  polarization.
00072 //     G4double aLept0E,                 // The Lepton  e- of interest Total energy.
00073 //     G4ThreeVector aParticl_01_Dir,    // The Lepton  e- of interest direction.
00074 //     G4double cos_aTetha_Angle            // The lepton of interest Scattering angle.
00075 //     )
00076 
00077 // ***********************************************************
00078 // ************  added by Karim   Polarization transfer to e- in PhotoelectricEffect.
00079 // ************ 
00080 // ***********************************************************
00081     G4double Gfactor   = aLept0E/electron_mass_c2+1;
00082     G4double Gfactor_2 = Gfactor * Gfactor;
00083 
00084     G4double BETA    = sqrt(1. - 1./(Gfactor_2));
00085 
00086     G4double Stokes_P3  = beamPol.z()   ;                       
00087     
00088     G4double m0_c2  = electron_mass_c2; 
00089     G4double Lept0E = aLept0E/m0_c2+1,   Lept0E2 = Lept0E * Lept0E ;
00090     G4double GammaE = aGammaE/m0_c2;
00091 
00092 
00093 //     G4double cosTheta = cos_aTetha_Angle;
00094 //     G4double sinTheta = sqrt(1- cos_aTetha_Angle * cos_aTetha_Angle);
00095     G4double cosTheta = std::sqrt(1. - sinTheta*sinTheta);
00096 
00097     G4double D_Lepton0 = (1./GammaE) * ((2./(GammaE*Lept0E*(1-BETA*cosTheta)))-1.);   
00098 
00099     G4double I_Lepton0 = 1.0+D_Lepton0;
00100 
00101     G4double A_Lepton0 = (Lept0E/(Lept0E+1))*(2.0/(GammaE*Lept0E) 
00102                                               + BETA*cosTheta 
00103                                               +(2.0/((GammaE*Lept0E2)*(1.0-BETA*cosTheta)))) / I_Lepton0 ;
00104     
00105     G4double B_Lepton0 = (Lept0E/(Lept0E+1.0)) * BETA * sinTheta * (2.0/(GammaE*Lept0E*(1-BETA*cosTheta))-1.0)/I_Lepton0;   
00106 
00107     G4double Stokes_S1 = (Stokes_P3 * B_Lepton0) ;
00108     G4double Stokes_S2 = 0.;
00109     G4double Stokes_S3 = (Stokes_P3 * A_Lepton0) ; 
00110 
00111     
00112     theFinalElectronPolarization.setX(Stokes_S1);
00113     theFinalElectronPolarization.setY(Stokes_S2);
00114     theFinalElectronPolarization.setZ(Stokes_S3);
00115 
00116     if((theFinalElectronPolarization.x()*theFinalElectronPolarization.x()
00117         + theFinalElectronPolarization.y()* theFinalElectronPolarization.y()
00118         + theFinalElectronPolarization.z()* theFinalElectronPolarization.z())>1)
00119         
00120     {
00121         cout<<"Warning: PhotoelectricEffect Problem in pol-transfer photon to lepton:Px2 + Py2 + Pz2 > 1"<<endl;
00122         cout<<"Polarization transfer forced to be total and similar as incoming Photo"<<endl;
00123         // *KL* Surprising if it arrives (never seen it up to now)
00124         theFinalElectronPolarization = beamPol; // suplement de securite
00125 //      cout<<"PhotoEffect okay :"
00126 //            <<"\t"<<(aLept0E-m0_c2)/aGammaE
00127 //            <<"\t"<<aGammaE
00128 //            <<"\t"<<aLept0E
00129 //            <<"\t"<<cos_aTetha_Angle
00130 //            <<"\t"<<beamPol
00131 //            <<"\t"<<theFinalElectronPolarization
00132 //            <<"\t"<<A_Lepton0
00133 //            <<endl;
00134     }
00135 }
00136 
00137 
00138 
00139 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00140 
00141 G4double G4PolarizedPEEffectCrossSection::XSection(const G4StokesVector & /*pol2*/,
00142                                                    const G4StokesVector & /*pol3*/)
00143 {
00144   cout<<"ERROR dummy routine G4PolarizedPEEffectCrossSection::XSection() called\n";
00145   return 0;
00146 } 
00147 
00148 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00149 
00150 G4StokesVector G4PolarizedPEEffectCrossSection::GetPol2()
00151 {
00152   return theFinalElectronPolarization;
00153 }
00154 
00155 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00156 
00157 G4StokesVector G4PolarizedPEEffectCrossSection::GetPol3()
00158 {
00159   return G4StokesVector();
00160 }
00161 
00162 
00163 
00164 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
00165 
00166 
00167 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......

Generated on Mon May 27 17:49:23 2013 for Geant4 by  doxygen 1.4.7