//---------------------------------------------------------------------------// //bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu// // // // // // MaGe Simulation // // // // This code implementation is the intellectual property of the // // MAJORANA and Gerda Collaborations. It is based on Geant4, an // // intellectual property of the RD44 GEANT4 collaboration. // // // // ********************* // // // // Neither the authors of this software system, nor their employing // // institutes, nor the agencies providing financial support for this // // work make any representation or warranty, express or implied, // // regarding this software system or assume any liability for its use. // // By copying, distributing or modifying the Program (or any work based // // on on the Program) you indicate your acceptance of this statement, // // and all its terms. // // // //bb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nubb0nu// //---------------------------------------------------------------------------// /* 11-02-2004 First Implementation, Luciano 11-12-2004 added materials: enriched Ge76, liquid N2, copper, water 12-3-2004 new materials and naming conventions, C.Tomei 12-16-2004 added Beryllium and Lithium, Luciano (used to by-pass a stupid tracking in bug) 02-15-2005 added Optical Properties of water, Luciano 03-16-2005 changed ssteel material composition and density, added Cr, Mn and Ni elements, dkollar (material name should be changed to StainlessSteel) 04-07-2005 added Quartz, Luciano 07-11-2005 added nickel metal and gold metal, Kevin 09-14-2005 added polyethylene, Kevin 06-09-2006 fixed typo with the density of gaseous nitrogen, Luciano 10-12-2006 added method to define materials on-the-fly from external files, Luciano 12-07-2006 added tungsten, Manuela 02-22-2007 added materials: Brass, ProportionalGas, N2Gas, GeLi,Acrylic, PVC, HD100, Monel,PE, and element Zinc, Dusan 03-15-2007 added silver, Manuela 04-03-2007 added materials: Bronze, GeB, and element Tin, Hardy 06-27-2007 added hydrogen material for testing purposes, Luciano 07-02-2007 added gaseous argon, Luciano 10-16-2007 added PET, Daniel 10-18-2007 added Polyimid aka. Kapton, Daniel 11-28-2007 added Polycarbonate(makrolon), Jens 12-07-2007 added Torlon, Jens 04-29-2008 added LNGS Rock, Luciano 07-02-2008 modified density of enriched Germanium, according to data from the existing detectors at LNGS (density about 5.5 g/cm3, while it is 5.32 g/cm3 for natural metallic Ge). 07-10-2008 modified again density of enriched Germanium, according to calculations performed by Jens (density = 1.04236*density_nat = 5.54 g/cm3) 17-02-2009 added Tantalum, Francis 14-05-2009 added Erbium and Krypton elements, Marc 22-07-2009 added borated polyethylene of 5% and 10% Boron content, Alan Poon 21-01-2010 added lithiated polyethylene (7.5% by weight), Alan Poon 08-03-2010 added PPO, DMP, generic epoxyresin, CU5, polypropylene, Matthias Laubenstein */ #include "MGGerdaLocalMaterialTable.hh" #include "G4Isotope.hh" #include "G4Element.hh" #include "G4Material.hh" #include "G4MaterialPropertiesTable.hh" #include "G4NistManager.hh" #include "G4SystemOfUnits.hh" MGGerdaLocalMaterialTable::MGGerdaLocalMaterialTable() {;} MGGerdaLocalMaterialTable::~MGGerdaLocalMaterialTable() {;} void MGGerdaLocalMaterialTable::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 abundance; G4double density; G4double fractionmass; G4int in; G4int iz; G4String name; G4int natoms; G4int ncomponents; G4int nIsotopes; G4int nIsotopes_Se; G4double pressure; G4String symbol; G4double temperature; G4double z; // atomic number // ================ // define elements // ================ G4Element* elH = new G4Element (name="Hydrogen",symbol="H",z=1.,a=1.00794*g/mole); G4Element* elLi = new G4Element (name="Lithium",symbol="Li",z=3.,a=6.941*g/mole); G4Element* elBe = new G4Element (name="Beryllium",symbol="Be",z=4.,a=9.012*g/mole); G4Element* elB = new G4Element (name="Boron",symbol="B",z=5.,a=10.811*g/mole); G4Element* elC = new G4Element (name="Carbon",symbol="C",z=6.,a=12.011*g/mole); G4Element* elN = new G4Element (name="Nitrogen",symbol="N",z=7.,a=14.01*g/mole); G4Element* elO = new G4Element (name="Oxygen",symbol="O",z=8.,a=16.00*g/mole); G4Element* elF = new G4Element (name="Fluorine",symbol="F",z=9.,a=19.00*g/mole); G4Element* elNa = new G4Element (name="Natrium",symbol="Na",z=11.,a=22.990*g/mole); G4Element* elMg = new G4Element (name="Magnesium",symbol="Mg",z=12.,a=24.305*g/mole); G4Element* elAl = new G4Element (name="Aluminum" ,symbol="Al",z=13.,a=26.98*g/mole); G4Element* elSi = new G4Element (name="Silicon",symbol="Si",z=14.,a=28.09*g/mole); G4Element* elP = new G4Element (name="Phosphorous",symbol="P",z=15.,a=30.974*g/mole); G4Element* elS = new G4Element (name="Sulfur",symbol="S",z=16.,a=32.066*g/mole); G4Element* elCl = new G4Element (name="Chlorine",symbol="Cl",z=17.,a=35.45*g/mole); G4Element* elAr = new G4Element (name="Argon",symbol="Ar",z=18.,a=39.95*g/mole); G4Element* elK = new G4Element (name="Potassium",symbol="K",z=19.,a=39.098*g/mole); G4Element* elCa = new G4Element (name="Calcium",symbol="Ca",z=20.,a=40.078*g/mole); G4Element* elTi = new G4Element (name="Titanium",symbol="Ti",z=22.,a=47.867*g/mole); G4Element* elV = new G4Element (name="Vanadium",symbol="V",z=23.,a=50.942*g/mole); G4Element* elCr = new G4Element (name="Chromium",symbol="Cr",z=24.,a=51.9961*g/mole); G4Element* elMn = new G4Element (name="Manganese",symbol="Mn",z=25.,a=54.93805*g/mole); G4Element* elFe = new G4Element (name="Iron",symbol="Fe",z=26.,a=55.845*g/mole); G4Element* elCo = new G4Element (name="Cobalt",symbol="Co",z=27.,a=58.9332*g/mole); G4Element* elNi = new G4Element (name="Nickel",symbol="Ni",z=28.,a=58.6934*g/mole); G4Element* elCu = new G4Element (name="Copper",symbol="Cu",z=29.,a=63.55*g/mole); G4Element* elZn = new G4Element (name="Zinc", symbol="Zn", z=30.,a=65.38*g/mole); G4Element* elSe = new G4Element (name="Selenium",symbol="Se",z=34.,a=78.96*g/mole); G4Element* elBr = new G4Element (name="Bromine",symbol="Br",z=35.,a=79.904*g/mole); G4Element* elKr = new G4Element (name="Krypton",symbol="Kr",z=36.,a=83.798*g/mole); G4Element* elSr = new G4Element (name="Strontium",symbol="Sr",z=38.,a=87.62*g/mole); G4Element* elY = new G4Element (name="Yttrium",symbol="Y",z=39.,a=88.90584*g/mole); G4Element* elZr = new G4Element (name="Zirconium",symbol="Zr",z=40.,a=91.224*g/mole); G4Element* elNb = new G4Element (name="Niobium",symbol="Nb",z=41.,a=92.91*g/mole); G4Element* elMo = new G4Element (name="Molybdenum",symbol="Mo",z=42.,a=95.94*g/mole); G4Element* elPd = new G4Element (name="Palladium", symbol="Pd", z=46, a=106.42*g/mole); G4Element* elAg = new G4Element (name="Silver", symbol="Ag", z=47, a=107.87*g/mole); G4Element* elCd = new G4Element (name="Cadmium",symbol="Cd",z=48.,a=112.41*g/mole); G4Element* elIn = new G4Element (name="Indium",symbol="In",z=49.,a=114.82*g/mole); G4Element* elSn = new G4Element (name="Tin", symbol="Sn", z=50.,a=118.71*g/mole); G4Element* elSb = new G4Element (name="Antimony",symbol="Sb",z=51.,a=121.76*g/mole); G4Element* elTe = new G4Element (name="Tellurium",symbol="Te",z=52.,a=127.60*g/mole); G4Element* elI = new G4Element (name="Iodine",symbol="I",z=53.,a=126.90*g/mole); G4Element* elCs = new G4Element (name="Cesium",symbol="Cs",z=55.,a=132.91*g/mole); G4Element* elBa = new G4Element (name="Barium",symbol="Ba",z=56.,a=137.33*g/mole); G4Element* elLa = new G4Element (name="Lanthanum",symbol="La",z=57.,a=138.90547*g/mole); G4Element* elCe = new G4Element (name="Cerium",symbol="Ce",z=58.,a=140.12*g/mole); G4Element* elPr = new G4Element (name="Paseodyimium",symbol="Pr",z=59.,a=140.90766*g/mole); G4Element* elNd = new G4Element (name="Neodymium",symbol="Nd",z=60.,a=144.24*g/mole); G4Element* elSm = new G4Element (name="Samarium",symbol="Sm",z=62.,a=150.36*g/mole); G4Element* elEu = new G4Element (name="Europium",symbol="Eu",z=63.,a=151.96*g/mole); G4Element* elGd = new G4Element (name="Gadolinium",symbol="Gd",z=64.,a=157.25*g/mole); G4Element* elTb = new G4Element (name="Terbium",symbol="Tb",z=65.,a=158.92535*g/mole); G4Element* elEr = new G4Element (name="Erbium",symbol="Er",z=68.,a=167.26*g/mole); G4Element* elTm = new G4Element (name="Thulium",symbol="Tm",z=69.,a=168.93422*g/mole); G4Element* elYb = new G4Element (name="Ytterbium",symbol="Yb",z=70.,a=173.04*g/mole); G4Element* elHf = new G4Element (name="Hafnium", symbol="Hf", z=72., a=178.49*g/mole); G4Element* elTa = new G4Element (name="Tantalum", symbol="Ta", z=73., a=180.94*g/mole); G4Element* elW = new G4Element (name="Tungsten", symbol="W", z=74, a=182.30*g/mole); G4Element* elOs = new G4Element (name="Osmium", symbol="Os", z=76, a=190.23*g/mole); G4Element* elPt = new G4Element (name="Platinum", symbol="Pt", z=78, a=195.08*g/mole); G4Element* elAu = new G4Element (name="Gold", symbol="Au", z=79, a=196.97*g/mole); G4Element* elPb = new G4Element (name="Lead",symbol="Pb",z=82.,a=207.2*g/mole); G4Element* elBi = new G4Element (name="Bismuth",symbol="Bi",z=83.,a=208.98040*g/mole); G4Element* elTh = new G4Element (name="Thorium",symbol="Th",z=90.,a=232.04*g/mole); // ================ // define isotopes // ================ 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* elGa = new G4Element (name="Gallium",symbol="Ga",ncomponents=2); elGa->AddIsotope(Ga69, abundance=0.61); elGa->AddIsotope(Ga71, abundance=0.39); G4Isotope* Ge70 = new G4Isotope(name="Ge70", 32, 70, 69.92*g/mole); G4Isotope* Ge72 = new G4Isotope(name="Ge72", 32, 72, 71.92*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); G4Element* elGeNat=new G4Element(name="naturalGermanium", symbol="GeNat",nIsotopes=5); elGeNat->AddIsotope(Ge70,abundance= 20.9*perCent); elGeNat->AddIsotope(Ge72,abundance= 27.5*perCent); elGeNat->AddIsotope(Ge73,abundance= 7.7*perCent); elGeNat->AddIsotope(Ge74,abundance= 36.3*perCent); elGeNat->AddIsotope(Ge76,abundance= 7.6*perCent); G4Element* elGeEnr = new G4Element(name="enrichedGermanium", symbol="GeEnr",nIsotopes=5); elGeEnr->AddIsotope(Ge70,abundance= 0.0*perCent); elGeEnr->AddIsotope(Ge72,abundance= 0.1*perCent); elGeEnr->AddIsotope(Ge73,abundance= 0.2*perCent); elGeEnr->AddIsotope(Ge74,abundance= 13.1*perCent); elGeEnr->AddIsotope(Ge76,abundance= 86.6*perCent); G4Isotope* Se76 = new G4Isotope(name="Se76", 34, 76, 75.919*g/mole); G4Isotope* Se77 = new G4Isotope(name="Se77", 34, 77, 76.920*g/mole); G4Isotope* Se78 = new G4Isotope(name="Se78", 34, 78, 77.917*g/mole); G4Isotope* Se80 = new G4Isotope(name="Se80", 34, 80, 79.917*g/mole); G4Isotope* Se82 = new G4Isotope(name="Se82", 34, 82, 81.917*g/mole); G4Element* elSeEnrHigh = new G4Element(name="highenrichedSelenium", symbol="SeEnrHigh",nIsotopes_Se=5); elSeEnrHigh->AddIsotope(Se76,abundance= 0.1*perCent); elSeEnrHigh->AddIsotope(Se77,abundance= 0.1*perCent); elSeEnrHigh->AddIsotope(Se78,abundance= 0.3*perCent); elSeEnrHigh->AddIsotope(Se80,abundance= 4.2*perCent); elSeEnrHigh->AddIsotope(Se82,abundance= 95.3*perCent); G4Element* elSeEnrLow = new G4Element(name="lowenrichedSelenium", symbol="SeEnrLow",nIsotopes_Se=5); elSeEnrLow->AddIsotope(Se76,abundance= 1.3*perCent); elSeEnrLow->AddIsotope(Se77,abundance= 1.0*perCent); elSeEnrLow->AddIsotope(Se78,abundance= 3.4*perCent); elSeEnrLow->AddIsotope(Se80,abundance= 10.2*perCent); elSeEnrLow->AddIsotope(Se82,abundance= 84.1*perCent); // ================ // define materials // ================ // Acrylic G4Material* Acrylic = new G4Material (name="Acrylic", density= 1.15*g/cm3, ncomponents=3); Acrylic->AddElement(elC,natoms=5); Acrylic->AddElement(elH,natoms=8); Acrylic->AddElement(elO,natoms=2); //austenitic steel for clamps (CF811/316) G4Material* auststeel = new G4Material (name="auststeel",density=8.15727*g/cm3,ncomponents=5); auststeel->AddElement(elMo,fractionmass= 0.0200); auststeel->AddElement(elCr,fractionmass= 0.1800); auststeel->AddElement(elNi,fractionmass= 0.0900); auststeel->AddElement(elFe,fractionmass= 0.7092); auststeel->AddElement(elC, fractionmass= 0.0008); //BaF2 G4Material* BaF2 = new G4Material (name="BaF2",density=4.89*g/cm3,ncomponents=2); BaF2->AddElement(elBa,natoms= 1); BaF2->AddElement(elF,natoms= 2); //Cu-Be alloy (8.129 g cm-3) according to (ASTM B194, B196, B197), there can be Co instead of Ni G4Material* becu = new G4Material (name="becu",density=8.129*g/cm3,ncomponents=3); becu->AddElement(elCu,fractionmass= 0.9775); becu->AddElement(elBe,fractionmass= 0.02); becu->AddElement(elNi,fractionmass= 0.0025); //Cu-Be alloy in pins G4Material* becup = new G4Material (name="becup",density=5.7982*g/cm3,ncomponents=3); becup->AddElement(elCu,fractionmass= 0.9775); becup->AddElement(elBe,fractionmass= 0.02); becup->AddElement(elNi,fractionmass= 0.0025); // cable 1 G4Material* cable1 = new G4Material (name="cable1", density=1.1636*g/cm3, ncomponents=3); cable1->AddElement(elF,fractionmass= 34.8*perCent); cable1->AddElement(elC,fractionmass= 11.0*perCent); cable1->AddElement(elCu,fractionmass= 54.2*perCent); // cable 2 G4Material* cable2 = new G4Material (name="cable2", density=1.4986*g/cm3, ncomponents=3); cable2->AddElement(elF,fractionmass= 34.8*perCent); cable2->AddElement(elC,fractionmass= 11.0*perCent); cable2->AddElement(elCu,fractionmass= 54.2*perCent); // cable 3 G4Material* cable3 = new G4Material (name="cable3", density=1.758946*g/cm3, ncomponents=3); cable3->AddElement(elF,fractionmass= 64.8*perCent); cable3->AddElement(elC,fractionmass= 20.5*perCent); cable3->AddElement(elCu,fractionmass= 14.7*perCent); // cable 4 G4Material* cable4 = new G4Material (name="cable4", density=3.67989*g/cm3, ncomponents=3); cable4->AddElement(elF,fractionmass= 64.8*perCent); cable4->AddElement(elC,fractionmass= 20.5*perCent); cable4->AddElement(elCu,fractionmass= 14.7*perCent); // cable 5 G4Material* cable5 = new G4Material (name="cable5", density=1.853278*g/cm3, ncomponents=3); cable5->AddElement(elF,fractionmass= 64.8*perCent); cable5->AddElement(elC,fractionmass= 20.5*perCent); cable5->AddElement(elCu,fractionmass= 14.7*perCent); // cable 6 G4Material* cable6 = new G4Material (name="cable6", density=2.11092*g/cm3, ncomponents=3); cable6->AddElement(elF,fractionmass= 64.8*perCent); cable6->AddElement(elC,fractionmass= 20.5*perCent); cable6->AddElement(elCu,fractionmass= 14.7*perCent); // cable 7 G4Material* cable7 = new G4Material (name="cable7", density=2.08218*g/cm3, ncomponents=3); cable7->AddElement(elF,fractionmass= 64.8*perCent); cable7->AddElement(elC,fractionmass= 20.5*perCent); cable7->AddElement(elCu,fractionmass= 14.7*perCent); // cable 8 G4Material* cable8 = new G4Material (name="cable8", density=1.0149*g/cm3, ncomponents=3); cable8->AddElement(elF,fractionmass= 64.8*perCent); cable8->AddElement(elC,fractionmass= 20.5*perCent); cable8->AddElement(elCu,fractionmass= 14.7*perCent); // cable 9 G4Material* cable9 = new G4Material (name="cable9", density=1.38307*g/cm3, ncomponents=4); cable9->AddElement(elF,fractionmass= 39.35*perCent); cable9->AddElement(elC,fractionmass= 12.42*perCent); cable9->AddElement(elCu,fractionmass= 5.84*perCent); cable9->AddElement(elAg,fractionmass= 42.39*perCent); // cable 12, thin AXON cable LEGEND, 3.75 mg/cm, TUM G4Material* cable12 = new G4Material (name="cable12", density=1.464498*g/cm3, ncomponents=3); cable12->AddElement(elF,fractionmass= 34.8*perCent); cable12->AddElement(elC,fractionmass= 11.0*perCent); cable12->AddElement(elCu,fractionmass= 54.2*perCent); // cable 15a G4Material* cable15a = new G4Material (name="cable15a", density=1.287760*g/cm3, ncomponents=3); cable15a->AddElement(elF,fractionmass= 34.8*perCent); cable15a->AddElement(elC,fractionmass= 11.0*perCent); cable15a->AddElement(elCu,fractionmass= 54.2*perCent); // cable 15b G4Material* cable15b = new G4Material (name="cable15b", density=1.301907*g/cm3, ncomponents=3); cable15b->AddElement(elF,fractionmass= 34.8*perCent); cable15b->AddElement(elC,fractionmass= 11.0*perCent); cable15b->AddElement(elCu,fractionmass= 54.2*perCent); //CaF2 G4Material* CaF2 = new G4Material (name="CaF2",density=3.18*g/cm3,ncomponents=2); CaF2->AddElement(elCa,natoms= 1); CaF2->AddElement(elF,natoms= 2); //CaF2 grains, DAMA, 105.2 g G4Material* CaF2a = new G4Material (name="CaF2a",density=1.9234*g/cm3,ncomponents=2); CaF2a->AddElement(elCa,natoms= 1); CaF2a->AddElement(elF,natoms= 2); //CaF2 grains, DAMA, 102.2 g & 101.6 g G4Material* CaF2b = new G4Material (name="CaF2b",density=1.8633*g/cm3,ncomponents=2); CaF2b->AddElement(elCa,natoms= 1); CaF2b->AddElement(elF,natoms= 2); //camera orcaflash4.0, CYGNUS, on GeMPI G4Material* camera = new G4Material (name="Camera",density=2.341933*g/cm3,ncomponents=5); camera->AddElement(elSi,fractionmass= 0.01); camera->AddElement(elCr,fractionmass= 0.20); camera->AddElement(elMn,fractionmass= 0.02); camera->AddElement(elFe,fractionmass= 0.67); camera->AddElement(elNi,fractionmass= 0.10); //camera_01 PRIME-BSI EXPRESS, Teledyne, CYGNUS, on GeMPI G4Material* camera_01 = new G4Material (name="Camera_01",density=1.284731*g/cm3,ncomponents=5); camera_01->AddElement(elSi,fractionmass= 0.01); camera_01->AddElement(elCr,fractionmass= 0.20); camera_01->AddElement(elMn,fractionmass= 0.02); camera_01->AddElement(elFe,fractionmass= 0.67); camera_01->AddElement(elNi,fractionmass= 0.10); //meteorite, CM chondrite G4Material* CMch = new G4Material (name="CMch",density=2.20*g/cm3,ncomponents=13); CMch->AddElement(elH,fractionmass= 0.014); CMch->AddElement(elC,fractionmass= 0.022); CMch->AddElement(elN,fractionmass= 0.0015); CMch->AddElement(elO,fractionmass= 0.432); CMch->AddElement(elNa,fractionmass= 0.002); CMch->AddElement(elFe,fractionmass= 0.210); CMch->AddElement(elSi,fractionmass= 0.129); CMch->AddElement(elMg,fractionmass= 0.117); CMch->AddElement(elCr,fractionmass= 0.003); CMch->AddElement(elS,fractionmass= 0.033); CMch->AddElement(elCa,fractionmass= 0.0127); CMch->AddElement(elAl,fractionmass= 0.0118); CMch->AddElement(elNi,fractionmass= 0.012); //meteorite, CM chondrite G4Material* CMchh = new G4Material (name="CMchh",density=2.38*g/cm3,ncomponents=13); CMchh->AddElement(elH,fractionmass= 0.014); CMchh->AddElement(elC,fractionmass= 0.022); CMchh->AddElement(elN,fractionmass= 0.0015); CMchh->AddElement(elO,fractionmass= 0.432); CMchh->AddElement(elNa,fractionmass= 0.002); CMchh->AddElement(elFe,fractionmass= 0.210); CMchh->AddElement(elSi,fractionmass= 0.129); CMchh->AddElement(elMg,fractionmass= 0.117); CMchh->AddElement(elCr,fractionmass= 0.003); CMchh->AddElement(elS,fractionmass= 0.033); CMchh->AddElement(elCa,fractionmass= 0.0127); CMchh->AddElement(elAl,fractionmass= 0.0118); CMchh->AddElement(elNi,fractionmass= 0.012); //coins G4Material* coins = new G4Material (name="coins",density=6.7*g/cm3,ncomponents=3); coins->AddElement(elCr,fractionmass= 0.02); coins->AddElement(elFe,fractionmass= 0.94); coins->AddElement(elNi,fractionmass= 0.04); //concrete G4Material* concrete = new G4Material (name="Concrete",density=2.3*g/cm3,ncomponents=6); concrete->AddElement(elSi,fractionmass= 0.227915); concrete->AddElement(elO,fractionmass= 0.60541); concrete->AddElement(elH,fractionmass= 0.09972); concrete->AddElement(elCa,fractionmass= 0.04986); concrete->AddElement(elAl,fractionmass= 0.014245); concrete->AddElement(elFe,fractionmass= 0.00285); //concrete1 G4Material* concrete1 = new G4Material (name="Concrete1",density=2.2*g/cm3,ncomponents=6); concrete1->AddElement(elSi,fractionmass= 0.227915); concrete1->AddElement(elO,fractionmass= 0.60541); concrete1->AddElement(elH,fractionmass= 0.09972); concrete1->AddElement(elCa,fractionmass= 0.04986); concrete1->AddElement(elAl,fractionmass= 0.014245); concrete1->AddElement(elFe,fractionmass= 0.00285); //concrete2 G4Material* concrete2 = new G4Material (name="Concrete2",density=2.1*g/cm3,ncomponents=6); concrete2->AddElement(elSi,fractionmass= 0.227915); concrete2->AddElement(elO,fractionmass= 0.60541); concrete2->AddElement(elH,fractionmass= 0.09972); concrete2->AddElement(elCa,fractionmass= 0.04986); concrete2->AddElement(elAl,fractionmass= 0.014245); concrete2->AddElement(elFe,fractionmass= 0.00285); //CsI G4Material* CsI = new G4Material (name="CsI",density=4.51*g/cm3,ncomponents=2); CsI->AddElement(elCs,natoms= 1); CsI->AddElement(elI,natoms= 1); //Er2O3, before purification, KINR, on GePV G4Material* Er2O3 = new G4Material (name="Er2O3",density=3.298*g/cm3,ncomponents=2); Er2O3->AddElement(elEr,natoms= 2); Er2O3->AddElement(elO,natoms= 3); //Er2O3, after purification, KINR, on GeCris G4Material* Er2O3after = new G4Material (name="Er2O3after",density=1.037634*g/cm3,ncomponents=2); Er2O3after->AddElement(elEr,natoms= 2); Er2O3after->AddElement(elO,natoms= 3); //meteorite, Eucrite G4Material* Euc = new G4Material (name="Euc",density=2.7*g/cm3,ncomponents=7); Euc->AddElement(elO,fractionmass= 0.438); Euc->AddElement(elFe,fractionmass= 0.158); Euc->AddElement(elSi,fractionmass= 0.223); Euc->AddElement(elMg,fractionmass= 0.038); Euc->AddElement(elCa,fractionmass= 0.074); Euc->AddElement(elAl,fractionmass= 0.064); Euc->AddElement(elTi,fractionmass= 0.005); //meteorite, Eucrite mod (before 2.7*g/cm3, now, (2.87 +- 0.07) g/cm3 G4Material* Eucmod = new G4Material (name="Eucmod",density=2.87*g/cm3,ncomponents=7); Eucmod->AddElement(elO,fractionmass= 0.438); Eucmod->AddElement(elFe,fractionmass= 0.158); Eucmod->AddElement(elSi,fractionmass= 0.223); Eucmod->AddElement(elMg,fractionmass= 0.038); Eucmod->AddElement(elCa,fractionmass= 0.074); Eucmod->AddElement(elAl,fractionmass= 0.064); Eucmod->AddElement(elTi,fractionmass= 0.005); //EuF2 G4Material* EuF2 = new G4Material (name="EuF2",density=6.495*g/cm3,ncomponents=2); EuF2->AddElement(elEu,natoms= 1); EuF2->AddElement(elF,natoms= 2); //stainless steel (flexible tube DS50) G4Material* flexsteel = new G4Material (name="FlexSteel",density=1.04*g/cm3,ncomponents=5); flexsteel->AddElement(elSi,fractionmass= 0.01); flexsteel->AddElement(elCr,fractionmass= 0.20); flexsteel->AddElement(elMn,fractionmass= 0.02); flexsteel->AddElement(elFe,fractionmass= 0.67); flexsteel->AddElement(elNi,fractionmass= 0.10); //stainless steel 1 (various tubing CUORE-US) G4Material* flexsteel1 = new G4Material (name="FlexSteel1",density=1.3050575*g/cm3,ncomponents=5); flexsteel1->AddElement(elSi,fractionmass= 0.01); flexsteel1->AddElement(elCr,fractionmass= 0.20); flexsteel1->AddElement(elMn,fractionmass= 0.02); flexsteel1->AddElement(elFe,fractionmass= 0.67); flexsteel1->AddElement(elNi,fractionmass= 0.10); //stainless steel 2 (various tubing CUORE-US) G4Material* flexsteel2 = new G4Material (name="FlexSteel2",density=0.4872761*g/cm3,ncomponents=5); flexsteel2->AddElement(elSi,fractionmass= 0.01); flexsteel2->AddElement(elCr,fractionmass= 0.20); flexsteel2->AddElement(elMn,fractionmass= 0.02); flexsteel2->AddElement(elFe,fractionmass= 0.67); flexsteel2->AddElement(elNi,fractionmass= 0.10); //stainless steel 3 (steel braided flexible tube with inner PTFE lining, SABRE) G4Material* flexsteel3 = new G4Material (name="FlexSteel3",density=1.789628*g/cm3,ncomponents=5); flexsteel3->AddElement(elSi,fractionmass= 0.01); flexsteel3->AddElement(elCr,fractionmass= 0.20); flexsteel3->AddElement(elMn,fractionmass= 0.02); flexsteel3->AddElement(elFe,fractionmass= 0.67); flexsteel3->AddElement(elNi,fractionmass= 0.10); //stainless forged stainless steel G4Material* fsteel = new G4Material (name="fsteel",density=8.12*g/cm3,ncomponents=5); fsteel->AddElement(elSi,fractionmass= 0.01); fsteel->AddElement(elCr,fractionmass= 0.20); fsteel->AddElement(elMn,fractionmass= 0.02); fsteel->AddElement(elFe,fractionmass= 0.67); fsteel->AddElement(elNi,fractionmass= 0.10); // Germanium/Lithium mixture of 1 per cent assumed (in dead layer) G4Material* GeLi = new G4Material(name="Germanium/Lithium", density = 5.323*g/cm3, ncomponents=2); GeLi -> AddElement(elGeNat, fractionmass=0.99); GeLi -> AddElement(elLi, fractionmass=0.01); //glass (Pyrex) G4Material* glass = new G4Material (name="glass",density=1.48*g/cm3,ncomponents=6); glass->AddElement(elB,fractionmass= 0.040064); glass->AddElement(elO,fractionmass= 0.539562); glass->AddElement(elNa,fractionmass= 0.028191); glass->AddElement(elAl,fractionmass= 0.011644); glass->AddElement(elSi,fractionmass= 0.377220); glass->AddElement(elK,fractionmass= 0.003321); //powder-like soda lime borosilicate glass (50.0 SiO2: 30.0 B2O3: 10.0 CaO: 10.0 Na2O mol%)bubbles, 3M, types k&s k1 G4Material* glassbubbles3M = new G4Material (name="glassbubbles3M",density=0.067892*g/cm3,ncomponents=5); glassbubbles3M->AddElement(elB,fractionmass= 0.103393); glassbubbles3M->AddElement(elO,fractionmass= 0.535566); glassbubbles3M->AddElement(elNa,fractionmass= 0.073290); glassbubbles3M->AddElement(elCa,fractionmass= 0.063882); glassbubbles3M->AddElement(elSi,fractionmass= 0.223870); //glass_fibre, Cryopillar, Cosinus G4Material* glass_fibre = new G4Material (name="glass_fibre",density=1.781591*g/cm3,ncomponents=4); glass_fibre->AddElement(elO,fractionmass= 0.4598); glass_fibre->AddElement(elNa,fractionmass= 0.096442); glass_fibre->AddElement(elSi,fractionmass= 0.336553); glass_fibre->AddElement(elCa,fractionmass= 0.107205); //glass, KM3 sphere G4Material* glassKM3 = new G4Material (name="glassKM3",density=2.70*g/cm3,ncomponents=6); glassKM3->AddElement(elB,fractionmass= 0.040064); glassKM3->AddElement(elO,fractionmass= 0.539562); glassKM3->AddElement(elNa,fractionmass= 0.028191); glassKM3->AddElement(elAl,fractionmass= 0.011644); glassKM3->AddElement(elSi,fractionmass= 0.377220); glassKM3->AddElement(elK,fractionmass= 0.003321); //glass_plate G4Material* glass_plate = new G4Material (name="glass_plate",density=2.4165*g/cm3,ncomponents=4); glass_plate->AddElement(elO,fractionmass= 0.4598); glass_plate->AddElement(elNa,fractionmass= 0.096442); glass_plate->AddElement(elSi,fractionmass= 0.336553); glass_plate->AddElement(elCa,fractionmass= 0.107205); //stainless steel for shell of Glenair Micro-D, DS-20K G4Material* glenair = new G4Material (name="glenair",density=4.129481*g/cm3,ncomponents=5); glenair->AddElement(elSi,fractionmass= 0.018); glenair->AddElement(elCr,fractionmass= 0.004); glenair->AddElement(elMn,fractionmass= 0.010); glenair->AddElement(elFe,fractionmass= 0.962); glenair->AddElement(elC,fractionmass= 0.006); //meteorite, H chondrite G4Material* Hch = new G4Material (name="Hch",density=3.40*g/cm3,ncomponents=8); Hch->AddElement(elO,fractionmass= 0.357); Hch->AddElement(elFe,fractionmass= 0.275); Hch->AddElement(elSi,fractionmass= 0.169); Hch->AddElement(elMg,fractionmass= 0.140); Hch->AddElement(elS,fractionmass= 0.020); Hch->AddElement(elCa,fractionmass= 0.012); Hch->AddElement(elAl,fractionmass= 0.011); Hch->AddElement(elNi,fractionmass= 0.016); //meteorite, H chondrite, high density G4Material* Hchh = new G4Material (name="Hchh",density=3.58*g/cm3,ncomponents=8); Hchh->AddElement(elO,fractionmass= 0.357); Hchh->AddElement(elFe,fractionmass= 0.275); Hchh->AddElement(elSi,fractionmass= 0.169); Hchh->AddElement(elMg,fractionmass= 0.140); Hchh->AddElement(elS,fractionmass= 0.020); Hchh->AddElement(elCa,fractionmass= 0.012); Hchh->AddElement(elAl,fractionmass= 0.011); Hchh->AddElement(elNi,fractionmass= 0.016); //meteorite, H chondrite, low density G4Material* Hchl = new G4Material (name="Hchl",density=3.22*g/cm3,ncomponents=8); Hchl->AddElement(elO,fractionmass= 0.357); Hchl->AddElement(elFe,fractionmass= 0.275); Hchl->AddElement(elSi,fractionmass= 0.169); Hchl->AddElement(elMg,fractionmass= 0.140); Hchl->AddElement(elS,fractionmass= 0.020); Hchl->AddElement(elCa,fractionmass= 0.012); Hchl->AddElement(elAl,fractionmass= 0.011); Hchl->AddElement(elNi,fractionmass= 0.016); //meteorite, Howardite G4Material* How = new G4Material (name="How",density=2.925*g/cm3,ncomponents=9); How->AddElement(elO,fractionmass= 0.435); How->AddElement(elFe,fractionmass= 0.140); How->AddElement(elSi,fractionmass= 0.239); How->AddElement(elMg,fractionmass= 0.088); How->AddElement(elCa,fractionmass= 0.046); How->AddElement(elAl,fractionmass= 0.043); How->AddElement(elTi,fractionmass= 0.002); How->AddElement(elNa,fractionmass= 0.003); How->AddElement(elMn,fractionmass= 0.004); //meteorite, (BGIC0001/2018, from asteroid 2018 LA), Botswana, Howardite G4Material* How_bgic = new G4Material (name="How_bgic",density=2.85*g/cm3,ncomponents=10); How_bgic->AddElement(elO,fractionmass= 0.4429); How_bgic->AddElement(elFe,fractionmass= 0.1395); How_bgic->AddElement(elSi,fractionmass= 0.235); How_bgic->AddElement(elMg,fractionmass= 0.0575); How_bgic->AddElement(elCa,fractionmass= 0.0635); How_bgic->AddElement(elAl,fractionmass= 0.0575); How_bgic->AddElement(elTi,fractionmass= 0.000358); How_bgic->AddElement(elNi,fractionmass= 0.000028); How_bgic->AddElement(elMn,fractionmass= 0.0041); How_bgic->AddElement(elK,fractionmass= 0.00036); //stainless steel INVAR, per Darkside G4Material* invar = new G4Material (name="Invar",density=8.04*g/cm3,ncomponents=3); invar->AddElement(elMn,fractionmass= 0.003); invar->AddElement(elFe,fractionmass= 0.64); invar->AddElement(elNi,fractionmass= 0.357); //meteorite, iron, 1AB - MG G4Material* ironmet = new G4Material (name="ironmet",density=7.0*g/cm3,ncomponents=4); ironmet->AddElement(elFe,fractionmass= 0.924); ironmet->AddElement(elNi,fractionmass= 0.0755); ironmet->AddElement(elGa,fractionmass= 0.0001); ironmet->AddElement(elGeNat,fractionmass= 0.0004); //Kovar (8.36 g cm-3), Nilo K (8.16 g cm-3) G4Material* kovar = new G4Material (name="kovar",density=8.16*g/cm3,ncomponents=4); kovar->AddElement(elSi,fractionmass= 0.01); kovar->AddElement(elCo,fractionmass= 0.17); kovar->AddElement(elFe,fractionmass= 0.53); kovar->AddElement(elNi,fractionmass= 0.29); //meteorite, L chondrite (3.35 g cm-3 average, 3.29 g cm-3 before) G4Material* Lch = new G4Material (name="Lch",density=3.35*g/cm3,ncomponents=8); Lch->AddElement(elO,fractionmass= 0.380); Lch->AddElement(elFe,fractionmass= 0.220); Lch->AddElement(elSi,fractionmass= 0.190); Lch->AddElement(elMg,fractionmass= 0.150); Lch->AddElement(elS,fractionmass= 0.022); Lch->AddElement(elCa,fractionmass= 0.013); Lch->AddElement(elAl,fractionmass= 0.013); Lch->AddElement(elNi,fractionmass= 0.012); //meteorite, L chondrite, high density G4Material* Lchh = new G4Material (name="Lchh",density=3.540668*g/cm3,ncomponents=8); Lchh->AddElement(elO,fractionmass= 0.380); Lchh->AddElement(elFe,fractionmass= 0.220); Lchh->AddElement(elSi,fractionmass= 0.190); Lchh->AddElement(elMg,fractionmass= 0.150); Lchh->AddElement(elS,fractionmass= 0.022); Lchh->AddElement(elCa,fractionmass= 0.013); Lchh->AddElement(elAl,fractionmass= 0.013); Lchh->AddElement(elNi,fractionmass= 0.012); //meteorite, LL chondrite (3.21 g cm-3 average, 3.138 g cm-3 Stubenberg) G4Material* LLch = new G4Material (name="LLch",density=3.21*g/cm3,ncomponents=13); LLch->AddElement(elO,fractionmass= 0.400); LLch->AddElement(elFe,fractionmass= 0.185); LLch->AddElement(elSi,fractionmass= 0.189); LLch->AddElement(elMg,fractionmass= 0.153); LLch->AddElement(elS,fractionmass= 0.023); LLch->AddElement(elCa,fractionmass= 0.013); LLch->AddElement(elAl,fractionmass= 0.012); LLch->AddElement(elNi,fractionmass= 0.010); LLch->AddElement(elNa,fractionmass= 0.007); LLch->AddElement(elCr,fractionmass= 0.004); LLch->AddElement(elMn,fractionmass= 0.002); LLch->AddElement(elP,fractionmass= 0.001); LLch->AddElement(elC,fractionmass= 0.001); //stainless steel AISI 316LN, per CUORE G4Material* lnsteel = new G4Material (name="LNSteel",density=8.0*g/cm3,ncomponents=10); lnsteel->AddElement(elSi,fractionmass= 0.01); lnsteel->AddElement(elCr,fractionmass= 0.17); lnsteel->AddElement(elMn,fractionmass= 0.02); lnsteel->AddElement(elFe,fractionmass= 0.65195); lnsteel->AddElement(elNi,fractionmass= 0.12); lnsteel->AddElement(elMo,fractionmass= 0.025); lnsteel->AddElement(elP,fractionmass= 0.00045); lnsteel->AddElement(elS,fractionmass= 0.0003); lnsteel->AddElement(elC,fractionmass= 0.0003); lnsteel->AddElement(elN,fractionmass= 0.002); //Magnox alloy (1.8 g cm-3) 99.2 % Mg, 0.8 % Al for GeBaradello endcap G4Material* magnox = new G4Material (name="magnox",density=1.8*g/cm3,ncomponents=2); magnox->AddElement(elMg,fractionmass= 0.992); magnox->AddElement(elAl,fractionmass= 0.008); //stainless steel mesh transfer lines DS50 G4Material* meshsteel = new G4Material (name="MeshSteel",density=5.1565*g/cm3,ncomponents=5); meshsteel->AddElement(elSi,fractionmass= 0.01); meshsteel->AddElement(elCr,fractionmass= 0.20); meshsteel->AddElement(elMn,fractionmass= 0.02); meshsteel->AddElement(elFe,fractionmass= 0.67); meshsteel->AddElement(elNi,fractionmass= 0.10); //step motor, DEC, HM6GT-F00-1A, n. 20111122, OPERA G4Material* motor = new G4Material (name="motor",density=3.6961*g/cm3,ncomponents=5); motor->AddElement(elSi,fractionmass= 0.01); motor->AddElement(elCr,fractionmass= 0.20); motor->AddElement(elMn,fractionmass= 0.02); motor->AddElement(elFe,fractionmass= 0.67); motor->AddElement(elNi,fractionmass= 0.10); //meteorite, MS-179, Almahata Sitta, EL 3-5 breccia, nominal average density 3.55 g/cm3 G4Material* MS179 = new G4Material (name="MS179",density=3.1233*g/cm3,ncomponents=11); MS179->AddElement(elO,fractionmass= 0.367); MS179->AddElement(elFe,fractionmass= 0.210); MS179->AddElement(elNi,fractionmass= 0.012); MS179->AddElement(elSi,fractionmass= 0.235); MS179->AddElement(elMg,fractionmass= 0.124); MS179->AddElement(elCa,fractionmass= 0.010); MS179->AddElement(elAl,fractionmass= 0.010); MS179->AddElement(elS,fractionmass= 0.02); MS179->AddElement(elNa,fractionmass= 0.006); MS179->AddElement(elC,fractionmass= 0.004); MS179->AddElement(elMn,fractionmass= 0.002); //NaI G4Material* NaI = new G4Material (name="NaI",density=3.67*g/cm3,ncomponents=2); NaI->AddElement(elNa,natoms= 1); NaI->AddElement(elI,natoms= 1); //Cu-Ni alloy in pins G4Material* nicup = new G4Material (name="nicup",density=5.22347*g/cm3,ncomponents=3); nicup->AddElement(elCu,fractionmass= 0.9775); nicup->AddElement(elBe,fractionmass= 0.0025); nicup->AddElement(elNi,fractionmass= 0.02); // Nitrogen Gas density = 0.0011467*g/cm3; G4Material* N2Gas = new G4Material(name="NitrogenGas", density, 1); N2Gas -> AddElement(elN, natoms=2); //meteorite, unknown chondrite, Oman, MiB G4Material* Oman = new G4Material (name="Oman",density=3.00*g/cm3,ncomponents=8); Oman->AddElement(elO,fractionmass= 0.380); Oman->AddElement(elFe,fractionmass= 0.220); Oman->AddElement(elSi,fractionmass= 0.190); Oman->AddElement(elMg,fractionmass= 0.150); Oman->AddElement(elS,fractionmass= 0.022); Oman->AddElement(elCa,fractionmass= 0.013); Oman->AddElement(elAl,fractionmass= 0.013); Oman->AddElement(elNi,fractionmass= 0.012); //stainless steel for shell of Omnetics, 37 pins, Micro-D, DS-20K G4Material* omnetics37 = new G4Material (name="omnetics37",density=4.312517*g/cm3,ncomponents=5); omnetics37->AddElement(elSi,fractionmass= 0.018); omnetics37->AddElement(elCr,fractionmass= 0.004); omnetics37->AddElement(elMn,fractionmass= 0.010); omnetics37->AddElement(elFe,fractionmass= 0.962); omnetics37->AddElement(elC,fractionmass= 0.006); //stainless steel for shell of Omnetics, rectangular mount, 37 pins, Micro-D, DS-20K G4Material* omnetics37r = new G4Material (name="omnetics37r",density=5.504219*g/cm3,ncomponents=5); omnetics37r->AddElement(elSi,fractionmass= 0.018); omnetics37r->AddElement(elCr,fractionmass= 0.004); omnetics37r->AddElement(elMn,fractionmass= 0.010); omnetics37r->AddElement(elFe,fractionmass= 0.962); omnetics37r->AddElement(elC,fractionmass= 0.006); // polyethylen terephthalat (PET, Mylar) G4Material* PET = new G4Material (name = "PET", density = 1370 * kg/m3, ncomponents = 3); PET -> AddElement(elC, natoms = 10); PET -> AddElement(elH, natoms = 8); PET -> AddElement(elO, natoms = 4); //Quartzel, SilPower HD, 99.95% SiO2 fibre G4Material* Quartzel = new G4Material (name="Quartzel",density=1.11529*g/cm3,ncomponents=2); Quartzel->AddElement(elSi,natoms= 1); Quartzel->AddElement(elO,natoms= 2); //old ship steel G4Material* shipsteel = new G4Material (name="ShipSteel",density=7.9*g/cm3,ncomponents=3); shipsteel->AddElement(elCr,fractionmass= 0.02); shipsteel->AddElement(elFe,fractionmass= 0.94); shipsteel->AddElement(elNi,fractionmass= 0.04); //SiO2 G4Material* SiO2 = new G4Material (name="SiO2",density=2.00*g/cm3,ncomponents=2); SiO2->AddElement(elSi,natoms= 1); SiO2->AddElement(elO,natoms= 2); //stainless steel for springs (61SiCr7, DIN 60SiCr7, 1.7108) G4Material* springsteel = new G4Material (name="springsteel",density=7.80*g/cm3,ncomponents=5); springsteel->AddElement(elSi,fractionmass= 0.018); springsteel->AddElement(elCr,fractionmass= 0.004); springsteel->AddElement(elMn,fractionmass= 0.010); springsteel->AddElement(elFe,fractionmass= 0.962); springsteel->AddElement(elC,fractionmass= 0.006); //SrI2(Eu), real density 4.58 g/cm3 G4Material* SrI = new G4Material (name="SrI",density=4.30*g/cm3,ncomponents=3); SrI->AddElement(elSr,natoms= 988); SrI->AddElement(elEu,natoms= 12); SrI->AddElement(elI,natoms= 2000); //stainless steel G4Material* ssteel = new G4Material (name="Steel",density=7.892295*g/cm3,ncomponents=5); ssteel->AddElement(elSi,fractionmass= 0.01); ssteel->AddElement(elCr,fractionmass= 0.20); ssteel->AddElement(elMn,fractionmass= 0.02); ssteel->AddElement(elFe,fractionmass= 0.67); ssteel->AddElement(elNi,fractionmass= 0.10); //CRISTEK custom steel shell, DS-20k G4Material* ssteel_cristek = new G4Material (name="Steel_cristek",density=3.253561*g/cm3,ncomponents=5); ssteel_cristek->AddElement(elSi,fractionmass= 0.01); ssteel_cristek->AddElement(elCr,fractionmass= 0.20); ssteel_cristek->AddElement(elMn,fractionmass= 0.02); ssteel_cristek->AddElement(elFe,fractionmass= 0.67); ssteel_cristek->AddElement(elNi,fractionmass= 0.10); //stainless steel, Jingke Xu G4Material* ssteel_mod = new G4Material (name="Steel_mod",density=7.946556*g/cm3,ncomponents=5); // stainless steel, AISI 316, 1.4301, sample n. 27, Xenon1t, GeMPI3 ssteel_mod->AddElement(elSi,fractionmass= 0.01); ssteel_mod->AddElement(elCr,fractionmass= 0.20); ssteel_mod->AddElement(elMn,fractionmass= 0.02); ssteel_mod->AddElement(elFe,fractionmass= 0.67); ssteel_mod->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M8, Xenon1t, lowest layer in box G4Material* ssteel_mod1 = new G4Material (name="Steel_mod1",density=5.891992*g/cm3,ncomponents=5); ssteel_mod1->AddElement(elSi,fractionmass= 0.01); ssteel_mod1->AddElement(elCr,fractionmass= 0.20); ssteel_mod1->AddElement(elMn,fractionmass= 0.02); ssteel_mod1->AddElement(elFe,fractionmass= 0.67); ssteel_mod1->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M8, hexagonal head, M8 1.25 x 40 mm, Xenon1t, lowest & uppermost layer in box G4Material* ssteel_mod1a = new G4Material (name="Steel_mod1a",density=4.355740*g/cm3,ncomponents=5); ssteel_mod1a->AddElement(elSi,fractionmass= 0.01); ssteel_mod1a->AddElement(elCr,fractionmass= 0.20); ssteel_mod1a->AddElement(elMn,fractionmass= 0.02); ssteel_mod1a->AddElement(elFe,fractionmass= 0.67); ssteel_mod1a->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M5, round head, M5 0.8 x 20 mm, Xenon1t, lowest & uppermost layer in box G4Material* ssteel_mod1b = new G4Material (name="Steel_mod1b",density=2.703669*g/cm3,ncomponents=5); ssteel_mod1b->AddElement(elSi,fractionmass= 0.01); ssteel_mod1b->AddElement(elCr,fractionmass= 0.20); ssteel_mod1b->AddElement(elMn,fractionmass= 0.02); ssteel_mod1b->AddElement(elFe,fractionmass= 0.67); ssteel_mod1b->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M6, round head, M6 1.0 x 30 mm, Xenon1t, lowest & uppermost layer in box G4Material* ssteel_mod1c = new G4Material (name="Steel_mod1c",density=3.147026*g/cm3,ncomponents=5); ssteel_mod1c->AddElement(elSi,fractionmass= 0.01); ssteel_mod1c->AddElement(elCr,fractionmass= 0.20); ssteel_mod1c->AddElement(elMn,fractionmass= 0.02); ssteel_mod1c->AddElement(elFe,fractionmass= 0.67); ssteel_mod1c->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M8, Xenon1t, middle layer in box G4Material* ssteel_mod2 = new G4Material (name="Steel_mod2",density=4.918717*g/cm3,ncomponents=5); ssteel_mod2->AddElement(elSi,fractionmass= 0.01); ssteel_mod2->AddElement(elCr,fractionmass= 0.20); ssteel_mod2->AddElement(elMn,fractionmass= 0.02); ssteel_mod2->AddElement(elFe,fractionmass= 0.67); ssteel_mod2->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M8, hexagonal head, M8 1.25 x 40 mm, Xenon1t, middle layer in box G4Material* ssteel_mod2a = new G4Material (name="Steel_mod2a",density=3.150771*g/cm3,ncomponents=5); ssteel_mod2a->AddElement(elSi,fractionmass= 0.01); ssteel_mod2a->AddElement(elCr,fractionmass= 0.20); ssteel_mod2a->AddElement(elMn,fractionmass= 0.02); ssteel_mod2a->AddElement(elFe,fractionmass= 0.67); ssteel_mod2a->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M5, round head, M5 0.8 x 20 mm, Xenon1t, middle layer in box G4Material* ssteel_mod2b = new G4Material (name="Steel_mod2b",density=3.326283*g/cm3,ncomponents=5); ssteel_mod2b->AddElement(elSi,fractionmass= 0.01); ssteel_mod2b->AddElement(elCr,fractionmass= 0.20); ssteel_mod2b->AddElement(elMn,fractionmass= 0.02); ssteel_mod2b->AddElement(elFe,fractionmass= 0.67); ssteel_mod2b->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M6, round head, M6 1.0 x 30 mm, Xenon1t, middle layer in box G4Material* ssteel_mod2c = new G4Material (name="Steel_mod2c",density=4.285839*g/cm3,ncomponents=5); ssteel_mod2c->AddElement(elSi,fractionmass= 0.01); ssteel_mod2c->AddElement(elCr,fractionmass= 0.20); ssteel_mod2c->AddElement(elMn,fractionmass= 0.02); ssteel_mod2c->AddElement(elFe,fractionmass= 0.67); ssteel_mod2c->AddElement(elNi,fractionmass= 0.10); //silver plated stainless steel bolts, M8, Xenon1t, uppermost layer in box G4Material* ssteel_mod3 = new G4Material (name="Steel_mod3",density=5.933780*g/cm3,ncomponents=5); ssteel_mod3->AddElement(elSi,fractionmass= 0.01); ssteel_mod3->AddElement(elCr,fractionmass= 0.20); ssteel_mod3->AddElement(elMn,fractionmass= 0.02); ssteel_mod3->AddElement(elFe,fractionmass= 0.67); ssteel_mod3->AddElement(elNi,fractionmass= 0.10); //silver stainless steel bolts, M8, Xenon1t, 18 bolts on top in box G4Material* ssteel_mod4 = new G4Material (name="Steel_mod4",density=7.9*g/cm3,ncomponents=5); ssteel_mod4->AddElement(elSi,fractionmass= 0.01); ssteel_mod4->AddElement(elCr,fractionmass= 0.20); ssteel_mod4->AddElement(elMn,fractionmass= 0.02); ssteel_mod4->AddElement(elFe,fractionmass= 0.67); ssteel_mod4->AddElement(elNi,fractionmass= 0.10); //stainless steel, Acciaierie Terni, 1st sample G4Material* ssteel_terni_01 = new G4Material (name="Steel_terni_01",density=7.7*g/cm3,ncomponents=5); ssteel_terni_01->AddElement(elSi,fractionmass= 0.01); ssteel_terni_01->AddElement(elCr,fractionmass= 0.20); ssteel_terni_01->AddElement(elMn,fractionmass= 0.02); ssteel_terni_01->AddElement(elFe,fractionmass= 0.67); ssteel_terni_01->AddElement(elNi,fractionmass= 0.10); //stainless steel welding rods, XenonNT, 2614.9 g, GeMPI4 G4Material* steel_rods = new G4Material (name="Steel_rods",density=5.67174*g/cm3,ncomponents=5); steel_rods->AddElement(elSi,fractionmass= 0.01); steel_rods->AddElement(elCr,fractionmass= 0.20); steel_rods->AddElement(elMn,fractionmass= 0.02); steel_rods->AddElement(elFe,fractionmass= 0.67); steel_rods->AddElement(elNi,fractionmass= 0.10); //stainless steel pcs for levelmeter 2, TPC, Xenon G4Material* steeltpc = new G4Material (name="steeltpc",density=7.1049*g/cm3,ncomponents=5); steeltpc->AddElement(elSi,fractionmass= 0.018); steeltpc->AddElement(elCr,fractionmass= 0.004); steeltpc->AddElement(elMn,fractionmass= 0.010); steeltpc->AddElement(elFe,fractionmass= 0.962); steeltpc->AddElement(elC,fractionmass= 0.006); //tektite G4Material* tektite = new G4Material (name="tektite",density=2.45*g/cm3,ncomponents=6); tektite->AddElement(elB,fractionmass= 0.040064); tektite->AddElement(elO,fractionmass= 0.539562); tektite->AddElement(elNa,fractionmass= 0.028191); tektite->AddElement(elAl,fractionmass= 0.011644); tektite->AddElement(elSi,fractionmass= 0.377220); tektite->AddElement(elK,fractionmass= 0.003321); //Cu-Ti alloy in clamps, DS20K G4Material* ticu = new G4Material (name="ticu",density=8.7*g/cm3,ncomponents=2); ticu->AddElement(elCu,fractionmass= 0.97); ticu->AddElement(elTi,fractionmass= 0.03); // Vacuum G4Material* vacuum = new G4Material("Vacuum", z=1, a = 1.008*g/mole, density = 1.e-25*g/cm3, kStateGas, temperature = 2.73*kelvin, pressure = 3.e-18*pascal); //no ordine alfabetico //Sm2O3, before purification, KINR, on GePaolo G4Material* Sm2O3 = new G4Material (name="Sm2O3",density=4.3855*g/cm3,ncomponents=2); Sm2O3->AddElement(elSm,natoms= 2); Sm2O3->AddElement(elO,natoms= 3); //Sm2O3, after purification, KINR, on GeCris G4Material* Sm2O3after = new G4Material (name="Sm2O3after",density=1.087863*g/cm3,ncomponents=2); Sm2O3after->AddElement(elSm,natoms= 2); Sm2O3after->AddElement(elO,natoms= 3); //Yb2(C2O4)3, ytterbium oxalate, after purification, LENS, on GeMPI G4Material* Yb2C6O12 = new G4Material (name="Yb2C6O12",density=0.785547*g/cm3,ncomponents=3); Yb2C6O12->AddElement(elYb,natoms= 2); Yb2C6O12->AddElement(elC,natoms= 6); Yb2C6O12->AddElement(elO,natoms= 12); //Yb2O3, before purification, KINR, on GePaolo G4Material* Yb2O3 = new G4Material (name="Yb2O3",density=2.3112*g/cm3,ncomponents=2); Yb2O3->AddElement(elYb,natoms= 2); Yb2O3->AddElement(elO,natoms= 3); //Yb2O3, after purification, KINR, on GeCris G4Material* Yb2O3after = new G4Material (name="Yb2O3after",density=1.179419*g/cm3,ncomponents=2); Yb2O3after->AddElement(elYb,natoms= 2); Yb2O3after->AddElement(elO,natoms= 3); //Gd2O3, before purification, KINR, on GePaolo G4Material* Gd2O3b = new G4Material (name="Gd2O3b",density=1.90*g/cm3,ncomponents=2); Gd2O3b->AddElement(elGd,natoms= 2); Gd2O3b->AddElement(elO,natoms= 3); //Gd2O3, after purification, KINR, on GePaolo G4Material* Gd2O3a = new G4Material (name="Gd2O3a",density=1.13*g/cm3,ncomponents=2); Gd2O3a->AddElement(elGd,natoms= 2); Gd2O3a->AddElement(elO,natoms= 3); //Gd2O3, sample "A", SOM 04-JUN-2018, before purification, KINR, on GeCris G4Material* Gd2O3c = new G4Material (name="Gd2O3c",density=1.889613*g/cm3,ncomponents=2); Gd2O3c->AddElement(elGd,natoms= 2); Gd2O3c->AddElement(elO,natoms= 3); //Gd2O3, sample "B", SOM 27-JUL-2018, before purification, KINR, on GeCris G4Material* Gd2O3d = new G4Material (name="Gd2O3d",density=1.644811*g/cm3,ncomponents=2); Gd2O3d->AddElement(elGd,natoms= 2); Gd2O3d->AddElement(elO,natoms= 3); //Gd2O3, XenonNT, SOM 22-OCT-2021, three batches, on GeMPI G4Material* Gd2O3e = new G4Material (name="Gd2O3e",density=1.5847*g/cm3,ncomponents=2); Gd2O3e->AddElement(elGd,natoms= 2); Gd2O3e->AddElement(elO,natoms= 3); //Gd2O3, after purification, STS, on GeMPI4 G4Material* Gd2O3_sts = new G4Material (name="Gd2O3_sts",density=1.330597*g/cm3,ncomponents=2); Gd2O3_sts->AddElement(elGd,natoms= 2); Gd2O3_sts->AddElement(elO,natoms= 3); //Gadolinium Sulphate octohydrate min 99.99%, LOT: GDSU19-04, TREIBACHER INDUSTRIE AG, XenonNT, Univ. Mainz, GeMPI2, SOM 08-MAR-2019; G4Material* Gd2O12S3H16O8_mpik = new G4Material (name="Gd2O12S3H16O8_mpik",density=1.757697*g/cm3,ncomponents=4); Gd2O12S3H16O8_mpik->AddElement(elGd,natoms= 2); Gd2O12S3H16O8_mpik->AddElement(elO,natoms= 20); Gd2O12S3H16O8_mpik->AddElement(elS,natoms= 3); Gd2O12S3H16O8_mpik->AddElement(elH,natoms= 16); //Eu2O3, before purification, KINR, on GeMPI3 G4Material* Eu2O3 = new G4Material (name="Eu2O3",density=2.261375*g/cm3,ncomponents=2); Eu2O3->AddElement(elEu,natoms= 2); Eu2O3->AddElement(elO,natoms= 3); //BaWO4, crystal, density 6.34287 g/cm3 G4Material* BaWO4 = new G4Material (name="BaWO4",density=6.34287*g/cm3,ncomponents=3); BaWO4->AddElement(elBa,natoms= 1); BaWO4->AddElement(elW,natoms= 1); BaWO4->AddElement(elO,natoms= 4); //PbWO4, crystal, DAMA, density 8.3530105 g/cm3 (official 8.28 g/cm3) G4Material* PbWO4 = new G4Material (name="PbWO4",density=8.3530105*g/cm3,ncomponents=3); PbWO4->AddElement(elPb,natoms= 1); PbWO4->AddElement(elW,natoms= 1); PbWO4->AddElement(elO,natoms= 4); //Tm3Al5O12 crystal, Univ. San Pietroburgo, density 6.34287 g/cm3 G4Material* Tm3Al5O12 = new G4Material (name="Tm3Al5O12",density=6.34287*g/cm3,ncomponents=3); Tm3Al5O12->AddElement(elTm,natoms= 3); Tm3Al5O12->AddElement(elAl,natoms= 5); Tm3Al5O12->AddElement(elO,natoms= 12); //Gd2(SO4)3, purified, SuperK, INFN-Bari //Kanto, lot 70214Z, 3647.5 g, GeMPI4; //Shin-etsu, lot RGD-OSF-005, 3509.4 g, GeMPI2; G4Material* Gd2S3O12 = new G4Material (name="Gd2S3O12",density=0.853869*g/cm3,ncomponents=3); Gd2S3O12->AddElement(elGd,natoms= 2); Gd2S3O12->AddElement(elO,natoms= 12); Gd2S3O12->AddElement(elS,natoms= 3); //Gd2(SO4)3, Treibacher, XenonNT, 996.14 g, GeMPI2; //Gd2(SO4)3, NYC University, XenonNT, 996.71 g, GeMPI2; G4Material* Gd2S3O12_xnt = new G4Material (name="Gd2S3O12_xnt",density=1.471676*g/cm3,ncomponents=3); Gd2S3O12_xnt->AddElement(elGd,natoms= 2); Gd2S3O12_xnt->AddElement(elO,natoms= 12); Gd2S3O12_xnt->AddElement(elS,natoms= 3); //LiNO3, salt used for Li diffusion in Ge crystals G4Material* LiNO3 = new G4Material (name="LiNO3",density=0.8138*g/cm3,ncomponents=3); LiNO3->AddElement(elLi,natoms= 1); LiNO3->AddElement(elN,natoms= 1); LiNO3->AddElement(elO,natoms= 3); //NaI, Astro grade, 99.999+%, ultra low potassium, // HT-NAI101XG, three times 2 kg, batch: 0000019212, COSINUS, SIGMA-ALDRICH G4Material* NaI_AG = new G4Material (name="NaI_AG",density=1.538959*g/cm3,ncomponents=2); NaI_AG->AddElement(elNa,natoms= 1); NaI_AG->AddElement(elI,natoms= 1); //WO3, powder, Alfa Aesar, Luc. G4Material* WO3 = new G4Material (name="WO3",density=2.5919*g/cm3,ncomponents=2); WO3->AddElement(elW,natoms= 1); WO3->AddElement(elO,natoms= 4); //MoO3, powder, Lumineu, Bordeaux G4Material* MoO3p = new G4Material (name="MoO3p",density=1.391995*g/cm3,ncomponents=2); MoO3p->AddElement(elMo,natoms= 1); MoO3p->AddElement(elO,natoms= 3); //CaWO4, powder, TUM, CRESST G4Material* CaWO4p = new G4Material (name="CaWO4p",density=2.979586*g/cm3,ncomponents=3); CaWO4p->AddElement(elCa,natoms= 1); CaWO4p->AddElement(elW,natoms= 1); CaWO4p->AddElement(elO,natoms= 4); //CaWO4, real density 5.9 - 6.12 g/cm3 G4Material* CaWO4 = new G4Material (name="CaWO4",density=5.8*g/cm3,ncomponents=3); CaWO4->AddElement(elCa,natoms= 1); CaWO4->AddElement(elW,natoms= 1); CaWO4->AddElement(elO,natoms= 4); //CaWO4, real density 5.9 - 6.12 g/cm3 G4Material* CaWO4a = new G4Material (name="CaWO4a",density=6.0*g/cm3,ncomponents=3); CaWO4a->AddElement(elCa,natoms= 1); CaWO4a->AddElement(elW,natoms= 1); CaWO4a->AddElement(elO,natoms= 4); //CdWO4, real density 7.9 g/cm3 G4Material* CdWO4 = new G4Material (name="CdWO4",density=7.9*g/cm3,ncomponents=3); // (name="CdWO4",density=8.0*g/cm3,ncomponents=3); CdWO4->AddElement(elCd,natoms= 1); CdWO4->AddElement(elW,natoms= 1); CdWO4->AddElement(elO,natoms= 4); //CdWO4, scraps, density 6.859 g/cm3 G4Material* CdWO4scraps = new G4Material (name="CdWO4scraps",density=6.859*g/cm3,ncomponents=3); CdWO4scraps->AddElement(elCd,natoms= 1); CdWO4scraps->AddElement(elW,natoms= 1); CdWO4scraps->AddElement(elO,natoms= 4); //CeF3, real density 6.16 g/cm3 G4Material* CeF3 = new G4Material (name="CeF3",density=6.125*g/cm3,ncomponents=2); CeF3->AddElement(elCe,natoms= 1); CeF3->AddElement(elF,natoms= 3); //CdO, from 116CdWO4 scraps, DAMA, KINR G4Material* CdO = new G4Material (name="CdO",density=1.411*g/cm3,ncomponents=2); CdO->AddElement(elCd,natoms= 1); CdO->AddElement(elO,natoms= 1); //CeO2, jeweller's rouge, Darkside-50 G4Material* CeO2 = new G4Material (name="CeO2",density=1.38*g/cm3,ncomponents=2); CeO2->AddElement(elCe,natoms= 1); CeO2->AddElement(elO,natoms= 2); //CeO2, real density 7.65 g/cm3, Stanford, before purification G4Material* CeO2bp = new G4Material (name="CeO2bp",density=2.518*g/cm3,ncomponents=2); CeO2bp->AddElement(elCe,natoms= 1); CeO2bp->AddElement(elO,natoms= 2); //CeO2, real density 7.65 g/cm3, Stanford, after purification G4Material* CeO2ap = new G4Material (name="CeO2ap",density=1.7011*g/cm3,ncomponents=2); CeO2ap->AddElement(elCe,natoms= 1); CeO2ap->AddElement(elO,natoms= 2); //CeO2, Stanford, after purification, long measurement, KINR, GeCris G4Material* CeO2long = new G4Material (name="CeO2long",density=2.558*g/cm3,ncomponents=2); CeO2long->AddElement(elCe,natoms= 1); CeO2long->AddElement(elO,natoms= 2); //CeO2, Stanford, after second purification, 626.94 g, long measurement, KINR, GeCris G4Material* CeO2long2 = new G4Material (name="CeO2long2",density=2.19*g/cm3,ncomponents=2); CeO2long2->AddElement(elCe,natoms= 1); CeO2long2->AddElement(elO,natoms= 2); //CeO2, SOX, density 4.0 g/cm^3 G4Material* CeO2_40 = new G4Material (name="CeO2_40",density=4.0*g/cm3,ncomponents=2); CeO2_40->AddElement(elCe,natoms= 1); CeO2_40->AddElement(elO,natoms= 2); //CeO2, SOX, density 4.5 g/cm^3 G4Material* CeO2_45 = new G4Material (name="CeO2_45",density=4.5*g/cm3,ncomponents=2); CeO2_45->AddElement(elCe,natoms= 1); CeO2_45->AddElement(elO,natoms= 2); //CeO2, SOX, density 5.0 g/cm^3 G4Material* CeO2_50 = new G4Material (name="CeO2_50",density=5.0*g/cm3,ncomponents=2); CeO2_50->AddElement(elCe,natoms= 1); CeO2_50->AddElement(elO,natoms= 2); //Gd(THD)3, Tris(2,2,6,6-tetramethyl-3,5-heptanedionato)gadolinium(III), Linear Formula Gd(C11H19O2)3, PS box D 50 mm G4Material* GdTHD3small = new G4Material (name="GdTHD3small",density=0.570525*g/cm3,ncomponents=4); GdTHD3small->AddElement(elGd,natoms= 1); GdTHD3small->AddElement(elO,natoms= 33); GdTHD3small->AddElement(elC,natoms= 6); GdTHD3small->AddElement(elH,natoms= 57); //Gd(THD)3, Tris(2,2,6,6-tetramethyl-3,5-heptanedionato)gadolinium(III), Linear Formula Gd(C11H19O2)3, PS box D 70 mm G4Material* GdTHD3big = new G4Material (name="GdTHD3big",density=0.575371*g/cm3,ncomponents=4); GdTHD3big->AddElement(elGd,natoms= 1); GdTHD3big->AddElement(elO,natoms= 33); GdTHD3big->AddElement(elC,natoms= 6); GdTHD3big->AddElement(elH,natoms= 57); //Gd(THD)3, Tris(2,2,6,6-tetramethyl-3,5-heptanedionato)gadolinium(III), Linear Formula Gd(C11H19O2)3, Marinelli 438G, DoubleCHOOZ G4Material* GdTHD3bm = new G4Material (name="GdTHD3bm",density=0.833645*g/cm3,ncomponents=4); GdTHD3bm->AddElement(elGd,natoms= 1); GdTHD3bm->AddElement(elO,natoms= 33); GdTHD3bm->AddElement(elC,natoms= 6); GdTHD3bm->AddElement(elH,natoms= 57); //Gd2O3, real density 7.407 g/cm3, Stanford, before purification G4Material* Gd2O3bp = new G4Material (name="Gd2O3bp",density=2.05426*g/cm3,ncomponents=2); Gd2O3bp->AddElement(elGd,natoms= 2); Gd2O3bp->AddElement(elO,natoms= 3); //Gd2O3, real density 7.407 g/cm3, Stanford, after purification G4Material* Gd2O3ap = new G4Material (name="Gd2O3ap",density=1.50081*g/cm3,ncomponents=2); Gd2O3ap->AddElement(elGd,natoms= 2); Gd2O3ap->AddElement(elO,natoms= 3); // chromium G4Material* metalCr = new G4Material (name="MetalChromium",density=7.140*g/cm3,ncomponents=1); metalCr->AddElement(elCr,natoms=1); // terbium G4Material* metalTb = new G4Material (name="MetalTerbium",density=8.23*g/cm3,ncomponents=1); metalTb->AddElement(elTb,natoms=1); // terbium(III,IV) oxide, theoretical density as solid: 7.3 g/cm3, here in powder form G4Material* Tb4O7 = new G4Material (name="Terbiumoxide",density=1.591549*g/cm3,ncomponents=2); Tb4O7->AddElement(elTb,natoms=4); Tb4O7->AddElement(elO,natoms=7); // indium(III) oxide, theoretical density as solid: 7.18 g/cm3 G4Material* In2O3 = new G4Material (name="Indiumoxide",density=7.29412*g/cm3,ncomponents=2); In2O3->AddElement(elIn,natoms=2); In2O3->AddElement(elO,natoms=3); // copper G4Material* metalCu = new G4Material (name="MetalCopper",density=8.960*g/cm3,ncomponents=1); metalCu->AddElement(elCu,natoms=1); // copper from 3D printer) G4Material* metalCu3D = new G4Material (name="MetalCopper3D",density=7.577559*g/cm3,ncomponents=1); metalCu3D->AddElement(elCu,natoms=1); // Niobium G4Material* metalNb = new G4Material (name="MetalNiobium",density=8.57*g/cm3,ncomponents=1); metalNb->AddElement(elNb,natoms=1); // Vanadium ingots G4Material* metalV = new G4Material (name="MetalVanadium",density=6.11*g/cm3,ncomponents=1); metalV->AddElement(elV,natoms=1); // Vanadium ingots G4Material* metalV1 = new G4Material (name="MetalVanadium1",density=6.0163*g/cm3,ncomponents=1); metalV1->AddElement(elV,natoms=1); // Vanadium ingots G4Material* metalV2 = new G4Material (name="MetalVanadium2",density=5.888151*g/cm3,ncomponents=1); metalV2->AddElement(elV,natoms=1); // Vanadium granules in small Marinelli (141G) G4Material* metalVg = new G4Material (name="MetalVanadiumg",density=1.961*g/cm3,ncomponents=1); metalVg->AddElement(elV,natoms=1); // Vanadium disks purified, #1 ยง(n. 5, average 76 g) G4Material* metalVd1 = new G4Material (name="MetalVanadiumD1",density=5.80502*g/cm3,ncomponents=1); metalVd1->AddElement(elV,natoms=1); // Vanadium disks purified, #2 (n. 2, average 73 g) G4Material* metalVd2 = new G4Material (name="MetalVanadiumD2",density=5.61477*g/cm3,ncomponents=1); metalVd2->AddElement(elV,natoms=1); // Vanadium disks purified, #3 (n. 2, average 112 g) G4Material* metalVd3 = new G4Material (name="MetalVanadiumD3",density=5.88475*g/cm3,ncomponents=1); metalVd3->AddElement(elV,natoms=1); // Vanadium disks purified, #4 (n. 1, average 66.55 g) G4Material* metalVd4 = new G4Material (name="MetalVanadiumD4",density=5.6945*g/cm3,ncomponents=1); metalVd4->AddElement(elV,natoms=1); // Vanadium disks purified, #1 (73.80 g) G4Material* metalVd01 = new G4Material (name="MetalVanadiumDisk01",density=5.814532*g/cm3,ncomponents=1); metalVd01->AddElement(elV,natoms=1); // Vanadium disks purified, #2 (75.15 g) G4Material* metalVd02 = new G4Material (name="MetalVanadiumDisk02",density=6.071901*g/cm3,ncomponents=1); metalVd02->AddElement(elV,natoms=1); // Vanadium disks purified, #3 (76.45 g) G4Material* metalVd03 = new G4Material (name="MetalVanadiumDisk03",density=5.978545*g/cm3,ncomponents=1); metalVd03->AddElement(elV,natoms=1); // Vanadium disks purified, #4 (66.55 g) G4Material* metalVd04 = new G4Material (name="MetalVanadiumDisk04",density=5.846168*g/cm3,ncomponents=1); metalVd04->AddElement(elV,natoms=1); // Vanadium disks purified, #5 (76.02 g) G4Material* metalVd05 = new G4Material (name="MetalVanadiumDisk05",density=5.900746*g/cm3,ncomponents=1); metalVd05->AddElement(elV,natoms=1); // Vanadium disks purified, #6 (75.95 g) G4Material* metalVd06 = new G4Material (name="MetalVanadiumDisk06",density=5.984672*g/cm3,ncomponents=1); metalVd06->AddElement(elV,natoms=1); // Vanadium disks purified, #7 (115.60 g) G4Material* metalVd07 = new G4Material (name="MetalVanadiumDisk07",density=6.005607*g/cm3,ncomponents=1); metalVd07->AddElement(elV,natoms=1); // Vanadium disks purified, #8 (111.50 g) G4Material* metalVd08 = new G4Material (name="MetalVanadiumDisk08",density=6.056584*g/cm3,ncomponents=1); metalVd08->AddElement(elV,natoms=1); // Vanadium disks purified, #9 (75.14 g) G4Material* metalVd09 = new G4Material (name="MetalVanadiumDisk09",density=5.890763*g/cm3,ncomponents=1); metalVd09->AddElement(elV,natoms=1); // Vanadium disks purified, #10 (72.31 g) G4Material* metalVd10 = new G4Material (name="MetalVanadiumDisk10",density=5.725767*g/cm3,ncomponents=1); metalVd10->AddElement(elV,natoms=1); // Zirconium disk Prof. Pavel P. Povinec, 129.94 g, D 40.1 mm , h 5.18 mm G4Material* metalZrPPP = new G4Material (name="MetalZirconiumDiskPPP",density=6.501109*g/cm3,ncomponents=1); metalZrPPP->AddElement(elZr,natoms=1); // Zirconium disks purified, #1 (89.30 g) G4Material* metalZrd01 = new G4Material (name="MetalZirconiumDisk01",density=6.458434*g/cm3,ncomponents=1); metalZrd01->AddElement(elZr,natoms=1); // Zirconium disks purified, #2 (88.38 g) G4Material* metalZrd02 = new G4Material (name="MetalZirconiumDisk02",density=6.346214*g/cm3,ncomponents=1); metalZrd02->AddElement(elZr,natoms=1); // Zirconium disks purified, #3 (86.28 g) G4Material* metalZrd03 = new G4Material (name="MetalZirconiumDisk03",density=6.254992*g/cm3,ncomponents=1); metalZrd03->AddElement(elZr,natoms=1); // Zirconium disks purified, #4 (88.95 g) G4Material* metalZrd04 = new G4Material (name="MetalZirconiumDisk04",density=6.449669*g/cm3,ncomponents=1); metalZrd04->AddElement(elZr,natoms=1); // Zirconium disks purified, #5 (90.10 g) G4Material* metalZrd05 = new G4Material (name="MetalZirconiumDisk05",density=6.468689*g/cm3,ncomponents=1); metalZrd05->AddElement(elZr,natoms=1); // Zirconium disks purified, #6 (90.22 g) G4Material* metalZrd06 = new G4Material (name="MetalZirconiumDisk06",density=6.492829*g/cm3,ncomponents=1); metalZrd06->AddElement(elZr,natoms=1); // Zirconium disks purified, #7 (88.80 g) G4Material* metalZrd07 = new G4Material (name="MetalZirconiumDisk07",density=6.454224*g/cm3,ncomponents=1); metalZrd07->AddElement(elZr,natoms=1); // Zirconium disks purified, #8 (114.32 g) G4Material* metalZrd08 = new G4Material (name="MetalZirconiumDisk08",density=6.374529*g/cm3,ncomponents=1); metalZrd08->AddElement(elZr,natoms=1); // Palladium plates Pd1 in PS box, TUD G4Material* metalPd1 = new G4Material (name="MetalPalladium1",density=7.216082*g/cm3,ncomponents=1); metalPd1->AddElement(elPd,natoms=1); // Palladium plates Pd2 in PS box, TUD G4Material* metalPd2 = new G4Material (name="MetalPalladium2",density=7.594563*g/cm3,ncomponents=1); metalPd2->AddElement(elPd,natoms=1); // Selenium granules in PS box, Luc G4Material* metalSeg = new G4Material (name="MetalSeleniumg",density=2.488*g/cm3,ncomponents=1); metalSeg->AddElement(elSe,natoms=1); // Selenium granules in very small Marinelli (541G), Luc G4Material* metalSeg541 = new G4Material (name="MetalSeleniumg541",density=2.949*g/cm3,ncomponents=1); metalSeg541->AddElement(elSe,natoms=1); // Selenium granules in small Marinelli (141G), Luc G4Material* metalSeg141 = new G4Material (name="MetalSeleniumg141",density=2.646129*g/cm3,ncomponents=1); metalSeg141->AddElement(elSe,natoms=1); // Zinc granules in PS box, Luc G4Material* metalZng = new G4Material (name="MetalZincg",density=2.984*g/cm3,ncomponents=1); metalZng->AddElement(elZn,natoms=1); // Zinc granules, purified, in big Marinelli (441G), Luc G4Material* metalZng441 = new G4Material (name="MetalZincg441",density=3.62778*g/cm3,ncomponents=1); metalZng441->AddElement(elZn,natoms=1); // Zincselenide in PS box, Luc G4Material* ZnSe = new G4Material (name="ZnSe",density=2.7534*g/cm3,ncomponents=2); ZnSe->AddElement(elZn,natoms=1); ZnSe->AddElement(elSe,natoms=1); // Zincselenide powder in PS box, Luc G4Material* ZnSep = new G4Material (name="ZnSep",density=2.1007*g/cm3,ncomponents=2); ZnSep->AddElement(elZn,natoms=1); ZnSep->AddElement(elSe,natoms=1); // Zincselenide powder from Novosibirsk in PS box, Luc G4Material* ZnSen = new G4Material (name="ZnSen",density=3.5238*g/cm3,ncomponents=2); ZnSen->AddElement(elZn,natoms=1); ZnSen->AddElement(elSe,natoms=1); // Platinum G4Material* metalPt = new G4Material (name="MetalPlatinum",density=21.45*g/cm3,ncomponents=1); metalPt->AddElement(elPt,natoms=1); // Hafnium_mod (foil placed around crystal in detector GS1) G4Material* metalHf_mod = new G4Material (name="GS1_Hafnium",density=13.31*g/cm3,ncomponents=1); metalHf_mod->AddElement(elHf,natoms=1); // Platinum_mod (foil placed around crystal in detector GS1) G4Material* metalPt_mod = new G4Material (name="GS1_Platinum",density=20.33*g/cm3,ncomponents=1); metalPt_mod->AddElement(elPt,natoms=1); // Osmium (nominal density CRC 2003: 22.57 g/cm3) G4Material* metalOs = new G4Material (name="MetalOsmium",density=23.52*g/cm3,ncomponents=1); metalOs->AddElement(elOs,natoms=1); // Osmium (nominal density CRC 2003: 22.57 g/cm3) G4Material* metalOs_nominal = new G4Material (name="MetalOsmiumNominal",density=22.57*g/cm3,ncomponents=1); metalOs_nominal->AddElement(elOs,natoms=1); // magnets DAMA (Nd2Fe14B) // magnet sample #1, 3 pcs, l 2.4 cm, D 0.4 cm, 2.2699 g; horizontally; density = 7.526371 g cm^-3; G4Material* magnet_1_dama = new G4Material (name="magnet_1_dama", density=7.526371*g/cm3, ncomponents=3); magnet_1_dama->AddElement(elNd,natoms= 2); magnet_1_dama->AddElement(elFe,natoms= 14); magnet_1_dama->AddElement(elB,natoms= 1); // magnet sample #2, 3 pcs, l 3.0 cm, D 0.3 cm, 1.5717 g; horizontally; density = 7.411669 g cm^-3; G4Material* magnet_2_dama = new G4Material (name="magnet_2_dama", density=7.411669*g/cm3, ncomponents=3); magnet_2_dama->AddElement(elNd,natoms= 2); magnet_2_dama->AddElement(elFe,natoms= 14); magnet_2_dama->AddElement(elB,natoms= 1); // magnet sample #3, 3 pcs, l 0.6 cm, D 0.7 cm, 1.7451 g; vertically; density = 7.557586 g cm^-3; G4Material* magnet_3_dama = new G4Material (name="magnet_3_dama", density=7.557586*g/cm3, ncomponents=3); magnet_3_dama->AddElement(elNd,natoms= 2); magnet_3_dama->AddElement(elFe,natoms= 14); magnet_3_dama->AddElement(elB,natoms= 1); // magnet sample #4, 3 pcs, l 0.45 cm, D 0.5 cm, 0.6568 g; vertically; density = 7.433455 g cm^-3; G4Material* magnet_4_dama = new G4Material (name="magnet_4_dama", density=7.433455*g/cm3, ncomponents=3); magnet_4_dama->AddElement(elNd,natoms= 2); magnet_4_dama->AddElement(elFe,natoms= 14); magnet_4_dama->AddElement(elB,natoms= 1); // magnet sample #5, 2 spheres, D 0.5 cm each, 0.9834 g; density = 7.512623 g cm^-3; G4Material* magnet_5_dama = new G4Material (name="magnet_5_dama", density=7.512623*g/cm3, ncomponents=3); magnet_5_dama->AddElement(elNd,natoms= 2); magnet_5_dama->AddElement(elFe,natoms= 14); magnet_5_dama->AddElement(elB,natoms= 1); // magnet sample #6, 6 spheres, D 0.3 cm each, 0.6212 g; horizontally; density = 7.323485 g cm^-3; G4Material* magnet_6_dama = new G4Material (name="magnet_6_dama", density=7.323485*g/cm3, ncomponents=3); magnet_6_dama->AddElement(elNd,natoms= 2); magnet_6_dama->AddElement(elFe,natoms= 14); magnet_6_dama->AddElement(elB,natoms= 1); // magnet sample #7, 4 pcs, l 0.8 cm, D 0.2 cm, 0.1983 g; horizontally; density = 7.890106 g cm^-3; G4Material* magnet_7_dama = new G4Material (name="magnet_7_dama", density=7.890106*g/cm3, ncomponents=3); magnet_7_dama->AddElement(elNd,natoms= 2); magnet_7_dama->AddElement(elFe,natoms= 14); magnet_7_dama->AddElement(elB,natoms= 1); //water G4Material* water = new G4Material (name="Water", density=1.00*g/cm3, ncomponents=2); water->AddElement(elH,natoms=2); water->AddElement(elO,natoms=1); //sea water G4Material* seawater = new G4Material (name="SeaWater", density=1.02*g/cm3, ncomponents=2); seawater->AddElement(elH,natoms=2); seawater->AddElement(elO,natoms=1); //treated sea water G4Material* tseawater = new G4Material (name="TSeaWater", density=1.09*g/cm3, ncomponents=2); tseawater->AddElement(elH,natoms=2); tseawater->AddElement(elO,natoms=1); //Nd2O3 G4Material* nd2o3 = new G4Material (name="nd2o3", density=0.835*g/cm3, ncomponents=2); nd2o3->AddElement(elNd,natoms=2); nd2o3->AddElement(elO,natoms=3); //Nd2O3, Barabash, sample #1, purified G4Material* nd2o3a = new G4Material (name="nd2o3a", density=0.97954*g/cm3, ncomponents=2); nd2o3a->AddElement(elNd,natoms=2); nd2o3a->AddElement(elO,natoms=3); //Nd2O3, Barabash, sample #1, not purified G4Material* nd2o3b = new G4Material (name="nd2o3b", density=1.0674*g/cm3, ncomponents=2); nd2o3b->AddElement(elNd,natoms=2); nd2o3b->AddElement(elO,natoms=3); // lead G4Material* metalPb = new G4Material (name="MetalLead", density=11.340*g/cm3, ncomponents=1); metalPb->AddElement(elPb,natoms=1); // Roman lead G4Material* metalRPb = new G4Material (name="MetalRlead", density=11.16127*g/cm3, ncomponents=1); metalRPb->AddElement(elPb,natoms=1); // lead & antimony G4Material* metalPbSb = new G4Material (name="MetalLeadSb", density=11.052*g/cm3, ncomponents=2); metalPbSb->AddElement(elPb,fractionmass= 0.98); metalPbSb->AddElement(elSb,fractionmass= 0.02); // lead free solder, SnCu G4Material* metalSnCu = new G4Material (name="metalSnCu", density=2.5347*g/cm3, ncomponents=2); metalSnCu->AddElement(elSn,fractionmass= 0.993); metalSnCu->AddElement(elCu,fractionmass= 0.007); // Mumetal, ASTM 753 Alloy 4 G4Material* mumetal = new G4Material (name="mumetal", density=8.7*g/cm3, ncomponents=5); mumetal->AddElement(elNi,fractionmass= 0.8000); mumetal->AddElement(elFe,fractionmass= 0.1450); mumetal->AddElement(elMo,fractionmass= 0.0500); mumetal->AddElement(elSi,fractionmass= 0.0048); mumetal->AddElement(elCu,fractionmass= 0.0002); // lead free solder, Heraeus Microbond, 91121 P, SAC 105, Sn96.5/Ag3.0/Cu0.5, Xenon, TPC G4Material* solder_xenon = new G4Material (name="solder_xenon", density=4.50704*g/cm3, ncomponents=3); solder_xenon->AddElement(elSn,fractionmass= 0.965); solder_xenon->AddElement(elAg,fractionmass= 0.030); solder_xenon->AddElement(elCu,fractionmass= 0.005); // near eutectic lead free solder, Caplinq, Lot #LNCA17A2307BF, 0.5 mm diameter, Bi57/Sn42/Ag1, DS-20K G4Material* solder_ds20k = new G4Material (name="solder_ds20k", density=4.836295*g/cm3, ncomponents=3); solder_ds20k->AddElement(elBi,fractionmass= 0.57); solder_ds20k->AddElement(elSn,fractionmass= 0.42); solder_ds20k->AddElement(elAg,fractionmass= 0.01); // STANNOL, S-Sn62Pb36Ag2, type HS10, Lot 201778/46, 0.5 mm diameter, LUNA-400 G4Material* solder_luna400 = new G4Material (name="solder_luna400", density=4.371456*g/cm3, ncomponents=3); solder_luna400->AddElement(elSn,fractionmass= 0.62); solder_luna400->AddElement(elPb,fractionmass= 0.36); solder_luna400->AddElement(elAg,fractionmass= 0.02); // thorium G4Material* metalTh = new G4Material (name="MetalThorium", density=11.72*g/cm3, ncomponents=1); metalTh->AddElement(elTh,natoms=1); // tungsten G4Material* metalW = new G4Material (name="MetalTungsten", density=19.250*g/cm3, ncomponents=1); metalW->AddElement(elW,natoms=1); // LUMA, Sweden, tungsten wire, gold plated, d 25 micron, Lot 20rE988-U9, plated weight 1.24, base weight 1.17, Au 85% G4Material* AuWWire = new G4Material (name="AuWWire", density=19.17*g/cm3, ncomponents=2); AuWWire->AddElement(elAu,fractionmass= 0.05); AuWWire->AddElement(elW,fractionmass= 0.95); // tungsten powder CeLAND G4Material* powderW = new G4Material (name="PowderTungsten", density=18.0*g/cm3, ncomponents=1); powderW->AddElement(elW,natoms=1); // HP tungsten powder in Saclay bottle, CeLAND G4Material* powderbHPW = new G4Material (name="PowderBHPTungsten", density=4.067925*g/cm3, ncomponents=1); powderbHPW->AddElement(elW,natoms=1); // tungsten powder in Saclay bottle, CeLAND G4Material* powderbW = new G4Material (name="PowderBTungsten", density=5.637015*g/cm3, ncomponents=1); powderbW->AddElement(elW,natoms=1); // iron powder in Saclay bottle, CeLAND G4Material* powderbFe = new G4Material (name="PowderBIron", density=3.403911*g/cm3, ncomponents=1); powderbFe->AddElement(elFe,natoms=1); // nickel powder in Saclay bottle, CeLAND G4Material* powderbNi = new G4Material (name="PowderBNickel", density=2.679352*g/cm3, ncomponents=1); powderbNi->AddElement(elNi,natoms=1); // copper powder in Saclay bottle, CeLAND G4Material* powderbCu = new G4Material (name="PowderBCopper", density=2.88196*g/cm3, ncomponents=1); powderbCu->AddElement(elCu,natoms=1); // copper rods for TPC in GeMPI2, Xenon G4Material* TPCCu = new G4Material (name="TPCCu", density=5.542*g/cm3, ncomponents=1); TPCCu->AddElement(elCu,natoms=1); // tungsten alloys Plansee, CeLAND G4Material* alloyW = new G4Material (name="AlloyTungsten", density=18.327*g/cm3, ncomponents=3); alloyW->AddElement(elW,fractionmass= 0.96); alloyW->AddElement(elNi,fractionmass= 0.02); alloyW->AddElement(elCu,fractionmass= 0.02); // thoriated tungsten wire (W 99/ Th 1), Cuore G4Material* wireWTh = new G4Material (name="Tungstenwire", density=5.391*g/cm3, ncomponents=2); wireWTh->AddElement(elW,fractionmass= 0.99); wireWTh->AddElement(elTh,fractionmass= 0.01); // copper wire TUM, Gerda (Jozsef Janicsko) G4Material* wireCu = new G4Material (name="Copperwire", density=6.7845743*g/cm3, ncomponents=1); wireCu->AddElement(elCu,fractionmass= 1.00); // getter for Hamamatsu PMT R11410-21, Xenon G4Material* getter = new G4Material (name="getter", density=2.6956*g/cm3, ncomponents=1); getter->AddElement(elBa,natoms=1); // 100 getters for Hamamatsu PMT R11410-21, Xenon ,sample 13, 2016/2017 G4Material* getter02 = new G4Material (name="getter02", density=2.01288*g/cm3, ncomponents=1); getter02->AddElement(elBa,natoms=1); // tantalum G4Material* metalTa = new G4Material (name="MetalTantalum", density=16.69*g/cm3, ncomponents=1); metalTa->AddElement(elTa,natoms=1); // tellurium G4Material* metalTe = new G4Material (name="MetalTellurium", density=6.24*g/cm3, ncomponents=1); metalTe->AddElement(elTa,natoms=1); //aluminium G4Material* metalAl = new G4Material (name="MetalAluminium", density=2.702*g/cm3, ncomponents=1); // (name="MetalAluminium", density=0.52*g/cm3, ncomponents=1); metalAl->AddElement(elAl,natoms=1); //superinsulation, DS50 //aluminium G4Material* metalAl_a = new G4Material (name="MetalAluminium_a", density=0.046*g/cm3, ncomponents=1); metalAl_a->AddElement(elAl,natoms=1); //superinsulation, DS50, 146400-003, embossed //aluminium G4Material* metalAl_b = new G4Material (name="MetalAluminium_b", density=0.934*g/cm3, ncomponents=1); metalAl_b->AddElement(elAl,natoms=1); //iron G4Material* metalFe = new G4Material (name="MetalIron", density=7.874*g/cm3, ncomponents=1); metalFe->AddElement(elFe,natoms=1); //silicon G4Material* metalSi = new G4Material (name="MetalSilicon", density=2.330*g/cm3, ncomponents=1); metalSi->AddElement(elSi,natoms=1); //Asics for Cobra (silicon) G4Material* ASICSi = new G4Material (name="ASICSilicon", density=1.11403*g/cm3, ncomponents=1); ASICSi->AddElement(elSi,natoms=1); //100 FETs for MPIK (silicon) G4Material* FETSi = new G4Material (name="FETSilicon", density=1.12293*g/cm3, ncomponents=1); FETSi->AddElement(elSi,natoms=1); //MPPC (Multi Pixel Photon Counter) HPD (APD, Si-PM) S13371-6050CQ-02/NG, Hamamatsu for XenonNT (silicon) G4Material* HPDSi = new G4Material (name="HPDSilicon", density=2.291128*g/cm3, ncomponents=1); HPDSi->AddElement(elSi,natoms=1); //Silicone paste for GeMPI3 (oxysilane polymer) G4Material* silicone = new G4Material (name="Silicone", density=1.23255*g/cm3, ncomponents=5); silicone->AddElement(elSi,fractionmass=0.1267); silicone->AddElement(elO,fractionmass= 0.2172); silicone->AddElement(elN,fractionmass= 0.0633); silicone->AddElement(elH,fractionmass= 0.1041); silicone->AddElement(elC,fractionmass= 0.4887); //molecular sieve, 3A, 2-5 mm beads, COSINUS G4Material* molsieve3A = new G4Material (name="molecular_sieve_3A", density=0.747003*g/cm3, ncomponents=6); molsieve3A->AddElement(elH,natoms= 18); molsieve3A->AddElement(elO,natoms= 429); molsieve3A->AddElement(elNa,natoms= 80); molsieve3A->AddElement(elK,natoms= 120); molsieve3A->AddElement(elAl,natoms= 200); molsieve3A->AddElement(elSi,natoms= 200); //silicon PMs, LEGEND, TUM G4Material* sipms_legend = new G4Material (name="sipms_legend", density=2.469136*g/cm3, ncomponents=1); sipms_legend->AddElement(elSi,natoms=1); //silicon PMs, XENONnT, MPI-K-HD G4Material* sipms_mpik_xenon = new G4Material (name="sipms_mpik_xenon", density=2.5*g/cm3, ncomponents=1); sipms_mpik_xenon->AddElement(elSi,natoms=1); //timing belt Brecoflex BrecoProtect, Xenon, metal enforced polyurethane G4Material* breco = new G4Material (name="breco", density=1.356815*g/cm3, ncomponents=4); breco->AddElement(elO,fractionmass= 0.1622); breco->AddElement(elN,fractionmass= 0.0946); breco->AddElement(elH,fractionmass= 0.0540); breco->AddElement(elC,fractionmass= 0.6892); //timing belt MegaDyne MXC, 73.6, 2012, OPERA, metal enforced polyurethane G4Material* megadyne = new G4Material (name="megadyne", density=1.371079*g/cm3, ncomponents=4); megadyne->AddElement(elO,fractionmass= 0.1622); megadyne->AddElement(elN,fractionmass= 0.0946); megadyne->AddElement(elH,fractionmass= 0.0540); megadyne->AddElement(elC,fractionmass= 0.6892); //timing belt Bando 160XL BF6, OPERA, metal enforced polyurethane G4Material* bando = new G4Material (name="bando", density=1.141535*g/cm3, ncomponents=4); bando->AddElement(elO,fractionmass= 0.1622); bando->AddElement(elN,fractionmass= 0.0946); bando->AddElement(elH,fractionmass= 0.0540); bando->AddElement(elC,fractionmass= 0.6892); //pins Xenon G4Material* pinsb = new G4Material (name="pinsb", density=1.2841*g/cm3, ncomponents=1); pinsb->AddElement(elAl,natoms=1); //pins cresst G4Material* pinsc = new G4Material (name="pinsc", density=2.549*g/cm3, ncomponents=1); pinsc->AddElement(elAl,natoms=1); //pins DS50 G4Material* pinsd = new G4Material (name="pinsd", density=2.3884*g/cm3, ncomponents=1); pinsd->AddElement(elAl,natoms=1); //pins and receipts, #1, Cu-Be, Gerda (MIB - Cattadori) G4Material* pinse = new G4Material (name="pinse",density=2.4788*g/cm3,ncomponents=3); pinse->AddElement(elCu,fractionmass= 0.9775); pinse->AddElement(elBe,fractionmass= 0.02); pinse->AddElement(elNi,fractionmass= 0.0025); G4Material* rece = new G4Material (name="rece",density=2.929*g/cm3,ncomponents=3); rece->AddElement(elCu,fractionmass= 0.9775); rece->AddElement(elBe,fractionmass= 0.02); rece->AddElement(elNi,fractionmass= 0.0025); //pins and receipts, #2, Cu-Be, Gerda (MIB - Cattadori) G4Material* pinsf = new G4Material (name="pinsf",density=3.8970*g/cm3,ncomponents=3); pinsf->AddElement(elCu,fractionmass= 0.9775); pinsf->AddElement(elBe,fractionmass= 0.02); pinsf->AddElement(elNi,fractionmass= 0.0025); G4Material* recf = new G4Material (name="recf",density=2.1155*g/cm3,ncomponents=3); recf->AddElement(elCu,fractionmass= 0.9775); recf->AddElement(elBe,fractionmass= 0.02); recf->AddElement(elNi,fractionmass= 0.0025); //pins and receipts, #3, Cu-Be, Gerda (TUM, Janicsko) G4Material* pinsg = new G4Material (name="pinsg",density=2.1967*g/cm3,ncomponents=3); pinsg->AddElement(elCu,fractionmass= 0.9775); pinsg->AddElement(elBe,fractionmass= 0.02); pinsg->AddElement(elNi,fractionmass= 0.0025); //pins and receipts, #4, bronze, Gerda (MPIK, Schwingenheuer) G4Material* pinsh = new G4Material (name="pinsh",density=3.2926*g/cm3,ncomponents=2); pinsh->AddElement(elCu,fractionmass= 0.95); pinsh->AddElement(elSn,fractionmass= 0.05); G4Material* rech = new G4Material (name="rech",density=1.7384*g/cm3,ncomponents=2); rech->AddElement(elCu,fractionmass= 0.95); rech->AddElement(elSn,fractionmass= 0.05); //pins and receipts, #5, copper alloy tin coated, Gerda (MPIK, Schwingenheuer) G4Material* pinsi = new G4Material (name="pinsi",density=1.59155*g/cm3,ncomponents=3); pinsi->AddElement(elCu,fractionmass= 0.9775); pinsi->AddElement(elBe,fractionmass= 0.02); pinsi->AddElement(elNi,fractionmass= 0.0025); //pins and receipts (sockets), Xenon G4Material* pinsj = new G4Material (name="pinsj",density=2.277409*g/cm3,ncomponents=3); pinsj->AddElement(elCu,fractionmass= 0.9775); pinsj->AddElement(elBe,fractionmass= 0.02); pinsj->AddElement(elNi,fractionmass= 0.0025); //pins and receipts (sockets), Gerda (MPIK, Schwingenheuer) G4Material* pinsk = new G4Material (name="pinsk",density=2.716531*g/cm3,ncomponents=3); pinsk->AddElement(elCu,fractionmass= 0.9775); pinsk->AddElement(elBe,fractionmass= 0.02); pinsk->AddElement(elNi,fractionmass= 0.0025); //pins MPIK G4Material* pinsl = new G4Material (name="pinsl", density=2.322764*g/cm3, ncomponents=3); pinsl->AddElement(elCu,fractionmass= 0.9775); pinsl->AddElement(elBe,fractionmass= 0.02); pinsl->AddElement(elNi,fractionmass= 0.0025); //pins Cu-Be, Gerda (TUM, Krause) G4Material* pinsm = new G4Material (name="pinsm",density=2.678446*g/cm3,ncomponents=3); pinsm->AddElement(elCu,fractionmass= 0.9775); pinsm->AddElement(elBe,fractionmass= 0.02); pinsm->AddElement(elNi,fractionmass= 0.0025); //receipts, Cu-Be, Gerda (TUM, Krause) G4Material* pinsn = new G4Material (name="pinsn",density=3.148043*g/cm3,ncomponents=3); pinsn->AddElement(elCu,fractionmass= 0.9775); pinsn->AddElement(elBe,fractionmass= 0.02); pinsn->AddElement(elNi,fractionmass= 0.0025); //receipts, Cu-Be, CC4, 2021, LEGEND G4Material* pinso = new G4Material (name="pinso",density=0.4388505*g/cm3,ncomponents=3); pinso->AddElement(elCu,fractionmass= 0.9775); pinso->AddElement(elBe,fractionmass= 0.02); pinso->AddElement(elNi,fractionmass= 0.0025); // Densimet D185, CeLAND G4Material* D185 = new G4Material (name="D185", density=18.5*g/cm3, ncomponents=3); D185->AddElement(elW,fractionmass= 0.97); D185->AddElement(elNi,fractionmass= 0.015); D185->AddElement(elFe,fractionmass= 0.015); // pure W powder, CeLAND G4Material* w_powder = new G4Material (name="w_powder", density=5.3289*g/cm3, ncomponents=1); w_powder->AddElement(elW,natoms= 1); // Cuflon, GERDA G4Material* cuflon = new G4Material (name="cuflon", density=2.3567*g/cm3, ncomponents=3); cuflon->AddElement(elF,fractionmass= 0.49); cuflon->AddElement(elC,fractionmass= 0.16); cuflon->AddElement(elCu,fractionmass= 0.35); // composite pcs (Cu + Teflon + PEEK etc) for levelmeter, batch 1, TPC, Xenon G4Material* levelmeter1 = new G4Material (name="levelmeter1", density=3.117512*g/cm3, ncomponents=3); levelmeter1->AddElement(elF,fractionmass= 0.49); levelmeter1->AddElement(elC,fractionmass= 0.16); levelmeter1->AddElement(elCu,fractionmass= 0.35); // composite pcs (Cu + PEEK) for levelmeter, batch 2, TPC, Xenon G4Material* levelmeter2 = new G4Material (name="levelmeter2", density=0.904075*g/cm3, ncomponents=3); levelmeter2->AddElement(elF,fractionmass= 0.49); levelmeter2->AddElement(elC,fractionmass= 0.16); levelmeter2->AddElement(elCu,fractionmass= 0.35); // organic Cu adsorbent (Cu3(BTC)2(H2O)x(EtOH)y, BTC = benzene-1,3,5-tricarboxylate (C9H6O6 = carboxylic acid) G4Material* CuBTC = new G4Material (name="cubtc", density=0.510868*g/cm3, ncomponents=4); CuBTC -> AddElement(elC, natoms =18); CuBTC -> AddElement(elH, natoms = 6); CuBTC -> AddElement(elO, natoms = 12); CuBTC -> AddElement(elCu, natoms = 3); // Glyptal 1201, red insulating coating, alkyd varnish (C19H13O7 = Alkyd) G4Material* Glyptal = new G4Material (name="Glyptal", density=0.951*g/cm3, ncomponents=3); Glyptal -> AddElement(elC, natoms =19); Glyptal -> AddElement(elH, natoms = 13); Glyptal -> AddElement(elO, natoms = 7); // Allectra, 312-KAP-MAN-014-5M, KAPTON MANGANIN WIRE G4Material* kapman = new G4Material (name="kapman", density=3.048884*g/cm3, ncomponents=7); kapman->AddElement(elH,fractionmass= 0.14*perCent); kapman->AddElement(elC,fractionmass= 3.82*perCent); kapman->AddElement(elO,fractionmass= 0.52*perCent); kapman->AddElement(elN,fractionmass= 0.22*perCent); kapman->AddElement(elCu,fractionmass= 80.24*perCent); kapman->AddElement(elMn,fractionmass= 11.53*perCent); kapman->AddElement(elNi,fractionmass= 3.53*perCent); // Allectra, 301-KAPM-035 UHV KAPTON WIRE G4Material* kapm = new G4Material (name="kapm", density=3.5876*g/cm3, ncomponents=5); kapm->AddElement(elH,fractionmass= 0.84*perCent); kapm->AddElement(elC,fractionmass= 22.45*perCent); kapm->AddElement(elO,fractionmass= 3.00*perCent); kapm->AddElement(elN,fractionmass= 1.31*perCent); kapm->AddElement(elCu,fractionmass= 72.4*perCent); // Allectra, 311-KAP-RIB10, Kapton ribbon cable G4Material* ribbon = new G4Material (name="ribbon", density=0.828391*g/cm3, ncomponents=5); ribbon->AddElement(elH,fractionmass= 2.13*perCent); ribbon->AddElement(elC,fractionmass= 56.79*perCent); ribbon->AddElement(elO,fractionmass= 7.57*perCent); ribbon->AddElement(elN,fractionmass= 3.32*perCent); ribbon->AddElement(elCu,fractionmass= 30.19*perCent); // PCB, Cu lined Kapton (35 micron Cu on both sides), empty TPC, Xenon G4Material* pcbtpcx = new G4Material (name="pcbtpcx", density=0.901558*g/cm3, ncomponents=5); pcbtpcx->AddElement(elH,fractionmass= 2.13*perCent); pcbtpcx->AddElement(elC,fractionmass= 56.79*perCent); pcbtpcx->AddElement(elO,fractionmass= 7.57*perCent); pcbtpcx->AddElement(elN,fractionmass= 3.32*perCent); pcbtpcx->AddElement(elCu,fractionmass= 30.19*perCent); // modified PCB, GERDA G4Material* pcbmod = new G4Material (name="pcbmod", density=3.147*g/cm3, ncomponents=3); pcbmod->AddElement(elF,fractionmass= 0.49); pcbmod->AddElement(elC,fractionmass= 0.16); pcbmod->AddElement(elCu,fractionmass= 0.35); // modified PCB (CC4), no components, GERDA G4Material* pcbcc4 = new G4Material (name="pcbcc4", density=3.821995*g/cm3, ncomponents=3); pcbcc4->AddElement(elF,fractionmass= 0.49); pcbcc4->AddElement(elC,fractionmass= 0.16); pcbcc4->AddElement(elCu,fractionmass= 0.35); // modified PCB (CC4), with components, GERDA G4Material* pcbcc4w = new G4Material (name="pcbcc4w", density=3.46*g/cm3, ncomponents=3); pcbcc4w->AddElement(elF,fractionmass= 0.49); pcbcc4w->AddElement(elC,fractionmass= 0.16); pcbcc4w->AddElement(elCu,fractionmass= 0.35); // modified PCB (CC4), no components, 2021, GERDA G4Material* pcbcc4_21 = new G4Material (name="pcbcc4_21", density=5.542816*g/cm3, ncomponents=3); pcbcc4_21->AddElement(elF,fractionmass= 0.49); pcbcc4_21->AddElement(elC,fractionmass= 0.16); pcbcc4_21->AddElement(elCu,fractionmass= 0.35); // modified PCB (CC4), with components, 2021, GERDA G4Material* pcbcc4w_21 = new G4Material (name="pcbcc4w_21", density=3.46*g/cm3, ncomponents=3); pcbcc4w_21->AddElement(elF,fractionmass= 0.49); pcbcc4w_21->AddElement(elC,fractionmass= 0.16); pcbcc4w_21->AddElement(elCu,fractionmass= 0.35); // modified PCB with pins and components, GERDA G4Material* pcbmodc = new G4Material (name="pcbmodc", density=3.613*g/cm3, ncomponents=3); pcbmodc->AddElement(elF,fractionmass= 0.49); pcbmodc->AddElement(elC,fractionmass= 0.16); pcbmodc->AddElement(elCu,fractionmass= 0.35); // PCB CC3, MiB, with pins and components, GERDA Phase 2 G4Material* pcbcc3 = new G4Material (name="pcbcc3", density=3.249799*g/cm3, ncomponents=3); pcbcc3->AddElement(elF,fractionmass= 0.49); pcbcc3->AddElement(elC,fractionmass= 0.16); pcbcc3->AddElement(elCu,fractionmass= 0.35); // Amphenol and Molex RF MCX connectors, DS50, GePV G4Material* RFconnds = new G4Material (name="RFconnds", density=2.376*g/cm3, ncomponents=3); RFconnds->AddElement(elF,fractionmass= 0.06); RFconnds->AddElement(elC,fractionmass= 0.02); RFconnds->AddElement(elCu,fractionmass= 0.92); // Amphenol RF MCX connectors G4Material* RFconn = new G4Material (name="RFconn", density=2.999189*g/cm3, ncomponents=3); RFconn->AddElement(elF,fractionmass= 0.06); RFconn->AddElement(elC,fractionmass= 0.02); RFconn->AddElement(elCu,fractionmass= 0.92); // Amphenol female MCX connectors for PCB board mounting G4Material* mcxf = new G4Material (name="mcxf", density=2.25982*g/cm3, ncomponents=2); mcxf->AddElement(elNi,fractionmass= 0.35); mcxf->AddElement(elCu,fractionmass= 0.65); // crimp contact, Xenon G4Material* crimp = new G4Material (name="crimp", density=2.207*g/cm3, ncomponents=2); crimp->AddElement(elNi,fractionmass= 0.35); crimp->AddElement(elCu,fractionmass= 0.65); // PCB with components, DS50 G4Material* pcbmodds50 = new G4Material (name="pcbmodds50", density=2.08*g/cm3, ncomponents=3); pcbmodds50->AddElement(elF,fractionmass= 0.49); pcbmodds50->AddElement(elC,fractionmass= 0.16); pcbmodds50->AddElement(elCu,fractionmass= 0.35); // Composition of SILVER HALIDES IN PHOTOGRAPHIC EMULSION G4Material* AgBrI = new G4Material (name="AgBrI", density=1.2274*g/cm3, ncomponents=3); AgBrI->AddElement(elBr,fractionmass= 0.422895); AgBrI->AddElement(elAg,fractionmass= 0.573748); AgBrI->AddElement(elI,fractionmass= 0.003357); // Composition of SILVER HALIDES IN PHOTOGRAPHIC EMULSION, potassium free G4Material* AgBrInoK = new G4Material (name="AgBrInoK", density=1.145206*g/cm3, ncomponents=3); AgBrInoK->AddElement(elBr,fractionmass= 0.422895); AgBrInoK->AddElement(elAg,fractionmass= 0.573748); AgBrInoK->AddElement(elI,fractionmass= 0.003357); // irradiated AgNO3, NEWS G4Material* AgNO3 = new G4Material (name="AgNO3", density=2.210306*g/cm3, ncomponents=3); AgNO3->AddElement(elN,natoms= 1); AgNO3->AddElement(elAg,natoms= 1); AgNO3->AddElement(elO,natoms= 3); // Composition of TISSUE, SOFT (ICRU FOUR-COMPONENT), Density (g/cm3) = 1.00000E+00, Mean Excitation Energy (eV) = 74.900000 G4Material* softtissue4c = new G4Material (name="softtissue4c", density=0.395364*g/cm3, ncomponents=4); softtissue4c->AddElement(elH,fractionmass= 10.1172*perCent); softtissue4c->AddElement(elC,fractionmass= 11.1000*perCent); softtissue4c->AddElement(elO,fractionmass= 76.1828*perCent); softtissue4c->AddElement(elN,fractionmass= 2.6000*perCent); // Rogers Duroid 5880 (gold plated PTFE) G4Material* duroid = new G4Material (name="duroid", density=2.79862*g/cm3, ncomponents=3); duroid->AddElement(elF,fractionmass= 0.5776); duroid->AddElement(elC,fractionmass= 0.1824); duroid->AddElement(elAu,fractionmass= 0.24); // perfluoroalkoxy alkane (PFA) G4Material* PFA = new G4Material (name="PFA", density=2.150*g/cm3, ncomponents=3); PFA->AddElement(elF,fractionmass= 0.714); PFA->AddElement(elC,fractionmass= 0.226); PFA->AddElement(elO,fractionmass= 0.060); // fluorinated ethylene propylene (FEP) G4Material* FEP = new G4Material (name="FEP", density=2.150*g/cm3, ncomponents=2); FEP->AddElement(elF,fractionmass= 0.76); FEP->AddElement(elC,fractionmass= 0.24); // fluorinated ethylene propylene (FEP), pellets, Dupont, 9494X, GERDA, TUM G4Material* FEP_dup = new G4Material (name="FEP_dup", density=1.208937*g/cm3, ncomponents=2); FEP_dup->AddElement(elF,fractionmass= 0.76); FEP_dup->AddElement(elC,fractionmass= 0.24); // fluorinated ethylene propylene (FEP), pellets, Daikin, NEOFLON NP20, GERDA, TUM G4Material* FEP_dai = new G4Material (name="FEP_dai", density=1.227443*g/cm3, ncomponents=2); FEP_dai->AddElement(elF,fractionmass= 0.76); FEP_dai->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), pellets for AXON cables, LEGEND G4Material* ptfe_pellets_axon = new G4Material (name="PTFE_pellets_axon", density=1.341304*g/cm3, ncomponents=2); ptfe_pellets_axon->AddElement(elF,fractionmass= 0.76); ptfe_pellets_axon->AddElement(elC,fractionmass= 0.24); // Ethylene propylene diene monomer (EPDM), big o-rings, MOSCA-B G4Material* EPDM_01 = new G4Material (name="EPDM_01", density=1.192*g/cm3, ncomponents=2); EPDM_01->AddElement(elH,fractionmass= 0.115291); EPDM_01->AddElement(elC,fractionmass= 0.884709); // Polymethylpentene (PMP), boxes, Nalgene, SNOlab G4Material* PMP = new G4Material (name="PMP", density=0.833*g/cm3, ncomponents=2); PMP->AddElement(elH, natoms = 12); PMP->AddElement(elC, natoms = 6); // Polyethylene naphtalate (PEN), granulate, Grade TN-8065S, Lot Nr. KS01-40127, GERDA, MPIM G4Material* PEN_01 = new G4Material (name="PEN_01", density=0.812862*g/cm3, ncomponents=3); PEN_01->AddElement(elH,fractionmass= 0.047619); PEN_01->AddElement(elC,fractionmass= 0.8); PEN_01->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), granulate, Grade TN-8050SC, Lot Nr. KS01-50622, GERDA, MPIM G4Material* PEN_02 = new G4Material (name="PEN_02", density=0.832034*g/cm3, ncomponents=3); PEN_02->AddElement(elH,fractionmass= 0.047619); PEN_02->AddElement(elC,fractionmass= 0.8); PEN_02->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), granulate, Grade 8065, GERDA-LEGEND200, MPIM G4Material* PEN_03 = new G4Material (name="PEN_03", density=0.880850*g/cm3, ncomponents=3); PEN_03->AddElement(elH,fractionmass= 0.047619); PEN_03->AddElement(elC,fractionmass= 0.8); PEN_03->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), pcs, made of PEN Grade 8065, GERDA-LEGEND200, MPIM G4Material* PEN_pcs = new G4Material (name="PEN_pcs", density=0.378501*g/cm3, ncomponents=3); PEN_pcs->AddElement(elH,fractionmass= 0.047619); PEN_pcs->AddElement(elC,fractionmass= 0.8); PEN_pcs->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), 131 sheets, GERDA-LEGEND200, MPIM/Univ. Lancaster/TU Dortmund, Univ. Michigan G4Material* PEN_04 = new G4Material (name="PEN_04", density=1.240369*g/cm3, ncomponents=3); PEN_04->AddElement(elH,fractionmass= 0.047619); PEN_04->AddElement(elC,fractionmass= 0.8); PEN_04->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), roll, DS20K G4Material* PEN_05a = new G4Material (name="PEN_05a", density=0.250263*g/cm3, ncomponents=3); PEN_05a->AddElement(elH,fractionmass= 0.047619); PEN_05a->AddElement(elC,fractionmass= 0.8); PEN_05a->AddElement(elO,fractionmass= 0.152381); G4Material* PEN_05b = new G4Material (name="PEN_05b", density=0.385564*g/cm3, ncomponents=3); PEN_05b->AddElement(elH,fractionmass= 0.047619); PEN_05b->AddElement(elC,fractionmass= 0.8); PEN_05b->AddElement(elO,fractionmass= 0.152381); G4Material* PEN_05c = new G4Material (name="PEN_05c", density=0.444964*g/cm3, ncomponents=3); PEN_05c->AddElement(elH,fractionmass= 0.047619); PEN_05c->AddElement(elC,fractionmass= 0.8); PEN_05c->AddElement(elO,fractionmass= 0.152381); // Polyethylene naphtalate (PEN), holders, GERDA-LEGEND200, MPIM/Univ. Lancaster/TU Dortmund, Univ. Michigan G4Material* PEN_06 = new G4Material (name="PEN_06", density=1.363886*g/cm3, ncomponents=3); PEN_06->AddElement(elH,fractionmass= 0.047619); PEN_06->AddElement(elC,fractionmass= 0.8); PEN_06->AddElement(elO,fractionmass= 0.152381); // Polyethyer ether ketone (PEEK), screws, Xenon TPC G4Material* PEEK_screws = new G4Material (name="PEEK_screws", density=0.433075*g/cm3, ncomponents=3); PEEK_screws->AddElement(elH,fractionmass= 0.041667); PEEK_screws->AddElement(elC,fractionmass= 0.791666); PEEK_screws->AddElement(elO,fractionmass= 0.166667); // Polyethyer ether ketone (PEEK), screws on GeDSG, new sample 2019, XenonNT TPC G4Material* PEEK_screws_nt = new G4Material (name="PEEK_screws_nt", density=0.475682*g/cm3, ncomponents=3); PEEK_screws_nt->AddElement(elH,fractionmass= 0.041667); PEEK_screws_nt->AddElement(elC,fractionmass= 0.791666); PEEK_screws_nt->AddElement(elO,fractionmass= 0.166667); // PCTFE-Polymonochlorotrifluoroethylene, ORTEC-Pd G4Material* pctfe = new G4Material (name="PCTFE", density=2.0*g/cm3, ncomponents=3); pctfe->AddElement(elC, natoms = 2); pctfe->AddElement(elCl, natoms = 1); pctfe->AddElement(elF, natoms = 3); // teflon (PTFE) G4Material* teflon = new G4Material (name="Teflon", density=2.165*g/cm3, ncomponents=2); teflon->AddElement(elF,fractionmass= 0.76); teflon->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), shrinking tube, #1, DS50 G4Material* teflon1 = new G4Material (name="Teflon_1", density=0.24348*g/cm3, ncomponents=2); teflon1->AddElement(elF,fractionmass= 0.76); teflon1->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), shrinking tube, #2, DS50 G4Material* teflon2 = new G4Material (name="Teflon_2", density=0.42383*g/cm3, ncomponents=2); teflon2->AddElement(elF,fractionmass= 0.76); teflon2->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), shrinking tube, #3, DS50 G4Material* teflon3 = new G4Material (name="Teflon_3", density=0.293425*g/cm3, ncomponents=2); teflon3->AddElement(elF,fractionmass= 0.76); teflon3->AddElement(elC,fractionmass= 0.24); // expanded teflon (ePTFE), Darkside G4Material* eteflon = new G4Material (name="eTeflon", density=0.6*g/cm3, ncomponents=2); eteflon->AddElement(elF,fractionmass= 0.76); eteflon->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), Arlon, PTFE woven laminate G4Material* teflon4 = new G4Material (name="Teflon_4", density=2.44925*g/cm3, ncomponents=2); teflon4->AddElement(elF,fractionmass= 0.76); teflon4->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), Rogers, PTFE laminate G4Material* teflon5 = new G4Material (name="Teflon_5", density=2.12914*g/cm3, ncomponents=2); teflon5->AddElement(elF,fractionmass= 0.76); teflon5->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), PTFE rings G4Material* teflon6 = new G4Material (name="Teflon_6", density=2.11376*g/cm3, ncomponents=2); teflon6->AddElement(elF,fractionmass= 0.76); teflon6->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), PTFE rings G4Material* teflon7 = new G4Material (name="Teflon_7", density=2.293565*g/cm3, ncomponents=2); teflon7->AddElement(elF,fractionmass= 0.76); teflon7->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), PTFE_Ref2, XenonNT G4Material* teflon8 = new G4Material (name="Teflon_8", density=2.206754*g/cm3, ncomponents=2); teflon8->AddElement(elF,fractionmass= 0.76); teflon8->AddElement(elC,fractionmass= 0.24); // teflon (PTFE), PTFE_HVFT, XenonNT G4Material* teflon9 = new G4Material (name="Teflon_9", density=2.25963*g/cm3, ncomponents=2); teflon9->AddElement(elF,fractionmass= 0.76); teflon9->AddElement(elC,fractionmass= 0.24); //superinsulation (RUAG COOLCAT 2 NW, aluminized PCT + polycyclohexylenedimethylene terephthalate // (polyester, PCT), in 3.123 g of superinsulation, 0.419 g aluminum (13.42 %) and 2.704 g PCT (86.58 %), Cuore G4Material* si_cuore = new G4Material (name = "si_cuore", density = 0.179413 * g/cm3, ncomponents = 4); si_cuore -> AddElement(elAl,fractionmass= 0.1342); si_cuore -> AddElement(elC, fractionmass= 0.6070); si_cuore -> AddElement(elH, fractionmass= 0.0571); si_cuore -> AddElement(elO, fractionmass= 0.2017); // teflon_mod (PTFEmod), before 0.430 // 0.430 g cm-3, GERDA PTFE cable (tube) // 2.269 g cm-3, PTFE pieces of material for coil, GERDA // 2.082 g cm-3, red PTFE wire, DS50 // 1.719595 g cm-3, ORTEC, Teflon buttons G4Material* teflon_mod = new G4Material (name="teflon_mod", density=1.719505*g/cm3, ncomponents=2); teflon_mod->AddElement(elF,fractionmass= 0.76); teflon_mod->AddElement(elC,fractionmass= 0.24); // PTFE body + Ag plated brass contacts, DS50, n-type HPGe detector, Jeff Martoff G4Material* ptfemod = new G4Material (name="ptfemod", density=3.114*g/cm3, ncomponents=5); ptfemod->AddElement(elF,fractionmass= 0.4636); ptfemod->AddElement(elC,fractionmass= 0.1464); ptfemod->AddElement(elCu,fractionmass= 0.2379); ptfemod -> AddElement(elZn, fractionmass=0.1443); ptfemod -> AddElement(elPb, fractionmass=0.0078); // operational amplifiers Analog AD8011ANZ, Xenon, Marcello Messina G4Material* opamp_a = new G4Material (name="opamp_a", density=1.778511*g/cm3, ncomponents=5); opamp_a->AddElement(elF,fractionmass= 0.4636); opamp_a->AddElement(elC,fractionmass= 0.1464); opamp_a->AddElement(elCu,fractionmass= 0.2379); opamp_a -> AddElement(elZn, fractionmass=0.1443); opamp_a -> AddElement(elPb, fractionmass=0.0078); // 400 pcs LMH6654 Single and Dual Low Power, 250 MHz, Low Noise Amplifiers,LEGEND, Stefano Riboldi G4Material* amp_smd_legend = new G4Material (name="amp_smd_legend", density=1.162264*g/cm3, ncomponents=5); amp_smd_legend->AddElement(elF,fractionmass= 0.4636); amp_smd_legend->AddElement(elC,fractionmass= 0.1464); amp_smd_legend->AddElement(elCu,fractionmass= 0.2379); amp_smd_legend -> AddElement(elZn, fractionmass=0.1443); amp_smd_legend -> AddElement(elPb, fractionmass=0.0078); // cellulose taps, vials RENOIR (cellulose) G4Material* renoirtaps = new G4Material (name="renoirtaps", density=0.18223*g/cm3, ncomponents=4); renoirtaps->AddElement(elH,fractionmass= 2.929*perCent); renoirtaps->AddElement(elN,fractionmass= 11.715*perCent); renoirtaps->AddElement(elO,fractionmass= 56.904*perCent); renoirtaps->AddElement(elC,fractionmass= 28.452*perCent); // Agar-Agar, food for flies RENOIR G4Material* agar = new G4Material (name="agar", density=0.95*g/cm3, ncomponents=3); agar->AddElement(elH,natoms = 24); agar->AddElement(elO,natoms = 9); agar->AddElement(elC,natoms = 14); // cardboard (cellulose) G4Material* cardboard = new G4Material (name="cardboard", density=0.66*g/cm3, ncomponents=4); cardboard->AddElement(elH,fractionmass= 2.929*perCent); cardboard->AddElement(elN,fractionmass= 11.715*perCent); cardboard->AddElement(elO,fractionmass= 56.904*perCent); cardboard->AddElement(elC,fractionmass= 28.452*perCent); // seaweed (cellulose) G4Material* seaweed = new G4Material (name="seaweed", density=0.96244*g/cm3, ncomponents=4); seaweed->AddElement(elH,fractionmass= 2.929*perCent); seaweed->AddElement(elN,fractionmass= 11.715*perCent); seaweed->AddElement(elO,fractionmass= 56.904*perCent); seaweed->AddElement(elC,fractionmass= 28.452*perCent); // MCE filter, cellulose nitrate + cellulose diacetate (C6H7011N3 (80%)+ C10H14O7(20%)) G4Material* mce = new G4Material (name="MCE", density=1.588*g/cm3, ncomponents=4); mce->AddElement(elH,fractionmass= 2.929*perCent); mce->AddElement(elN,fractionmass= 11.715*perCent); mce->AddElement(elO,fractionmass= 56.904*perCent); mce->AddElement(elC,fractionmass= 28.452*perCent); // food G4Material* food = new G4Material (name="food", density=0.679*g/cm3, ncomponents=2); food->AddElement(elF,fractionmass= 0.76); food->AddElement(elC,fractionmass= 0.24); // filter G4Material* filter = new G4Material (name = "filter", density = 0.273*g/cm3, ncomponents = 2); filter -> AddElement(elC, natoms = 2); filter -> AddElement(elH, natoms = 4); // gasket G4Material* gasket = new G4Material (name="gasket", density=1.851*g/cm3, ncomponents=3); gasket->AddElement(elH,fractionmass= 0.031); gasket->AddElement(elF,fractionmass= 0.594); gasket->AddElement(elC,fractionmass= 0.375); // tefzel (ETFE) G4Material* tefzel = new G4Material (name="Tefzel", density=1.70*g/cm3, ncomponents=3); tefzel->AddElement(elH,fractionmass= 0.031); tefzel->AddElement(elF,fractionmass= 0.594); tefzel->AddElement(elC,fractionmass= 0.375); // vespel G4Material* vespel = new G4Material (name="Vespel", density=1.430*g/cm3, ncomponents=2); vespel->AddElement(elH,fractionmass= 6.7*perCent); vespel->AddElement(elN,fractionmass= 93.3*perCent); // nylon // density =1.15*g/cm3; // http://www.boedeker.com G4Material* nylon = new G4Material (name="Nylon", density=1.15*g/cm3, ncomponents=4); nylon->AddElement(elH,fractionmass= 10.0*perCent); nylon->AddElement(elN,fractionmass= 10.0*perCent); nylon->AddElement(elO,fractionmass= 15.0*perCent); nylon->AddElement(elC,fractionmass= 65.0*perCent); // nylon cable ties, DS50 G4Material* nylon1 = new G4Material (name="Nylon1", density=0.76394*g/cm3, ncomponents=4); nylon1->AddElement(elH,fractionmass= 10.0*perCent); nylon1->AddElement(elN,fractionmass= 10.0*perCent); nylon1->AddElement(elO,fractionmass= 15.0*perCent); nylon1->AddElement(elC,fractionmass= 65.0*perCent); // nylon cable ties, DS50, GePV G4Material* nylon2 = new G4Material (name="Nylon2", density=0.550154*g/cm3, ncomponents=4); nylon2->AddElement(elH,fractionmass= 10.0*perCent); nylon2->AddElement(elN,fractionmass= 10.0*perCent); nylon2->AddElement(elO,fractionmass= 15.0*perCent); nylon2->AddElement(elC,fractionmass= 65.0*perCent); // Cu cable, in kapton, Gerda G4Material* cable11 = new G4Material (name="cable11", density=4.71411*g/cm3, ncomponents=5); cable11->AddElement(elH,fractionmass= 2.13*perCent); cable11->AddElement(elC,fractionmass= 56.79*perCent); cable11->AddElement(elO,fractionmass= 7.57*perCent); cable11->AddElement(elN,fractionmass= 3.32*perCent); cable11->AddElement(elCu,fractionmass= 30.19*perCent); // Cu cable, in kapton, XENON_nT G4Material* cable13 = new G4Material (name="cable13", density=0.459733*g/cm3, ncomponents=5); cable13->AddElement(elH,fractionmass= 2.13*perCent); cable13->AddElement(elC,fractionmass= 56.79*perCent); cable13->AddElement(elO,fractionmass= 7.57*perCent); cable13->AddElement(elN,fractionmass= 3.32*perCent); cable13->AddElement(elCu,fractionmass= 30.19*perCent); // cable 14a, cable a XENONnT G4Material* cable14a = new G4Material (name="cable14a", density=0.723052*g/cm3, ncomponents=3); cable14a->AddElement(elF,fractionmass= 34.8*perCent); cable14a->AddElement(elC,fractionmass= 11.0*perCent); cable14a->AddElement(elCu,fractionmass= 54.2*perCent); // cable 14, cable b XENONnT G4Material* cable14b = new G4Material (name="cable14b", density=0.933442*g/cm3, ncomponents=3); cable14b->AddElement(elF,fractionmass= 34.8*perCent); cable14b->AddElement(elC,fractionmass= 11.0*perCent); cable14b->AddElement(elCu,fractionmass= 54.2*perCent); // cable 14, cable c XENONnT G4Material* cable14c = new G4Material (name="cable14c", density=0.719440*g/cm3, ncomponents=3); cable14c->AddElement(elF,fractionmass= 34.8*perCent); cable14c->AddElement(elC,fractionmass= 11.0*perCent); cable14c->AddElement(elCu,fractionmass= 54.2*perCent); // cable 14, cable d XENONnT G4Material* cable14d = new G4Material (name="cable14d", density=2.703084*g/cm3, ncomponents=3); cable14d->AddElement(elF,fractionmass= 34.8*perCent); cable14d->AddElement(elC,fractionmass= 11.0*perCent); cable14d->AddElement(elCu,fractionmass= 54.2*perCent); // silver brazing #1, between ceramic body and Kovar flange, 670 plates, PMT Hamamatsu R11410, Xenon G4Material* agbrazing1 = new G4Material (name="agbrazing1", density=9.6828*g/cm3, ncomponents=1); agbrazing1->AddElement(elAg,fractionmass= 100*perCent); // silver brazing #2, between ceramic body and Kovar leads, PMT Hamamatsu R11410, Xenon G4Material* agbrazing2 = new G4Material (name="agbrazing2", density=3.1619*g/cm3, ncomponents=1); agbrazing2->AddElement(elAg,fractionmass= 100*perCent); // silver brazing #2, P02N12 between ceramic body and Kovar leads, PMT Hamamatsu R11410, Xenon G4Material* agbrazing22 = new G4Material (name="agbrazing22", density=2.694464*g/cm3, ncomponents=1); agbrazing22->AddElement(elAg,fractionmass= 100*perCent); // silver wire, sample 11, PMT Hamamatsu R11410, Xenon G4Material* agwire = new G4Material (name="agwire", density=2.077667*g/cm3, ncomponents=1); agwire->AddElement(elAg,fractionmass= 100*perCent); // ABS (acrilonitrile-butadiene-stirene: (C8H8 C4H6 C3H3N)n) G4Material* abs = new G4Material (name="ABS", density=1.2*g/cm3, ncomponents=3); abs->AddElement(elH,fractionmass= 8.1*perCent); abs->AddElement(elN,fractionmass= 6.6*perCent); abs->AddElement(elC,fractionmass= 85.3*perCent); // Parylene C (white powder) G4Material* parylenec = new G4Material (name="ParyleneC", density=0.590*g/cm3, ncomponents=3); parylenec->AddElement(elH,fractionmass= 5.09*perCent); parylenec->AddElement(elCl,fractionmass= 25.58*perCent); parylenec->AddElement(elC,fractionmass= 69.33*perCent); // Parylene C (evaporated and recondensed) G4Material* parylenec_cond = new G4Material (name="ParyleneC_cond", density=1.04306*g/cm3, ncomponents=3); parylenec_cond->AddElement(elH,fractionmass= 5.09*perCent); parylenec_cond->AddElement(elCl,fractionmass= 25.58*perCent); parylenec_cond->AddElement(elC,fractionmass= 69.33*perCent); //boron-ethanolamine G4Material* borEtha = new G4Material (name="BoronEthanolamine", density=1.435*g/cm3, ncomponents=5); borEtha->AddElement(elH,fractionmass= 12.4*perCent); borEtha->AddElement(elC,fractionmass= 31.3*perCent); borEtha->AddElement(elB,fractionmass= 11.8*perCent); borEtha->AddElement(elN,fractionmass= 18.3*perCent); borEtha->AddElement(elO,fractionmass= 26.1*perCent); //trimethylborate G4Material* TMB = new G4Material (name="TMB", density=0.932*g/cm3, ncomponents=4); TMB->AddElement(elH,natoms=9); TMB->AddElement(elC,natoms=3); TMB->AddElement(elB,natoms=1); TMB->AddElement(elO,natoms=3); //molybdenum oxide (Luc.), 1.3388 g/cm3 //molybdenum oxide, Pierre De Marcillac (Lucineux), Alfa Aesar, 5N, 1.310054 g/cm3 //molybdenum oxide, Pierre De Marcillac (Lucineux), Bordeaux, Industrial, 4N, 1.458282 g/cm3 G4Material* moo3 = new G4Material // (name="moo3", density=1.3388*g/cm3, ncomponents=2); // (name="moo3", density=1.310054*g/cm3, ncomponents=2); (name="moo3", density=1.458282*g/cm3, ncomponents=2); moo3->AddElement(elMo,natoms=1); moo3->AddElement(elO,natoms=3); //boron carbide G4Material* b4c = new G4Material (name="b4c", density=2.52*g/cm3, ncomponents=2); b4c->AddElement(elC,natoms=1); b4c->AddElement(elB,natoms=4); //thorium oxide G4Material* tho2 = new G4Material (name="tho2", density=1.91265*g/cm3, ncomponents=2); tho2->AddElement(elTh,natoms=1); tho2->AddElement(elO,natoms=2); //boron-polyethylen G4Material* borPoly = new G4Material (name="BoronPolyethylene", density=0.942*g/cm3, ncomponents=3); borPoly->AddElement(elH,fractionmass= 11.5*perCent); borPoly->AddElement(elC,fractionmass= 73.5*perCent); borPoly->AddElement(elB,fractionmass= 15.0*perCent); //boron-polyethylene 5%-B by mass G4Material* borPoly5 = new G4Material (name="BoronPolyethylene5", density=0.94*g/cm3, ncomponents=3); borPoly5->AddElement(elH,fractionmass= 13.5*perCent); borPoly5->AddElement(elC,fractionmass= 81.5*perCent); borPoly5->AddElement(elB,fractionmass= 5.0*perCent); //boron-polyethylene 10%-B by mass G4Material* borPoly10 = new G4Material (name="BoronPolyethylene10", density=0.94*g/cm3, ncomponents=3); borPoly10->AddElement(elH,fractionmass= 12.8*perCent); borPoly10->AddElement(elC,fractionmass= 77.2*perCent); borPoly10->AddElement(elB,fractionmass= 10.0*perCent); //lithium-polyethylene 7.5%-Li by mass G4Material* LiPoly = new G4Material (name="LithiumPolyethylene", density=1.06*g/cm3, ncomponents=3); LiPoly->AddElement(elH,fractionmass= 13.2*perCent); LiPoly->AddElement(elC,fractionmass= 79.3*perCent); LiPoly->AddElement(elLi,fractionmass= 7.5*perCent); // styropor G4Material* styropor = new G4Material (name="Styropor", density=20*kg/m3, ncomponents=2); styropor->AddElement(elH,fractionmass= 12.8*perCent); styropor->AddElement(elC,fractionmass= 87.2*perCent); //plastic scintillator G4Material* plasticScint = new G4Material (name="PlasticScintillator", density= 1.032*g/cm3, ncomponents=2); plasticScint->AddElement(elC,natoms=9); plasticScint->AddElement(elH,natoms=10); //liquid scintillator (Ultima Gold LLT & LS-221) G4Material* liquidScint = new G4Material (name="LiquidScintillator", density= 0.983*g/cm3, ncomponents=4); liquidScint->AddElement(elH,fractionmass= 9.7*perCent); liquidScint->AddElement(elC,fractionmass= 76.4*perCent); liquidScint->AddElement(elN,fractionmass= 0.1*perCent); liquidScint->AddElement(elO,fractionmass= 13.8*perCent); //plastic capacitor G4Material* plasticCap = new G4Material (name="PlasticCapacitor", density= 0.8*g/cm3, ncomponents=2); plasticCap->AddElement(elC,natoms=9); plasticCap->AddElement(elH,natoms=10); // Delrin (Polyoxymethylene), DS50, VT G4Material* delrin = new G4Material (name = "delrin", density = 1.424*g/cm3, ncomponents = 3); delrin -> AddElement(elC, fractionmass= 40.00*perCent); delrin -> AddElement(elH, fractionmass= 6.67*perCent); delrin -> AddElement(elO, fractionmass= 53.33*perCent); // polyethylene (shrinking tube #1), DS50 G4Material* shrink1 = new G4Material (name = "shrink1", density = 0.116097*g/cm3, ncomponents = 2); shrink1 -> AddElement(elC, fractionmass= 85.71*perCent); shrink1 -> AddElement(elH, fractionmass= 14.29*perCent); // polyethylene (shrinking tube #2), DS50 G4Material* shrink2 = new G4Material (name = "shrink2", density = 0.244967*g/cm3, ncomponents = 2); shrink2 -> AddElement(elC, fractionmass= 85.71*perCent); shrink2 -> AddElement(elH, fractionmass= 14.29*perCent); // polyethylene (scotch tape) G4Material* scotch = new G4Material (name = "scotch", density = 0.86*g/cm3, ncomponents = 2); scotch -> AddElement(elC, natoms = 2); scotch -> AddElement(elH, natoms = 4); // polyethylene G4Material* polyethylene = new G4Material (name = "Polyethylene", density = 0.92*g/cm3, ncomponents = 2); polyethylene -> AddElement(elC, natoms = 2); polyethylene -> AddElement(elH, natoms = 4); // polyethylene COOP CUORE, lower density G4Material* PECOOPLD = new G4Material (name = "PECOOPLD", density = 0.755*g/cm3, ncomponents = 2); PECOOPLD -> AddElement(elC, natoms = 2); PECOOPLD -> AddElement(elH, natoms = 4); // polyethylene COOP CUORE, lower density, exposure measurements G4Material* PECOOPEXP = new G4Material (name = "PECOOPEXP", density = 0.494*g/cm3, ncomponents = 2); PECOOPEXP -> AddElement(elC, natoms = 2); PECOOPEXP -> AddElement(elH, natoms = 4); // polyethylene TYVEK KINR-CUORE, apparent density G4Material* tyvekkinr = new G4Material (name = "tyvekkinr", density = 0.397*g/cm3, ncomponents = 2); tyvekkinr -> AddElement(elC, natoms = 2); tyvekkinr -> AddElement(elH, natoms = 4); // polyethylene TYVEK Darkside, apparent density G4Material* tyvekds = new G4Material (name = "tyvekds", density = 0.327345*g/cm3, ncomponents = 2); tyvekds -> AddElement(elC, natoms = 2); tyvekds -> AddElement(elH, natoms = 4); // polyethylene COOP CUORE, lowest density G4Material* PECOOPLSTD = new G4Material (name = "PECOOPLSTD", density = 0.314*g/cm3, ncomponents = 2); PECOOPLSTD -> AddElement(elC, natoms = 2); PECOOPLSTD -> AddElement(elH, natoms = 4); // polypropylene G4Material* polypropylene = new G4Material (name = "Polypropylene", density = 0.90*g/cm3, ncomponents = 2); polypropylene -> AddElement(elC, natoms = 3); polypropylene -> AddElement(elH, natoms = 6); // polyethylen terephthalat (PET, Mylar) - KINR G4Material* PETkinr = new G4Material (name = "PETkinr", density = 0.1 * g/cm3, ncomponents = 3); PETkinr -> AddElement(elC, natoms = 10); PETkinr -> AddElement(elH, natoms = 8); PETkinr -> AddElement(elO, natoms = 4); // polyetherimide pellets (Ultem) G4Material* ultem_pellets = new G4Material (name = "ultem_pellets", density = 0.870863*g/cm3, ncomponents = 4); ultem_pellets -> AddElement(elC, natoms = 37); ultem_pellets -> AddElement(elH, natoms = 24); ultem_pellets -> AddElement(elO, natoms = 6); ultem_pellets -> AddElement(elN, natoms = 2); // polyetherimide plates (Ultem) G4Material* ultem_plates = new G4Material (name = "ultem_plates", density = 1.272401*g/cm3, ncomponents = 4); ultem_plates -> AddElement(elC, natoms = 37); ultem_plates -> AddElement(elH, natoms = 24); ultem_plates -> AddElement(elO, natoms = 6); ultem_plates -> AddElement(elN, natoms = 2); // polyetherimide pellets (Siltem) G4Material* siltem_pellets = new G4Material (name = "siltem_pellets", density = 0.758314*g/cm3, ncomponents = 5); siltem_pellets -> AddElement(elC, natoms = 37); siltem_pellets -> AddElement(elH, natoms = 24); siltem_pellets -> AddElement(elO, natoms = 6); siltem_pellets -> AddElement(elN, natoms = 2); siltem_pellets -> AddElement(elSi, natoms = 4); // poly(4,4'-oxydiphenylene-pyromellitimide) (Kapton) G4Material* kapton = new G4Material (name = "kapton", density = 1530 * kg/m3, ncomponents = 4); kapton -> AddElement(elC, natoms = 20); kapton -> AddElement(elH, natoms = 9); kapton -> AddElement(elO, natoms = 2); kapton -> AddElement(elN, natoms = 1); // Pyralux (Kapton) G4Material* pyralux = new G4Material (name = "pyralux", density = 1773.717 * kg/m3, ncomponents = 4); pyralux -> AddElement(elC, natoms = 20); pyralux -> AddElement(elH, natoms = 9); pyralux -> AddElement(elO, natoms = 2); pyralux -> AddElement(elN, natoms = 1); // Pyralux (Kapton), PMT bases Sabre with no components G4Material* pyralux1 = new G4Material (name = "pyralux1", density = 1844.097 * kg/m3, ncomponents = 4); pyralux1 -> AddElement(elC, natoms = 20); pyralux1 -> AddElement(elH, natoms = 9); pyralux1 -> AddElement(elO, natoms = 2); pyralux1 -> AddElement(elN, natoms = 1); // Pyralux (Kapton) + Cu, PMT bases DS20K with no components G4Material* pyralux2 = new G4Material (name="pyralux2", density=3.245587*g/cm3, ncomponents=5); pyralux2->AddElement(elH,fractionmass= 2.13*perCent); pyralux2->AddElement(elC,fractionmass= 56.79*perCent); pyralux2->AddElement(elO,fractionmass= 7.57*perCent); pyralux2->AddElement(elN,fractionmass= 3.32*perCent); pyralux2->AddElement(elCu,fractionmass= 30.19*perCent); // Pyralux (Kapton) + Cu, PMT bases DAMA without components G4Material* pyralux3 = new G4Material (name="pyralux3", density=2.58416*g/cm3, ncomponents=5); pyralux3->AddElement(elH,fractionmass= 2.13*perCent); pyralux3->AddElement(elC,fractionmass= 56.79*perCent); pyralux3->AddElement(elO,fractionmass= 7.57*perCent); pyralux3->AddElement(elN,fractionmass= 3.32*perCent); pyralux3->AddElement(elCu,fractionmass= 30.19*perCent); // Pyralux (Kapton) + Cu, PMT bases DS20K with no components, T-shaped and regular G4Material* pyralux4 = new G4Material (name="pyralux4", density=3.086434*g/cm3, ncomponents=5); pyralux4->AddElement(elH,fractionmass= 2.13*perCent); pyralux4->AddElement(elC,fractionmass= 56.79*perCent); pyralux4->AddElement(elO,fractionmass= 7.57*perCent); pyralux4->AddElement(elN,fractionmass= 3.32*perCent); pyralux4->AddElement(elCu,fractionmass= 30.19*perCent); // polycyclohexylenedimethylene terephthalate (polyester, PCT) G4Material* PCT = new G4Material (name = "PCT", density = 1270 * kg/m3, ncomponents = 3); PCT -> AddElement(elC, natoms = 16); PCT -> AddElement(elH, natoms = 18); PCT -> AddElement(elO, natoms = 4); // polyvinyl alcohol (PVA) G4Material* PVA = new G4Material (name = "PVA", density = 1270 * kg/m3, ncomponents = 3); PVA -> AddElement(elC, natoms = 2); PVA -> AddElement(elH, natoms = 4); PVA -> AddElement(elO, natoms = 1); // Lumirror, Toray Plastics, Inc., type 188 micrometer, E6SR (reflective polyester film) G4Material* E6SR = new G4Material (name = "E6SR", density = 0.423929 * g/cm3, ncomponents = 3); E6SR -> AddElement(elC, natoms = 16); E6SR -> AddElement(elH, natoms = 18); E6SR -> AddElement(elO, natoms = 4); // white thread (polyester tape), wrapped around the nose of a marinelli 438G, DS20K G4Material* PCT_thread = new G4Material (name = "PCT_thread", density = 0.592461 * g/cm3, ncomponents = 3); PCT_thread -> AddElement(elC, natoms = 16); PCT_thread -> AddElement(elH, natoms = 18); PCT_thread -> AddElement(elO, natoms = 4); // Kapton heating bands G4Material* kapton1 = new G4Material (name = "kapton1", density = 5.656 *g/cm3, ncomponents = 4); kapton1 -> AddElement(elC, natoms = 20); kapton1 -> AddElement(elH, natoms = 9); kapton1 -> AddElement(elO, natoms = 2); kapton1 -> AddElement(elN, natoms = 1); // Pecorino cheese G4Material* cheese = new G4Material (name = "cheese", density = 0.764251*g/cm3, ncomponents = 6); cheese -> AddElement(elC, natoms = 1952); cheese -> AddElement(elH, natoms = 3556); cheese -> AddElement(elO, natoms = 258); cheese -> AddElement(elN, natoms = 64); cheese -> AddElement(elCl, natoms = 2); cheese -> AddElement(elNa, natoms = 2); // Arlon 55, Darkside-20K, nominal density 1.38 g/cm3, real density 0.93 g/cm3 G4Material* arlon55 = new G4Material (name = "Arlon55", density = 930 * kg/m3, ncomponents = 4); arlon55 -> AddElement(elC, natoms = 41); arlon55 -> AddElement(elH, natoms = 49); arlon55 -> AddElement(elO, natoms = 7); arlon55 -> AddElement(elN, natoms = 3); // Arlon PCBs, bonded, Darkside-20K, nominal density 1.38 g/cm3, apparent density 1.3325 g/cm3 G4Material* arlon_pcb_bonded = new G4Material (name = "Arlon_pcb_bonded", density = 1.3325 * g/cm3, ncomponents = 4); arlon_pcb_bonded -> AddElement(elC, natoms = 41); arlon_pcb_bonded -> AddElement(elH, natoms = 49); arlon_pcb_bonded -> AddElement(elO, natoms = 7); arlon_pcb_bonded -> AddElement(elN, natoms = 3); // silver epoxy resin (40 % epoxy, 60 % Ag) (Loctite Ablestik 56C CAT 9, 5 % hardener), Gerda TUM G4Material* epoxy_ag_01 = new G4Material (name = "Epoxy_ag_01", density = 3048 * kg/m3, ncomponents = 5); epoxy_ag_01 -> AddElement(elC, natoms = 41); epoxy_ag_01 -> AddElement(elH, natoms = 49); epoxy_ag_01 -> AddElement(elO, natoms = 7); epoxy_ag_01 -> AddElement(elN, natoms = 3); epoxy_ag_01 -> AddElement(elAg, natoms = 35); // silver epoxy resin (40 % epoxy, 60 % Ag) (Loctite Ablestik 56C CAT 9, 2.5 % hardener), Gerda TUM G4Material* epoxy_ag_02 = new G4Material (name = "Epoxy_ag_02", density = 2935 * kg/m3, ncomponents = 5); epoxy_ag_02 -> AddElement(elC, natoms = 41); epoxy_ag_02 -> AddElement(elH, natoms = 49); epoxy_ag_02 -> AddElement(elO, natoms = 7); epoxy_ag_02 -> AddElement(elN, natoms = 3); epoxy_ag_02 -> AddElement(elAg, natoms = 33); // epoxy resin Araldit 2012B, CUPID-0 G4Material* epoxy_cupid0 = new G4Material (name = "Epoxy_cupid0", density = 1.200149 * g/cm3, ncomponents = 4); epoxy_cupid0 -> AddElement(elC, natoms = 41); epoxy_cupid0 -> AddElement(elH, natoms = 49); epoxy_cupid0 -> AddElement(elO, natoms = 7); epoxy_cupid0 -> AddElement(elN, natoms = 3); // epoxy resin (generic) G4Material* epoxy = new G4Material (name = "Epoxy", density = 1129 * kg/m3, ncomponents = 4); epoxy -> AddElement(elC, natoms = 41); epoxy -> AddElement(elH, natoms = 49); epoxy -> AddElement(elO, natoms = 7); epoxy -> AddElement(elN, natoms = 3); // modified epoxy resin (generic) G4Material* epoxy_mod = new G4Material (name = "Epoxy_mod", density = 711 * kg/m3, ncomponents = 4); epoxy_mod -> AddElement(elC, natoms = 41); epoxy_mod -> AddElement(elH, natoms = 49); epoxy_mod -> AddElement(elO, natoms = 7); epoxy_mod -> AddElement(elN, natoms = 3); // modified epoxy resin (capacitors Darkside, GeMi) G4Material* epoxy_mod1 = new G4Material (name = "Epoxy_mod1", density = 1515 * kg/m3, ncomponents = 4); epoxy_mod1 -> AddElement(elC, natoms = 41); epoxy_mod1 -> AddElement(elH, natoms = 49); epoxy_mod1 -> AddElement(elO, natoms = 7); epoxy_mod1 -> AddElement(elN, natoms = 3); // modified epoxy resin (capacitors Darkside, GeMi) G4Material* epoxy_mod2 = new G4Material (name = "Epoxy_mod2", density = 1420 * kg/m3, ncomponents = 4); // (name = "Epoxy_mod2", density = 1320 * kg/m3, ncomponents = 4); epoxy_mod2 -> AddElement(elC, natoms = 41); epoxy_mod2 -> AddElement(elH, natoms = 49); epoxy_mod2 -> AddElement(elO, natoms = 7); epoxy_mod2 -> AddElement(elN, natoms = 3); // modified epoxy resin (Stycast 2850FT blue, Gerda, GeMi) G4Material* epoxy_mod3 = new G4Material (name = "Epoxy_mod3", density = 2198 * kg/m3, ncomponents = 4); epoxy_mod3 -> AddElement(elC, natoms = 41); epoxy_mod3 -> AddElement(elH, natoms = 49); epoxy_mod3 -> AddElement(elO, natoms = 7); epoxy_mod3 -> AddElement(elN, natoms = 3); // modified epoxy resin (PCB for fisheye lens, GeCris) G4Material* epoxy_mod4 = new G4Material (name = "Epoxy_mod4", density = 3.7 * g/cm3, ncomponents = 4); epoxy_mod4 -> AddElement(elC, natoms = 41); epoxy_mod4 -> AddElement(elH, natoms = 49); epoxy_mod4 -> AddElement(elO, natoms = 7); epoxy_mod4 -> AddElement(elN, natoms = 3); // modified epoxy resin (fr2 bases for PMT Darkside, GeMi;clear epoxy, Gerda, Janiscko, TUM, GeMPI4) G4Material* epoxy_mod5 = new G4Material (name = "Epoxy_mod5", density = 1.1617 * g/cm3, ncomponents = 4); epoxy_mod5 -> AddElement(elC, natoms = 41); epoxy_mod5 -> AddElement(elH, natoms = 49); epoxy_mod5 -> AddElement(elO, natoms = 7); epoxy_mod5 -> AddElement(elN, natoms = 3); // modified epoxy resin (fr2 bases for electronics lab) G4Material* epoxy_mod6 = new G4Material (name = "Epoxy_mod6", density = 1.311979 * g/cm3, ncomponents = 4); epoxy_mod6 -> AddElement(elC, natoms = 41); epoxy_mod6 -> AddElement(elH, natoms = 49); epoxy_mod6 -> AddElement(elO, natoms = 7); epoxy_mod6 -> AddElement(elN, natoms = 3); // NORPLEX Micarta, NP315 phenolic laminate with graphite, black G4Material* epoxy_mod7 = new G4Material (name = "Epoxy_mod7", density = 1.4 * g/cm3, ncomponents = 4); epoxy_mod7 -> AddElement(elC, natoms = 41); epoxy_mod7 -> AddElement(elH, natoms = 49); epoxy_mod7 -> AddElement(elO, natoms = 7); epoxy_mod7 -> AddElement(elN, natoms = 3); // Current In., C770 ESD (Electro-Static Dissipative material), G10/FR4 (glass/epoxy) G4Material* epoxy_mod8 = new G4Material (name = "Epoxy_mod8", density = 1.9 * g/cm3, ncomponents = 4); epoxy_mod8 -> AddElement(elC, natoms = 41); epoxy_mod8 -> AddElement(elH, natoms = 49); epoxy_mod8 -> AddElement(elO, natoms = 7); epoxy_mod8 -> AddElement(elN, natoms = 3); // boards LUNA-400, put in truncated pyramid-like shape, unknown composition G4Material* epoxy_luna400 = new G4Material (name = "epoxy_luna400", density = 2.252633 * g/cm3, ncomponents = 4); epoxy_luna400 -> AddElement(elC, natoms = 41); epoxy_luna400 -> AddElement(elH, natoms = 49); epoxy_luna400 -> AddElement(elO, natoms = 7); epoxy_luna400 -> AddElement(elN, natoms = 3); // resin, yellowish colour, MPP CMAT DRESIN C, Xenon G4Material* epoxy_mod9 = new G4Material (name = "Epoxy_mod9", density = 1.584668 * g/cm3, ncomponents = 4); epoxy_mod9 -> AddElement(elC, natoms = 41); epoxy_mod9 -> AddElement(elH, natoms = 49); epoxy_mod9 -> AddElement(elO, natoms = 7); epoxy_mod9 -> AddElement(elN, natoms = 3); // TeO2 powder CUORE, 2kg in big Marinelli density = 1.2775*g/cm3; G4Material* CTeO2 = new G4Material(name="CTeO2",density,ncomponents=2); CTeO2->AddElement(elTe,natoms=1); CTeO2->AddElement(elO,natoms=2); // 2,5 - diphenyl oxazole (PPO) G4Material* PPO = new G4Material (name = "PPO", density = 1421 * kg/m3, ncomponents = 4); PPO -> AddElement(elC, natoms = 15); PPO -> AddElement(elH, natoms = 11); PPO -> AddElement(elO, natoms = 1); PPO -> AddElement(elN, natoms = 1); // dimethyl benzene-1,2-dicarboxylate (Dimethyl phtalate, DMP) G4Material* dmp = new G4Material (name = "DMP", density = 1073 * kg/m3, ncomponents = 3); dmp->AddElement(elH, natoms = 10); dmp->AddElement(elC, natoms = 10); dmp->AddElement(elO, natoms = 4); G4Material* makrolon = new G4Material (name = "polycarbonate", density = 1200 * kg/m3, ncomponents = 3); makrolon->AddElement(elH,fractionmass= 5.5491*perCent); makrolon->AddElement(elC,fractionmass= 75.5751*perCent); makrolon->AddElement(elO,fractionmass= 18.8758*perCent); // Cirlex G4Material* cirlex = new G4Material (name = "cirlex", density = 1.320 * g/cm3, ncomponents = 4); cirlex->AddElement(elH , natoms=10 ); cirlex->AddElement(elC , natoms=22 ); cirlex->AddElement(elN , natoms=2 ); cirlex->AddElement(elO , natoms=5 ); // Cirlex1, boards for DS50 G4Material* cirlex1 = new G4Material (name = "cirlex1", density = 1.506248 * g/cm3, ncomponents = 4); cirlex1->AddElement(elH , natoms=10 ); cirlex1->AddElement(elC , natoms=22 ); cirlex1->AddElement(elN , natoms=2 ); cirlex1->AddElement(elO , natoms=5 ); // Cirlex2, small boards for DS50 G4Material* cirlex2 = new G4Material (name = "cirlex2", density = 2.285223 * g/cm3, ncomponents = 4); cirlex2->AddElement(elH , natoms=10 ); cirlex2->AddElement(elC , natoms=22 ); cirlex2->AddElement(elN , natoms=2 ); cirlex2->AddElement(elO , natoms=5 ); // Cirlex3, boards for LUNA-400 G4Material* cirlex3 = new G4Material (name = "cirlex3", density = 1.537715 * g/cm3, ncomponents = 4); cirlex3->AddElement(elH , natoms=10 ); cirlex3->AddElement(elC , natoms=22 ); cirlex3->AddElement(elN , natoms=2 ); cirlex3->AddElement(elO , natoms=5 ); // Polyamid-imide G4Material* torlon = new G4Material (name = "torlon", density = 1.420 * g/cm3, ncomponents = 4); torlon->AddElement(elH , natoms=4 ); torlon->AddElement(elC , natoms=8 ); torlon->AddElement(elN , natoms=2 ); torlon->AddElement(elO , natoms=3 ); // Polyamid-imide, ORTEC-Pd G4Material* torlon_ortec = new G4Material (name = "torlon_ortec", density = 1.63508 * g/cm3, ncomponents = 4); torlon_ortec->AddElement(elH , natoms=4 ); torlon_ortec->AddElement(elC , natoms=8 ); torlon_ortec->AddElement(elN , natoms=2 ); torlon_ortec->AddElement(elO , natoms=3 ); // ZIF connectors, large (GePV), for GERDA G4Material* zif1 = new G4Material (name = "zif1", density = 1.5267 * g/cm3, ncomponents = 4); zif1->AddElement(elH,fractionmass= 10.0*perCent); zif1->AddElement(elN,fractionmass= 10.0*perCent); zif1->AddElement(elO,fractionmass= 15.0*perCent); zif1->AddElement(elC,fractionmass= 65.0*perCent); // ZIF connectors, small, horizontal (GePaolo), for GERDA G4Material* zif2 = new G4Material (name = "zif2", density = 1.64096 * g/cm3, ncomponents = 4); zif2->AddElement(elH,fractionmass= 10.0*perCent); zif2->AddElement(elN,fractionmass= 10.0*perCent); zif2->AddElement(elO,fractionmass= 15.0*perCent); zif2->AddElement(elC,fractionmass= 65.0*perCent); // ZIF connectors, small, vertical (GeDSG), for GERDA G4Material* zif3 = new G4Material (name = "zif3", density = 1.38069 * g/cm3, ncomponents = 4); zif3->AddElement(elH,fractionmass= 10.0*perCent); zif3->AddElement(elN,fractionmass= 10.0*perCent); zif3->AddElement(elO,fractionmass= 15.0*perCent); zif3->AddElement(elC,fractionmass= 65.0*perCent); // Korrodin CU5 bolts (CU5 former CuNi60, Ni (1.0 - 1.6 %), Si (0.4 - 0.7 %), rest Cu; put apparent density) G4Material* CU5 = new G4Material (name="CU5",density=3672*kg/m3,ncomponents=3); CU5->AddElement(elSi,fractionmass= 0.55*perCent); CU5->AddElement(elNi,fractionmass= 1.3*perCent); CU5->AddElement(elCu,fractionmass= 98.15*perCent); // Sapphire disks for GeDSG G4Material* sapphire = new G4Material (name="Sapphire",density=3972*kg/m3,ncomponents=2); sapphire->AddElement(elAl, natoms=2 ); sapphire->AddElement(elO, natoms=3 ); // Sapphire "Verneuil-Birne", MPIK-HD G4Material* verneuil = new G4Material (name="verneuil",density=3773*kg/m3,ncomponents=2); verneuil->AddElement(elAl, natoms=2 ); verneuil->AddElement(elO, natoms=3 ); // Sapphire Canberra, MPIK-HD G4Material* canberra = new G4Material (name="canberra",density=4045.7*kg/m3,ncomponents=2); canberra->AddElement(elAl, natoms=2 ); canberra->AddElement(elO, natoms=3 ); // Aluminium nitride, INFN-BO, DS20t G4Material* aln = new G4Material (name="aln",density=3347*kg/m3,ncomponents=2); aln->AddElement(elAl, natoms=2 ); aln->AddElement(elN, natoms=3 ); // n25 MICROCHIP TECHNOLOGY, HV3418PG-G-ND, IC 64BIT SRL PARALLEL 180V 80QFP, DS20K, SOM 11-OCT-2019 G4Material* IC_01 = new G4Material (name = "ic_01", density = 1.621981 * g/cm3, ncomponents = 4); IC_01->AddElement(elH,fractionmass= 10.0*perCent); IC_01->AddElement(elN,fractionmass= 10.0*perCent); IC_01->AddElement(elO,fractionmass= 15.0*perCent); IC_01->AddElement(elC,fractionmass= 65.0*perCent); // 11 Electrolytic capacitors, DS50 for GeDSG G4Material* electrolytic = new G4Material (name="electrolytic",density=830*kg/m3,ncomponents=2); electrolytic->AddElement(elAl, natoms=2 ); electrolytic->AddElement(elO, natoms=3 ); // 12 Electrolytic capacitors, DS50 for GeDSG G4Material* electrolytic1 = new G4Material (name="electrolytic1",density=4149*kg/m3,ncomponents=2); electrolytic1->AddElement(elAl, natoms=2 ); electrolytic1->AddElement(elO, natoms=3 ); // Susumu thick film resistors #4, DS50, GeDSG G4Material* electrolytic2a = new G4Material (name="electrolytic2a",density=2198.5*kg/m3,ncomponents=2); electrolytic2a->AddElement(elAl, natoms=2 ); electrolytic2a->AddElement(elO, natoms=3 ); // Susumu thick film resistors #5, DS50, GeDSG G4Material* electrolytic2b = new G4Material (name="electrolytic2b",density=1901.4*kg/m3,ncomponents=2); electrolytic2b->AddElement(elAl, natoms=2 ); electrolytic2b->AddElement(elO, natoms=3 ); // electronics pieces for four PMT bases, DS50, GeDSG G4Material* electrolytic2c = new G4Material (name="electrolytic2c",density=2.533*g/cm3,ncomponents=2); electrolytic2c->AddElement(elAl, natoms=2 ); electrolytic2c->AddElement(elO, natoms=3 ); // 120 resistors , Xenon, GeDSG G4Material* electrolytic2d = new G4Material (name="electrolytic2d",density=0.512821*g/cm3,ncomponents=2); electrolytic2d->AddElement(elAl, natoms=2 ); electrolytic2d->AddElement(elO, natoms=3 ); // 150 ceramic thick film resistors, lead bearing, Vishay, Mfg CRCW0603470MJPTAHR, GERDA, MIB (Cattadori), GeDSG G4Material* electrolytic2e = new G4Material (name="electrolytic2e",density=1.50501*g/cm3,ncomponents=2); electrolytic2e->AddElement(elAl, natoms=2 ); electrolytic2e->AddElement(elO, natoms=3 ); // 50 resistors, Ohmite, model HVC0402E5006KET, format 0402, GERDA, MIB (Cattadori), GeDSG G4Material* electrolytic2f = new G4Material (name="electrolytic2f",density=0.278136*g/cm3,ncomponents=2); electrolytic2f->AddElement(elAl, natoms=2 ); electrolytic2f->AddElement(elO, natoms=3 ); // 4973 resistors, Vishay Draloric, 51 Ohm, model CRCW80551ROFKEA, thick film, Xenon, GsOr G4Material* electrolytic2g = new G4Material (name="electrolytic2g",density=2.011515*g/cm3,ncomponents=2); electrolytic2g->AddElement(elAl, natoms=2 ); electrolytic2g->AddElement(elO, natoms=3 ); // 4784 resistors, Vishay Draloric, 10 MOhm, model CRCW80510MOFKEA, thick film, Xenon, GsOr G4Material* electrolytic2h = new G4Material (name="electrolytic2h",density=1.248366*g/cm3,ncomponents=2); electrolytic2h->AddElement(elAl, natoms=2 ); electrolytic2h->AddElement(elO, natoms=3 ); // 236 resistors, Kapton, Xenon, GeMi G4Material* electrolytic2i = new G4Material (name="electrolytic2i",density=1.533058*g/cm3,ncomponents=2); electrolytic2i->AddElement(elAl, natoms=2 ); electrolytic2i->AddElement(elO, natoms=3 ); // 4977 resistors, Stackpole Electronics, Inc., 15 MOhm, model RMCF0805JT15M0, thick film, Xenon, GeCris G4Material* electrolytic2j = new G4Material (name="electrolytic2j",density=1.394325*g/cm3,ncomponents=2); electrolytic2j->AddElement(elAl, natoms=2 ); electrolytic2j->AddElement(elO, natoms=3 ); // 4963 resistors, Vishay Draloric, 1 kOhm, model CRCW8051k0OFKT, thick film, Xenon, GeDSG G4Material* electrolytic2k = new G4Material (name="electrolytic2k",density=1.677413*g/cm3,ncomponents=2); electrolytic2k->AddElement(elAl, natoms=2 ); electrolytic2k->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #1, 1.86 g in std vial type, GeDSG G4Material* electrolytic2l = new G4Material (name="electrolytic2l",density=2.436446*g/cm3,ncomponents=2); electrolytic2l->AddElement(elAl, natoms=2 ); electrolytic2l->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #2, 6.32 g in new std vial type 2, GeDSG G4Material* electrolytic2m = new G4Material // (name="electrolytic2m",density=2.146214*g/cm3,ncomponents=2); (name="electrolytic2m",density=2.156402*g/cm3,ncomponents=2); electrolytic2m->AddElement(elAl, natoms=2 ); electrolytic2m->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #3, 6.38 g in new std vial type 2, GeDSG G4Material* electrolytic2n = new G4Material (name="electrolytic2n",density=2.166589*g/cm3,ncomponents=2); electrolytic2n->AddElement(elAl, natoms=2 ); electrolytic2n->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #4, 2.05 g in std vial type, GeDSG G4Material* electrolytic2o = new G4Material (name="electrolytic2o",density=2.500134*g/cm3,ncomponents=2); electrolytic2o->AddElement(elAl, natoms=2 ); electrolytic2o->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #5, 6.32 g in std vial type, GeDSG G4Material* electrolytic2p = new G4Material (name="electrolytic2p",density=1.677413*g/cm3,ncomponents=2); electrolytic2p->AddElement(elAl, natoms=2 ); electrolytic2p->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #6, 6.32 g in std vial type, GeDSG G4Material* electrolytic2q = new G4Material (name="electrolytic2q",density=1.677413*g/cm3,ncomponents=2); electrolytic2q->AddElement(elAl, natoms=2 ); electrolytic2q->AddElement(elO, natoms=3 ); // SMD components for PMT bases, Xenon, sample #7, 4.87 g in new std vial type 2, GeDSG G4Material* electrolytic2r = new G4Material (name="electrolytic2r",density=1.747329*g/cm3,ncomponents=2); electrolytic2r->AddElement(elAl, natoms=2 ); electrolytic2r->AddElement(elO, natoms=3 ); // MLCC capacitors Accu-p, DS20k, , GeDSG G4Material* electrolytic2s = new G4Material (name="electrolytic2s",density=2.628422*g/cm3,ncomponents=2); electrolytic2s->AddElement(elAl, natoms=2 ); electrolytic2s->AddElement(elO, natoms=3 ); // XENON resistors, 310 pcs, OHMITE (1P) P/N: HVF2512T500FE, 5 GOHM, FLIP CHIP RES 1% (4L) COO:CH G4Material* electrolytic2t = new G4Material (name="electrolytic2t",density=2.038763*g/cm3,ncomponents=2); electrolytic2t->AddElement(elAl, natoms=2 ); electrolytic2t->AddElement(elO, natoms=3 ); // n. 200 ANALOG DEVICE, used n. 180 LTC6820HMS#PBF, IC INTERFACE SPECIALIZED 16MSOP,DS-20k, SOM 08-OCT-2019 G4Material* electrolytic2u = new G4Material (name="electrolytic2u",density=1.134318*g/cm3,ncomponents=2); electrolytic2u->AddElement(elAl, natoms=2 ); electrolytic2u->AddElement(elO, natoms=3 ); // components on CC4, LEGEND, 2021 G4Material* electrolytic2v = new G4Material (name="electrolytic2v",density=0.777778*g/cm3,ncomponents=2); electrolytic2v->AddElement(elAl, natoms=2 ); electrolytic2v->AddElement(elO, natoms=3 ); // n.525, Vishay/Beyschlag, MCT06030C1005FP500, thin film, SMD, 0.1W, 10 Mohm, 1%, 0603, 50ppm, Auto, DS-20K, SOM 24-MAR-2022 G4Material* electrolytic2w = new G4Material (name="electrolytic2w",density=1.744827*g/cm3,ncomponents=2); electrolytic2w->AddElement(elAl, natoms=2 ); electrolytic2w->AddElement(elO, natoms=3 ); // 40 Tantalum capacitors, DS50 for GeDSG G4Material* electrolytic3 = new G4Material (name="electrolytic3",density=1.401*g/cm3,ncomponents=2); electrolytic3->AddElement(elTa , natoms=2 ); electrolytic3->AddElement(elO , natoms=5 ); // 430 Tantalum capacitors, F981C475MMA, LEGEND for GeCris G4Material* electrolytic3a = new G4Material (name="electrolytic3a",density=1.684594*g/cm3,ncomponents=2); electrolytic3a->AddElement(elTa , natoms=2 ); electrolytic3a->AddElement(elO , natoms=5 ); // n. 1, 500 pcs, 47 Ohm, 0.3 g, (19.5 x 7.15 x 2.9) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4a = new G4Material (name="electrolytic4a",density=0.741964*g/cm3,ncomponents=2); electrolytic4a->AddElement(elAl, natoms=2 ); electrolytic4a->AddElement(elO, natoms=3 ); // n. 2, 250 pcs, 100 Ohm, 0.2 g, (9 x 9 x 1.8) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4b = new G4Material (name="electrolytic4b",density=1.371742*g/cm3,ncomponents=2); electrolytic4b->AddElement(elAl, natoms=2 ); electrolytic4b->AddElement(elO, natoms=3 ); // n. 3, 500 pcs, 22 Ohm, 0.3 g, (15.9 x 12 x 1.26) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4c = new G4Material (name="electrolytic4c",density=1.247876*g/cm3,ncomponents=2); electrolytic4c->AddElement(elAl, natoms=2 ); electrolytic4c->AddElement(elO, natoms=3 ); // n. 4, 500 pcs, 820 Ohm, 0.3 g, (11.5 x 13.3 x 2.4) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4d = new G4Material (name="electrolytic4d",density=0.817261*g/cm3,ncomponents=2); electrolytic4d->AddElement(elAl, natoms=2 ); electrolytic4d->AddElement(elO, natoms=3 ); // n. 5, 300 pcs, 1.5 kOhm, 0.2 g, (14 x 12 x 2.3) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4e = new G4Material (name="electrolytic4e",density=0.517598*g/cm3,ncomponents=2); electrolytic4e->AddElement(elAl, natoms=2 ); electrolytic4e->AddElement(elO, natoms=3 ); // n. 6, 250 pcs, 1000 Ohm, 0.2 g, (12 x 13 x 1.7) mm, LEGEND, UniMi (Riboldi), GeMPI2 G4Material* electrolytic4f = new G4Material (name="electrolytic4f",density=0.754148*g/cm3,ncomponents=2); electrolytic4f->AddElement(elAl, natoms=2 ); electrolytic4f->AddElement(elO, natoms=3 ); // MLG series, 0402, SMD inductor, 220 nHy, 100 mA, Qty: 4400, stock No: 741-3797P (PN: MLG1005SR22JT), on reel, DS20K, GeMPI2 G4Material* electrolytic4g = new G4Material (name="electrolytic4g",density=0.774483*g/cm3,ncomponents=2); electrolytic4g->AddElement(elAl, natoms=2 ); electrolytic4g->AddElement(elO, natoms=3 ); // hybrid HRC0301, 2 kMEG Ohm, 0.5 pF capacitor, ORTEC-Pd G4Material* electrolytic5a = new G4Material (name="electrolytic5a",density=3.464567*g/cm3,ncomponents=2); electrolytic5a->AddElement(elAl, natoms=2 ); electrolytic5a->AddElement(elO, natoms=3 ); // CdZnTe crystals COBRA G4Material* czt = new G4Material (name="czt",density=5642*kg/m3,ncomponents=3); czt->AddElement(elCd , natoms=9 ); czt->AddElement(elZn , natoms=1 ); czt->AddElement(elTe , natoms=0 ); // YVO4 crystal STS G4Material* yvo = new G4Material (name="yvo",density=4.283429*g/cm3,ncomponents=3); yvo->AddElement(elY , natoms=1 ); yvo->AddElement(elV , natoms=1 ); yvo->AddElement(elO , natoms=4 ); // YCl3 powder for CLYC, Milano Bicocca, Carla Cattadori G4Material* ycl3 = new G4Material (name="ycl3",density=1.708218*g/cm3,ncomponents=2); ycl3->AddElement(elY , natoms=1 ); ycl3->AddElement(elCl , natoms=3 ); // CLYC (Cs2LiYCl6:Ce) crystal, Milano Bicocca, Carla Cattadori, nominal density 3.31 g/cm3 G4Material* CLYCcry = new G4Material (name="CLYCcry",density=1.728737*g/cm3,ncomponents=5); CLYCcry->AddElement(elCs , natoms=2000 ); CLYCcry->AddElement(elLi , natoms=1000 ); CLYCcry->AddElement(elY , natoms=1000 ); CLYCcry->AddElement(elCl , natoms=6000 ); CLYCcry->AddElement(elCe , natoms=3 ); // CLYC (Cs2LiYCl6:Ce) powder, Milano Bicocca, Carla Cattadori G4Material* CLYCpow = new G4Material (name="CLYCpow",density=1.976787*g/cm3,ncomponents=5); CLYCpow->AddElement(elCs , natoms=2000 ); CLYCpow->AddElement(elLi , natoms=1000 ); CLYCpow->AddElement(elY , natoms=1000 ); CLYCpow->AddElement(elCl , natoms=6000 ); CLYCpow->AddElement(elCe , natoms=3 ); // Cesium Hafnium Chloride CLYC (Cs2HfCl2) crystal, Prof. Burger, Fisk University, nominal density 3.67 g/cm3 G4Material* CHCcry = new G4Material (name="CHCcry",density=3.456169*g/cm3,ncomponents=3); CHCcry->AddElement(elCs , natoms=2000 ); CHCcry->AddElement(elHf , natoms=1000 ); CHCcry->AddElement(elCl , natoms=2000 ); // Cesium Hafnium Chloride CLYC (Cs2HfCl2) crystal, Prof. Burger, Fisk University, DAMA, Nagorny, nominal density 3.67 g/cm3 G4Material* CHCcry_02 = new G4Material (name="CHCcry_02",density=3.888802*g/cm3,ncomponents=3); CHCcry_02->AddElement(elCs , natoms=2000 ); CHCcry_02->AddElement(elHf , natoms=1000 ); CHCcry_02->AddElement(elCl , natoms=2000 ); // Cesium Chloride, CsCl, raw material for growth of Cs2HfCl2) crystal, Prof. Burger, Fisk University, nominal density 3.67 g/cm3 G4Material* cscl = new G4Material (name="cscl",density=3.323238*g/cm3,ncomponents=2); cscl->AddElement(elCs , natoms=1 ); cscl->AddElement(elCl , natoms=1 ); // Hafnium Chloride HfCl3, raw material for growth of Cs2HfCl2) crystal, Prof. Burger, Fisk University, nominal density 3.67 g/cm3 G4Material* hfcl3 = new G4Material (name="hfcl3",density=2.021113*g/cm3,ncomponents=2); hfcl3->AddElement(elCl , natoms=3 ); hfcl3->AddElement(elHf , natoms=1 ); // Yttrium oxide Y2O3, purified, in two new standard vials type 2, GeDSG G4Material* y2o3 = new G4Material (name="Yttriumoxide",density=0.945875*g/cm3,ncomponents=2); y2o3->AddElement(elY , natoms=2 ); y2o3->AddElement(elO , natoms=3 ); // Yttrium oxide Y2O3, purified, in flat box D 65 mm, BEGe G4Material* y2o3_box = new G4Material (name="Yttriumoxide_box",density=1.36921*g/cm3,ncomponents=2); y2o3_box->AddElement(elY , natoms=2 ); y2o3_box->AddElement(elO , natoms=3 ); // Cesium Zirconium Chloride CLYC (Cs2ZrCl6) crystal, cylindric, BEGe, Serge Nagorny, DAMA G4Material* CZCcyl = new G4Material (name="CZCcyl",density=3.26849*g/cm3,ncomponents=3); CZCcyl->AddElement(elCs , natoms=2 ); CZCcyl->AddElement(elZr , natoms=1 ); CZCcyl->AddElement(elCl , natoms=6 ); // Cesium Zirconium Chloride CLYC (Cs2ZrCl6) crystal, truncated cone, BEGe, Serge Nagorny, DAMA G4Material* CZCcone = new G4Material (name="CZCcone",density=3.710017*g/cm3,ncomponents=3); CZCcone->AddElement(elCs , natoms=2 ); CZCcone->AddElement(elZr , natoms=1 ); CZCcone->AddElement(elCl , natoms=6 ); //metallic beryllium G4Material* metalBe = new G4Material (name="MetalBeryllium",density=1.848*g/cm3,ncomponents=1); metalBe->AddElement(elBe,natoms=1); //graphene, 191B, CMG Italia, DS50 G4Material* graphene = new G4Material (name="graphene",density=0.00903762*g/cm3,ncomponents=1); graphene->AddElement(elC,natoms=1); //synthetic charcoal, Temple University, DS50 G4Material* synchar = new G4Material (name="synchar",density=0.13824*g/cm3,ncomponents=1); synchar->AddElement(elC,natoms=1); //charcoal Bluecher 102688, MPI-K-HD, Xenon G4Material* Bluecher102688 = new G4Material (name="Bluecher102688",density=0.6011*g/cm3,ncomponents=1); Bluecher102688->AddElement(elC,natoms=1); //carbon fibre //charcoal G4Material* carbonf = new G4Material (name="CarbonFibre",density=0.44*g/cm3,ncomponents=1); carbonf->AddElement(elC,natoms=1); //carbon fibre, ORTEC, DS50 //charcoal G4Material* carbonf1 = new G4Material (name="CarbonFibre1",density=0.52838*g/cm3,ncomponents=1); carbonf1->AddElement(elC,natoms=1); //carbon felt, Luc. //charcoal G4Material* carbonfelt = new G4Material (name="CarbonFelt",density=0.119897*g/cm3,ncomponents=1); carbonfelt->AddElement(elC,natoms=1); // paper 1, cellulose, MiB, Previtali G4Material* paper1 = new G4Material (name="paper1", density=0.509386*g/cm3, ncomponents=4); paper1->AddElement(elH,fractionmass= 2.929*perCent); paper1->AddElement(elN,fractionmass= 11.715*perCent); paper1->AddElement(elO,fractionmass= 56.904*perCent); paper1->AddElement(elC,fractionmass= 28.452*perCent); // paper 2, cellulose, MiB, Previtali G4Material* paper2 = new G4Material (name="paper2", density=0.617533*g/cm3, ncomponents=4); paper2->AddElement(elH,fractionmass= 2.929*perCent); paper2->AddElement(elN,fractionmass= 11.715*perCent); paper2->AddElement(elO,fractionmass= 56.904*perCent); paper2->AddElement(elC,fractionmass= 28.452*perCent); // ICRU sphere, 30 cm diameter G4Material* icru = new G4Material (name="icru", density=1.0*g/cm3, ncomponents=4); icru->AddElement(elH,fractionmass= 10.1*perCent); icru->AddElement(elN,fractionmass= 2.6*perCent); icru->AddElement(elO,fractionmass= 76.2*perCent); icru->AddElement(elC,fractionmass= 11.1*perCent); // lithium carbonate Alfa Aesar, Luc. (Li2CO3) G4Material* Li2CO3 = new G4Material (name="Li2CO3",density=0.68349*g/cm3,ncomponents=3); Li2CO3->AddElement(elLi,2); Li2CO3->AddElement(elC,1); Li2CO3->AddElement(elO,3); // lithium carbonate Sigma-Aldrich, Luc. (Li2CO3) G4Material* Li2CO3SA = new G4Material (name="Li2CO3SA",density=0.781187*g/cm3,ncomponents=3); Li2CO3SA->AddElement(elLi,2); Li2CO3SA->AddElement(elC,1); Li2CO3SA->AddElement(elO,3); // lithium carbonate, Bordeaux, Lucineu (Li2CO3) G4Material* Li2CO3bdx = new G4Material (name="Li2CO3bdx",density=0.87455*g/cm3,ncomponents=3); Li2CO3bdx->AddElement(elLi,2); Li2CO3bdx->AddElement(elC,1); Li2CO3bdx->AddElement(elO,3); // lithium carbonate, Novosibirsk, Lucineu (Li2CO3) G4Material* Li2CO3nsbk = new G4Material (name="Li2CO3nsbk",density=0.715281*g/cm3,ncomponents=3); Li2CO3nsbk->AddElement(elLi,2); Li2CO3nsbk->AddElement(elC,1); Li2CO3nsbk->AddElement(elO,3); // lithium carbonate, Stanford Laboratories, 99.999%, Luc. (Li2CO3) G4Material* Li2CO3sfl = new G4Material (name="Li2CO3sfl",density=1.35985*g/cm3,ncomponents=3); Li2CO3sfl->AddElement(elLi,2); Li2CO3sfl->AddElement(elC,1); Li2CO3sfl->AddElement(elO,3); // lithium carbonate, Cina, CUPID (Li2CO3) G4Material* Li2CO3cncupid = new G4Material (name="Li2CO3cncupid",density=0.741512*g/cm3,ncomponents=3); Li2CO3cncupid->AddElement(elLi,2); Li2CO3cncupid->AddElement(elC,1); Li2CO3cncupid->AddElement(elO,3); G4Isotope* Li6 = new G4Isotope(name="Li6",iz= 3,in=6,a=6.*g/mole); G4Isotope* Li7 = new G4Isotope(name="Li7",iz= 3,in=7,a=7.*g/mole); G4Element* enr7Li = new G4Element (name="enr7Li",symbol="eLi7",ncomponents=2); enr7Li->AddIsotope(Li6, abundance=0.0006); enr7Li->AddIsotope(Li7, abundance=0.9994); G4Element* enr6Li = new G4Element (name="enr6Li",symbol="eLi6",ncomponents=2); enr6Li->AddIsotope(Li6, abundance=0.95); enr6Li->AddIsotope(Li7, abundance=0.05); // lithium carbonate, enriched in 7Li, Lot N1701/DAMRI, (Li2CO3) G4Material* enr7Li2CO3 = new G4Material (name="enr7Li2CO3",density=0.647697*g/cm3,ncomponents=3); enr7Li2CO3->AddElement(enr7Li,2); enr7Li2CO3->AddElement(elC,1); enr7Li2CO3->AddElement(elO,3); // lithium carbonate, enriched in 6Li, Lot I-18129, (Li2CO3) G4Material* enr6Li2CO3 = new G4Material (name="enr6Li2CO3",density=0.647697*g/cm3,ncomponents=3); enr6Li2CO3->AddElement(enr6Li,2); enr6Li2CO3->AddElement(elC,1); enr6Li2CO3->AddElement(elO,3); // lithium molybdate, Luc. (Li2MoO4) G4Material* Li2MoO4 = new G4Material (name="Li2MoO4",density=1.373935*g/cm3,ncomponents=3); Li2MoO4->AddElement(elLi,2); Li2MoO4->AddElement(elMo,1); Li2MoO4->AddElement(elO,4); //metallic lithium G4Material* metalLi = new G4Material (name="MetalLithium",density=0.535*g/cm3,ncomponents=1); metalLi->AddElement(elLi,natoms=1); //metallic molybdenum G4Material* metalMo = new G4Material (name="Mo",density=10.2*g/cm3,ncomponents=1); metalMo->AddElement(elMo,natoms=1); //quartz (used in 3pmt DS50) G4Material* quartz = new G4Material (name="Quartz",density=2.2*g/cm3,ncomponents=2); quartz->AddElement(elSi,1); quartz->AddElement(elO,2); //quartz, window for VT camera system, DS50 G4Material* quartzDS = new G4Material (name="QuartzDS",density=3.9612*g/cm3,ncomponents=2); quartzDS->AddElement(elSi,1); quartzDS->AddElement(elO,2); //quartz, broken quartz tube for ZnSe synthesis, Luc. G4Material* quartzLuc = new G4Material (name="QuartzLuc",density=0.9853*g/cm3,ncomponents=2); quartzLuc->AddElement(elSi,1); quartzLuc->AddElement(elO,2); //ceramic G4Material* ceramic = new G4Material //(name="ceramic",density=3.8683*g/cm3,ncomponents=2); //(name="ceramic",density=3.9304*g/cm3,ncomponents=2); // ceramic disks (AH03b (2x) and AH04a (5x)), DS50 //(name="ceramic",density=3.35*g/cm3,ncomponents=2); // 8 types of MLCC (samples #1 & #6-12), DS50 //(name="ceramic",density=3.23*g/cm3,ncomponents=2); // ceramic supports for PMT R11065 Hamamatsu, red, 19 pcs //(name="ceramic",density=3.8393*g/cm3,ncomponents=2); // ceramic Martoff, DS50 //(name="ceramic",density=3.8683*g/cm3,ncomponents=2); // ceramic for n-type HPGe, ORTEC, Martoff, DS50 //(name="ceramic",density=3.917*g/cm3,ncomponents=2); // ceramic disks for stem body of PMTs Hamamatsu 11410, Xenon //(name="ceramic",density=3.984*g/cm3,ncomponents=2); // ceramic cylinders for stem body of PMTs Hamamatsu 11410, Xenon //(name="ceramic",density=3.685636*g/cm3,ncomponents=2); // ceramic cylinders for HPGe detector, Conus //(name="ceramic",density=3.327962*g/cm3,ncomponents=2); // ceramic PCB for Si detectors Lakeshore DT-471-CU, VIP, GePaolo (name="ceramic",density=3.327962*g/cm3,ncomponents=2); ceramic->AddElement(elAl,2); ceramic->AddElement(elO,3); //new stem for Hamamatsu PMT 11410, Xenon (Kovar 8.16 g/cm3, ceramic 3.917 g/cm3) G4Material* stem = new G4Material (name="stem",density=7.12578*g/cm3,ncomponents=6); stem->AddElement(elAl,fractionmass= 0.04); stem->AddElement(elO,fractionmass= 0.03); stem->AddElement(elSi,fractionmass= 0.01); stem->AddElement(elCo,fractionmass= 0.16); stem->AddElement(elFe,fractionmass= 0.49); stem->AddElement(elNi,fractionmass= 0.27); //ceramic heater Luc. G4Material* ceramicLuc = new G4Material (name="ceramicLuc",density=2.1755*g/cm3,ncomponents=2); ceramicLuc->AddElement(elAl,2); ceramicLuc->AddElement(elO,3); //ceramic generic for CCD cam Cygnus G4Material* ceramicCCD = new G4Material (name="ceramicCCD",density=1.777855*g/cm3,ncomponents=2); ceramicCCD->AddElement(elAl,2); ceramicCCD->AddElement(elO,3); //resistors 1 GOhm, Xenon G4Material* res1go = new G4Material (name="res1go",density=2.127*g/cm3,ncomponents=2); res1go->AddElement(elAl,2); res1go->AddElement(elO,3); // titanium G4Material* metalTi = new G4Material(name = "MetalTitanium", density = 4.51*g/cm3, ncomponents = 1); metalTi -> AddElement(elTi, natoms = 1); // nickel G4Material* metalNi = new G4Material(name = "MetalNickel", density = 8.800*g/cm3, ncomponents = 1); metalNi -> AddElement(elNi, natoms = 1); // gold G4Material* metalAu = new G4Material(name = "MetalGold", density = 19.320*g/cm3, ncomponents = 1); metalAu -> AddElement(elAu, natoms = 1); // indium G4Material* metalIn = new G4Material(name = "MetalIndium", density = 7.31*g/cm3, ncomponents = 1); metalIn -> AddElement(elIn, natoms = 1); // indium wire for DS50 on spool G4Material* indiumspool = new G4Material(name = "indiumspool", density = 2.546*g/cm3, ncomponents = 1); indiumspool -> AddElement(elIn, natoms = 1); // silver G4Material* metalAg = new G4Material(name = "MetalSilver", density = 10.49*g/cm3, ncomponents = 1); metalAg -> AddElement(elAg, natoms = 1); // titanium sponge, DarkSide, JSC Solikamsk Magnesium Works G4Material* Tisponge = new G4Material(name = "Tisponge", density = 1.149206*g/cm3, ncomponents = 1); Tisponge -> AddElement(elTi, natoms = 1); // PMT photocathode aluminium G4Material* CathodeMetalAluminium = new G4Material (name="CathodeMetalAluminium", density=2.7*g/cm3, ncomponents=1); CathodeMetalAluminium->AddElement(elAl, 1); // Brass density = 8.5*g/cm3; G4Material* Brass = new G4Material(name="Brass", density, ncomponents=3); Brass -> AddElement(elCu, fractionmass=0.61); Brass -> AddElement(elZn, fractionmass=0.37); Brass -> AddElement(elPb, fractionmass=0.02); // Brass set 1 of nuts, washers and bolts for CUORE density = 3.6101*g/cm3; G4Material* Brassbolts1 = new G4Material(name="Brassbolts1", density, ncomponents=3); Brassbolts1 -> AddElement(elCu, fractionmass=0.61); Brassbolts1 -> AddElement(elZn, fractionmass=0.37); Brassbolts1 -> AddElement(elPb, fractionmass=0.02); // Brass set 2 of nuts, washers and bolts for CUORE density = 2.921522*g/cm3; G4Material* Brassbolts2 = new G4Material(name="Brassbolts2", density, ncomponents=3); Brassbolts2 -> AddElement(elCu, fractionmass=0.61); Brassbolts2 -> AddElement(elZn, fractionmass=0.37); Brassbolts2 -> AddElement(elPb, fractionmass=0.02); // Brass springs for cryostat of Andrea Giuliani in Canfranc density = 0.816868*g/cm3; G4Material* Brasssprings = new G4Material(name="Brasssprings", density, ncomponents=3); Brasssprings -> AddElement(elCu, fractionmass=0.61); Brasssprings -> AddElement(elZn, fractionmass=0.37); Brasssprings -> AddElement(elPb, fractionmass=0.02); // Brass contacts on CCD of disassembled camera for Cygnus density = 0.340414*g/cm3; G4Material* Brasscontacts = new G4Material(name="Brasscontacts", density, ncomponents=3); Brasscontacts -> AddElement(elCu, fractionmass=0.61); Brasscontacts -> AddElement(elZn, fractionmass=0.37); Brasscontacts -> AddElement(elPb, fractionmass=0.02); // FPD sensor for CUORE density = 5.697609*g/cm3; G4Material* FPD = new G4Material(name="FPD", density, ncomponents=3); FPD -> AddElement(elCu, fractionmass=0.61); FPD -> AddElement(elZn, fractionmass=0.37); FPD -> AddElement(elPb, fractionmass=0.02); // CT sensor for CUORE density = 3.55638*g/cm3; G4Material* CT = new G4Material(name="CT", density, ncomponents=3); CT -> AddElement(elCu, fractionmass=0.61); CT -> AddElement(elZn, fractionmass=0.37); CT -> AddElement(elPb, fractionmass=0.02); // CMN sensor for CUORE density = 4.786537*g/cm3; G4Material* CMN = new G4Material(name="CMN", density, ncomponents=3); CMN -> AddElement(elCu, fractionmass=0.61); CMN -> AddElement(elZn, fractionmass=0.37); CMN -> AddElement(elPb, fractionmass=0.02); // Brass set 1 of screws, TUM - GERDA (Janicsko) density = 1.9412*g/cm3; G4Material* Brassscrews1 = new G4Material(name="Brassscrews1", density, ncomponents=3); Brassscrews1 -> AddElement(elCu, fractionmass=0.61); Brassscrews1 -> AddElement(elZn, fractionmass=0.37); Brassscrews1 -> AddElement(elPb, fractionmass=0.02); // Brass set 2 of screws, TUM - GERDA (Janicsko) density = 2.02*g/cm3; G4Material* Brassscrews2 = new G4Material(name="Brassscrews2", density, ncomponents=3); Brassscrews2 -> AddElement(elCu, fractionmass=0.61); Brassscrews2 -> AddElement(elZn, fractionmass=0.37); Brassscrews2 -> AddElement(elPb, fractionmass=0.02); // Proportional Gas density = 0.001563*g/cm3; G4Material* ProportionalGas = new G4Material(name="ProportionalGas", density, ncomponents=3); ProportionalGas -> AddElement(elAr, fractionmass=0.90); ProportionalGas -> AddElement(elC, fractionmass=0.0748694); ProportionalGas -> AddElement(elH, fractionmass=0.0251306); // composition corrected (original : C-> 0.02 H-> 0.08 ) // Argon Gas (@ 1 bar, 293 K) density = 0.001663*g/cm3; G4Material* ArGas = new G4Material(name="ArgonGas", density, 1); ArGas -> AddElement(elAr, natoms=1); // Acrylic Juno, Milano Bicocca G4Material* Acrylic_juno = new G4Material (name="Acrylic_juno", density= 1.18534*g/cm3, ncomponents=3); Acrylic_juno->AddElement(elC,natoms=5); Acrylic_juno->AddElement(elH,natoms=8); Acrylic_juno->AddElement(elO,natoms=2); // Acrylic CYGNUS, GSSI G4Material* Acrylic_cygno = new G4Material (name="Acrylic_cygno", density= 1.25*g/cm3, ncomponents=3); Acrylic_cygno->AddElement(elC,natoms=5); Acrylic_cygno->AddElement(elH,natoms=8); Acrylic_cygno->AddElement(elO,natoms=2); // Acrylic SuperKamiokande, Napoli G4Material* Acrylic_sk = new G4Material (name="Acrylic_sk", density= 1.17178*g/cm3, ncomponents=3); Acrylic_sk->AddElement(elC,natoms=5); Acrylic_sk->AddElement(elH,natoms=8); Acrylic_sk->AddElement(elO,natoms=2); // Acrylic DS-20K G4Material* Acrylic_ds20k = new G4Material (name="Acrylic_ds20k", density= 1.17707*g/cm3, ncomponents=3); Acrylic_ds20k->AddElement(elC,natoms=5); Acrylic_ds20k->AddElement(elH,natoms=8); Acrylic_ds20k->AddElement(elO,natoms=2); // Acrylic Darkside, forks G4Material* Acrylic_fork = new G4Material (name="Acrylic_fork", density= 0.364533*g/cm3, ncomponents=3); Acrylic_fork->AddElement(elC,natoms=5); Acrylic_fork->AddElement(elH,natoms=8); Acrylic_fork->AddElement(elO,natoms=2); // Acrylic Darkside, cages G4Material* Acrylic_cage = new G4Material (name="Acrylic_cage", density= 0.1964*g/cm3, ncomponents=3); Acrylic_cage->AddElement(elC,natoms=5); Acrylic_cage->AddElement(elH,natoms=8); Acrylic_cage->AddElement(elO,natoms=2); // PolyVinylChlorid G4Material* PVC = new G4Material (name="PVC", density= 1.38*g/cm3, ncomponents=3); PVC->AddElement(elC,natoms=2); PVC->AddElement(elH,natoms=3); PVC->AddElement(elCl,natoms=1); // PolyVinylChlorid, Darkside G4Material* PVCds = new G4Material (name="PVCds", density= 1.4178*g/cm3, ncomponents=3); PVCds->AddElement(elC,natoms=2); PVCds->AddElement(elH,natoms=3); PVCds->AddElement(elCl,natoms=1); // PolyVinylChlorid, Darkside, Virginia Tech G4Material* PVCvt = new G4Material (name="PVCvt", density= 1.2747*g/cm3, ncomponents=3); PVCvt->AddElement(elC,natoms=2); PVCvt->AddElement(elH,natoms=3); PVCvt->AddElement(elCl,natoms=1); // gloves made of PolyVinylChlorid, Stella G4Material* PVCg = new G4Material (name="PVCg", density= 0.3045*g/cm3, ncomponents=3); PVCg->AddElement(elC,natoms=2); PVCg->AddElement(elH,natoms=3); PVCg->AddElement(elCl,natoms=1); // Kynar, Polyvinylidene fluoride (PVDF) G4Material* PVDF = new G4Material (name="PVDF", density= 1.78*g/cm3, ncomponents=3); PVDF->AddElement(elC,natoms=2); PVDF->AddElement(elH,natoms=2); PVDF->AddElement(elF,natoms=2); // HD1000 UHMW Polyethylene density = 0.93*g/cm3; G4Material* HD1000 = new G4Material(name="HD1000", density, ncomponents=2); HD1000 -> AddElement(elH, natoms=4); HD1000 -> AddElement(elC, natoms=2); // HDPE Polyethylene, T2K density = 0.965*g/cm3; G4Material* HDPE = new G4Material(name="HDPE", density, ncomponents=2); HDPE -> AddElement(elH, natoms=4); HDPE -> AddElement(elC, natoms=2); // MONEL 400 G4Material* MONEL = new G4Material (name="MONEL",density=8.83*g/cm3,ncomponents=5); MONEL->AddElement(elNi,fractionmass= 0.64); MONEL->AddElement(elCu,fractionmass= 0.31); MONEL->AddElement(elMn,fractionmass= 0.02); MONEL->AddElement(elFe,fractionmass= 0.025); MONEL->AddElement(elSi,fractionmass= 0.005); // Polyethylene density = 0.945*g/cm3; G4Material* PE = new G4Material(name="PE", density, ncomponents=2); PE -> AddElement(elH, natoms=4); PE -> AddElement(elC, natoms=2); // Polystyrene density = 1.05*g/cm3; G4Material* PS = new G4Material(name="PS", density, ncomponents=2); PS -> AddElement(elH, natoms=8); PS -> AddElement(elC, natoms=8); // expanded Polystyrene (EPS) density = 0.023*g/cm3; G4Material* EPS = new G4Material(name="EPS", density, ncomponents=2); EPS -> AddElement(elH, natoms=8); EPS -> AddElement(elC, natoms=8); // Viton (Polyfluoroelastomer) density = 0.631708*g/cm3; G4Material* viton_cyg = new G4Material(name="viton_cyg", density, ncomponents=3); viton_cyg -> AddElement(elH, fractionmass=0.009417); viton_cyg -> AddElement(elC, fractionmass=0.280555); viton_cyg -> AddElement(elF, fractionmass=0.710028); // Viton (Polyfluoroelastomer) for CCD Hamamatsu, CYGNUS density = 1.8*g/cm3; G4Material* viton = new G4Material(name="viton", density, ncomponents=3); viton -> AddElement(elH, fractionmass=0.009417); viton -> AddElement(elC, fractionmass=0.280555); viton -> AddElement(elF, fractionmass=0.710028); /* Phos Bronze C 51000, CDA 510 Phosphor Bronze , ASTM B 139 */ density = 8.86*g/cm3; G4Material* PBronze = new G4Material(name="PBronze", density, ncomponents=3); PBronze->AddElement(elCu, fractionmass=0.948); PBronze->AddElement(elSn, fractionmass=0.05); PBronze->AddElement(elP, fractionmass=0.002); /* ***Bronze*** mixture of 60% Cu and 40% Sn is assumed */ density = 8.8*g/cm3; G4Material* Bronze = new G4Material(name="Bronze", density, ncomponents=2); Bronze->AddElement(elCu, natoms=14); Bronze->AddElement(elSn, natoms=5); /* ***silver plated bronze clamps - CRESST*** mixture of 60% Cu and 40% Sn is assumed */ density = 2.978*g/cm3; G4Material* Clamps = new G4Material(name="Clamps", density, ncomponents=2); Clamps->AddElement(elCu, natoms=14); Clamps->AddElement(elSn, natoms=5); /* *** Germanium/Boron mixture of 1 per cent assumed (in dead layer *** */ density = 5.18*g/cm3; G4Material* GeB = new G4Material(name="Germanium/Boron", density, ncomponents=2); GeB -> AddElement(elGeNat, fractionmass=0.99); GeB -> AddElement(elB, fractionmass=0.01); // gaseous hydrogen density=0.0899*mg/cm3; G4Material *hydrogen = new G4Material(name="Hydrogen",density,ncomponents=1, kStateGas,300.00*kelvin,1.0*atmosphere); hydrogen->AddElement(elH,natoms=1); //LNGS rock from Wulandari et al. Astrop Phys 22 (2004) 313 density = 2.71*g/cm3; G4Material* lngsRock = new G4Material(name="LNGSRock",density,ncomponents=7); lngsRock->AddElement(elC,fractionmass=11.88*perCent); lngsRock->AddElement(elO,fractionmass=48.92*perCent); lngsRock->AddElement(elMn,fractionmass=5.58*perCent); lngsRock->AddElement(elAl,fractionmass=1.03*perCent); lngsRock->AddElement(elSi,fractionmass=1.27*perCent); lngsRock->AddElement(elK,fractionmass=1.03*perCent); lngsRock->AddElement(elCa,fractionmass=30.29*perCent); //monazite-ce sand from IRMM, density = 3.694323*g/cm3; G4Material* monazitece = new G4Material(name="monazitece",density,ncomponents=10); monazitece->AddElement(elCe,natoms=43); monazitece->AddElement(elLa,natoms=33); monazitece->AddElement(elTh,natoms=11); monazitece->AddElement(elNd,natoms=9); monazitece->AddElement(elCa,natoms=4); monazitece->AddElement(elPr,natoms=3); monazitece->AddElement(elSm,natoms=1); monazitece->AddElement(elP,natoms=94); monazitece->AddElement(elSi,natoms=2); monazitece->AddElement(elO,natoms=400); //rock density = 2.5*g/cm3; G4Material* Rock = new G4Material(name="Rock",density,ncomponents=7); Rock->AddElement(elC,fractionmass=11.88*perCent); Rock->AddElement(elO,fractionmass=48.92*perCent); Rock->AddElement(elMn,fractionmass=5.58*perCent); Rock->AddElement(elAl,fractionmass=1.03*perCent); Rock->AddElement(elSi,fractionmass=1.27*perCent); Rock->AddElement(elK,fractionmass=1.03*perCent); Rock->AddElement(elCa,fractionmass=30.29*perCent); //rock1 density = 1.7*g/cm3; G4Material* Rock1 = new G4Material(name="Rock1",density,ncomponents=7); Rock1->AddElement(elC,fractionmass=11.88*perCent); Rock1->AddElement(elO,fractionmass=48.92*perCent); Rock1->AddElement(elMn,fractionmass=5.58*perCent); Rock1->AddElement(elAl,fractionmass=1.03*perCent); Rock1->AddElement(elSi,fractionmass=1.27*perCent); Rock1->AddElement(elK,fractionmass=1.03*perCent); Rock1->AddElement(elCa,fractionmass=30.29*perCent); //rock wool density = 0.0465465*g/cm3; G4Material* Rockwool = new G4Material(name="Rockwool",density,ncomponents=7); Rockwool->AddElement(elC,fractionmass=11.88*perCent); Rockwool->AddElement(elO,fractionmass=48.92*perCent); Rockwool->AddElement(elMn,fractionmass=5.58*perCent); Rockwool->AddElement(elAl,fractionmass=1.03*perCent); Rockwool->AddElement(elSi,fractionmass=1.27*perCent); Rockwool->AddElement(elK,fractionmass=1.03*perCent); Rockwool->AddElement(elCa,fractionmass=30.29*perCent); //Peruvian soil from Nist, SRM 4355 //density = 0.91*g/cm3; Nist, SRM 4355 //density = 1.21*g/cm3; ?? //density = 1.12*g/cm3; ?? //density = 0.94122*g/cm3; IAEA 412, Pacific Ocean Sediment //density = 1.4865*g/cm3; IAEA TEL-2014-05, Sediment //density = 0.912959*g/cm3; NIST, SRM 4357, Ocean Sediment density = 0.912959*g/cm3; G4Material* peruvSoil = new G4Material(name="peruvSoil",density,ncomponents=8); peruvSoil->AddElement(elSi,fractionmass=33*perCent); peruvSoil->AddElement(elMn,fractionmass=1.6*perCent); peruvSoil->AddElement(elCa,fractionmass=2.2*perCent); peruvSoil->AddElement(elFe,fractionmass=4.5*perCent); peruvSoil->AddElement(elAl,fractionmass=8.2*perCent); peruvSoil->AddElement(elK,fractionmass=1.9*perCent); peruvSoil->AddElement(elNa,fractionmass=1.9*perCent); peruvSoil->AddElement(elO,fractionmass=46.7*perCent); //soil from IAEA TEL-2012-03 //density = 0.999403*g/cm3; IAEA CU-2006-03 //density = 0.73787*g/cm3; IAEA TEL-2012-03 //density = 0.72441*g/cm3; IAEA 410 Bikini Atol Sediment //density = 0.750568*g/cm3; IAEA TEL-2015-05, Syrian soil //density = 0.342995*g/cm3; IAEA 465, Baltic Sea Sediment, bottle #571 density = 0.999403*g/cm3; G4Material* iaeaSoil = new G4Material(name="iaeaSoil",density,ncomponents=8); iaeaSoil->AddElement(elSi,fractionmass=33*perCent); iaeaSoil->AddElement(elMn,fractionmass=1.6*perCent); iaeaSoil->AddElement(elCa,fractionmass=2.2*perCent); iaeaSoil->AddElement(elFe,fractionmass=4.5*perCent); iaeaSoil->AddElement(elAl,fractionmass=8.2*perCent); iaeaSoil->AddElement(elK,fractionmass=1.9*perCent); iaeaSoil->AddElement(elNa,fractionmass=1.9*perCent); iaeaSoil->AddElement(elO,fractionmass=46.7*perCent); //residue from pipeline, sample #1472, NYU1 //density = 1.843789*g/cm3; (BEFORE change to type 2 new std vial) //density = 1.222527*g/cm3; (AFTER change to type 2 new std vial) density = 1.222527*g/cm3; G4Material* nyu1 = new G4Material(name="nyu1",density,ncomponents=6); nyu1->AddElement(elSi,fractionmass=31.3*perCent); nyu1->AddElement(elFe,fractionmass=3.4*perCent); nyu1->AddElement(elS,fractionmass=0.6*perCent); nyu1->AddElement(elH,fractionmass=3.7*perCent); nyu1->AddElement(elC,fractionmass=6.1*perCent); nyu1->AddElement(elO,fractionmass=54.9*perCent); //residue from pipeline, sample #1467, NYU2 //density = 1.638924*g/cm3; (BEFORE change to type 2 new std vial) //density = 1.086690*g/cm3; (AFTER change to type 2 new std vial) density = 1.086690*g/cm3; G4Material* nyu2 = new G4Material(name="nyu2",density,ncomponents=6); nyu2->AddElement(elSi,fractionmass=31.3*perCent); nyu2->AddElement(elFe,fractionmass=3.4*perCent); nyu2->AddElement(elS,fractionmass=0.6*perCent); nyu2->AddElement(elH,fractionmass=3.7*perCent); nyu2->AddElement(elC,fractionmass=6.1*perCent); nyu2->AddElement(elO,fractionmass=54.9*perCent); //residue from pipeline, sample #1468, NYU3 //density = 1.485275*g/cm3; (BEFORE change to type 2 new std vial) //density = 0.984813*g/cm3; (AFTER change to type 2 new std vial) density = 0.984813*g/cm3; G4Material* nyu3 = new G4Material(name="nyu3",density,ncomponents=6); nyu3->AddElement(elSi,fractionmass=31.3*perCent); nyu3->AddElement(elFe,fractionmass=3.4*perCent); nyu3->AddElement(elS,fractionmass=0.6*perCent); nyu3->AddElement(elH,fractionmass=3.7*perCent); nyu3->AddElement(elC,fractionmass=6.1*perCent); nyu3->AddElement(elO,fractionmass=54.9*perCent); //hay from IAEA TEL-2012-03 //density = 0.454097*g/cm3; //rice from IAEA TEL-2015 //density = 0.862712*g/cm3; //clover from IAEA TEL-2016-04 //density = 0.433768*g/cm3; density = 0.433768*g/cm3; G4Material* iaeaHay = new G4Material(name="iaeaHay",density,ncomponents=3); iaeaHay->AddElement(elH,fractionmass=6.2162*perCent); iaeaHay->AddElement(elC,fractionmass=44.4462*perCent); iaeaHay->AddElement(elO,fractionmass=49.3376*perCent); //clover from IAEA TEL-2016-04 //density = 0.412585*g/cm3; density = 0.412585*g/cm3; G4Material* iaeaHay32 = new G4Material(name="iaeaHay32",density,ncomponents=3); iaeaHay32->AddElement(elH,fractionmass=6.2162*perCent); iaeaHay32->AddElement(elC,fractionmass=44.4462*perCent); iaeaHay32->AddElement(elO,fractionmass=49.3376*perCent); //spruce needles from IAEA TEL-2016-03 //density = 0.1592056*g/cm3; density = 0.1592056*g/cm3; G4Material* iaeaspruce = new G4Material(name="iaeaspruce",density,ncomponents=3); iaeaspruce->AddElement(elH,fractionmass=6.2162*perCent); iaeaspruce->AddElement(elC,fractionmass=44.4462*perCent); iaeaspruce->AddElement(elO,fractionmass=49.3376*perCent); //gelatin 405 for photographic emulsion (OPERA) density = 0.624181*g/cm3; G4Material* gelatin405 = new G4Material(name="gelatin405",density,ncomponents=5); gelatin405->AddElement(elH,fractionmass=5.7171*perCent); gelatin405->AddElement(elC,fractionmass=34.5824*perCent); gelatin405->AddElement(elN,fractionmass=12.2063*perCent); gelatin405->AddElement(elO,fractionmass=47.2812*perCent); gelatin405->AddElement(elS,fractionmass=0.2130*perCent); //gelatin 406 for photographic emulsion (OPERA) density = 0.587013*g/cm3; G4Material* gelatin406 = new G4Material(name="gelatin406",density,ncomponents=5); gelatin406->AddElement(elH,fractionmass=5.7171*perCent); gelatin406->AddElement(elC,fractionmass=34.5824*perCent); gelatin406->AddElement(elN,fractionmass=12.2063*perCent); gelatin406->AddElement(elO,fractionmass=47.2812*perCent); gelatin406->AddElement(elS,fractionmass=0.2130*perCent); //gelatin P6406 for photographic emulsion (OPERA) density = 0.601903*g/cm3; G4Material* gelatin6406 = new G4Material(name="gelatin6406",density,ncomponents=5); gelatin6406->AddElement(elH,fractionmass=5.7171*perCent); gelatin6406->AddElement(elC,fractionmass=34.5824*perCent); gelatin6406->AddElement(elN,fractionmass=12.2063*perCent); gelatin6406->AddElement(elO,fractionmass=47.2812*perCent); gelatin6406->AddElement(elS,fractionmass=0.2130*perCent); //fish gelatin P6407 for photographic emulsion (OPERA) density = 0.636121*g/cm3; G4Material* gelatin6407 = new G4Material(name="gelatin6407",density,ncomponents=5); gelatin6407->AddElement(elH,fractionmass=5.7171*perCent); gelatin6407->AddElement(elC,fractionmass=34.5824*perCent); gelatin6407->AddElement(elN,fractionmass=12.2063*perCent); gelatin6407->AddElement(elO,fractionmass=47.2812*perCent); gelatin6407->AddElement(elS,fractionmass=0.2130*perCent); //pure gel + NaBr, NIT 04-SEP-2015, FAN 042-04X (NEWS) (NIST material database composition) density = 0.9972875*g/cm3; G4Material* gelnabr = new G4Material(name="gelnabr",density,ncomponents=8); gelnabr->AddElement(elH,fractionmass=1.41*perCent); gelnabr->AddElement(elC,fractionmass=7.2261*perCent); gelnabr->AddElement(elN,fractionmass=1.932*perCent); gelnabr->AddElement(elO,fractionmass=6.6101*perCent); gelnabr->AddElement(elS,fractionmass=0.189*perCent); gelnabr->AddElement(elBr,fractionmass= 34.9103*perCent); gelnabr->AddElement(elAg,fractionmass= 47.4105*perCent); gelnabr->AddElement(elI,fractionmass= 0.3120*perCent); // 1N HCl G4Material* HCl1N = new G4Material (name="HCl1N", density= 1.0159*g/cm3, ncomponents=3); HCl1N->AddElement(elO,fractionmass=90.81*perCent); HCl1N->AddElement(elH,fractionmass=5.77*perCent); HCl1N->AddElement(elCl,fractionmass=3.42*perCent); // 0.1N HCl G4Material* HCl01N = new G4Material (name="HCl01N", density= 1.001*g/cm3, ncomponents=3); HCl01N->AddElement(elO,fractionmass=93.78*perCent); HCl01N->AddElement(elH,fractionmass=5.87*perCent); HCl01N->AddElement(elCl,fractionmass=0.35*perCent); //Luciano, debugging only density = 6.0*g/cm3; G4Material* TeO2 = new G4Material(name="TeO2",density,ncomponents=2); TeO2->AddElement(elTe,natoms=1); TeO2->AddElement(elO,natoms=2); density = 6.0*g/cm3; G4Material* O2Te = new G4Material(name="O2Te",density,ncomponents=2); O2Te->AddElement(elO,natoms=2); O2Te->AddElement(elTe,natoms=1); } void MGGerdaLocalMaterialTable::AddNewMaterial() { }