| Geant4 Cross Reference |
1 // 1 //
2 // ******************************************* 2 // ********************************************************************
3 // * License and Disclaimer 3 // * License and Disclaimer *
4 // * 4 // * *
5 // * The Geant4 software is copyright of th 5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided 6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License 7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/ 8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. 9 // * include a list of copyright holders. *
10 // * 10 // * *
11 // * Neither the authors of this software syst 11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing fin 12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warran 13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assum 14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file 15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitatio 16 // * for the full disclaimer and the limitation of liability. *
17 // * 17 // * *
18 // * This code implementation is the result 18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboratio 19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distri 20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you ag 21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publicati 22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Sof 23 // * acceptance of all terms of the Geant4 Software license. *
24 // ******************************************* 24 // ********************************************************************
25 // 25 //
26 // G4SPSRandomGenerator << 26 ///////////////////////////////////////////////////////////////////////////////
>> 27 //
>> 28 // MODULE: G4SPSRandomGenerator.hh
>> 29 //
>> 30 // Version: 1.0
>> 31 // Date: 5/02/04
>> 32 // Author: Fan Lei
>> 33 // Organisation: QinetiQ ltd.
>> 34 // Customer: ESA/ESTEC
>> 35 //
>> 36 ///////////////////////////////////////////////////////////////////////////////
>> 37 //
>> 38 // CHANGE HISTORY
>> 39 // --------------
>> 40 //
>> 41 // 26/10/2004 F Lei
>> 42 // Created separated the theta, phi generators for position distributions.
>> 43 //
>> 44 // Version 1.0, 05/02/2004, Fan Lei, Created.
>> 45 // Based on the G4GeneralParticleSource class in Geant4 v6.0
>> 46 //
>> 47 ///////////////////////////////////////////////////////////////////////////////
27 // 48 //
28 // Class Description: 49 // Class Description:
29 // 50 //
30 // Special random number generator used by G4G 51 // Special random number generator used by G4GeneralParticleSource to allow
31 // biasing applied at the lowest level for all 52 // biasing applied at the lowest level for all distributions.
32 // This is a shared class between threads. << 53 //
33 // Only one thread should use the set-methods << 54 ///////////////////////////////////////////////////////////////////////////////
34 // Note that this is exactly what is achieved << 55 //
35 // If you use the set methods to set defaults << 56 // MEMBER FUNCTIONS
36 // application take care that only one thread << 57 // ----------------
37 // In addition take care of calling these meth << 58 //
38 // Do not use the setters during the event loo << 59 // G4SPSRandomGenerator ()
39 << 60 // Constructor: Initializes variables
40 // Author: Fan Lei, QinetiQ ltd. << 61 //
41 // Customer: ESA/ESTEC << 62 // ~G4SPSRandomGenerator ()
42 // History: << 63 // Destructor:
43 // - 05/02/2004, Fan Lei - Created. << 64 //
44 // Based on the G4GeneralParticleSource cl << 65 // void SetXBias(G4ThreeVector)
45 // - 06/06/2014, Andrea Dotti << 66 // Allows the user to re-distribute the random
46 // Added a mutex to protect access to shar << 67 // numbers used to generate x co-ordinates.
47 // Getters and Setters are mutex'd but not << 68 //
48 // because it is assumed these are called << 69 // void SetYBias(G4ThreeVector)
49 // during which the status of this class i << 70 // Allows the user to re-distribute the random
50 // ------------------------------------------- << 71 // numbers used to generate y co-ordinates.
51 #ifndef G4SPSRandomGenerator_hh << 72 //
52 #define G4SPSRandomGenerator_hh 1 << 73 // void SetZBias(G4ThreeVector)
>> 74 // Allows the user to re-distribute the random
>> 75 // numbers used to generate z co-ordinates.
>> 76 //
>> 77 // void SetThetaBias(G4ThreeVector)
>> 78 // Allows the user to re-distribute the random
>> 79 // numbers used to generate values of theta.
>> 80 //
>> 81 // void SetPhiBias(G4ThreeVector)
>> 82 // Allows the user to re-distribute the random
>> 83 // numbers used to generate values of phi.
>> 84 //
>> 85 // void SetPosThetaBias(G4ThreeVector)
>> 86 // Allows the user to re-distribute the random
>> 87 // numbers used to generate values of theta for position distribution.
>> 88 //
>> 89 // void SetPosPhiBias(G4ThreeVector)
>> 90 // Allows the user to re-distribute the random
>> 91 // numbers used to generate values of phi for position distribution.
>> 92 //
>> 93 // void SetEnergyBias(G4ThreeVector)
>> 94 // Allows the user to re-distribute the random
>> 95 // numbers used to generate the energies.
>> 96 //
>> 97 // G4double GenRandX()
>> 98 // Generates the random number for x, with or without biasing.
>> 99 //
>> 100 // G4double GenRandY()
>> 101 // Generates the random number for y, with or without biasing.
>> 102 //
>> 103 // G4double GenRandZ()
>> 104 // Generates the random number for z, with or without biasing.
>> 105 //
>> 106 // G4double GenRandTheta()
>> 107 // Generates the random number for theta, with or without biasing.
>> 108 //
>> 109 // G4double GenRandPhi()
>> 110 // Generates the random number for phi, with or without biasing.
>> 111 //
>> 112 // G4double GenRandEnergy()
>> 113 // Generates the random number for energy, with or without biasing.
>> 114 //
>> 115 // G4double GenRandPosTheta()
>> 116 // Generates the random number for theta, with or without biasing for position distribution.
>> 117 //
>> 118 // G4double GenRandPosPhi()
>> 119 // Generates the random number for phi, with or without biasing for position distribution.
>> 120 //
>> 121 // inline G4double GetBiasWeight()
>> 122 // Returns the weight change after biasing
>> 123 //
>> 124 // void ReSetHist(G4String);
>> 125 // Re-sets the histogram for user defined distribution
>> 126 //
>> 127 // void SetVerbosity(G4int)
>> 128 // Sets the verbosity level.
>> 129 //
>> 130 ///////////////////////////////////////////////////////////////////////////////
>> 131 //
>> 132 #ifndef G4SPSRandomGenerator_h
>> 133 #define G4SPSRandomGenerator_h 1
53 134
54 #include "G4PhysicsFreeVector.hh" << 135 #include "G4PhysicsOrderedFreeVector.hh"
55 #include "G4DataInterpolation.hh" 136 #include "G4DataInterpolation.hh"
56 #include "G4ThreeVector.hh" <<
57 #include "G4Threading.hh" <<
58 #include "G4Cache.hh" <<
59 137
60 class G4SPSRandomGenerator 138 class G4SPSRandomGenerator
61 { 139 {
62 public: << 140 public:
>> 141 G4SPSRandomGenerator ();
>> 142 ~G4SPSRandomGenerator ();
>> 143
>> 144 // static G4SPSRandomGenerator* getInstance ();
>> 145
>> 146 // Biasing Methods
>> 147 void SetXBias(G4ThreeVector);
>> 148 void SetYBias(G4ThreeVector);
>> 149 void SetZBias(G4ThreeVector);
>> 150 void SetThetaBias(G4ThreeVector);
>> 151 void SetPhiBias(G4ThreeVector);
>> 152 void SetEnergyBias(G4ThreeVector);
>> 153 void SetPosThetaBias(G4ThreeVector);
>> 154 void SetPosPhiBias(G4ThreeVector);
>> 155 G4double GenRandX();
>> 156 G4double GenRandY();
>> 157 G4double GenRandZ();
>> 158 G4double GenRandTheta();
>> 159 G4double GenRandPhi();
>> 160 G4double GenRandEnergy();
>> 161 G4double GenRandPosTheta();
>> 162 G4double GenRandPosPhi();
63 163
64 G4SPSRandomGenerator(); << 164 inline void SetIntensityWeight(G4double weight) {bweights[8]=weight;};
65 // Constructor: initializes variables <<
66 165
67 ~G4SPSRandomGenerator(); << 166 inline G4double GetBiasWeight()
68 // Destructor << 167 { return bweights[0]*bweights[1]*bweights[2]*bweights[3]*bweights[4]*bweights[5]*bweights[6]*bweights[7]*bweights[8];};
>> 168
>> 169 // method to re-set the histograms
>> 170 void ReSetHist(G4String);
69 171
70 // Biasing Methods << 172 // Set the verbosity level.
71 << 173 void SetVerbosity(G4int a) {verbosityLevel = a; } ;
72 void SetXBias(const G4ThreeVector&); <<
73 // Allows the user to re-distribute the <<
74 // numbers used to generate x co-ordinat <<
75 <<
76 void SetYBias(const G4ThreeVector&); <<
77 // Allows the user to re-distribute the <<
78 // numbers used to generate y co-ordinat <<
79 <<
80 void SetZBias(const G4ThreeVector&); <<
81 // Allows the user to re-distribute the <<
82 // numbers used to generate z co-ordinat <<
83 <<
84 void SetThetaBias(const G4ThreeVector&); <<
85 // Allows the user to re-distribute the <<
86 // numbers used to generate values of th <<
87 <<
88 void SetPhiBias(const G4ThreeVector&); <<
89 // Allows the user to re-distribute the <<
90 // numbers used to generate values of ph <<
91 <<
92 void SetEnergyBias(const G4ThreeVector&); <<
93 // Allows the user to re-distribute the ra <<
94 // numbers used to generate the energies <<
95 <<
96 void SetPosThetaBias(const G4ThreeVector&) <<
97 // Allows the user to re-distribute the <<
98 // numbers used to generate values of th <<
99 <<
100 void SetPosPhiBias(const G4ThreeVector&); <<
101 // Allows the user to re-distribute the <<
102 // numbers used to generate values of ph <<
103 <<
104 G4double GenRandX(); <<
105 // Generates the random number for x, wi <<
106 <<
107 G4double GenRandY(); <<
108 // Generates the random number for y, wi <<
109 <<
110 G4double GenRandZ(); <<
111 // Generates the random number for z, wi <<
112 <<
113 G4double GenRandTheta(); <<
114 // Generates the random number for theta <<
115 <<
116 G4double GenRandPhi(); <<
117 // Generates the random number for phi, <<
118 <<
119 G4double GenRandEnergy(); <<
120 // Generates the random number for energ <<
121 <<
122 G4double GenRandPosTheta(); <<
123 // Generates the random number for theta <<
124 // for position distribution <<
125 <<
126 G4double GenRandPosPhi(); <<
127 // Generates the random number for phi, <<
128 // for position distribution <<
129 <<
130 void SetIntensityWeight(G4double weight); <<
131 <<
132 G4double GetBiasWeight() const ; <<
133 // Returns the weight change after biasi <<
134 <<
135 // method to re-set the histograms <<
136 void ReSetHist(const G4String&); <<
137 // Resets the histogram for user defined <<
138 <<
139 void SetVerbosity(G4int a); <<
140 // Sets the verbosity level <<
141 <<
142 private: <<
143 <<
144 // Encapsulate in a struct to guarantee th <<
145 // initial state is set via constructor <<
146 // <<
147 struct a_check <<
148 { <<
149 G4bool val; <<
150 a_check() { val = false; } <<
151 }; <<
152 <<
153 // See .cc for an explanation of this in m <<
154 // <<
155 G4Cache<a_check> local_IPDFXBias; <<
156 G4bool XBias, IPDFXBias; <<
157 G4PhysicsFreeVector XBiasH; <<
158 G4PhysicsFreeVector IPDFXBiasH; <<
159 G4Cache<a_check> local_IPDFYBias; <<
160 G4bool YBias, IPDFYBias; <<
161 G4PhysicsFreeVector YBiasH; <<
162 G4PhysicsFreeVector IPDFYBiasH; <<
163 G4Cache<a_check> local_IPDFZBias; <<
164 G4bool ZBias, IPDFZBias; <<
165 G4PhysicsFreeVector ZBiasH; <<
166 G4PhysicsFreeVector IPDFZBiasH; <<
167 G4Cache<a_check> local_IPDFThetaBias; <<
168 G4bool ThetaBias, IPDFThetaBias; <<
169 G4PhysicsFreeVector ThetaBiasH; <<
170 G4PhysicsFreeVector IPDFThetaBiasH; <<
171 G4Cache<a_check> local_IPDFPhiBias; <<
172 G4bool PhiBias, IPDFPhiBias; <<
173 G4PhysicsFreeVector PhiBiasH; <<
174 G4PhysicsFreeVector IPDFPhiBiasH; <<
175 G4Cache<a_check> local_IPDFEnergyBias; <<
176 G4bool EnergyBias, IPDFEnergyBias; <<
177 G4PhysicsFreeVector EnergyBiasH; <<
178 G4PhysicsFreeVector IPDFEnergyBiasH; <<
179 G4Cache<a_check> local_IPDFPosThetaBias; <<
180 G4bool PosThetaBias, IPDFPosThetaBias; <<
181 G4PhysicsFreeVector PosThetaBiasH; <<
182 G4PhysicsFreeVector IPDFPosThetaBiasH; <<
183 G4Cache<a_check> local_IPDFPosPhiBias; <<
184 G4bool PosPhiBias, IPDFPosPhiBias; <<
185 G4PhysicsFreeVector PosPhiBiasH; <<
186 G4PhysicsFreeVector IPDFPosPhiBiasH; <<
187 <<
188 struct bweights_t <<
189 { <<
190 G4double w[9]; <<
191 bweights_t(); <<
192 G4double& operator[] (const int i); <<
193 }; <<
194 G4Cache<bweights_t> bweights; <<
195 // record x,y,z,theta,phi,energy,posThet <<
196 174
197 G4int verbosityLevel; << 175 private:
198 // Verbosity << 176
199 << 177 // static G4SPSRandomGenerator *instance;
200 G4Mutex mutex; << 178
201 // Protect shared resources << 179 G4bool XBias, IPDFXBias;
>> 180 G4PhysicsOrderedFreeVector XBiasH;
>> 181 G4PhysicsOrderedFreeVector IPDFXBiasH;
>> 182 G4bool YBias, IPDFYBias;
>> 183 G4PhysicsOrderedFreeVector YBiasH;
>> 184 G4PhysicsOrderedFreeVector IPDFYBiasH;
>> 185 G4bool ZBias, IPDFZBias;
>> 186 G4PhysicsOrderedFreeVector ZBiasH;
>> 187 G4PhysicsOrderedFreeVector IPDFZBiasH;
>> 188 G4bool ThetaBias, IPDFThetaBias;
>> 189 G4PhysicsOrderedFreeVector ThetaBiasH;
>> 190 G4PhysicsOrderedFreeVector IPDFThetaBiasH;
>> 191 G4bool PhiBias, IPDFPhiBias;
>> 192 G4PhysicsOrderedFreeVector PhiBiasH;
>> 193 G4PhysicsOrderedFreeVector IPDFPhiBiasH;
>> 194 G4bool EnergyBias, IPDFEnergyBias;
>> 195 G4PhysicsOrderedFreeVector EnergyBiasH;
>> 196 G4PhysicsOrderedFreeVector IPDFEnergyBiasH;
>> 197 G4bool PosThetaBias, IPDFPosThetaBias;
>> 198 G4PhysicsOrderedFreeVector PosThetaBiasH;
>> 199 G4PhysicsOrderedFreeVector IPDFPosThetaBiasH;
>> 200 G4bool PosPhiBias, IPDFPosPhiBias;
>> 201 G4PhysicsOrderedFreeVector PosPhiBiasH;
>> 202 G4PhysicsOrderedFreeVector IPDFPosPhiBiasH;
>> 203
>> 204 G4double bweights[9]; //record x,y,z,theta,phi,energy,posThet,posPhi,intensity weights
>> 205
>> 206 // Verbosity
>> 207 G4int verbosityLevel;
>> 208
>> 209 G4PhysicsOrderedFreeVector ZeroPhysVector ; // for re-set only
>> 210
202 }; 211 };
203 212
>> 213
204 #endif 214 #endif
>> 215
>> 216
>> 217
>> 218
205 219