LArSoft  v07_13_02
Liquid Argon Software toolkit - http://larsoft.org/
Public Member Functions | List of all members
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 ( )

Definition at line 58 of file MyMaterials.cc.

References Z.

Referenced by BOOST_PYTHON_MODULE().

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 }
Z
Float_t Z
Definition: plot.C:39
The documentation for this class was generated from the following files: