Geant4 10.7.0
Toolkit for the simulation of the passage of particles through matter
Loading...
Searching...
No Matches
Public Member Functions | List of all members
G4eeToTwoGammaModel Class Reference

#include <G4eeToTwoGammaModel.hh>

+ Inheritance diagram for G4eeToTwoGammaModel:

Public Member Functions

 G4eeToTwoGammaModel (const G4ParticleDefinition *p=0, const G4String &nam="eplus2gg")
 
virtual ~G4eeToTwoGammaModel ()
 
virtual void Initialise (const G4ParticleDefinition *, const G4DataVector &) override
 
virtual G4double ComputeCrossSectionPerElectron (G4double kinEnergy)
 
virtual G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0., G4double maxEnergy=DBL_MAX) override
 
virtual G4double CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX) override
 
virtual void SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin, G4double maxEnergy) override
 
G4eeToTwoGammaModel & operator= (const G4eeToTwoGammaModel &right)=delete
 
 G4eeToTwoGammaModel (const G4eeToTwoGammaModel &)=delete
 
- Public Member Functions inherited from G4VEmModel
 G4VEmModel (const G4String &nam)
 
virtual ~G4VEmModel ()
 
virtual void Initialise (const G4ParticleDefinition *, const G4DataVector &)=0
 
virtual void SampleSecondaries (std::vector< G4DynamicParticle * > *, const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double tmin=0.0, G4double tmax=DBL_MAX)=0
 
virtual void InitialiseLocal (const G4ParticleDefinition *, G4VEmModel *masterModel)
 
virtual void InitialiseForMaterial (const G4ParticleDefinition *, const G4Material *)
 
virtual void InitialiseForElement (const G4ParticleDefinition *, G4int Z)
 
virtual G4double ComputeDEDXPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
 
