{ TGeoMedium *medium[44]; char medname[44][15] = { "Helium","Scintilator","NaI","N2","Magnet","Air","N2inside","DrG","Lead","Sil", "N2track","HRW","HPW","HKW","CuFoil","vaccuum","RPCCbox","ICGas","Air","Scinfield", "H20Field","Leadfield","MPCGas","G10","MylPlane","Alu","MPCInGas","Uniformdist","MPCSensGas","MPCInSensGas", "TarNi","Plexifield","Ironfield","LH2field","PVC","Kapton","Glass","SF6","Isobutan","C2F4H2", "RPCGas","C","PP","PE" }; double medprop[44][8] = { // nmat, isvol,ifield,fieldm,tmaxfd,dmaxms,deemax,epsil,stmin { 0. , 1. , 6.0 , 45.0 , 10.0 , 0.5 , 0.1 , 1.0 }, // Helium (#1) { 1. , 1. , 6.0 , 45.0 , 0.5 , 0.05 , 0.01 , 0.1 }, // Scintilator (#2) { 1. , 1. , 6.0 , 45.0 , 2.0 , 0.2 , 0.01 , 0.5 }, // NaI (#3) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.2 , 0.01 , 5.0 }, // N2 (#4) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.2 , 0.01 , 10.0 }, // Magnet Coil (#5) { 0. , 1. , 6.0 , 45.0 , 2.0 , 0.2 , 0.1 , 2.0 }, // Air (#6) { 0. , 1. , 6.0 , 10.0 , 10.0 , 0.5 , 0.1 , 5.0 }, // N2 for inside stuff (#7) { 1. , 1. , 6.0 , 1.0 , 0.2 , 0.1 , 0.01 , 0.1 }, // Gas fot CDC (#8) { 0. , 1. , 6.0 , 45.0 , 2.0 , 0.2 , 0.001 , 0.5 }, // Lead (#9) { 1. , 1. , 6.0 , 45.0 , 0.2 , 0.1 , 0.001 , 0.1 }, // Silicon (#10) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // N2 for tracking (#11) { 0. , 0. , 0.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // Ni-CR (#12) { 0. , 0. , 0.0 , 45.0 , 30.0 , 0.5 , 0.01 , 0.1 }, // Cu-Be (#13) { 0. , 0. , 0.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // Tungsten (#14) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // Cu for Foil (#15) { 0. , 0. , 0.0 , 45.0 , 30.0 , 1.0 , 1.0 , 20.0 }, // Vaccum inside pipe (#16) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // Carbon (#17) { 1. , 1. , 6.0 , 45.0 , 30.0 , 0.2 , 0.5 , 1.0 }, // Gas for IC (#18) { 0. , 1. , 6.0 , 10.0 , 0.2 , 0.1 , 0.1 , 1.0 }, // Air+Field (#19) { 1. , 1. , 6.0 , 1.0 , 0.2 , 0.01 , 0.01 , 0.1 }, // Scintillatir + Field (#20) { 1. , 1. , 6.0 , 1.0 , 0.2 , 0.01 , 0.01 , 0.1 }, // Water + Field (#21) { 0. , 1. , 6.0 , 10.0 , 1.0 , 0.2 , 0.1 , 0.5 }, // Lead in field (#22) { 1. , 1. , 6.0 , 1.0 , 0.5 , 0.1 , 0.01 , 0.1 }, // Pure MPC Sens Gas (#23) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.1 , 0.2 }, // MPC Frame (#24) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // MPC Mylar Plane (#25) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.0005 , 0.0018}, // MPC Aluminum Plane (#26) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.05 , 0.1 }, // MPC Insensitive gas (#27) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.0005 , 0.002}, // Uniform ditr. (#28) { 1. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.01 , 0.1 }, // MPC Sensitive gas (#29) { 0. , 1. , 6.0 , 1.0 , 0.5 , 0.1 , 0.05 , 0.1 }, // Pure MPC Insensitive gas (#30) { 0. , 1. , 6.0 , 45.0 , 30.0 , 0.5 , 0.001 , 0.1 }, // Ni for target (#31) { 1. , 1. , 6.0 , 1.0 , 0.2 , 0.01 , 0.01 , 0.1 }, // Plex + field (#32) { 0. , 1. , 6.0 , 1.0 , 0.1 , 0.1 , 0.03 , 0.05 }, // Iron + field (#33) { 0. , 1. , 6.0 , 10.0 , 0.2 , 0.1 , 0.1 , 1.0 }, // liquid hydrongen + field (#34) { 0. , 1. , 6.0 , 45.0 , 0.5 , 0.05 , 0.01 , 0.1 }, // PVC (#35) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // Kapton (#36) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // glass (#37) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.2 , 0.5 , 0.005}, // sf6(#38) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // Isobutan (#39) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // C2f4h2 (#40) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // Rpcgas (#41) { 0. , 1. , 6.0 , 1.0 , 0.2 , 0.5 , 0.001 , 0.005}, // C (#42) { 0. , 1. , 6.0 , 10.0 , 0.2 , 0.1 , 0.1 , 0.005}, // PP (#43) { 0. , 1. , 6.0 , 10.0 , 0.2 , 0.1 , 0.1 , 0.005}}; // PE(#44) TGeoElement *hydrogen = new TGeoElement( "hydrogen","hydrogen" , 1.01 , 1.000 ); TGeoElement *carbon = new TGeoElement( "carbon","carbon" , 12.0107 , 6.000 ); TGeoElement *nitrogen = new TGeoElement( "nitrogen","nitrogen" , 14.01 , 7.000 ) ; TGeoElement *oxygen = new TGeoElement( "oxygen","oxygen" , 15.9994 , 8.000 ); TGeoElement *flourid = new TGeoElement( "flourid","flourid" , 18.998 , 9.000 ); TGeoElement *ch4 = new TGeoElement( "ch4","ch4" , 18. , 10. ); TGeoElement *sodium = new TGeoElement( "sodium","sodium" , 22.9890 , 11.000 ); TGeoElement *silicon = new TGeoElement( "silicon","silicon" , 28.0855 , 15.994 ); TGeoElement *sulfur = new TGeoElement( "sulfur","sulfur" , 32.065 , 16.0 ); TGeoElement *argon = new TGeoElement( "argon","argon" , 39.95 , 16.05 ); TGeoMaterial *mat1 = new TGeoMaterial("HYDROGEN",1.01,1,0.7080000E-01); TGeoMaterial *mat2 = new TGeoMaterial("DEUTERIUM",2.01,1,0.162); TGeoMaterial *mat3 = new TGeoMaterial( "Helium" , 4.0 , 2.0 , 0.0001785 , 529900.0 , 335500.0 ); TGeoMaterial *mat4 = new TGeoMaterial("LITHIUM",6.94,3,0.534); TGeoMaterial *mat5 = new TGeoMaterial("BERILLIUM",9.01,4,1.848); TGeoMaterial *mat6 = new TGeoMaterial("CARBON",12.01,6,2.265); TGeoMaterial *mat7 = new TGeoMaterial( "Nitrogen (Gas)" , 14.01 , 7.0 , 0.00125 , 30392.0 , 67890.0 ); TGeoMaterial *mat8 = new TGeoMaterial("NEON",20.18,10,1.207); TGeoMaterial *mat9 = new TGeoMaterial("ALUMINIUM",26.98,13,2.7); TGeoMaterial *mat10 = new TGeoMaterial("IRON",55.85,26,7.87); TGeoMaterial *mat11 = new TGeoMaterial("COPPER",63.54,29,8.96); TGeoMaterial *mat12 = new TGeoMaterial("TUNGSTEN",183.85,74,19.3); TGeoMaterial *mat13 = new TGeoMaterial("LEAD",207.19,82,11.35); TGeoMaterial *mat14 = new TGeoMaterial("URANIUM",238.03,92,18.95); TGeoMaterial *mat15 = new TGeoMaterial("AIR",14.61,7.,0.1205000E-02); TGeoMaterial *mat16 = new TGeoMaterial("VACUUM",0,0,0); TGeoMaterial *mat50=new TGeoMaterial( "Silicon" , 28.09 , 14.0 , 2.33 , 9.36 , 45.49 ); TGeoMaterial *mat51=new TGeoMaterial( "DC Gas" , 14.670 , 7. , 0.002 , 18300.0 , 45800.0 ); TGeoMaterial *mat57=new TGeoMaterial( "Scintillator" , 11.16 , 5. , 1.032 , 42.4 , 79.50 ); TGeoMaterial *mat61=new TGeoMaterial( "Water" , 14.322 , 7.0 , 1.0 , 36.1. , 78.8 ); TGeoMaterial *mat62=new TGeoMaterial( "Plexi Glas" , 12.397 , 6. , 1.18 , 34.4 , 70.8 ); TGeoMaterial *mat71=new TGeoMaterial( "(Ar+CH4 +Iso)" , 39.2809 , 16. , 0.00189339 , 14914.0 , 68790.0 ); TGeoMaterial *mat72=new TGeoMaterial( "G10 Plate" , 18.14 , 9. , 1.7 , 19.4. , 56.72 ); TGeoMaterial *mat73=new TGeoMaterial( "Mylar" , 12.88 , 6. , 1.39 , 28.7 , 56.1 ); TGeoMaterial *mat74=new TGeoMaterial( "MPC_ GAS+W (plane)" , 35.31 , 16.00 , 0.5781 , 38.09 , 176.04 ); TGeoMaterial *mat75=new TGeoMaterial( "MPC: Gas+AL+Mylar" , 30.68 , 14.00 , 0.00799 , 2965.0 , 12590.0 ); TGeoMaterial *mat76=new TGeoMaterial( "MPC: Gas+W (full)" , 30.74 , 14.00 , 0.00233 , 10160.0 , 43140.0 ); TGeoMaterial *mat77=new TGeoMaterial( "Nickel" , 58. , 28. , 8.9 , 12.6820/8.9 ,9.6 ); TGeoMaterial *mat758=new TGeoMaterial( "C" , 12.01 , 6.0 , 2.265 , 18.8 , 666. ); TGeoMaterial *mat921=new TGeoMaterial( "PVC" , 24.783 , 12.0 , 1.05 , 24.2 , 104. ); TGeoMixture *mat78=new TGeoMixture( "Gas (Ar + CH4)" , 2 , 0.00168 ); mat78->AddElement( argon , 0.9 ); mat78->AddElement( ch4 , 0.1 ); TGeoMixture *mat79=new TGeoMixture( "PE" , 2 , 0.9075 ); mat79->AddElement( carbon , 3. ); mat79->AddElement( hydrogen , 6. ); TGeoMixture *mat80=new TGeoMixture( "PE" , 2 , 0.9075 ); mat80->AddElement( carbon , 0.85629 ); mat80->AddElement( hydrogen , 0.143711 ); TGeoMixture *mat666=new TGeoMixture( "Kapton" , 4 , 1.42 ); mat666->AddElement( carbon , 22. ); mat666->AddElement( hydrogen , 10. ); mat666->AddElement( nitrogen , 2. ); mat666->AddElement( oxygen , 5. ); TGeoMixture *mat753=new TGeoMixture( "Glass" , 2 , 2.5 ); mat753->AddElement( silicon , 0.85629 ); mat753->AddElement( oxygen , 0.143711 ); TGeoMixture *mat754=new TGeoMixture( "SF6" , 2 , 0.00618 ); mat754->AddElement( sulfur , 1. ); mat754->AddElement( flourid , 6. ); TGeoMixture *mat755=new TGeoMixture( "Isobutan" , 2 , 0.0027 ); mat755->AddElement( carbon , 4. ); mat755->AddElement( hydrogen , 10. ); TGeoMixture *mat756=new TGeoMixture( "C2F4H2" , 3 , 0.0044 ); mat756->AddElement( carbon , 2. ); mat756->AddElement( hydrogen , 2. ); mat756->AddElement( flourid , 4. ); TGeoMixture *mat757=new TGeoMixture( "RPCGas" , 2 , 0.0044 ); mat757->AddElement(new TGeoElement( "rpc1","rpc1" , 12.75 , 6.25 ), 0.85 ); mat757->AddElement(new TGeoElement( "rpc2","rpc2" , 20.86 , 10.0 ), 0.1 ); mat757->AddElement(new TGeoElement( "rpc3","rpc3" , 4.15 , 2.42 ), 0.05 ); medium[0] = new TGeoMedium(medname[0],1,mat3,medprop[0]); medium[1] = new TGeoMedium(medname[1],2,mat57,medprop[1]); medium[3] = new TGeoMedium(medname[3],4,mat7,medprop[3]); medium[4] = new TGeoMedium(medname[4],5,mat10,medprop[4]); medium[5] = new TGeoMedium(medname[5],6,mat15,medprop[5]); medium[6] = new TGeoMedium(medname[6],7,mat7,medprop[6]); medium[7] = new TGeoMedium(medname[7],8,mat71,medprop[7]); medium[8] = new TGeoMedium(medname[8],9,mat13,medprop[8]); medium[9] = new TGeoMedium(medname[9],10,mat50,medprop[9]); medium[10] = new TGeoMedium(medname[10],11,mat7,medprop[10]); medium[11] = new TGeoMedium(medname[11],12,mat10,medprop[11]); medium[12] = new TGeoMedium(medname[12],13,mat11,medprop[12]); medium[13] = new TGeoMedium(medname[13],14,mat12,medprop[13]); medium[14] = new TGeoMedium(medname[14],15,mat11,medprop[14]); medium[15] = new TGeoMedium(medname[15],16,mat16,medprop[15]); medium[16] = new TGeoMedium(medname[16],17,mat6,medprop[16]); medium[17] = new TGeoMedium(medname[17],18,mat78,medprop[17]); medium[18] = new TGeoMedium(medname[18],19,mat15,medprop[18]); medium[19] = new TGeoMedium(medname[19],20,mat57,medprop[19]); medium[20] = new TGeoMedium(medname[20],21,mat61,medprop[19]); medium[21] = new TGeoMedium(medname[21],22,mat13,medprop[21]); medium[22] = new TGeoMedium(medname[22],23,mat71,medprop[22]); medium[23] = new TGeoMedium(medname[23],24,mat72,medprop[23]); medium[24] = new TGeoMedium(medname[24],25,mat73,medprop[24]); medium[25] = new TGeoMedium(medname[25],26,mat9,medprop[25]); medium[26] = new TGeoMedium(medname[26],27,mat75,medprop[26]); medium[27] = new TGeoMedium(medname[27],28,mat74,medprop[27]); medium[28] = new TGeoMedium(medname[28],29,mat76,medprop[28]); medium[29] = new TGeoMedium(medname[29],30,mat71,medprop[29]); medium[30] = new TGeoMedium(medname[30],31,mat77,medprop[30]); medium[31] = new TGeoMedium(medname[31],32,mat62,medprop[31]); medium[32] = new TGeoMedium(medname[32],33,mat1,medprop[32]); medium[33] = new TGeoMedium(medname[33],34,mat10,medprop[33]); medium[34] = new TGeoMedium(medname[34],35,mat921,medprop[34]); medium[35] = new TGeoMedium(medname[35],36,mat666,medprop[35]); medium[36] = new TGeoMedium(medname[36],37,mat753,medprop[36]); medium[37] = new TGeoMedium(medname[37],38,mat754,medprop[37]); medium[38] = new TGeoMedium(medname[38],39,mat755,medprop[38]); medium[39] = new TGeoMedium(medname[39],40,mat756,medprop[39]); medium[40] = new TGeoMedium(medname[40],41,mat757,medprop[40]); medium[41] = new TGeoMedium(medname[41],42,mat757,medprop[41]); medium[42] = new TGeoMedium(medname[42],43,mat79,medprop[42]); medium[43] = new TGeoMedium(medname[43],44,mat80,medprop[43]); }