source: trunk/source/processes/electromagnetic/polarisation/src/G4PolarizedPEEffectCrossSection.cc@ 819

Last change on this file since 819 was 819, checked in by garnier, 17 years ago

import all except CVS
File size: 6.4 KB
Line 
1//
2// ********************************************************************
3// * License and Disclaimer *
4// * *
5// * The Geant4 software is copyright of the Copyright Holders of *
6// * the Geant4 Collaboration. It is provided under the terms and *
7// * conditions of the Geant4 Software License, included in the file *
8// * LICENSE and available at http://cern.ch/geant4/license . These *
9// * include a list of copyright holders. *
10// * *
11// * Neither the authors of this software system, nor their employing *
12// * institutes,nor the agencies providing financial support for this *
13// * work make any representation or warranty, express or implied, *
14// * regarding this software system or assume any liability for its *
15// * use. Please see the license in the file LICENSE and URL above *
16// * for the full disclaimer and the limitation of liability. *
17// * *
18// * This code implementation is the result of the scientific and *
19// * technical work of the GEANT4 collaboration. *
20// * By using, copying, modifying or distributing the software (or *
21// * any work based on the software) you agree to acknowledge its *
22// * use in resulting scientific publications, and indicate your *
23// * acceptance of all terms of the Geant4 Software license. *
24// ********************************************************************
25//
26// $Id: G4PolarizedPEEffectCrossSection.cc,v 1.1 2007/11/01 17:29:09 schaelic Exp $
27// GEANT4 tag $Name: $
28//
29// GEANT4 Class file
30//
31//
32// File name: G4PolarizedPEEffectCrossSection
33//
34// Author: Karim Laihem
35//
36// Creation date: 15.03.2007
37//
38// Modifications:
39// 19-03-07 Modified to fit in g4.8.2 framework (A.Schaelicke)
40//
41// Class Description:
42//
43#include "G4PolarizedPEEffectCrossSection.hh"
44//#include "Randomize.hh"
45
46using namespace std;
47
48//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
49G4PolarizedPEEffectCrossSection::G4PolarizedPEEffectCrossSection()
50 {
51 cout<<"G4PolarizedPEEffectCrossSection() init\n";
52
53}
54
55//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
56G4PolarizedPEEffectCrossSection::~G4PolarizedPEEffectCrossSection()
57{}
58
59//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
60void G4PolarizedPEEffectCrossSection::Initialize(G4double aGammaE,
61 G4double aLept0E,
62 G4double sinTheta,
63 const G4StokesVector & beamPol,
64 const G4StokesVector & /*p1*/,
65 G4int /*flag*/)
66{
67 // cout<<"G4PolarizedPEEffectCrossSection::Initialize()\n";
68
69// G4StokesVector PolarizedPhotoElectricEffect::Transfer_G4StokesVector(
70// G4double aGammaE, // Incoming Primary Gamma Energy.
71// G4ThreeVector aGammaDir, // Incoming Primary Gamma Direction.
72// G4StokesVector beamPol, // Incoming Primary Gamma polarization.
73// G4double aLept0E, // The Lepton e- of interest Total energy.
74// G4ThreeVector aParticl_01_Dir, // The Lepton e- of interest direction.
75// G4double cos_aTetha_Angle // The lepton of interest Scattering angle.
76// )
77
78// ***********************************************************
79// ************ added by Karim Polarization transfer to e- in PhotoelectricEffect.
80// ************
81// ***********************************************************
82 G4double Gfactor = aLept0E/electron_mass_c2+1;
83 G4double Gfactor_2 = Gfactor * Gfactor;
84
85 G4double BETA = sqrt(1. - 1./(Gfactor_2));
86
87 G4double Stokes_P3 = beamPol.z() ;
88
89 G4double m0_c2 = electron_mass_c2;
90 G4double Lept0E = aLept0E/m0_c2+1, Lept0E2 = Lept0E * Lept0E ;
91 G4double GammaE = aGammaE/m0_c2;
92
93
94// G4double cosTheta = cos_aTetha_Angle;
95// G4double sinTheta = sqrt(1- cos_aTetha_Angle * cos_aTetha_Angle);
96 G4double cosTheta = std::sqrt(1. - sinTheta*sinTheta);
97
98 G4double D_Lepton0 = (1./GammaE) * ((2./(GammaE*Lept0E*(1-BETA*cosTheta)))-1.);
99
100 G4double I_Lepton0 = 1.0+D_Lepton0;
101
102 G4double A_Lepton0 = (Lept0E/(Lept0E+1))*(2.0/(GammaE*Lept0E)
103 + BETA*cosTheta
104 +(2.0/((GammaE*Lept0E2)*(1.0-BETA*cosTheta)))) / I_Lepton0 ;
105
106 G4double B_Lepton0 = (Lept0E/(Lept0E+1.0)) * BETA * sinTheta * (2.0/(GammaE*Lept0E*(1-BETA*cosTheta))-1.0)/I_Lepton0;
107
108 G4double Stokes_S1 = (Stokes_P3 * B_Lepton0) ;
109 G4double Stokes_S2 = 0.;
110 G4double Stokes_S3 = (Stokes_P3 * A_Lepton0) ;
111
112
113 theFinalElectronPolarization.setX(Stokes_S1);
114 theFinalElectronPolarization.setY(Stokes_S2);
115 theFinalElectronPolarization.setZ(Stokes_S3);
116
117 if((theFinalElectronPolarization.x()*theFinalElectronPolarization.x()
118 + theFinalElectronPolarization.y()* theFinalElectronPolarization.y()
119 + theFinalElectronPolarization.z()* theFinalElectronPolarization.z())>1)
120
121 {
122 cout<<"Warning: PhotoelectricEffect Problem in pol-transfer photon to lepton:Px2 + Py2 + Pz2 > 1"<<endl;
123 cout<<"Polarization transfer forced to be total and similar as incoming Photo"<<endl;
124 // *KL* Surprising if it arrives (never seen it up to now)
125 theFinalElectronPolarization = beamPol; // suplement de securite
126// cout<<"PhotoEffect okay :"
127// <<"\t"<<(aLept0E-m0_c2)/aGammaE
128// <<"\t"<<aGammaE
129// <<"\t"<<aLept0E
130// <<"\t"<<cos_aTetha_Angle
131// <<"\t"<<beamPol
132// <<"\t"<<theFinalElectronPolarization
133// <<"\t"<<A_Lepton0
134// <<endl;
135 }
136}
137
138
139
140//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
141
142G4double G4PolarizedPEEffectCrossSection::XSection(const G4StokesVector & /*pol2*/,
143 const G4StokesVector & /*pol3*/)
144{
145 cout<<"ERROR dummy routine G4PolarizedPEEffectCrossSection::XSection() called\n";
146 return 0;
147}
148
149//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
150
151G4StokesVector G4PolarizedPEEffectCrossSection::GetPol2()
152{
153 return theFinalElectronPolarization;
154}
155
156//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
157
158G4StokesVector G4PolarizedPEEffectCrossSection::GetPol3()
159{
160 return G4StokesVector();
161}
162
163
164
165//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
166
167
168//....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo......
Note: See TracBrowser for help on using the repository browser.