| 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 // G4DynamicParticle
27 //
28 // Class description:
29 //
30 // A G4DynamicParticle aggregates the informat
31 // of a G4Particle, such as energy, momentum,
32 // as well as the "particle definition", holdi
33 // It contains the purely dynamic aspects of a
34
35 // History:
36 // - 2 December 1995, G.Cosmo - first design,
37 // - 29 January 1996, M.Asai - first implement
38 // - 1996 - 2007, H.Kurashige - revisions.
39 // - 15 March 2019, M.Novak - log-kinetic en
40 // on demand if its stored
41 // -------------------------------------------
42 #ifndef G4DynamicParticle_hh
43 #define G4DynamicParticle_hh 1
44
45 #include "G4Allocator.hh"
46 #include "G4ElectronOccupancy.hh"
47 #include "G4Log.hh"
48 #include "G4LorentzVector.hh"
49 #include "G4ParticleDefinition.hh"
50 #include "G4ParticleMomentum.hh" // NOTE: mea
51 #include "G4ios.hh"
52 #include "globals.hh"
53
54 #include <CLHEP/Units/PhysicalConstants.h>
55 #include <CLHEP/Units/SystemOfUnits.h>
56
57 #include <cmath>
58
59 class G4PrimaryParticle;
60 class G4DecayProducts;
61
62 class G4DynamicParticle
63 {
64 public:
65 //- constructors
66 G4DynamicParticle();
67
68 G4DynamicParticle(const G4ParticleDefiniti
69 const G4ThreeVector& aMo
70 G4DynamicParticle(const G4ParticleDefiniti
71 const G4ThreeVector& aPa
72 G4DynamicParticle(const G4ParticleDefiniti
73 const G4LorentzVector& a
74 G4DynamicParticle(const G4ParticleDefiniti
75 const G4ThreeVector& aPa
76 G4DynamicParticle(const G4ParticleDefiniti
77 const G4ThreeVector& aMo
78 const G4double dynamical
79
80 G4DynamicParticle(const G4DynamicParticle&
81
82 //- destructor
83 ~G4DynamicParticle();
84
85 //- operators
86 G4DynamicParticle& operator=(const G4Dynam
87 G4bool operator==(const G4DynamicParticle&
88 G4bool operator!=(const G4DynamicParticle&
89
90 //- Move constructor & operator
91 G4DynamicParticle(G4DynamicParticle&& from
92 G4DynamicParticle& operator=(G4DynamicPart
93
94 //- new/delete operators are oberloded to
95 inline void* operator new(size_t);
96 inline void operator delete(void* aDynamic
97
98 //- Set/Get methods
99
100 // Returns the normalized direction of the
101 inline const G4ThreeVector& GetMomentumDir
102
103 // Sets the normalized direction of the mo
104 inline void SetMomentumDirection(const G4T
105
106 // Sets the normalized direction of the mo
107 inline void SetMomentumDirection(G4double
108
109 // Returns the current particle momentum v
110 inline G4ThreeVector GetMomentum() const;
111
112 // set the current particle momentum vecto
113 void SetMomentum(const G4ThreeVector& mome
114
115 // Returns the current particle energy-mom
116 inline G4LorentzVector Get4Momentum() cons
117
118 // Set the current particle energy-momentu
119 void Set4Momentum(const G4LorentzVector& m
120
121 // Returns the module of the momentum vect
122 inline G4double GetTotalMomentum() const;
123
124 // Returns the total energy of the particl
125 inline G4double GetTotalEnergy() const;
126
127 // Returns the kinetic energy of a particl
128 inline G4double GetKineticEnergy() const;
129
130 // Returns:
131 // - natural logarithm of the particle kin
132 // - LOG_EKIN_MIN otherwise
133 inline G4double GetLogKineticEnergy() cons
134
135 // Sets the kinetic energy of a particle
136 inline void SetKineticEnergy(G4double aEne
137
138 // Access Lorentz beta
139 inline G4double GetBeta() const;
140
141 // Returns the current particle proper tim
142 inline G4double GetProperTime() const;
143
144 // Set the current particle Proper Time
145 inline void SetProperTime(G4double);
146
147 // Set/Get polarization vector
148 inline const G4ThreeVector& GetPolarizatio
149 inline void SetPolarization(const G4ThreeV
150 inline void SetPolarization(G4double polX,
151
152 // Set/Get dynamical mass
153 // The dynamical mass is set to PDG mass i
154 inline G4double GetMass() const;
155 inline void SetMass(G4double mass);
156
157 // Set/Get dynamical charge
158 // The dynamical mass is set to PDG charge
159 inline G4double GetCharge() const;
160 inline void SetCharge(G4double charge);
161 inline void SetCharge(G4int chargeInUnitOf
162
163 // Set/Get dynamical spin
164 // The dynamical spin is set to PDG spin i
165 inline G4double GetSpin() const;
166 inline void SetSpin(G4double spin);
167 inline void SetSpin(G4int spinInUnitOfHalf
168
169 // Set/Get dynamical MagneticMoment
170 // The dynamical mass is set to PDG Magnet
171 inline G4double GetMagneticMoment() const;
172 inline void SetMagneticMoment(G4double mag
173
174 // Get electron occupancy
175 // ElectronOccupancy is valid only if the
176 inline const G4ElectronOccupancy* GetElect
177 inline G4int GetTotalOccupancy() const;
178 inline G4int GetOccupancy(G4int orbit) con
179 inline void AddElectron(G4int orbit, G4int
180 inline void RemoveElectron(G4int orbit, G4
181
182 // Set/Get particle definition
183 inline const G4ParticleDefinition* GetPart
184 void SetDefinition(const G4ParticleDefinit
185
186 // Following method of GetDefinition() rem
187 // because of backward compatiblity. May b
188 inline G4ParticleDefinition* GetDefinition
189
190 // Set/Get pre-assigned decay channel
191 inline const G4DecayProducts* GetPreAssign
192 inline void SetPreAssignedDecayProducts(G4
193
194 // Set/Get pre-assigned proper time when t
195 inline G4double GetPreAssignedDecayProperT
196 inline void SetPreAssignedDecayProperTime(
197
198 // Print out information
199 // - mode 0 : default )(minimum)
200 // - mode 1 : 0 + electron occupancy
201 void DumpInfo(G4int mode = 0) const;
202
203 // Set/Get controle flag for output messag
204 // - 0: Silent
205 // - 1: Warning message
206 // - 2: More
207 inline void SetVerboseLevel(G4int value);
208 inline G4int GetVerboseLevel() const;
209
210 inline void SetPrimaryParticle(G4PrimaryPa
211 inline void SetPDGcode(G4int c);
212
213 // Return the pointer to the corresponding
214 // if this particle is a primary particle
215 // pre-assigned decay product. Otherwise r
216 inline G4PrimaryParticle* GetPrimaryPartic
217
218 // Return the PDG code of this particle. I
219 // Geant4, its PDG code defined in G4Parti
220 // If it is unknown (i.e. PDG code in G4Pa
221 // PDG code defined in the corresponding p
222 // pre-assigned decay product will be retu
223 // Otherwise (e.g. for geantino) returns 0
224 inline G4int GetPDGcode() const;
225
226 protected:
227 void AllocateElectronOccupancy();
228 G4double GetElectronMass() const;
229
230 private:
231 inline void ComputeBeta() const;
232
233 // The normalized momentum vector
234 G4ThreeVector theMomentumDirection;
235
236 G4ThreeVector thePolarization;
237
238 // Contains the static information of this
239 const G4ParticleDefinition* theParticleDef
240
241 G4ElectronOccupancy* theElectronOccupancy
242
243 G4DecayProducts* thePreAssignedDecayProduc
244
245 // This void pointer is used by G4EventMan
246 // link between pre-assigned decay product
247 // primary particle
248 G4PrimaryParticle* primaryParticle = nullp
249
250 G4double theKineticEnergy = 0.0;
251
252 mutable G4double theLogKineticEnergy = DBL
253
254 mutable G4double theBeta = -1.0;
255
256 G4double theProperTime = 0.0;
257
258 G4double theDynamicalMass = 0.0;
259
260 G4double theDynamicalCharge = 0.0;
261
262 G4double theDynamicalSpin = 0.0;
263
264 G4double theDynamicalMagneticMoment = 0.0;
265
266 G4double thePreAssignedDecayTime = -1.0;
267
268 G4int verboseLevel = 1;
269
270 G4int thePDGcode = 0;
271 };
272
273 #include "G4DynamicParticle.icc"
274
275 #endif
276