Project

General

Profile

Feature #1551 » LiH.patch

Rogers, Chris, 18 September 2014 11:53

View differences:

src/legacy/Simulation/FillMaterials.cc 2014-09-17 15:51:25 +0000
84 84
  aerogel0_2->AddMaterial(water, 0.03);
85 85
  materials_list->addMaterial( aerogel0_2, name );
86 86

  
87
  // MICE Lithium Hydride
88
  // density of 0.69 g/cm3 is a myth/rumour; NO citation
89
  // Alan Bross is contact #1551
90
  density = 0.69*g/cm3;
91
  name = "MICE_LITHIUM_HYDRIDE";
92
  G4Material* miceLiH = new G4Material(name, density, nComp=2);
93
  G4Element* elMiceLi = new G4Element("MICE_LITHIUM", "Li", 2);
94
  G4Isotope* Li6 = new G4Isotope("Li6", 3, 6);
95
  G4Isotope* Li7 = new G4Isotope("Li7", 3, 7);
96
  elMiceLi->AddIsotope(Li6, 0.98);
97
  elMiceLi->AddIsotope(Li7, 0.02);
98
  miceLiH->AddElement(elMiceLi, 0.126797);
99
  G4Element* elH = man->FindOrBuildElement("H");
100
  miceLiH->AddElement(elH, 0.873203);
101
  materials_list->addMaterial( miceLiH, name );
102

  
87 103
  // BC600, here represented only as bisphenol (rho_bisphenol=1.20gcm-3)
88 104
  density = 1.18*g/cm3; // refractive index = 1.56 
89 105
  name = "BC600";
90 106
  G4Material* bc600 = new G4Material(name, density, nComp=3);
91 107
  G4Element* elC = man->FindOrBuildElement("C");
92 108
  bc600->AddElement(elC, 15);
93
  G4Element* elH = man->FindOrBuildElement("H");
94 109
  bc600->AddElement(elH, 16);
95 110
  G4Element* elO = man->FindOrBuildElement("O");
96 111
  bc600->AddElement(elO,  2);		//ME - Tamas thinks this should be 2, not 12!!!
(2-2/7)