MJGerdaLocalMaterialTable.cc

Go to the documentation of this file.

/*
// 11-02-2004 First Implementation, Luciano
// 11-12-2004 added materials: enriched Ge76,
                 liquid N2, copper, water
*/

#include "materials/MJGerdaLocalMaterialTable.hh"

#include "G4Isotope.hh"
#include "G4Element.hh"
#include "G4Material.hh"


MJGerdaLocalMaterialTable::MJGerdaLocalMaterialTable()
{;}

MJGerdaLocalMaterialTable::~MJGerdaLocalMaterialTable()
{;}

void MJGerdaLocalMaterialTable::CreateMaterials()
{
  //here one should define the material he/she wants to use 
  //in the geometry, if not read from the DB
  G4double a;  // atomic mass
  G4double z;  // atomic number
  G4double density,fractionmass,abundance;
  G4int iz,in;
  G4String name,symbol;
  G4int ncomponents;
  G4int natoms;
  G4double temperature, pressure;
  G4State  state;

 
  G4Element *N = new G4Element (name="Nitrogen",
                                symbol="N",z=7.,a=14.01*g/mole);
  G4Element* O  = new G4Element
    (name="Oxygen"  ,symbol="O" , z= 8., a=16.00*g/mole);
 
  G4Element* Si = new G4Element
    (name="Silicon",symbol="Si" , z= 14., a=28.09*g/mole);
  G4Element* Al = new G4Element
    (name="Aluminium"  ,symbol="Al" , z= 13., a=26.98*g/mole);
  G4Element* Fe = new G4Element
    (name="Iron"  ,symbol="Fe" , z= 26., a=55.85*g/mole);
  G4Element* H = new G4Element
    (name="Hydrogen",symbol="H" , z= 1., a=1.00794*g/mole);
  G4Element* Ca = new G4Element
    (name="Calcium",symbol="Ca" , z= 20., a=40.078*g/mole);
  G4Element* Cl = new G4Element
    (name="Chlorine",symbol="Cl", z=17., a=35.45*g/mole);
  G4Element* C  = new G4Element( "Carbon", "C",   6. , 12.011*g/mole);
  G4Element* Co = new G4Element( "Cobalt", "Co", 27. , 58.9332*g/mole);
 
  G4Isotope* Ga69 = new G4Isotope
    (name="Gallium69", iz= 31, in=69, a=69*g/mole);
  G4Isotope* Ga71 = new G4Isotope
    (name="Gallium71",iz=31,in=71,a=71*g/mole);
  G4Element* Ga = new G4Element
        (name="Gallium",symbol="Ga", ncomponents=2);
  Ga->AddIsotope(Ga69, abundance=0.61);
  Ga->AddIsotope(Ga71, abundance=0.39);

  G4Element* Cu = new G4Element
    (name="Copper"  ,symbol="Cu" , z= 29., a=63.55*g/mole);

    G4Isotope* Ge73 = new G4Isotope(name="Ge73",  32, 73, 73.0*g/mole);
    G4Isotope* Ge74 = new G4Isotope(name="Ge74",  32, 74, 74.0*g/mole);
    G4Isotope* Ge76 = new G4Isotope(name="Ge76",  32, 76, 76.0*g/mole);

    a = 75.71*g/mole;
    G4int nIsotopes=3;
    G4Element* Germanium=new G4Element(name="Germanium",symbol="Ge",nIsotopes);
    Germanium->AddIsotope(Ge73,0.01);
    Germanium->AddIsotope(Ge74,0.13);
    Germanium->AddIsotope(Ge76,0.86);

   G4Element* Pb = new G4Element
    (name="Lead",symbol="Pb" , z= 82., a=207.2*g/mole);
 
  //-- then materials
 
  density=1.290*mg/cm3;
  G4Material *air = new G4Material(name="Air",density,ncomponents=2,
                                    kStateGas,300.00*kelvin,1.0*atmosphere);
  air->AddElement(N,fractionmass=0.8);
  air->AddElement(O,fractionmass=0.2);

  //concrete
  G4Material* concrete = new G4Material
    (name="Concrete", density=2.3*g/cm3, ncomponents=6);
  concrete->AddElement(Si, 0.227915);
  concrete->AddElement(O, 0.60541);
  concrete->AddElement(H, 0.09972);
  concrete->AddElement(Ca, 0.04986);
  concrete->AddElement(Al, 0.014245);
  concrete->AddElement(Fe, 0.00285);
 
  //stainless steel
  G4Material* ssteel = new G4Material
    (name="Steel", density=7.7*g/cm3, ncomponents=3);
  ssteel->AddElement(C, 0.04);
  ssteel->AddElement(Fe, 0.88);
  ssteel->AddElement(Co, 0.08);

  // vacuum
    a= 14.000674*g/mole;
    density = 1.E-10*g/cm3;
    G4Material* vacuum = new G4Material(name="Vacuum", z=7., a, density);

  // liquid N2
    density=0.808*g/cm3;
    ncomponents=1;
    natoms=1;
    state=kStateLiquid;
    temperature=80.0*kelvin;
    pressure=1.0*bar;
    name="LiquidNitrogen";
    G4Material* liquidnitrogen = new G4Material(name,density,ncomponents,
                                     state,temperature,pressure);
    liquidnitrogen->AddElement(N,natoms);

  // enriched germanium
    density = 5.32*g/cm3;
    G4Material* ge = new G4Material(name="Ge", density, 1);
    G4double fraction=1.;
    ge->AddElement(Germanium, fraction);

  // copper
  G4Material* metalCu = new G4Material
    (name="MetalCopper", density=8.960*g/cm3, ncomponents=1);
  metalCu->AddElement(Cu, 1);

 //water
  G4Material* water = new G4Material
    (name="water", density=1.00*g/cm3, ncomponents=2);
  water->AddElement(H , 2);
  water->AddElement(O , 1);

  // lead
  G4Material* metalPb = new G4Material
    (name="MetalLead", density=11.340*g/cm3, ncomponents=1);
  metalPb->AddElement(Pb, 1);
 
}

---