| Geant4 Cross Reference |
1 // 1
2 // *******************************************
3 // * License and Disclaimer
4 // *
5 // * The Geant4 software is copyright of th
6 // * the Geant4 Collaboration. It is provided
7 // * conditions of the Geant4 Software License
8 // * LICENSE and available at http://cern.ch/
9 // * include a list of copyright holders.
10 // *
11 // * Neither the authors of this software syst
12 // * institutes,nor the agencies providing fin
13 // * work make any representation or warran
14 // * regarding this software system or assum
15 // * use. Please see the license in the file
16 // * for the full disclaimer and the limitatio
17 // *
18 // * This code implementation is the result
19 // * technical work of the GEANT4 collaboratio
20 // * By using, copying, modifying or distri
21 // * any work based on the software) you ag
22 // * use in resulting scientific publicati
23 // * acceptance of all terms of the Geant4 Sof
24 // *******************************************
25 //
26 // -------------------------------------------
27 //
28 // Geant4 Class header file
29 //
30 // File name: G4PolarizedComptonModel
31 //
32 // Author: Andreas Schaelicke
33 //
34 // Class Description:
35 // Implementation of polarized gamma Compton
36
37 #ifndef G4PolarizedComptonModel_h
38 #define G4PolarizedComptonModel_h 1
39
40 #include "globals.hh"
41 #include "G4KleinNishinaCompton.hh"
42 #include "G4StokesVector.hh"
43 #include "G4ThreeVector.hh"
44
45 class G4DynamicParticle;
46 class G4MaterialCutsCouple;
47 class G4ParticleDefinition;
48 class G4ParticleChangeForGamma;
49 class G4VPolarizedXS;
50
51 class G4PolarizedComptonModel : public G4Klein
52 {
53 public:
54 explicit G4PolarizedComptonModel(const G4Par
55 const G4Str
56
57 virtual ~G4PolarizedComptonModel() override;
58
59 virtual G4double ComputeCrossSectionPerAtom(
60
61
62
63 G4double ComputeAsymmetryPerAtom(G4double ga
64
65 virtual void SampleSecondaries(std::vector<G
66 const G4Mater
67 const G4Dynam
68 G4double maxE
69
70 // polarized routines
71 inline void SetTargetPolarization(const G4Th
72 inline void SetBeamPolarization(const G4Thre
73 inline const G4ThreeVector& GetTargetPolariz
74 inline const G4ThreeVector& GetBeamPolarizat
75 inline const G4ThreeVector& GetFinalGammaPol
76 inline const G4ThreeVector& GetFinalElectron
77
78 // hide assignment operator
79 G4PolarizedComptonModel& operator=(const G4P
80 delete;
81 G4PolarizedComptonModel(const G4PolarizedCom
82
83 private:
84 void PrintWarning(const G4DynamicParticle*,
85 G4double onecos, G4double
86
87 static constexpr G4int fLoopLim = 10000;
88
89 G4VPolarizedXS* fCrossSectionCalculator;
90 // incoming
91 G4StokesVector fBeamPolarization; // phot
92 G4StokesVector fTargetPolarization; // elec
93 // outgoing
94 G4StokesVector fFinalGammaPolarization;
95 G4StokesVector finalElectronPolarization;
96
97 G4int fVerboseLevel;
98 };
99
100 //....oooOO0OOooo........oooOO0OOooo........oo
101
102 inline void G4PolarizedComptonModel::SetTarget
103 const G4ThreeVector& pTarget)
104 {
105 fTargetPolarization = G4StokesVector(pTarget
106 }
107 inline void G4PolarizedComptonModel::SetBeamPo
108 const G4ThreeVector& pBeam)
109 {
110 fBeamPolarization = G4StokesVector(pBeam);
111 }
112 inline const G4ThreeVector& G4PolarizedCompton
113 const
114 {
115 return fTargetPolarization;
116 }
117 inline const G4ThreeVector& G4PolarizedCompton
118 {
119 return fBeamPolarization;
120 }
121 inline const G4ThreeVector& G4PolarizedCompton
122 const
123 {
124 return fFinalGammaPolarization;
125 }
126 inline const G4ThreeVector&
127 G4PolarizedComptonModel::GetFinalElectronPolar
128 {
129 return finalElectronPolarization;
130 }
131
132 #endif
133