| 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 Authors:
28 M. Omer and R. Hajima on 15 November 2019
29 contact:
30 omer.mohamed@jaea.go.jp and hajima.ryoichi@qst
31 Publication Information:
32 1- M. Omer, R. Hajima, Validating polarization
33 Carlo simulation, New J. Phys., vol. 21, 2019,
34 https://doi.org/10.1088/1367-2630/ab4d8a
35 */
36
37 #ifndef G4JAEAPolarizedElasticScatteringModel_
38 #define G4JAEAPolarizedElasticScatteringModel_
39
40 #include "G4VEmModel.hh"
41 #include "G4ParticleChangeForGamma.hh"
42 #include "G4PhysicsFreeVector.hh"
43 #include "G4ProductionCutsTable.hh"
44 #include "G4DataVector.hh"
45
46 class G4JAEAPolarizedElasticScatteringModel :
47 {
48 public:
49 explicit G4JAEAPolarizedElasticScatteringMod
50 virtual ~G4JAEAPolarizedElasticScatteringMod
51
52 void Initialise(const G4ParticleDefinition*,
53
54 void InitialiseLocal(const G4ParticleDefinit
55 G4VEmModel* masterModel) override;
56
57 void InitialiseForElement(const G4ParticleDe
58
59 G4double ComputeCrossSectionPerAtom(
60 const G4ParticleDefinition*,
61 G4double
62 G4double
63 G4double
64 G4double
65 G4double
66
67 void SampleSecondaries(std::vector<G4Dynamic
68 const G4MaterialCutsCouple*,
69 const G4DynamicParticle*,
70 G4double tmin,
71 G4double maxEnergy) override;
72
73 inline void SetLowEnergyThreshold(G4double v
74 inline void SetPolarizationSensitvity(G4bool
75
76 void SetDebugVerbosity(G4int val){verboseLev
77
78 G4JAEAPolarizedElasticScatteringModel & oper
79 G4JAEAPolarizedElasticScatteringModel(const
80
81 private:
82 void ReadData(std::size_t Z, const char* pat
83 G4double GeneratePolarizedPhi(G4double Sigma
84
85 static const G4int maxZ = 99;
86 static G4PhysicsFreeVector* dataCS[maxZ+1];
87 static G4DataVector* Polarized_ES_Data[maxZ+
88 G4double distribution[181];
89 G4double cdistribution[181];
90
91 G4ParticleChangeForGamma* fParticleChange;
92 G4double lowEnergyLimit;
93
94 G4int verboseLevel;
95
96 G4bool fLinearPolarizationSensitvity1;
97 G4bool fLinearPolarizationSensitvity2;
98 G4bool fCircularPolarizationSensitvity;
99 G4bool isInitialised;
100 };
101
102 //....oooOO0OOooo........oooOO0OOooo........oo
103
104 inline void G4JAEAPolarizedElasticScatteringMo
105 G4bool linear2,
106 G4bool circular)
107 {
108 fLinearPolarizationSensitvity1=linear1;
109 fLinearPolarizationSensitvity2=linear2;
110 fCircularPolarizationSensitvity=circular;
111 }
112
113 //....oooOO0OOooo........oooOO0OOooo........oo
114
115 #endif
116