virtual G4double CrossSectionPerVolume (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double GetPartialCrossSection (const G4Material *, G4int level, const G4ParticleDefinition *, G4double kineticEnergy)
 
virtual G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, G4double kinEnergy, G4double Z, G4double A=0., G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double ComputeCrossSectionPerShell (const G4ParticleDefinition *, G4int Z, G4int shellIdx, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
virtual G4double ChargeSquareRatio (const G4Track &)
 
virtual G4double GetChargeSquareRatio (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual G4double GetParticleCharge (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual void StartTracking (G4Track *)
 
virtual void CorrectionsAlongStep (const G4MaterialCutsCouple *, const G4DynamicParticle *, G4double &eloss, G4double &niel, G4double length)
 
virtual G4double Value (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy)
 
virtual G4double MinPrimaryEnergy (const G4Material *, const G4ParticleDefinition *, G4double cut=0.0)
 
virtual G4double MinEnergyCut (const G4ParticleDefinition *, const G4MaterialCutsCouple *)
 
virtual void SetupForMaterial (const G4ParticleDefinition *, const G4Material *, G4double kineticEnergy)
 
virtual void DefineForRegion (const G4Region *)
 
virtual void ModelDescription (std::ostream &outFile) const
 
void InitialiseElementSelectors (const G4ParticleDefinition *, const G4DataVector &)
 
std::vector< G4EmElementSelector * > * GetElementSelectors ()
 
void SetElementSelectors (std::vector< G4EmElementSelector * > *)
 
virtual G4double ComputeDEDX (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=DBL_MAX)
 
G4double CrossSection (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4double ComputeMeanFreePath (const G4ParticleDefinition *, G4double kineticEnergy, const G4Material *, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4double ComputeCrossSectionPerAtom (const G4ParticleDefinition *, const G4Element *, G4double kinEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
const G4Element * SelectRandomAtom (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
const G4Element * SelectTargetAtom (const G4MaterialCutsCouple *, const G4ParticleDefinition *, G4double kineticEnergy, G4double logKineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
const G4Element * SelectRandomAtom (const G4Material *, const G4ParticleDefinition *, G4double kineticEnergy, G4double cutEnergy=0.0, G4double maxEnergy=DBL_MAX)
 
G4int SelectRandomAtomNumber (const G4Material *)
 
G4int SelectIsotopeNumber (const G4Element *)
 
void SetParticleChange (G4VParticleChange *, G4VEmFluctuationModel *f=nullptr)
 
void SetCrossSectionTable (G4PhysicsTable *, G4bool isLocal)
 
G4ElementData * GetElementData ()
 
G4PhysicsTable * GetCrossSectionTable ()
 
G4VEmFluctuationModel * GetModelOfFluctuations ()
 
G4VEmAngularDistribution * GetAngularDistribution ()
 
G4VEmModel * GetTripletModel ()
 
void SetTripletModel (G4VEmModel *)
 
void SetAngularDistribution (G4VEmAngularDistribution *)
 
G4double HighEnergyLimit () const
 
G4double LowEnergyLimit () const
 
G4double HighEnergyActivationLimit () const
 
G4double LowEnergyActivationLimit () const
 
G4double PolarAngleLimit () const
 
G4double SecondaryThreshold () const
 
G4bool LPMFlag () const
 
G4bool DeexcitationFlag () const
 
G4bool ForceBuildTableFlag () const
 
G4bool UseAngularGeneratorFlag () const
 
void SetAngularGeneratorFlag (G4bool)
 
void SetHighEnergyLimit (G4double)
 
void SetLowEnergyLimit (G4double)
 
void SetActivationHighEnergyLimit (G4double)
 
void SetActivationLowEnergyLimit (G4double)
 
G4bool IsActive (G4double kinEnergy) const
 
void SetPolarAngleLimit (G4double)
 
void SetSecondaryThreshold (G4double)
 
void SetLPMFlag (G4bool val)
 
void SetDeexcitationFlag (G4bool val)
 
void SetForceBuildTable (G4bool val)
 
void SetFluctuationFlag (G4bool val)
 
void SetMasterThread (G4bool val)
 
G4bool IsMaster () const
 
void SetUseBaseMaterials (G4bool val)
 
G4bool UseBaseMaterials () const
 
G4double MaxSecondaryKinEnergy (const G4DynamicParticle *dynParticle)
 
const G4String & GetName () const
 
void SetCurrentCouple (const G4MaterialCutsCouple *)
 
const G4Element * GetCurrentElement () const
 
const G4Isotope * GetCurrentIsotope () const
 
G4bool IsLocked () const
 
void SetLocked (G4bool)
 
G4VEmModel & operator= (const G4VEmModel &right)=delete
 
 G4VEmModel (const G4VEmModel &)=delete
 

Additional Inherited Members

- Protected Member Functions inherited from G4VEmModel
G4ParticleChangeForLoss * GetParticleChangeForLoss ()
 
G4ParticleChangeForGamma * GetParticleChangeForGamma ()
 
virtual G4double MaxSecondaryEnergy (const G4ParticleDefinition *, G4double kineticEnergy)
 
const G4MaterialCutsCouple * CurrentCouple () const
 
void SetCurrentElement (const G4Element *)
 
- Protected Attributes inherited from G4VEmModel
G4ElementData * fElementData
 
G4VParticleChange * pParticleChange
 
G4PhysicsTable * xSectionTable
 
const G4Material * pBaseMaterial
 
const std::vector< G4double > * theDensityFactor
 
const std::vector< G4int > * theDensityIdx
 
size_t idxTable
 
G4bool lossFlucFlag
 
G4double inveplus
 
G4double pFactor
 

Detailed Description

Definition at line 59 of file G4eeToTwoGammaModel.hh.

Constructor & Destructor Documentation

◆ G4eeToTwoGammaModel() [1/2]

G4eeToTwoGammaModel::G4eeToTwoGammaModel ( const G4ParticleDefinition *  p = 0,
const G4String &  nam = "eplus2gg" 
)
explicit

Definition at line 90 of file G4eeToTwoGammaModel.cc.

92 : G4VEmModel(nam),
93 pi_rcl2(pi*classic_electr_radius*classic_electr_radius)
94{
95 theGamma = G4Gamma::Gamma();
96 fParticleChange = nullptr;
97}
G4Gamma::Gamma
static G4Gamma * Gamma()
Definition: G4Gamma.cc:85

◆ ~G4eeToTwoGammaModel()

G4eeToTwoGammaModel::~G4eeToTwoGammaModel ( )
virtual

Definition at line 101 of file G4eeToTwoGammaModel.cc.

102{}

◆ G4eeToTwoGammaModel() [2/2]

G4eeToTwoGammaModel::G4eeToTwoGammaModel ( const G4eeToTwoGammaModel &  )
delete

Member Function Documentation

◆ ComputeCrossSectionPerAtom()

G4double G4eeToTwoGammaModel::ComputeCrossSectionPerAtom ( const G4ParticleDefinition *  ,
G4double  kinEnergy,
G4double  Z,
G4double  A = 0.,
G4double  cutEnergy = 0.,
G4double  maxEnergy = DBL_MAX 
)
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 156 of file G4eeToTwoGammaModel.cc.

160{
161 // Calculates the cross section per atom of annihilation into two photons
162 return Z*ComputeCrossSectionPerElectron(kineticEnergy);
163}
G4eeToTwoGammaModel::ComputeCrossSectionPerElectron
virtual G4double ComputeCrossSectionPerElectron(G4double kinEnergy)
Definition: G4eeToTwoGammaModel.cc:136

◆ ComputeCrossSectionPerElectron()

G4double G4eeToTwoGammaModel::ComputeCrossSectionPerElectron ( G4double  kinEnergy)
virtual

Reimplemented in G4PolarizedAnnihilationModel.

Definition at line 136 of file G4eeToTwoGammaModel.cc.

137{
138 // Calculates the cross section per electron of annihilation into two photons
139 // from the Heilter formula.
140
141 G4double ekin = std::max(eV,kineticEnergy);
142
143 G4double tau = ekin/electron_mass_c2;
144 G4double gam = tau + 1.0;
145 G4double gamma2= gam*gam;
146 G4double bg2 = tau * (tau+2.0);
147 G4double bg = sqrt(bg2);
148
149 G4double cross = pi_rcl2*((gamma2+4*gam+1.)*G4Log(gam+bg) - (gam+3.)*bg)
150 / (bg2*(gam+1.));
151 return cross;
152}
G4Log
G4double G4Log(G4double x)
Definition: G4Log.hh:226
G4double
double G4double
Definition: G4Types.hh:83

Referenced by ComputeCrossSectionPerAtom(), G4PolarizedAnnihilationModel::ComputeCrossSectionPerElectron(), and CrossSectionPerVolume().

◆ CrossSectionPerVolume()

G4double G4eeToTwoGammaModel::CrossSectionPerVolume ( const G4Material *  material,
const G4ParticleDefinition *  ,
G4double  kineticEnergy,
G4double  cutEnergy = 0.0,
G4double  maxEnergy = DBL_MAX 
)
overridevirtual

Reimplemented from G4VEmModel.

Definition at line 167 of file G4eeToTwoGammaModel.cc.

172{
173 // Calculates the cross section per volume of annihilation into two photons
174 return material->GetElectronDensity()*ComputeCrossSectionPerElectron(kineticEnergy);
175}
G4Material::GetElectronDensity
G4double GetElectronDensity() const
Definition: G4Material.hh:215

◆ Initialise()

void G4eeToTwoGammaModel::Initialise ( const G4ParticleDefinition *  ,
const G4DataVector &   
)
overridevirtual

Implements G4VEmModel.

Reimplemented in G4PolarizedAnnihilationModel.

Definition at line 106 of file G4eeToTwoGammaModel.cc.

108{
109 if(IsMaster()) {
110 G4int verbose = G4EmParameters::Instance()->Verbose();
111 // redo initialisation for each new run
112 fSampleAtomicPDF = false;
113 const auto& materialTable = G4Material::GetMaterialTable();
114 for (const auto& material: *materialTable) {
115 const G4double meanEnergyPerIonPair = material->GetIonisation()->GetMeanEnergyPerIonPair();
116 if (meanEnergyPerIonPair > 0.) {
117 fSampleAtomicPDF = true;
118 if(verbose > 0) {
119 G4cout << "### G4eeToTwoGammaModel: for " << material->GetName() << " mean energy per ion pair is "
120 << meanEnergyPerIonPair/CLHEP::eV << " eV" << G4endl;
121 }
122 }
123 }
124 }
125 // If no materials have meanEnergyPerIonPair set. This is probably the usual
126 // case, since most applications are not senstive to the slight
127 // non-collinearity of gammas in eeToTwoGamma. Do not issue any warning.
128
129 if(fParticleChange) { return; }
130 fParticleChange = GetParticleChangeForGamma();
131}
G4int
int G4int
Definition: G4Types.hh:85
G4endl
#define G4endl
Definition: G4ios.hh:57
G4cout
G4GLOB_DLL std::ostream G4cout
G4EmParameters::Instance
static G4EmParameters * Instance()
Definition: G4EmParameters.cc:67
G4EmParameters::Verbose
G4int Verbose() const
Definition: G4EmParameters.cc:975
G4Material::GetMaterialTable
static G4MaterialTable * GetMaterialTable()
Definition: G4Material.cc:637
G4VEmModel::GetParticleChangeForGamma
G4ParticleChangeForGamma * GetParticleChangeForGamma()
Definition: G4VEmModel.cc:133
G4VEmModel::IsMaster
G4bool IsMaster() const
Definition: G4VEmModel.hh:736

Referenced by G4PolarizedAnnihilationModel::Initialise().

◆ operator=()

G4eeToTwoGammaModel & G4eeToTwoGammaModel::operator= ( const G4eeToTwoGammaModel &  right)
delete

◆ SampleSecondaries()

void G4eeToTwoGammaModel::SampleSecondaries ( std::vector< G4DynamicParticle * > *  vdp,
const G4MaterialCutsCouple *  pCutsCouple,
const G4DynamicParticle *  dp,
G4double  tmin,
G4double  maxEnergy 
)
overridevirtual

!! likely problematic direction to be checked

Implements G4VEmModel.

Reimplemented in G4PolarizedAnnihilationModel.

Definition at line 182 of file G4eeToTwoGammaModel.cc.

187{
188 G4double posiKinEnergy = dp->GetKineticEnergy();
189 G4DynamicParticle *aGamma1, *aGamma2;
190
191 CLHEP::HepRandomEngine* rndmEngine = G4Random::getTheEngine();
192
193 // Case at rest
194 if(posiKinEnergy == 0.0) {
195
196 const G4double eGamma = electron_mass_c2;
197
198 // In rest frame of positronium gammas are back to back
199 const G4ThreeVector& dir1 = G4RandomDirection();
200 const G4ThreeVector& dir2 = -dir1;
201 aGamma1 = new G4DynamicParticle(G4Gamma::Gamma(),dir1,eGamma);
202 aGamma2 = new G4DynamicParticle(G4Gamma::Gamma(),dir2,eGamma);
203
204 // In rest frame the gammas are polarised perpendicular to each other - see
205 // Pryce and Ward, Nature No 4065 (1947) p.435.
206 // Snyder et al, Physical Review 73 (1948) p.440.
207 G4ThreeVector pol1 = (G4RandomDirection().cross(dir1)).unit();
208 G4ThreeVector pol2 = (pol1.cross(dir2)).unit();
209
210 // But the positronium is moving...
211 // A positron in matter slows down and combines with an atomic electron to
212 // make a neutral “atom” called positronium, about half the size of a normal
213 // atom. I expect that when the energy of the positron is small enough,
214 // less than the binding energy of positronium (6.8 eV), it is
215 // energetically favourable for an electron from the outer orbitals of a
216 // nearby atom or molecule to transfer and bind to the positron, as in an
217 // ionic bond, leaving behind a mildly ionised nearby atom/molecule. I
218 // would expect the positronium to come away with a kinetic energy of a
219 // few eV on average. In its para (spin 0) state it annihilates into two
220 // photons, which in the rest frame of the positronium are collinear
221 // (back-to-back) due to momentum conservation. Because of the motion of the
222 // positronium, photons will be not quite back-to-back in the laboratory.
223
224 // The positroniuim acquires an energy of order its binding energy and
225 // doesn't have time to thermalise. Nevertheless, here we approximate its
226 // energy distribution by a Maxwell-Boltzman with mean energy <KE>. In terms
227 // of a more familiar concept of temperature, and the law of equipartition
228 // of energy of translational motion, <KE>=3kT/2. Each component of velocity
229 // has a distribution exp(-mv^2/2kT), which is a Gaussian of mean zero
230 // and variance kT/m=2<KE>/3m, where m is the positronium mass.
231
232 // We take <KE> = material->GetIonisation()->GetMeanEnergyPerIonPair().
233
234 if(fSampleAtomicPDF) {
235 const G4Material* material = pCutsCouple->GetMaterial();
236 const G4double meanEnergyPerIonPair = material->GetIonisation()->GetMeanEnergyPerIonPair();
237 const G4double& meanKE = meanEnergyPerIonPair; // Just an alias
238 if (meanKE > 0.) { // Positronium haas motion
239 // Mass of positronium
240 const G4double mass = 2.*electron_mass_c2;
241 // Mean <KE>=3kT/2, as described above
242 // const G4double T = 2.*meanKE/(3.*k_Boltzmann);
243 // Component velocities: Gaussian, variance kT/m=2<KE>/3m.
244 const G4double sigmav = std::sqrt(2.*meanKE/(3.*mass));
245 // This is in units where c=1
246 const G4double vx = G4RandGauss::shoot(0.,sigmav);
247 const G4double vy = G4RandGauss::shoot(0.,sigmav);
248 const G4double vz = G4RandGauss::shoot(0.,sigmav);
249 const G4ThreeVector v(vx,vy,vz); // In unit where c=1
250 const G4ThreeVector& beta = v; // so beta=v/c=v
251
252 aGamma1->Set4Momentum(aGamma1->Get4Momentum().boost(beta));
253 aGamma2->Set4Momentum(aGamma2->Get4Momentum().boost(beta));
254
255 // Rotate polarisation vectors
256 const G4ThreeVector& newDir1 = aGamma1->GetMomentumDirection();
257 const G4ThreeVector& newDir2 = aGamma2->GetMomentumDirection();
258 const G4ThreeVector& axis1 = dir1.cross(newDir1); // No need to be unit
259 const G4ThreeVector& axis2 = dir2.cross(newDir2); // No need to be unit
260 const G4double& angle1 = std::acos(dir1*newDir1);
261 const G4double& angle2 = std::acos(dir2*newDir2);
262 if (axis1 != G4ThreeVector()) pol1.rotate(axis1,angle1);
263 if (axis2 != G4ThreeVector()) pol2.rotate(axis2,angle2);
264 }
265 }
266 aGamma1->SetPolarization(pol1.x(),pol1.y(),pol1.z());
267 aGamma2->SetPolarization(pol2.x(),pol2.y(),pol2.z());
268
269 } else { // Positron interacts in flight
270
271 G4ThreeVector posiDirection = dp->GetMomentumDirection();
272
273 G4double tau = posiKinEnergy/electron_mass_c2;
274 G4double gam = tau + 1.0;
275 G4double tau2 = tau + 2.0;
276 G4double sqgrate = sqrt(tau/tau2)*0.5;
277 G4double sqg2m1 = sqrt(tau*tau2);
278
279 // limits of the energy sampling
280 G4double epsilmin = 0.5 - sqgrate;
281 G4double epsilmax = 0.5 + sqgrate;
282 G4double epsilqot = epsilmax/epsilmin;
283
284 //
285 // sample the energy rate of the created gammas
286 //
287 G4double epsil, greject;
288
289 do {
290 epsil = epsilmin*G4Exp(G4Log(epsilqot)*rndmEngine->flat());
291 greject = 1. - epsil + (2.*gam*epsil-1.)/(epsil*tau2*tau2);
292 // Loop checking, 03-Aug-2015, Vladimir Ivanchenko
293 } while( greject < rndmEngine->flat());
294
295 //
296 // scattered Gamma angles. ( Z - axis along the parent positron)
297 //
298
299 G4double cost = (epsil*tau2-1.)/(epsil*sqg2m1);
300 if(std::abs(cost) > 1.0) {
301 G4cout << "### G4eeToTwoGammaModel WARNING cost= " << cost
302 << " positron Ekin(MeV)= " << posiKinEnergy
303 << " gamma epsil= " << epsil
304 << G4endl;
305 if(cost > 1.0) cost = 1.0;
306 else cost = -1.0;
307 }
308 G4double sint = sqrt((1.+cost)*(1.-cost));
309 G4double phi = twopi * rndmEngine->flat();
310
311 //
312 // kinematic of the created pair
313 //
314
315 G4double totalEnergy = posiKinEnergy + 2.0*electron_mass_c2;
316 G4double phot1Energy = epsil*totalEnergy;
317
318 G4ThreeVector phot1Direction(sint*cos(phi), sint*sin(phi), cost);
319 phot1Direction.rotateUz(posiDirection);
320 aGamma1 = new G4DynamicParticle (theGamma,phot1Direction, phot1Energy);
321 phi = twopi * rndmEngine->flat();
322 G4double cosphi = cos(phi);
323 G4double sinphi = sin(phi);
324 G4ThreeVector pol(cosphi, sinphi, 0.0);
325 pol.rotateUz(phot1Direction);
326 aGamma1->SetPolarization(pol.x(),pol.y(),pol.z());
327
328 G4double phot2Energy =(1.-epsil)*totalEnergy;
329 G4double posiP= sqrt(posiKinEnergy*(posiKinEnergy+2.*electron_mass_c2));
330 G4ThreeVector dir = posiDirection*posiP - phot1Direction*phot1Energy;
331 G4ThreeVector phot2Direction = dir.unit();
332
333 // create G4DynamicParticle object for the particle2
334 aGamma2 = new G4DynamicParticle (theGamma, phot2Direction, phot2Energy);
335
336 //!!! likely problematic direction to be checked
337 pol.set(-sinphi, cosphi, 0.0);
338 pol.rotateUz(phot1Direction);
339 cost = pol*phot2Direction;
340 pol -= cost*phot2Direction;
341 pol = pol.unit();
342 aGamma2->SetPolarization(pol.x(),pol.y(),pol.z());
343 /*
344 G4cout << "Annihilation on fly: e0= " << posiKinEnergy
345 << " m= " << electron_mass_c2
346 << " e1= " << phot1Energy
347 << " e2= " << phot2Energy << " dir= " << dir
348 << " -> " << phot1Direction << " "
349 << phot2Direction << G4endl;
350 */
351 }
352
353 vdp->push_back(aGamma1);
354 vdp->push_back(aGamma2);
355
356 // kill primary positron
357 fParticleChange->SetProposedKineticEnergy(0.0);
358 fParticleChange->ProposeTrackStatus(fStopAndKill);
359}
G4Exp
G4double G4Exp(G4double initial_x)
Exponential Function double precision.
Definition: G4Exp.hh:179
G4RandomDirection
G4ThreeVector G4RandomDirection()
Definition: G4RandomDirection.hh:58
G4ThreeVector
CLHEP::Hep3Vector G4ThreeVector
Definition: G4ThreeVector.hh:36
fStopAndKill
@ fStopAndKill
Definition: G4TrackStatus.hh:46
CLHEP::Hep3Vector::z
double z() const
CLHEP::Hep3Vector::unit
Hep3Vector unit() const
CLHEP::Hep3Vector::x
double x() const
CLHEP::Hep3Vector::y
double y() const
CLHEP::Hep3Vector::cross
Hep3Vector cross(const Hep3Vector &) const
CLHEP::Hep3Vector::rotate
Hep3Vector & rotate(double, const Hep3Vector &)
Definition: ThreeVectorR.cc:24
CLHEP::HepLorentzVector::boost
HepLorentzVector & boost(double, double, double)
Definition: LorentzVector.cc:55
CLHEP::HepRandomEngine::flat
virtual double flat()=0
G4DynamicParticle::SetPolarization
void SetPolarization(const G4ThreeVector &)
G4DynamicParticle::GetMomentumDirection
const G4ThreeVector & GetMomentumDirection() const
G4DynamicParticle::Get4Momentum
G4LorentzVector Get4Momentum() const
G4DynamicParticle::GetKineticEnergy
G4double GetKineticEnergy() const
G4DynamicParticle::Set4Momentum
void Set4Momentum(const G4LorentzVector &momentum)
Definition: G4DynamicParticle.cc:410
G4IonisParamMat::GetMeanEnergyPerIonPair
G4double GetMeanEnergyPerIonPair() const
Definition: G4IonisParamMat.hh:158
G4MaterialCutsCouple::GetMaterial
const G4Material * GetMaterial() const
Definition: G4MaterialCutsCouple.hh:166
G4Material::GetIonisation
G4IonisParamMat * GetIonisation() const
Definition: G4Material.hh:224
G4ParticleChangeForGamma::SetProposedKineticEnergy
void SetProposedKineticEnergy(G4double proposedKinEnergy)
Definition: G4ParticleChangeForGamma.hh:119
G4VParticleChange::ProposeTrackStatus
void ProposeTrackStatus(G4TrackStatus status)
G4AblaRandom::flat
double flat()
Definition: G4AblaRandom.cc:48
The documentation for this class was generated from the following files: