Geant4  10.02.p01
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 93079 2015-10-02 14:43:41Z 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 }
58 
60 {
61  // Returns Z or effective Z=sum(pi*Zi) (if compound/mixture)
62  // of given material
63  //
64  G4double z = 0.;
65  G4int nofElements = mat->GetNumberOfElements();
66  if (nofElements > 1)
67  {
68  for (G4int i=0; i<nofElements; i++) {
69  G4double zOfElement = mat->GetElement(i)->GetZ();
70  G4double massFraction = mat->GetFractionVector()[i];
71  // cout << mat->GetElement(i)->GetName()
72  // <<" Z= "<<zOfElement << " , Fraction= "<<massFraction <<endl;
73  z += zOfElement*massFraction;
74  }
75  }
76  else {
77  z = mat->GetZ();
78  }
79  return z;
80 }
81 
83 {
84  // Returns A or effective A=sum(pi*Ai) (if compound/mixture)
85  // of given material
86  //
87  G4double a = 0.;
88  G4int nofElements = mat->GetNumberOfElements();
89  if (nofElements > 1) {
90  for (G4int i=0; i<nofElements; i++) {
91  G4double aOfElement = mat->GetElement(i)->GetA()/(g/mole);
92  G4double massFraction = mat->GetFractionVector()[i];
93  a += aOfElement*massFraction;
94  }
95  }
96  else {
97  a = mat->GetA()/(g/mole);
98  }
99  return a;
100 }
101 
103 {
104  G4cout<<"/********************************************/ " << G4endl;
105  G4cout<<" - GVFlashShowerParameterisation::Material - " << G4endl;
106  G4cout<<" Material : " << mat->GetName() << G4endl;
107  G4cout<<" Z = " << Z << G4endl;
108  G4cout<<" A = " << A << G4endl;
109  G4cout<<" X0 = " << X0/cm << " cm" << G4endl;
110  G4cout<<" Rm = " << Rm/cm << " cm" << G4endl;
111  G4cout<<" Ec = " << Ec/MeV << " MeV"<< G4endl;
112  G4cout<<"/********************************************/ " << G4endl;
113 }
114 
116 {
117  G4double Phi = twopi*G4UniformRand() ;
118  return Phi;
119 }
120 
122 {
123  return fGamma->Gamma(a, x);
124 }
static const double cm
Definition: G4SIunits.hh:118
G4double GetEffZ(const G4Material *material)
static const double MeV
Definition: G4SIunits.hh:211
G4double GetZ() const
Definition: G4Material.cc:625
G4double z
Definition: TRTMaterials.hh:39
const G4String & GetName() const
Definition: G4Material.hh:178
G4double GetZ() const
Definition: G4Element.hh:131
G4double a
Definition: TRTMaterials.hh:39
const G4Element * GetElement(G4int iel) const
Definition: G4Material.hh:202
G4double GetA() const
Definition: G4Element.hh:138
int G4int
Definition: G4Types.hh:78
double Gamma(double z)
Definition: Gamma.cc:42
Definition: Gamma.hh:46
G4double density
Definition: TRTMaterials.hh:39
#define G4UniformRand()
Definition: Randomize.hh:97
G4GLOB_DLL std::ostream G4cout
double A(double temperature)
static const double twopi
Definition: G4SIunits.hh:75
G4double GetA() const
Definition: G4Material.cc:638
static const double g
Definition: G4SIunits.hh:180
const G4double x[NPOINTSGL]
G4double gam(G4double x, G4double a) const
static const double mole
Definition: G4SIunits.hh:283
#define G4endl
Definition: G4ios.hh:61
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