Geant4  10.01.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 69000 2013-04-15 09:24:35Z 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 }
53 
55 {
56 }
57 
59 {
60  // Returns Z or effective Z=sum(pi*Zi) (if compound/mixture)
61  // of given material
62  //
63  G4double z = 0.;
64  G4int nofElements = mat->GetNumberOfElements();
65  if (nofElements > 1)
66  {
67  for (G4int i=0; i<nofElements; i++) {
68  G4double zOfElement = mat->GetElement(i)->GetZ();
69  G4double massFraction = mat->GetFractionVector()[i];
70  // cout << mat->GetElement(i)->GetName()
71  // <<" Z= "<<zOfElement << " , Fraction= "<<massFraction <<endl;
72  z += zOfElement*massFraction;
73  }
74  }
75  else {
76  z = mat->GetZ();
77  }
78  return z;
79 }
80 
82 {
83  // Returns A or effective A=sum(pi*Ai) (if compound/mixture)
84  // of given material
85  //
86  G4double a = 0.;
87  G4int nofElements = mat->GetNumberOfElements();
88  if (nofElements > 1) {
89  for (G4int i=0; i<nofElements; i++) {
90  G4double aOfElement = mat->GetElement(i)->GetA()/(g/mole);
91  G4double massFraction = mat->GetFractionVector()[i];
92  a += aOfElement*massFraction;
93  }
94  }
95  else {
96  a = mat->GetA()/(g/mole);
97  }
98  return a;
99 }
100 
102 {
103  G4cout<<"/********************************************/ " << G4endl;
104  G4cout<<" - GVFlashShowerParameterisation::Material - " << G4endl;
105  G4cout<<" Material : " << mat->GetName() << G4endl;
106  G4cout<<" Z = "<< Z << G4endl;
107  G4cout<<" A = "<< A << G4endl;
108  G4cout<<" X0 = "<<X0/cm <<" cm" << G4endl;
109  G4cout<<" Rm= "<<Rm/cm <<" cm" << G4endl;
110  G4cout<<" Ec = "<<Ec/MeV << " MeV"<< G4endl;
111  G4cout<<"/********************************************/ " << G4endl;
112 }
113 
115 {
116  G4double Phi = twopi*G4UniformRand() ;
117  return Phi;
118 }
119 
121 {
122  static MyGamma theG;
123  return theG.Gamma(a, x);
124 }
static const double cm
Definition: G4SIunits.hh:106
G4double GetEffZ(const G4Material *material)
static const double MeV
Definition: G4SIunits.hh:193
G4double GetZ() const
Definition: G4Material.cc:632
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:95
G4GLOB_DLL std::ostream G4cout
static const G4double A[nN]
G4double GetA() const
Definition: G4Material.cc:645
static const double g
Definition: G4SIunits.hh:162
G4double gam(G4double x, G4double a) const
static const double mole
Definition: G4SIunits.hh:265
#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