void kopio_mc_18() { // // This file has been generated automatically via the root // utility g2root from an interactive version of GEANT // (see ROOT class TGeoManager for an example of use) // gSystem->Load("libGeom"); TGeoRotation *rot; TGeoNode *Node, *Node1; TGeoManager *kopio_mc_18 = new TGeoManager("kopio_mc_18","kopio_mc_18.C"); //-----------List of Materials and Mixtures-------------- TGeoMaterial *mat1 = new TGeoMaterial("HYDROGEN",1.01,1,0.7080000E-01); mat1->SetUniqueID( 1); TGeoMaterial *mat2 = new TGeoMaterial("DEUTERIUM",2.01,1,0.162); mat2->SetUniqueID( 2); TGeoMaterial *mat3 = new TGeoMaterial("HELIUM",4,2,0.125); mat3->SetUniqueID( 3); TGeoMaterial *mat4 = new TGeoMaterial("LITHIUM",6.94,3,0.534); mat4->SetUniqueID( 4); TGeoMaterial *mat5 = new TGeoMaterial("BERILLIUM",9.01,4,1.848); mat5->SetUniqueID( 5); TGeoMaterial *mat6 = new TGeoMaterial("CARBON",12.01,6,2.265); mat6->SetUniqueID( 6); TGeoMaterial *mat7 = new TGeoMaterial("NITROGEN",14.01,7,0.808); mat7->SetUniqueID( 7); TGeoMaterial *mat8 = new TGeoMaterial("NEON",20.18,10,1.207); mat8->SetUniqueID( 8); TGeoMaterial *mat9 = new TGeoMaterial("ALUMINIUM",26.98,13,2.7); mat9->SetUniqueID( 9); TGeoMaterial *mat10 = new TGeoMaterial("IRON",55.85,26,7.87); mat10->SetUniqueID( 10); TGeoMaterial *mat11 = new TGeoMaterial("COPPER",63.54,29,8.96); mat11->SetUniqueID( 11); TGeoMaterial *mat12 = new TGeoMaterial("TUNGSTEN",183.85,74,19.3); mat12->SetUniqueID( 12); TGeoMaterial *mat13 = new TGeoMaterial("LEAD",207.19,82,11.35); mat13->SetUniqueID( 13); TGeoMaterial *mat14 = new TGeoMaterial("URANIUM",238.03,92,18.95); mat14->SetUniqueID( 14); TGeoMaterial *mat15 = new TGeoMaterial("AIR",14.61,7.3,0.1205000E-02); mat15->SetUniqueID( 15); TGeoMaterial *mat16 = new TGeoMaterial("VACUUM",0,0,0); mat16->SetUniqueID( 16); TGeoMixture *mat17 = new TGeoMixture("SCINTILLATOR",2, 1.03200 ); mat17->SetUniqueID( 17); mat17->DefineElement(0,12.01,6,0.9153266); mat17->DefineElement(1,1.01,1,0.8467343E-01); TGeoMixture *mat18 = new TGeoMixture("MYLAR",3, 1.39000 ); mat18->SetUniqueID( 18); mat18->DefineElement(0,12.01,6,0.624935); mat18->DefineElement(1,1.01,1,0.4204392E-01); mat18->DefineElement(2,16,8,0.3330211); TGeoMixture *mat19 = new TGeoMixture("LUCITE",3, 1.42000 ); mat19->SetUniqueID( 19); mat19->DefineElement(0,12.01,6,0.5997204); mat19->DefineElement(1,1.01,1,0.8069509E-01); mat19->DefineElement(2,16,8,0.3195845); TGeoMixture *mat20 = new TGeoMixture("AEROGEL",2, 0.110000 ); mat20->SetUniqueID( 20); mat20->DefineElement(0,28.086,14,0.4674456); mat20->DefineElement(1,15.999,8,0.5325544); TGeoMaterial *mat21 = new TGeoMaterial("CARBONFIBEREPOXY",12.01,6,2); mat21->SetUniqueID( 21); TGeoMixture *mat22 = new TGeoMixture("WSFIBER",2, 1.03200 ); mat22->SetUniqueID( 22); mat22->DefineElement(0,12.01,6,0.922427); mat22->DefineElement(1,1.01,1,0.7757296E-01); TGeoMixture *mat23 = new TGeoMixture("PAPER",4, 1.00000 ); mat23->SetUniqueID( 23); mat23->DefineElement(0,12.01,6,0.51); mat23->DefineElement(1,1.01,1,0.6000000E-01); mat23->DefineElement(2,16,8,0.42); mat23->DefineElement(3,14.01,7,0.1000000E-01); TGeoMixture *mat24 = new TGeoMixture("CLADDING",3, 1.20000 ); mat24->SetUniqueID( 24); mat24->DefineElement(0,12.01,6,0.5997204); mat24->DefineElement(1,1.01,1,0.8069509E-01); mat24->DefineElement(2,16,8,0.3195845); TGeoMaterial *mat25 = new TGeoMaterial("GOLD",196.9665,79,19.32); mat25->SetUniqueID( 25); TGeoMaterial *mat26 = new TGeoMaterial("ARGON",39.95,18,0.1782000E-02); mat26->SetUniqueID( 26); TGeoMixture *mat27 = new TGeoMixture("MWPCHAMBER",5, 0.285000E-02); mat27->SetUniqueID( 27); mat27->DefineElement(0,63.54,29,0.457); mat27->DefineElement(1,183.85,74,0.9000000E-02); mat27->DefineElement(2,39.95,18,0.312); mat27->DefineElement(3,12.01,6,0.177); mat27->DefineElement(4,1.01,1,0.4500000E-01); TGeoMixture *mat28 = new TGeoMixture("G_10",6, 1.70000 ); mat28->SetUniqueID( 28); mat28->DefineElement(0,1.01,1,0.6200000E-01); mat28->DefineElement(1,12.01,6,0.368); mat28->DefineElement(2,16,8,0.307); mat28->DefineElement(3,28.09,14,0.219); mat28->DefineElement(4,10.82,5,0.2200000E-01); mat28->DefineElement(5,23,11,0.2200000E-01); TGeoMixture *mat29 = new TGeoMixture("MAGICGAS",3, 0.195200E-02); mat29->SetUniqueID( 29); mat29->DefineElement(0,39.95,18,0.7332293); mat29->DefineElement(1,12.01,6,0.2204276); mat29->DefineElement(2,1.01,1,0.4634303E-01); TGeoMixture *mat30 = new TGeoMixture("GOLD-TUNGSTENWIRE",3, 19.3481 ); mat30->SetUniqueID( 30); mat30->DefineElement(0,183.85,74,0.94); mat30->DefineElement(1,186.2,75,0.3000000E-01); mat30->DefineElement(2,196.97,79,0.3000000E-01); TGeoMixture *mat31 = new TGeoMixture("ETHANE",2, 0.135600E-02); mat31->SetUniqueID( 31); mat31->DefineElement(0,12.01,6,0.7985373); mat31->DefineElement(1,1.01,1,0.2014628); TGeoMixture *mat32 = new TGeoMixture("ARGON-ETHANE(50,50)",2, 0.156900E-02); mat32->SetUniqueID( 32); mat32->DefineElement(0,39.95,18,0.571); mat32->DefineElement(1,9.794,4.993,0.429); TGeoMaterial *mat33 = new TGeoMaterial("LOW-DENSITY-AIR",14.61,7.3,0.1205000E-05); mat33->SetUniqueID( 33); TGeoMaterial *mat34 = new TGeoMaterial("LOW-DENSITY AL",26.98,13,0.9620000E-01); mat34->SetUniqueID( 34); TGeoMaterial *mat35 = new TGeoMaterial("TITANIUM",47.867,22,4.54); mat35->SetUniqueID( 35); TGeoMixture *mat36 = new TGeoMixture("TI DIOXIDE",2, 2.33000 ); mat36->SetUniqueID( 36); mat36->DefineElement(0,47.867,22,0.5993339); mat36->DefineElement(1,16,8,0.4006661); TGeoMaterial *mat37 = new TGeoMaterial("VERY-LOW-DENSITY-AIR",14.61,7.3,0.1205000E-09); mat37->SetUniqueID( 37); TGeoMixture *mat38 = new TGeoMixture("WATER",2, 1.00000 ); mat38->SetUniqueID( 38); mat38->DefineElement(0,16.01,8,0.8879645); mat38->DefineElement(1,1.01,1,0.1120355); //-----------List of Tracking Media-------------- TGeoMedium *med1 = new TGeoMedium("ALUMINUM",1,9,0,0,0,1,0.1000000E+11,0.1829598,0.1000000E-02,0.2231814E-01); TGeoMedium *med2 = new TGeoMedium("COPPER",2,11,0,0,0,1,0.1000000E+11,0.25,0.1000000E-02,0.1964564E-01); TGeoMedium *med3 = new TGeoMedium("LEAD",3,13,0,0,0,1,0.1000000E+11,0.25,0.1000000E-02,0.3291413E-01); TGeoMedium *med4 = new TGeoMedium("ATMOSPHERE (AIR)",4,15,0,0,0,1,0.1000000E+11,0.2488534,0.1,0.7694501); TGeoMedium *med5 = new TGeoMedium("VACUUM",5,16,0,0,0,1,0.1000000E+11,0,0.1,0); TGeoMedium *med6 = new TGeoMedium("SCINTILLATOR",6,17,1,0,0,1,0.1000000E+11,0.3081371E-01,0.1000000E-01,0.8135138E-02); TGeoMedium *med7 = new TGeoMedium("MYLAR",7,18,0,0,0,1,0.1000000E+11,0.2125741,0.1000000E-01,0.2031698E-01); TGeoMedium *med8 = new TGeoMedium("LUCITE",8,19,0,0,0,1,0.1000000E+11,0.2124136,0.1000000E-01,0.1899158E-01); TGeoMedium *med9 = new TGeoMedium("AEROGEL3",9,20,1,0,0,1,0.1000000E+11,0.1277086E-01,0.1000000E-01,0.1563874E-01); TGeoMedium *med10 = new TGeoMedium("CARBONFIBEREPOXY",10,21,0,0,0,1,0.1000000E+11,0.2066109,0.1000000E-01,0.1715193E-01); TGeoMedium *med11 = new TGeoMedium("WSFIBER",11,22,0,0,0,1,0.1000000E+11,0.2191564,0.1000000E-01,0.2132961E-01); TGeoMedium *med12 = new TGeoMedium("IRON",12,10,0,0,0,1,0.1000000E+11,0.25,0.1000000E-01,0.1948038E-01); TGeoMedium *med13 = new TGeoMedium("PAPER",13,23,0,0,0,1,0.1000000E+11,0.2179931,0.1000000E-01,0.2365506E-01); TGeoMedium *med14 = new TGeoMedium("CLADDING",14,24,0,0,0,1,0.1000000E+11,0.2154476,0.1000000E-01,0.2065926E-01); TGeoMedium *med15 = new TGeoMedium("GOLD",15,25,0,0,0,1,0.1000000E+11,0.25,0.1000000E-01,0.2455679E-01); TGeoMedium *med16 = new TGeoMedium("ARGON",16,26,1,0,0,1,0.1000000E+11,0.1910578E-02,0.1000000E-03,0.2564451E-01); TGeoMedium *med17 = new TGeoMedium("MWPCHAMBER",17,27,1,0,0,1,0.1000000E+11,0.2566757E-02,0.1000000E-03,0.2598033E-01); TGeoMedium *med18 = new TGeoMedium("G_10",18,28,0,0,0,1,0.1000000E+11,0.2048386,0.1000000E-02,0.1998152E-01); TGeoMedium *med19 = new TGeoMedium("MAGICGAS",19,29,1,0,0,1,0.1000000E+11,0.1843721E-02,0.1000000E-03,0.1891250E-01); TGeoMedium *med20 = new TGeoMedium("GOLD-TUNGSTENWIRE",20,30,0,0,0,1,0.1000000E+11,0.25,0.1000000E-03,0.2372238E-01); TGeoMedium *med21 = new TGeoMedium("ETHANE",21,31,1,0,0,1,0.1000000E+11,0.1099004E-02,0.1000000E-03,0.6877135E-02); TGeoMedium *med22 = new TGeoMedium("ARGON-ETHANE(50,50)",22,32,1,0,0,1,0.1000000E+11,0.1549533E-02,0.1000000E-03,0.1684663E-01); TGeoMedium *med23 = new TGeoMedium("LOW-DENSITY-AIR",23,33,0,1,20,1,0.1000000E+11,0.2499637,0.1000000E-01,24.3364); TGeoMedium *med24 = new TGeoMedium("LOW-DENSITY AL",24,34,0,0,0,1,0.1000000E+11,0.2373279,0.1,0.1184024); TGeoMedium *med25 = new TGeoMedium("BERYLIUM",25,5,0,0,0,1,0.1000000E+11,0.2163378,0.1,0.1548684E-01); TGeoMedium *med26 = new TGeoMedium("TITANIUM",26,35,0,0,0,1,0.1000000E+11,0.143963,0.5000000E-03,0.2353502E-01); TGeoMedium *med27 = new TGeoMedium("TUNGSTEN",27,12,0,0,0,1,0.1000000E+11,0.25,0.1000000E-01,0.2374649E-01); TGeoMedium *med28 = new TGeoMedium("TI DIOXIDE",28,36,0,0,0,1,0.1000000E+11,0.1825245,0.1000000E-01,0.2668352E-01); TGeoMedium *med29 = new TGeoMedium("VERY-LOW-DENSITY-AIR",29,37,0,0,0,1,0.1000000E+11,0.2499996,0.1,7.253259); TGeoMedium *med30 = new TGeoMedium("WATER",30,38,0,0,0,1,0.1000000E+11,0.216568,0.1000000E-01,0.2345920E-01); TGeoMedium *med31 = new TGeoMedium("B ALUMINUM",31,9,0,1,20,1,0.1000000E+11,0.1829598,0.1000000E-02,0.2231814E-01); TGeoMedium *med32 = new TGeoMedium("B LEAD",32,13,0,1,20,1,0.1000000E+11,0.25,0.1000000E-02,0.3291413E-01); TGeoMedium *med33 = new TGeoMedium("B SCINTILLATOR",33,17,1,1,20,1,0.1000000E+11,0.3081371E-01,0.1000000E-01,0.8135138E-02); TGeoMedium *med34 = new TGeoMedium("B WSFIBER",34,22,0,1,20,1,0.1000000E+11,0.2191564,0.1000000E-01,0.2132961E-01); //-----------List of Rotation matrices-------------- //-----------List of Volumes-------------- TGeoVolume *ENVL = gGeoManager->MakeBox("ENVL",med4,600,600,5000); TGeoVolume *PRER = gGeoManager->MakeBox("PRER",med5,600,600,48.3094); PRER->SetVisibility(0); TGeoVolume *EHOL = gGeoManager->MakeBox("EHOL",med4,92,12,48.3094); TGeoVolume *CHOL = gGeoManager->MakeBox("CHOL",med4,95.72,22.06,37.3828); TGeoVolume *CAL0 = gGeoManager->MakeBox("CAL0",med4,242.66,110.3,29.705); TGeoVolume *CAL9 = gGeoManager->MakeBox("CAL9",med4,71.695,22.06,29.705); TGeoVolume *CLW1 = gGeoManager->MakeBox("CLW1",med12,5.515,5.515,29.705); TGeoVolume *CLW2 = gGeoManager->MakeBox("CLW2",med12,5.515,5.515,29.705); TGeoVolume *CLE1 = gGeoManager->MakeBox("CLE1",med13,5.51,5.51,29.705); TGeoVolume *CLE2 = gGeoManager->MakeBox("CLE2",med13,5.51,5.51,29.705); TGeoVolume *CLA1 = gGeoManager->MakeBox("CLA1",med4,5.5,5.5,29.705); TGeoVolume *CLA2 = gGeoManager->MakeBox("CLA2",med4,5.5,5.5,29.705); TGeoVolume *CLA3 = gGeoManager->MakeBox("CLA3",med4,5.5,5.5,28.125); TGeoVolume *CLA4 = gGeoManager->MakeBox("CLA4",med4,5.5,5.5,28.125); TGeoVolume *CALP = gGeoManager->MakeBox("CALP",med3,5.5,5.5,0.1375000E-01); TGeoVolume *CALH = gGeoManager->MakeTube("CALH",med4,0,0.7000000E-01,0.1375000E-01); TGeoVolume *CALC = gGeoManager->MakeTube("CALC",med14,0,0.5600000E-01,0.1375000E-01); TGeoVolume *CALF = gGeoManager->MakeTube("CALF",med11,0,0.5000000E-01,0.1375000E-01); TGeoVolume *CLS1 = gGeoManager->MakeBox("CLS1",med6,5.5,5.5,0.7500000E-01); TGeoVolume *CLS2 = gGeoManager->MakeBox("CLS2",med6,5.5,5.5,0.7500000E-01); TGeoVolume *CLHL = gGeoManager->MakeTube("CLHL",med4,0,0.7000000E-01,0.7500000E-01); TGeoVolume *CLCL = gGeoManager->MakeTube("CLCL",med14,0,0.5600000E-01,0.7500000E-01); TGeoVolume *CLFB = gGeoManager->MakeTube("CLFB",med11,0,0.5000000E-01,0.7500000E-01); TGeoVolume *CALR = gGeoManager->MakeBox("CALR",med13,5.5,5.5,0.2500000E-02); TGeoVolume *CLRH = gGeoManager->MakeTube("CLRH",med4,0,0.7000000E-01,0.2500000E-02); TGeoVolume *CLRC = gGeoManager->MakeTube("CLRC",med14,0,0.5600000E-01,0.2500000E-02); TGeoVolume *CLRF = gGeoManager->MakeTube("CLRF",med11,0,0.5000000E-01,0.2500000E-02); TGeoVolume *CLL0 = gGeoManager->MakeBox("CLL0",med11,5.5,5.5,0.2202309E-01); TGeoVolume *CLL1 = gGeoManager->MakeBox("CLL1",med12,5.5,5.5,0.4000000E-01); //-----------List of Nodes-------------- gGeoManager->SetTopVolume(ENVL); ENVL->AddNode(PRER,1,new TGeoTranslation(0,0,1496.712)); PRER->AddNode(EHOL,1,gGeoIdentity); ENVL->AddNode(CHOL,1,new TGeoTranslation(0,0,1582.346)); ENVL->AddNode(CAL0,1,new TGeoTranslation(0,-132.36,1574.668)); TGeoVolume *CAL1 = CAL0->Divide("CAL1",1,44,-242.66,11.03); TGeoVolume *CAL2 = CAL1->Divide("CAL2",2,20,-110.3,11.03); CAL2->AddNode(CLW1,1,gGeoIdentity); CLW1->AddNode(CLE1,1,gGeoIdentity); CLE1->AddNode(CLA1,1,gGeoIdentity); CLA1->AddNode(CLA3,1,new TGeoTranslation(0,0,1.58)); TGeoVolume *CAL3 = CLA3->Divide("CAL3",3,300,-28.125,0.1875); CAL3->AddNode(CALP,1,new TGeoTranslation(0,0,-0.8000000E-01)); TGeoVolume *CLPX = CALP->Divide("CLPX",1,12,-5.5,0.9166667); TGeoVolume *CLPY = CLPX->Divide("CLPY",2,12,-5.5,0.9166667); CLPY->AddNode(CALH,1,gGeoIdentity); CALH->AddNode(CALC,1,gGeoIdentity); CALC->AddNode(CALF,1,gGeoIdentity); CAL3->AddNode(CALR,1,new TGeoTranslation(0,0,-0.6375000E-01)); TGeoVolume *CLRX = CALR->Divide("CLRX",1,12,-5.5,0.9166667); TGeoVolume *CLRY = CLRX->Divide("CLRY",2,12,-5.5,0.9166667); CLRY->AddNode(CLRH,1,gGeoIdentity); CLRH->AddNode(CLRC,1,gGeoIdentity); CLRC->AddNode(CLRF,1,gGeoIdentity); CAL3->AddNode(CLS1,1,new TGeoTranslation(0,0,0.1375000E-01)); TGeoVolume *CLX1 = CLS1->Divide("CLX1",1,12,-5.5,0.9166667); TGeoVolume *CLY1 = CLX1->Divide("CLY1",2,12,-5.5,0.9166667); CLY1->AddNode(CLHL,1,gGeoIdentity); CLHL->AddNode(CLCL,1,gGeoIdentity); CLCL->AddNode(CLFB,1,gGeoIdentity); CAL3->AddNode(CALR,2,new TGeoTranslation(0,0,0.9125001E-01)); CLA1->AddNode(CLL0,1,new TGeoTranslation(0,0,-28.045)); CLA1->AddNode(CLL1,1,new TGeoTranslation(0,0,-29.665)); ENVL->AddNode(CAL0,2,new TGeoTranslation(0,132.36,1574.668)); ENVL->AddNode(CAL9,1,new TGeoTranslation(-167.415,0,1574.668)); TGeoVolume *CAL4 = CAL9->Divide("CAL4",1,13,-71.695,11.03); TGeoVolume *CAL5 = CAL4->Divide("CAL5",2,4,-22.06,11.03); CAL5->AddNode(CLW2,1,gGeoIdentity); CLW2->AddNode(CLE2,1,gGeoIdentity); CLE2->AddNode(CLA2,1,gGeoIdentity); CLA2->AddNode(CLA4,1,new TGeoTranslation(0,0,1.58)); TGeoVolume *CAL6 = CLA4->Divide("CAL6",3,300,-28.125,0.1875); CAL6->AddNode(CALP,1,new TGeoTranslation(0,0,-0.8000000E-01)); CAL6->AddNode(CALR,1,new TGeoTranslation(0,0,-0.6375000E-01)); CAL6->AddNode(CLS2,1,new TGeoTranslation(0,0,0.1375000E-01)); TGeoVolume *CLX2 = CLS2->Divide("CLX2",1,12,-5.5,0.9166667); TGeoVolume *CLY2 = CLX2->Divide("CLY2",2,12,-5.5,0.9166667); CLY2->AddNode(CLHL,1,gGeoIdentity); CAL6->AddNode(CALR,2,new TGeoTranslation(0,0,0.9125001E-01)); CLA2->AddNode(CLL0,1,new TGeoTranslation(0,0,-28.045)); CLA2->AddNode(CLL1,1,new TGeoTranslation(0,0,-29.665)); ENVL->AddNode(CAL9,2,new TGeoTranslation(167.415,0,1574.668)); gGeoManager->CloseGeometry(); }