Geant4  10.03.p03
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
MyMaterials Class Reference

#include <MyMaterials.hh>

Public Member Functions

 MyMaterials ()
 
 ~MyMaterials ()
 
void Construct ()
 

Detailed Description

Definition at line 42 of file MyMaterials.hh.

Constructor & Destructor Documentation

MyMaterials::MyMaterials ( )

Definition at line 44 of file MyMaterials.cc.

46 {
47 }
MyMaterials::~MyMaterials ( )

Definition at line 51 of file MyMaterials.cc.

53 {
54 }

Member Function Documentation

void MyMaterials::Construct ( void  )

Definition at line 58 of file MyMaterials.cc.

60 {
61  G4double A, Z;
62 
63  // ------------------------------------------------------------------------
64  // Elements
65  // ------------------------------------------------------------------------
66  G4Element* elH = new G4Element("Hydrogen","H", Z=1., A=1.00794*g/mole);
67  G4Element* elC = new G4Element("Carbon", "C", Z=6., A= 12.011 *g/mole);
68  G4Element* elN = new G4Element("Nitrogen","N", Z=7., A= 14.00674*g/mole);
69  G4Element* elO = new G4Element("Oxygen", "O", Z=8., A= 15.9994*g/mole);
70 
71  // ------------------------------------------------------------------------
72  // Materials
73  // ------------------------------------------------------------------------
74  G4double density, massfraction;
75  G4int natoms, nel;
76 
77  // temperature of experimental hall is controlled at 20 degree.
78  const G4double expTemp= STP_Temperature+20.*kelvin;
79 
80  // vacuum
81  density= universe_mean_density;
82  G4Material* Vacuum= new G4Material("Vacuum", density, nel=2);
83  Vacuum-> AddElement(elN, .7);
84  Vacuum-> AddElement(elO, .3);
85 
86  // air
87  density= 1.2929e-03 *g/cm3; // at 20 degree
88  G4Material* Air= new G4Material("Air", density, nel=2,
89  kStateGas, expTemp);
90  G4double ttt= 75.47+23.20;
91  Air-> AddElement(elN, massfraction= 75.47/ttt);
92  Air-> AddElement(elO, massfraction= 23.20/ttt);
93 
94  // water
95  density= 1.000*g/cm3;
96  G4Material* H2O= new G4Material("Water", density, nel=2);
97  H2O-> AddElement(elH, natoms=2);
98  H2O-> AddElement(elO, natoms=1);
99 
100  // alminium
101  A= 26.98 *g/mole;
102  density= 2.70 *g/cm3;
103  G4Material* Al= new G4Material("Al", Z=13., A, density);
104 
105  // iron
106  A= 55.847 *g/mole;
107  density= 7.87 *g/cm3;
108  G4Material* Fe= new G4Material("Iron", Z=26., A, density);
109 
110  // lead
111  A= 207.2 *g/mole;
112  density= 11.35 *g/cm3;
113  G4Material* Pb= new G4Material("Lead", Z=82., A, density);
114 
115  // scintillator (Polystyene(C6H5CH=CH2))
116  density= 1.032 *g/cm3;
117  G4Material* Scinti= new G4Material("Scinti", density, nel=2);
118  Scinti-> AddElement(elC, natoms=8);
119  Scinti-> AddElement(elH, natoms=8);
120 
121 }
int universe_mean_density
Definition: hepunit.py:307
float STP_Temperature
Definition: hepunit.py:302
int G4int
Definition: G4Types.hh:78
function g(Y1, Y2, PT2)
Definition: hijing1.383.f:5205
double A(double temperature)
static constexpr double kelvin
Definition: G4SIunits.hh:281
static constexpr double cm3
Definition: G4SIunits.hh:121
double G4double
Definition: G4Types.hh:76
static constexpr double mole
Definition: G4SIunits.hh:286

Here is the call graph for this function:

Here is the caller graph for this function:


The documentation for this class was generated from the following files: