Geant4  10.03
GVFlashShowerParameterisation.cc
Go to the documentation of this file.
1 //
2 // ********************************************************************
3 // * License and Disclaimer *
4 // * *
5 // * The Geant4 software is copyright of the Copyright Holders of *
6 // * the Geant4 Collaboration. It is provided under the terms and *
7 // * conditions of the Geant4 Software License, included in the file *
8 // * LICENSE and available at http://cern.ch/geant4/license . These *
9 // * include a list of copyright holders. *
10 // * *
11 // * Neither the authors of this software system, nor their employing *
12 // * institutes,nor the agencies providing financial support for this *
13 // * work make any representation or warranty, express or implied, *
14 // * regarding this software system or assume any liability for its *
15 // * use. Please see the license in the file LICENSE and URL above *
16 // * for the full disclaimer and the limitation of liability. *
17 // * *
18 // * This code implementation is the result of the scientific and *
19 // * technical work of the GEANT4 collaboration. *
20 // * By using, copying, modifying or distributing the software (or *
21 // * any work based on the software) you agree to acknowledge its *
22 // * use in resulting scientific publications, and indicate your *
23 // * acceptance of all terms of the Geant4 Software license. *
24 // ********************************************************************
25 //
26 // $Id: GVFlashShowerParameterisation.cc 100944 2016-11-03 11:18:29Z gcosmo $
27 //
28 //
29 // ------------------------------------------------------------
30 // GEANT 4 class implementation
31 //
32 // ------- GVFlashShowerParameterisation -------
33 //
34 // Authors: Joanna Weng - 11.2005
35 // ------------------------------------------------------------
36 
37 #include <cmath>
38 
40 
41 #include "G4PhysicalConstants.hh"
42 #include "G4SystemOfUnits.hh"
43 #include "Randomize.hh"
44 #include "G4ios.hh"
45 #include "G4Material.hh"
46 #include "Gamma.hh" // @@@@
47 #include "G4MaterialTable.hh"
48 
50  : thePar(0), density(0.), A(0.), Z(0.), X0(0.), Ec(0.), Rm(0.), NSpot(0.)
51 {
52  fGamma = new MyGamma;
53 }
54 
56 {
57  delete fGamma;
58 }
59 
61 {
62  // Returns Z or effective Z=sum(pi*Zi) (if compound/mixture)
63  // of given material
64  //
65  G4double z = 0.;
66  G4int nofElements = mat->GetNumberOfElements();
67  if (nofElements > 1)
68  {
69  for (G4int i=0; i<nofElements; i++) {
70  G4double zOfElement = mat->GetElement(i)->GetZ();
71  G4double massFraction = mat->GetFractionVector()[i];
72  // cout << mat->GetElement(i)->GetName()
73  // <<" Z= "<<zOfElement << " , Fraction= "<<massFraction <<endl;
74  z += zOfElement*massFraction;
75  }
76  }
77  else {
78  z = mat->GetZ();
79  }
80  return z;
81 }
82 
84 {
85  // Returns A or effective A=sum(pi*Ai) (if compound/mixture)
86  // of given material
87  //
88  G4double a = 0.;
89  G4int nofElements = mat->GetNumberOfElements();
90  if (nofElements > 1) {
91  for (G4int i=0; i<nofElements; i++) {
92  G4double aOfElement = mat->GetElement(i)->GetA()/(g/mole);
93  G4double massFraction = mat->GetFractionVector()[i];
94  a += aOfElement*massFraction;
95  }
96  }
97  else {
98  a = mat->GetA()/(g/mole);
99  }
100  return a;
101 }
102 
104 {
105  G4cout<<"/********************************************/ " << G4endl;
106  G4cout<<" - GVFlashShowerParameterisation::Material - " << G4endl;
107  G4cout<<" Material : " << mat->GetName() << G4endl;
108  G4cout<<" Z = " << Z << G4endl;
109  G4cout<<" A = " << A << G4endl;
110  G4cout<<" X0 = " << X0/cm << " cm" << G4endl;
111  G4cout<<" Rm = " << Rm/cm << " cm" << G4endl;
112  G4cout<<" Ec = " << Ec/MeV << " MeV"<< G4endl;
113  G4cout<<"/********************************************/ " << G4endl;
114 }
115 
117 {
118  G4double Phi = twopi*G4UniformRand() ;
119  return Phi;
120 }
121 
123 {
124  return fGamma->Gamma(a, x);
125 }
G4double GetEffZ(const G4Material *material)
G4double GetZ() const
Definition: G4Material.cc:623
std::vector< ExP01TrackerHit * > a
Definition: ExP01Classes.hh:33
const G4String & GetName() const
Definition: G4Material.hh:178
G4double GetZ() const
Definition: G4Element.hh:131
static constexpr double g
Definition: G4SIunits.hh:183
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:202
G4double GetA() const
Definition: G4Element.hh:139
int G4int
Definition: G4Types.hh:78
double Gamma(double z)
Definition: Gamma.cc:42
static constexpr double twopi
Definition: G4SIunits.hh:76
Definition: Gamma.hh:46
#define G4UniformRand()
Definition: Randomize.hh:97
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
static constexpr double cm
Definition: G4SIunits.hh:119
G4double GetA() const
Definition: G4Material.cc:636
G4double gam(G4double x, G4double a) const
#define G4endl
Definition: G4ios.hh:61
static constexpr double MeV
Definition: G4SIunits.hh:214
G4double GetEffA(const G4Material *material)
size_t GetNumberOfElements() const
Definition: G4Material.hh:186
double G4double
Definition: G4Types.hh:76
const G4double * GetFractionVector() const
Definition: G4Material.hh:194
static constexpr double mole
Definition: G4SIunits.hh:286