void create_gfi_geo(const char* geoTag) { fGlobalTrans->SetTranslation(0.0,0.0,0.0); // ------- Load media from media file ----------------------------------- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader"); FairGeoInterface* geoFace = geoLoad->getGeoInterface(); TString geoPath = gSystem->Getenv("VMCWORKDIR"); TString medFile = geoPath + "/geometry/media_r3b.geo"; geoFace->setMediaFile(medFile); geoFace->readMedia(); gGeoMan = gGeoManager; // -------------------------------------------------------------------------- // ------- Geometry file name (output) ---------------------------------- TString geoFileName = geoPath + "/geometry/gfi_"; geoFileName = geoFileName + geoTag + ".geo.root"; // -------------------------------------------------------------------------- // ----------------- Get and create the required media ----------------- FairGeoMedia* geoMedia = geoFace->getMedia(); FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder(); FairGeoMedium* mAir = geoMedia->getMedium("Air"); if ( ! mAir ) Fatal("Main", "FairMedium Air not found"); geoBuild->createMedium(mAir); TGeoMedium* pMed2 = gGeoMan->GetMedium("Air"); if ( ! pMed2 ) Fatal("Main", "Medium Air not found"); FairGeoMedium* mVac = geoMedia->getMedium("vacuum"); if ( ! mVac ) Fatal("Main", "FairMedium vacuum not found"); geoBuild->createMedium(mVac); TGeoMedium* pMed1 = gGeoMan->GetMedium("vacuum"); if ( ! pMed1 ) Fatal("Main", "Medium vacuum not found"); FairGeoMedium* mGfi = geoMedia->getMedium("plasticForGFI"); if ( ! mGfi ) Fatal("Main", "FairMedium plasticForGFI not found"); geoBuild->createMedium(mGfi); TGeoMedium* pMed35 = gGeoMan->GetMedium("plasticForGFI"); if ( ! pMed35 ) Fatal("Main", "Medium plasticForGFI not found"); FairGeoMedium* mAl = geoMedia->getMedium("aluminium"); if ( ! mAl ) Fatal("Main", "FairMedium aluminium not found"); geoBuild->createMedium(mAl); TGeoMedium* pMed21 = gGeoMan->GetMedium("aluminium"); if ( ! pMed21 ) Fatal("Main", "Medium aluminium not found"); // -------------------------------------------------------------------------- // -------------- Create geometry and top volume ------------------------- gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom"); gGeoMan->SetName("GFIgeom"); TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); // -------------------------------------------------------------------------- // out-of-file geometry definition Double_t dx,dy,dz; Double_t a; Double_t thx, phx, thy, phy, thz, phz; Double_t z, density, radl, absl, w; Int_t nel, numed; // TRANSFORMATION MATRICES // Combi transformation: //GFI1 position //dx = -71.973310; //Justyna //dy = 0.000000; //Justyna //dz = 513.967775; //Justyna //dx = -72.164874; //Justyna new //dy = -0.010000; //Justyna new //dz = 513.910302; //Justyna new //dx = -73.381; //dE tracker //dy = 0.070; //dz = 513.421; //LABPOS(GFI1,-73.274339,0.069976,513.649524) dx = -73.274339; //dE tracker, correction due to wrong angle dy = 0.069976; dz = 513.649524; /* dx = 73.700000; dy = 0.000000; dz = 525.400000; */ // Rotation: thx = -106.700000; phx = 0.000000; // thx = 106.700000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = -16.700000; phz = 0.000000; TGeoRotation *pMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); // Combi transformation: dx = 0.000000; dy = 0.000000; dz = 0.000000; // Rotation: thx = 90.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 0.000000; phz = 0.000000; TGeoRotation *pMatrix7 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix6 = new TGeoCombiTrans("", dx,dy,dz,pMatrix7); // Combi transformation: dx = 0.000000; dy = 27.000000; dz = 0.000000; // Rotation: thx = 90.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 0.000000; phz = 0.000000; TGeoRotation *pMatrix9 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix8 = new TGeoCombiTrans("", dx,dy,dz,pMatrix9); // Combi transformation: dx = 0.000000; dy = -27.000000; dz = 0.000000; // Rotation: thx = 90.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 0.000000; phz = 0.000000; TGeoRotation *pMatrix11 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix10 = new TGeoCombiTrans("", dx,dy,dz,pMatrix11); // Combi transformation: dx = 27.000000; dy = 0.000000; dz = 0.000000; // Rotation: thx = 90.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 0.000000; phz = 0.000000; TGeoRotation *pMatrix13 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix12 = new TGeoCombiTrans("", dx,dy,dz,pMatrix13); // Combi transformation: dx = -27.000000; dy = 0.000000; dz = 0.000000; // Rotation: thx = 90.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 0.000000; phz = 0.000000; TGeoRotation *pMatrix15 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix14 = new TGeoCombiTrans("", dx,dy,dz,pMatrix15); // Combi transformation: //GFI2 position //dx = -145.270039; //Justyna //dy = 0.000000; //Justyna //dz = 730.318956; //Justyna //dx = -145.212570; //Justyna new //dy = -0.010000; //Justyna new //dz = 730.336197; //Justyna new //dx = -147.486; //dE tracker //dy = 0.070; //dz = 729.798; //LABPOS(GFI2,-147.135037,0.069976,729.680342) dx = -147.135037; //dE tracker, correction due to wrong angle dy = 0.069976; dz = 729.680342; /* dx = 141.800000; dy = 0.000000; dz = 727.300000; */ // Rotation: thx = -106.700000; phx = 0.000000; // thx = 106.700000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = -16.700000; phz = 0.000000; TGeoRotation *pMatrix5 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix4 = new TGeoCombiTrans("", dx,dy,dz,pMatrix5); // World definition TGeoVolume* pWorld = gGeoManager->GetTopVolume(); pWorld->SetVisLeaves(kTRUE); // SHAPES, VOLUMES AND GEOMETRICAL HIERARCHY // Shape: GFIBoxWorld type: TGeoBBox dx = 29.00000; dy = 29.00000; dz = 0.050000; TGeoShape *pGFIBoxWorld = new TGeoBBox("GFIBoxWorld", dx,dy,dz); // Volume: GFILogWorld TGeoVolume* pGFILogWorld = new TGeoVolume("GFILogWorld",pGFIBoxWorld, pMed2); pGFILogWorld->SetVisLeaves(kTRUE); // Global positioning TGeoCombiTrans *pGlobal1 = GetGlobalPosition(pMatrix2); TGeoCombiTrans *pGlobal2 = GetGlobalPosition(pMatrix4); TGeoVolumeAssembly *pw1 = new TGeoVolumeAssembly("GFI"); pWorld->AddNode(pw1, 0); pw1->AddNode( pGFILogWorld, 0, pGlobal1 ); pw1->AddNode( pGFILogWorld, 1, pGlobal2 ); // Shape: GFIBox type: TGeoBBox dx = 25.000000; dy = 25.000000; dz = 0.050000; TGeoShape *pGFIBox = new TGeoBBox("GFIBox", dx,dy,dz); // Volume: GFILog TGeoVolume* pGFILog = new TGeoVolume("GFILog",pGFIBox, pMed35); pGFILog->SetVisLeaves(kTRUE); pGFILogWorld->AddNode(pGFILog, 0, pMatrix6); // Shape: UpFrame type: TGeoBBox dx = 29.000000; dy = 2.000000; dz = 0.050000; TGeoShape *pUpFrame = new TGeoBBox("UpFrame", dx,dy,dz); // Volume: logicUpFrame TGeoVolume* plogicUpFrame = new TGeoVolume("logicUpFrame",pUpFrame, pMed21); plogicUpFrame->SetVisLeaves(kTRUE); pGFILogWorld->AddNode(plogicUpFrame, 0, pMatrix8); // Shape: DownFrame type: TGeoBBox dx = 29.000000; dy = 2.000000; dz = 0.050000; TGeoShape *pDownFrame = new TGeoBBox("DownFrame", dx,dy,dz); // Volume: logicDownFrame TGeoVolume* plogicDownFrame = new TGeoVolume("logicDownFrame",pDownFrame, pMed21); plogicDownFrame->SetVisLeaves(kTRUE); pGFILogWorld->AddNode(plogicDownFrame, 0, pMatrix10); // Shape: RightFrame type: TGeoBBox dx = 2.000000; dy = 25.000000; dz = 0.050000; TGeoShape *pRightFrame = new TGeoBBox("RightFrame", dx,dy,dz); // Volume: logicRightFrame TGeoVolume* plogicRightFrame = new TGeoVolume("logicRightFrame",pRightFrame, pMed21); plogicRightFrame->SetVisLeaves(kTRUE); pGFILogWorld->AddNode(plogicRightFrame, 0, pMatrix12); // Shape: LeftFrame type: TGeoBBox dx = 2.000000; dy = 25.000000; dz = 0.050000; TGeoShape *pLeftFrame = new TGeoBBox("LeftFrame", dx,dy,dz); // Volume: logicLeftFrame TGeoVolume* plogicLeftFrame = new TGeoVolume("logicLeftFrame",pLeftFrame, pMed21); plogicLeftFrame->SetVisLeaves(kTRUE); pGFILogWorld->AddNode(plogicLeftFrame, 0, pMatrix14); // Add the sensitive part // AddSensitiveVolume(pGFILog); // fNbOfSensitiveVol+=1; // --------------- Finish ----------------------------------------------- gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* geoFile = new TFile(geoFileName, "RECREATE"); top->Write(); geoFile->Close(); // -------------------------------------------------------------------------- }
void create_tof_geo(const char* geoTag) { fGlobalTrans->SetTranslation(0.0,0.0,0.0); // ------- Load media from media file ----------------------------------- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader"); FairGeoInterface* geoFace = geoLoad->getGeoInterface(); TString geoPath = gSystem->Getenv("VMCWORKDIR"); TString medFile = geoPath + "/geometry/media_r3b.geo"; geoFace->setMediaFile(medFile); geoFace->readMedia(); gGeoMan = gGeoManager; // -------------------------------------------------------------------------- // ------- Geometry file name (output) ---------------------------------- TString geoFileName = geoPath + "/geometry/tof_"; geoFileName = geoFileName + geoTag + ".geo.root"; // -------------------------------------------------------------------------- // ----------------- Get and create the required media ----------------- FairGeoMedia* geoMedia = geoFace->getMedia(); FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder(); FairGeoMedium* mAir = geoMedia->getMedium("Air"); if ( ! mAir ) Fatal("Main", "FairMedium Air not found"); geoBuild->createMedium(mAir); TGeoMedium* pMed2 = gGeoMan->GetMedium("Air"); if ( ! pMed2 ) Fatal("Main", "Medium Air not found"); FairGeoMedium* mTof = geoMedia->getMedium("plasticForTOF"); if ( ! mTof ) Fatal("Main", "FairMedium plasticForTOF not found"); geoBuild->createMedium(mTof); TGeoMedium* pMed34 = gGeoMan->GetMedium("plasticForTOF"); if ( ! pMed34 ) Fatal("Main", "Medium plasticForTOF not found"); // -------------------------------------------------------------------------- // -------------- Create geometry and top volume ------------------------- gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom"); gGeoMan->SetName("TOFgeom"); TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); // -------------------------------------------------------------------------- // out-of-file geometry definition Double_t dx,dy,dz; Double_t a; // Double_t thx, phx, thy, phy, thz, phz; Double_t z, density, w; Int_t nel, numed; // TRANSFORMATION MATRICES // Combi transformation: dx = -417.359574; //Justyna dy = 2.400000; //Justyna dz = 960.777114; //Justyna // dx = -421.33683; //Christoph // dy = 2.12; //Christoph // dz = 958.387337; //Christoph /* dx = -171.1; //position directrly (15cm) after DCH2 dy = 2.400000; dz = 548.95;*/ // dz = 0.; /* // Rotation: thx = -121.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = -31.000000; phz = 0.000000;*/ //this // TGeoRotation *pMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); //this TGeoRotation *gRot = new TGeoRotation(); gRot->RotateX(0.); gRot->RotateY(-31.000000); gRot->RotateZ(0.); TGeoCombiTrans* // pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); //this pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,gRot); //this /* PREVIOUS!!! // TRANSFORMATION MATRICES // Combi transformation: dx = 419.700000; dy = 0.000000; dz = 952.400000; // dz = 0.; // Rotation: thx = 121.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 31.000000; phz = 0.000000; TGeoRotation *pMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); */ //Top Volume TGeoVolume* pWorld = gGeoManager->GetTopVolume(); pWorld->SetVisLeaves(kTRUE); TGeoVolumeAssembly *ptof = new TGeoVolumeAssembly("TOF"); // SHAPES, VOLUMES AND GEOMETRICAL HIERARCHY // Shape: TOFBox type: TGeoBBox dx = 94.450000; dy = 73.450000; dz = 0.500000; TGeoShape *pTOFBox = new TGeoBBox("TOFBox", dx,dy,dz); // Volume: TOFLog TGeoVolume* pTOFLog = new TGeoVolume("TOFLog",pTOFBox, pMed34); pTOFLog->SetVisLeaves(kTRUE); TGeoCombiTrans *t0 = new TGeoCombiTrans("t0"); ptof->AddNode(pTOFLog, 0, t0); TGeoCombiTrans *pGlobal = GetGlobalPosition(pMatrix2); if (pGlobal) { pWorld->AddNode(ptof, 0, pGlobal); } else { pWorld->AddNode(ptof, 0, pMatrix2); } // AddSensitiveVolume(pTOFLog); // fNbOfSensitiveVol+=1; // --------------- Finish ----------------------------------------------- gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* geoFile = new TFile(geoFileName, "RECREATE"); top->Write(); geoFile->Close(); // -------------------------------------------------------------------------- }
void create_mtof_geo(const char* geoTag) { // ------- Load media from media file ----------------------------------- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader"); FairGeoInterface* geoFace = geoLoad->getGeoInterface(); TString geoPath = gSystem->Getenv("VMCWORKDIR"); TString medFile = geoPath + "/geometry/media_r3b.geo"; geoFace->setMediaFile(medFile); geoFace->readMedia(); gGeoMan = gGeoManager; // -------------------------------------------------------------------------- // ------- Geometry file name (output) ---------------------------------- TString geoFileName = geoPath + "/geometry/mtof_"; geoFileName = geoFileName + geoTag + ".geo.root"; // -------------------------------------------------------------------------- // ----------------- Get and create the required media ----------------- FairGeoMedia* geoMedia = geoFace->getMedia(); FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder(); FairGeoMedium* mAir = geoMedia->getMedium("Air"); if ( ! mAir ) Fatal("Main", "FairMedium Air not found"); geoBuild->createMedium(mAir); TGeoMedium* pMed2 = gGeoMan->GetMedium("Air"); if ( ! pMed2 ) Fatal("Main", "Medium Air not found"); FairGeoMedium* mTof = geoMedia->getMedium("plasticFormTOF"); if ( ! mTof ) Fatal("Main", "FairMedium plasticFormTOF not found"); geoBuild->createMedium(mTof); TGeoMedium* pMed34 = gGeoMan->GetMedium("plasticFormTOF"); if ( ! pMed34 ) Fatal("Main", "Medium plasticFormTOF not found"); // -------------------------------------------------------------------------- // -------------- Create geometry and top volume ------------------------- gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom"); gGeoMan->SetName("mTOFgeom"); TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); // -------------------------------------------------------------------------- // out-of-file geometry definition Double_t dx,dy,dz; Double_t a; Double_t thx, phx, thy, phy, thz, phz; Double_t z, density, w; Int_t nel, numed; // TRANSFORMATION MATRICES // Combi transformation: //NTF position //dx = -154.815998;//Justyna //dy = 0.000000; //Justyna //dz = 761.755160;//Justyna //dx = -157.536214;//Justyna new //dy = -0.010000; //Justyna new ////dz = 760.939056;//Justyna new //dz = 760.139056;//Justyna //dx = -155.709;//dE tracker //dy = 0.524; ////dz = 761.487; //dz = 760.687; //try -0.8 like Justyna //LABPOS(FTF,-155.824045,0.523976,761.870346) dx = -155.824045;//dE tracker, correction due to wrong angle dy = 0.523976; dz = 761.870346; // Rotation: thx = -106.700000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = -16.700000; phz = 0.000000; /* dx = -171.1; dy = 2.400000; dz = 548.95; // dz = 0.; // Rotation: thx = -121.000000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = -31.000000; phz = 0.000000;*/ TGeoRotation *pMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); /* // TRANSFORMATION MATRICES // Combi transformation: dx = 151.000000; dy = 0.000000; dz = 758.000000; // Rotation: thx = 106.700000; phx = 0.000000; thy = 90.000000; phy = 90.000000; thz = 16.700000; phz = 0.000000; TGeoRotation *pMatrix3 = new TGeoRotation("",thx,phx,thy,phy,thz,phz); TGeoCombiTrans* pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); */ // Shape: World type: TGeoBBox TGeoVolume* pWorld = gGeoManager->GetTopVolume(); pWorld->SetVisLeaves(kTRUE); // Create a global Mother Volume /* dx = 200.000000; dy = 200.000000; dz = 200.000000; TGeoShape *pBoxWorld = new TGeoBBox("mTofBoxWorld", dx,dy,dz); TGeoVolume* pWorld = new TGeoVolume("mTofBoxLogWorld",pBoxWorld, pMed2); pWorld->SetVisLeaves(kTRUE); TGeoCombiTrans *pGlobalc = GetGlobalPosition(); // add the sphere as Mother Volume pAWorld->AddNode(pWorld, 0, pGlobalc); */ // SHAPES, VOLUMES AND GEOMETRICAL HIERARCHY // Shape: mTOFBox type: TGeoBBox dx = 24.000000; dy = 24.000000; //dz = 0.250000; //wrong: should be 0.5->1cm total dz = 0.500000; /* dx = 94.450000; //TFW size dy = 73.450000; dz = 0.500000;*/ TGeoShape *pmTOFBox = new TGeoBBox("mTOFBox", dx,dy,dz); // Volume: mTOFLog TGeoVolume * pmTOFLog = new TGeoVolume("mTOFLog",pmTOFBox, pMed34); pmTOFLog->SetVisLeaves(kTRUE); TGeoVolumeAssembly *pmTof = new TGeoVolumeAssembly("mTOF"); TGeoCombiTrans *t0 = new TGeoCombiTrans("t0"); pmTof->AddNode(pmTOFLog, 0, t0); TGeoCombiTrans *pGlobal = GetGlobalPosition(pMatrix2); if (pGlobal){ pWorld->AddNode(pmTof, 0, pGlobal); }else{ pWorld->AddNode(pmTof, 0, pMatrix2); } // AddSensitiveVolume(pmTOFLog); // fNbOfSensitiveVol+=1; // --------------- Finish ----------------------------------------------- gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* geoFile = new TFile(geoFileName, "RECREATE"); top->Write(); geoFile->Close(); // -------------------------------------------------------------------------- }
void create_dch_geo(const char* geoTag) { fGlobalTrans->SetTranslation(0.0,0.0,0.0); // ------- Load media from media file ----------------------------------- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader"); FairGeoInterface* geoFace = geoLoad->getGeoInterface(); TString geoPath = gSystem->Getenv("VMCWORKDIR"); TString medFile = geoPath + "/geometry/media_r3b.geo"; geoFace->setMediaFile(medFile); geoFace->readMedia(); gGeoMan = gGeoManager; // -------------------------------------------------------------------------- // ------- Geometry file name (output) ---------------------------------- TString geoFileName = geoPath + "/geometry/dch_"; geoFileName = geoFileName + geoTag + ".geo.root"; // -------------------------------------------------------------------------- // ----------------- Get and create the required media ----------------- FairGeoMedia* geoMedia = geoFace->getMedia(); FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder(); FairGeoMedium* mAir = geoMedia->getMedium("Air"); if ( ! mAir ) Fatal("Main", "FairMedium Air not found"); geoBuild->createMedium(mAir); TGeoMedium* pMed2 = gGeoMan->GetMedium("Air"); if ( ! pMed2 ) Fatal("Main", "Medium Air not found"); FairGeoMedium* mAl = geoMedia->getMedium("aluminium"); if ( ! mAl ) Fatal("Main", "FairMedium aluminium not found"); geoBuild->createMedium(mAl); TGeoMedium* pMed21 = gGeoMan->GetMedium("aluminium"); if ( ! pMed21 ) Fatal("Main", "Medium aluminium not found"); FairGeoMedium* mMylar = geoMedia->getMedium("mylar"); if ( ! mMylar ) Fatal("Main", "FairMedium mylar not found"); geoBuild->createMedium(mMylar); TGeoMedium* pMed15 = gGeoMan->GetMedium("mylar"); if ( ! pMed15 ) Fatal("Main", "Medium mylar not found"); FairGeoMedium* mHe = geoMedia->getMedium("helium"); if ( ! mHe ) Fatal("Main", "FairMedium helium not found"); geoBuild->createMedium(mHe); TGeoMedium* pMed4 = gGeoMan->GetMedium("helium"); if ( ! pMed4 ) Fatal("Main", "Medium helium not found"); FairGeoMedium* mDCH = geoMedia->getMedium("DCHgas"); if ( ! mDCH ) Fatal("Main", "FairMedium DCHgas not found"); geoBuild->createMedium(mDCH); TGeoMedium* pMed33 = gGeoMan->GetMedium("DCHgas"); if ( ! pMed33 ) Fatal("Main", "Medium DCHgas not found"); // -------------------------------------------------------------------------- // -------------- Create geometry and top volume ------------------------- gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom"); gGeoMan->SetName("DCHgeom"); TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); // -------------------------------------------------------------------------- Double_t dx,dy,dz; Double_t theta, phi; Double_t thx, phx, thy, phy, thz, phz; Double_t alpha; // Define DCH Geometry TGeoVolume* topDCH = new TGeoVolumeAssembly("DCH"); //Active volume (envelope box enclosing the sense wires) Double_t actGasDx = 99.07/2.; // [cm] //99.07326 Double_t actGasDy = 77.60/2.; // [cm] //77.59584 Double_t actGasDz = 5.20/2.; // [cm] //distance of x and y planes centres is 4 cm. Distance between the first x halfplane and the second y halfplane is 4cm + 2*0.69*cos30 cm. TGeoShape* pActGasBox = new TGeoBBox("ActGasBox", actGasDx, actGasDy, actGasDz); TGeoVolume* pActGasDchLog = new TGeoVolume("ActGASBoxLog",pActGasBox, pMed33); pActGasDchLog->SetVisLeaves(kTRUE); pActGasDchLog->SetVisContainers(kTRUE); // Gas box Double_t gasDx = 106.4/2.; // [cm] Double_t gasDy = 83.4/2.; // [cm] Double_t gasDz = 7.50; // [cm] TGeoShape* pGasBox = new TGeoBBox("GasBox", gasDx, gasDy, gasDz); TGeoVolume* pGasDchLog = new TGeoVolume("GASBoxLog",pGasBox, pMed33); pGasDchLog->SetVisLeaves(kTRUE); pGasDchLog->SetVisContainers(kTRUE); // Al Frame Double_t alDx = 125.8/2.; // [cm] Double_t alDy = 103.4/2.; // [cm] //Double_t alDz = 4.06; // [cm] Double_t alDz = 7.50; // [cm] //modified TGeoShape* pAlBox = new TGeoBBox("AlBox", alDx, alDy, alDz); TGeoVolume* pAlDchLog = new TGeoVolume("ALBoxLog",pAlBox, pMed21); pAlDchLog->SetVisLeaves(kTRUE); // Mylar Entrance exit windows Double_t mylDx= gasDx; //[cm] Double_t mylDy= gasDy; //[cm] Double_t mylDz= 0.0006; //[cm] TGeoShape* pMylarBox = new TGeoBBox("MylarBox", mylDx, mylDy, mylDz); TGeoVolume* pMylDchLog = new TGeoVolume("MYLBoxLog",pMylarBox, pMed15); pMylDchLog->SetVisLeaves(kTRUE); // First assembly TGeoVolume *dch1 = new TGeoVolumeAssembly("DCH1"); TGeoRotation *rot = new TGeoRotation(); rot->RotateX(0.); rot->RotateY(0.); rot->RotateZ(0.); //Double_t tx = -3.5; //correct values? active-area vs Al frame offsets? //Double_t ty = -5.; Double_t tx = +3.5; //corrected values, active-area vs Al frame offsets. F.Wamers. Double_t ty = -3.5; Double_t tz = 0.; TGeoCombiTrans* pTransfo1 = new TGeoCombiTrans("", 0.,0.,0.,rot); TGeoCombiTrans* pTransfo2 = new TGeoCombiTrans("", tx,ty,tz,rot); TGeoCombiTrans* pTransfo3 = new TGeoCombiTrans("", -tx,-ty,-tz,rot); //new, in order to compensate Al vs Gas dch1->AddNode(pAlDchLog,0,pTransfo3); //1->3 pAlDchLog->AddNode(pGasDchLog,0,pTransfo2); //place gas in aluminum, shifted pGasDchLog->AddNode(pActGasDchLog, 0,pTransfo1); //place active in gas, centrally // Mylar Windows front+back //dch1->AddNode(pMylDchLog,0,new TGeoCombiTrans("", tx,ty,-alDz-mylDz,rot)); //dch1->AddNode(pMylDchLog,1,new TGeoCombiTrans("", tx,ty, alDz+mylDz,rot)); dch1->AddNode(pMylDchLog,0,new TGeoCombiTrans("", 0.,0.,-alDz-mylDz,rot)); dch1->AddNode(pMylDchLog,1,new TGeoCombiTrans("", 0.,0., alDz+mylDz,rot)); // Global Positioning //in agreement with the s318 tracker, those are supposed to be the centres of the active volumes!!! //Double_t pDch1x = -123.22 ; //Justyna //Double_t pDch1y = 3.6 ; //Justyna //Double_t pDch1z = 444.13 ; //Justyna //Double_t pDch2x = -167.0 ; //Justyna //Double_t pDch2y = 1.02 ; //Justyna //Double_t pDch2z = 535.1 ; //Justyna //use identical values as for the tracker config and the digitizer. F. Wamers. Double_t pDch1x = -123.219446 ; //Felix Double_t pDch1y = 3.597104 ; //Felix Double_t pDch1z = 444.126271 ; //Felix Double_t pDch2x = -167.015888 ; //Felix Double_t pDch2y = 1.016917 ; //Felix Double_t pDch2z = 535.093884 ; //Felix // Double_t pDch1x = -132.233355 ; //Christoph // Double_t pDch1y = 1.037475 ; //Christoph // Double_t pDch1z = 438.710168 ; //Christoph // Double_t pDch2x = -170.8653 ; //Christoph // Double_t pDch2y = 2.075139 ; //Christoph // Double_t pDch2z = 538.614091 ; //Christoph //The order of rotation matters!!! Rotate first z, and then y! Felix Double_t aDch1 = -31.0 ; Double_t aDch2 = -31.0 ; TGeoRotation *gRot1 = new TGeoRotation(); gRot1->RotateX(0.); gRot1->RotateZ(-8.880000); //Justyna //gRot1->RotateZ(+8.880000); //Felix // gRot1->RotateZ(-2.5); //Christoph gRot1->RotateY(aDch1); TGeoRotation *gRot2 = new TGeoRotation(); gRot2->RotateX(0.); gRot2->RotateZ(9.350000); //Justyna //gRot2->RotateZ(-9.350000); //Felix // gRot2->RotateZ(8.4); //Christoph gRot2->RotateY(aDch2); // Helium Bag definition Double_t heDx= alDx ; //[cm] Double_t heDy= alDy ; //[cm] Double_t heDz=(pDch2z-pDch1z)*0.953874/2.; //[cm] alpha=0.; //[degre] Double_t beta =0.; //[degre] phi =15.20; //[degre] TGeoShape* pHePara = new TGeoPara("HePara", heDx, heDy, heDz, alpha,beta,phi); TGeoVolume* pHeDchLog = new TGeoVolume("HeParaLog",pHePara, pMed4); pHeDchLog->SetVisLeaves(kTRUE); topDCH->AddNode(dch1,0,new TGeoCombiTrans("",pDch1x,pDch1y,pDch1z,gRot1) ); topDCH->AddNode(dch1,1,new TGeoCombiTrans("",pDch2x,pDch2y,pDch2z,gRot2) ); topDCH->AddNode(pHeDchLog,0,new TGeoCombiTrans("",(pDch1x+pDch2x)/2., pDch2y, (pDch1z+pDch2z)/2., gRot1) ); TGeoCombiTrans *temp1 = new TGeoCombiTrans(); top->AddNode(topDCH, 0, GetGlobalPosition(temp1)); // --------------- Finish ----------------------------------------------- gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* geoFile = new TFile(geoFileName, "RECREATE"); top->Write(); geoFile->Close(); // -------------------------------------------------------------------------- }
/** * Initialise ROOT geometry objects from GEAR objects * * @param geomName name of ROOT geometry object * @param dumpRoot dump automatically generated ROOT geometry file for further inspection */ void EUTelGeometryTelescopeGeoDescription::initializeTGeoDescription( std::string const & geomName, bool dumpRoot = false ) { if( _isGeoInitialized ) { streamlog_out( WARNING3 ) << "EUTelGeometryTelescopeGeoDescription: Geometry already initialized, using old initialization" << std::endl; return; } else { _geoManager = std::make_unique<TGeoManager>("Telescope", "v0.1"); _geoManager->SetBit(kCanDelete); } if( !_geoManager ) { streamlog_out( ERROR3 ) << "Can't instantiate ROOT TGeoManager " << std::endl; return; } // Create top world volume containing telescope geometry // Create air mixture // see http://pdg.lbl.gov/2013/AtomicNuclearProperties/HTML_PAGES/104.html double air_density = 1.2e-3; // g/cm^3 double air_radlen = 36.62; // g/cm^2 //Must be -ve to use this value rather than internal one calculated. TGeoMixture* pMatAir = new TGeoMixture("AIR",3,air_density); pMatAir->DefineElement(0, 14.007, 7., 0.755267 ); //Nitrogen pMatAir->DefineElement(1, 15.999, 8., 0.231781 ); //Oxygen pMatAir->DefineElement(2, 39.948, 18., 0.012827 ); //Argon pMatAir->DefineElement(3, 12.011, 6., 0.000124 ); //Carbon pMatAir->SetRadLen( air_radlen ); // Medium: medium_World_AIR TGeoMedium* pMedAir = new TGeoMedium("medium_World_AIR", 3, pMatAir ); // The World is the 10 x 10m x 10m box filled with air mixture Double_t dx,dy,dz; dx = 5000.000000; // [mm] dy = 5000.000000; // [mm] dz = 5000.000000; // [mm] TGeoShape *pBoxWorld = new TGeoBBox("Box_World", dx,dy,dz); // Volume: volume_World TGeoVolume* pvolumeWorld = new TGeoVolume("volume_World",pBoxWorld, pMedAir); pvolumeWorld->SetLineColor(4); pvolumeWorld->SetLineWidth(3); pvolumeWorld->SetVisLeaves(kTRUE); // Set top volume of geometry gGeoManager->SetTopVolume( pvolumeWorld ); IntVec::const_iterator itrPlaneId; for ( itrPlaneId = _sensorIDVec.begin(); itrPlaneId != _sensorIDVec.end(); ++itrPlaneId ) { translateSiPlane2TGeo(pvolumeWorld, *itrPlaneId ); } _geoManager->CloseGeometry(); _isGeoInitialized = true; // Dump ROOT TGeo object into file if ( dumpRoot ) _geoManager->Export( geomName.c_str() ); for(auto& mapEntry: _planePath) { auto pathName = mapEntry.second; auto sensorID = mapEntry.first; _geoManager->cd( pathName.c_str() ); _TGeoMatrixMap[sensorID] = _geoManager->GetCurrentNode()->GetMatrix(); } return; }
void EUTelGeometryTelescopeGeoDescription::translateSiPlane2TGeo(TGeoVolume* pvolumeWorld, int SensorId ) { double xc, yc, zc; // volume center position double alpha, beta, gamma; double rotRef1, rotRef2, rotRef3, rotRef4; std::stringstream strId; strId << SensorId; // Get sensor center position xc = siPlaneXPosition( SensorId ); yc = siPlaneYPosition( SensorId ); zc = siPlaneZPosition( SensorId ); // Get sensor orientation alpha = siPlaneXRotation( SensorId ); // in degrees ! beta = siPlaneYRotation( SensorId ); // gamma = siPlaneZRotation( SensorId ); // rotRef1 = siPlaneRotation1( SensorId ); rotRef2 = siPlaneRotation2( SensorId ); rotRef3 = siPlaneRotation3( SensorId ); rotRef4 = siPlaneRotation4( SensorId ); //We must check that the input is correct. Since this is a combination of initial rotations and reflections the determinate must be 1 or -1 float determinant = rotRef1*rotRef4 - rotRef2*rotRef3 ; if(determinant==1 or determinant==-1) { streamlog_out(DEBUG5) << "SensorID: " << SensorId << ". Determinant = " <<determinant <<" This is the correct determinate for this transformation." << std::endl; } else { streamlog_out(ERROR5) << "SensorID: " << SensorId << ". Determinant = " <<determinant << std::endl; throw(lcio::Exception("The initial rotation and reflection matrix does not have determinant of 1 or -1. Gear file input must be wrong.")); } //Create spatial TGeoTranslation object. std::string stTranslationName = "matrixTranslationSensor"; stTranslationName.append( strId.str() ); TGeoTranslation* pMatrixTrans = new TGeoTranslation( stTranslationName.c_str(), xc, yc, zc ); //ALL clsses deriving from TGeoMatrix are not owned by the ROOT geometry manager, invoking RegisterYourself() transfers //ownership and thus ROOT will clean up pMatrixTrans->RegisterYourself(); //Create TGeoRotation object. //Translations are of course just positional changes in the global frame. //Note that each subsequent rotation is using the new coordinate system of the last transformation all the way back to the global frame. //The way to think about this is that each rotation is the multiplication of the last rotation matrix by a new one. //The order is: //Integer Z rotation and reflections. //Z rotations specified by in degrees. //X rotations //Y rotations TGeoRotation* pMatrixRotRefCombined = new TGeoRotation(); //We have to ensure that we retain a right handed coordinate system, i.e. if we only flip the x or y axis, we have to also flip the z-axis. If we flip both we have to flip twice. double integerRotationsAndReflections[9]={rotRef1,rotRef2,0,rotRef3,rotRef4,0,0,0, determinant}; pMatrixRotRefCombined->SetMatrix(integerRotationsAndReflections); std::cout << "Rotating plane " << SensorId << " to gamma: " << gamma << std::endl; pMatrixRotRefCombined->RotateZ(gamma);//Z Rotation (degrees)//This will again rotate a vector around z axis usign the right hand rule. pMatrixRotRefCombined->RotateX(alpha);//X Rotations (degrees)//This will rotate a vector usign the right hand rule round the x-axis pMatrixRotRefCombined->RotateY(beta);//Y Rotations (degrees)//Same again for Y axis pMatrixRotRefCombined->RegisterYourself();//We must allow the matrix to be used by the TGeo manager. // Combined translation and orientation TGeoCombiTrans* combi = new TGeoCombiTrans( *pMatrixTrans, *pMatrixRotRefCombined ); //This is to print to screen the rotation and translation matrices used to transform from local to global frame. streamlog_out(MESSAGE9) << "THESE MATRICES ARE USED TO TAKE A POINT IN THE LOCAL FRAME AND MOVE IT TO THE GLOBAL FRAME." << std::endl; streamlog_out(MESSAGE9) << "SensorID: " << SensorId << " Rotation/Reflection matrix for this object." << std::endl; const double* rotationMatrix = combi->GetRotationMatrix(); streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[0]<<" "<<rotationMatrix[1]<<" "<<rotationMatrix[2]<< std::endl; streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[3]<<" "<<rotationMatrix[4]<<" "<<rotationMatrix[5]<< std::endl; streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[6]<<" "<<rotationMatrix[7]<<" "<<rotationMatrix[8]<< std::endl; //streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[0] << std::setw(10) <<rotationMatrix[1]<< std::setw(10) <<rotationMatrix[2]<< std::setw(10)<< std::endl<< std::endl; //streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[3] << std::setw(10) <<rotationMatrix[4]<< std::setw(10) <<rotationMatrix[5]<< std::setw(10)<< std::endl<< std::endl; //streamlog_out(MESSAGE9) << std::setw(10) <<rotationMatrix[6] << std::setw(10) <<rotationMatrix[7]<< std::setw(10) <<rotationMatrix[8]<< std::setw(10)<< std::endl<< std::endl; const double* translationMatrix = combi->GetTranslation(); streamlog_out(MESSAGE9) << "SensorID: " << SensorId << " Translation vector for this object." << std::endl; streamlog_out(MESSAGE9) << std::setw(10) <<translationMatrix[0] << std::setw(10) <<translationMatrix[1]<< std::setw(10) <<translationMatrix[2]<< std::setw(10)<< std::endl; combi->RegisterYourself(); // Construct object medium. Required for radiation length determination // assume SILICON, though all information except of radiation length is ignored double a = 28.085500; double z = 14.000000; double density = 2.330000; double radl = siPlaneRadLength( SensorId ); double absl = 45.753206; std::string stMatName = "materialSensor"; stMatName.append( strId.str() ); TGeoMaterial* pMat = new TGeoMaterial( stMatName.c_str(), a, z, density, -radl, absl ); pMat->SetIndex( 1 ); // Medium: medium_Sensor_SILICON int numed = 0; // medium number double par[8]; par[0] = 0.000000; // isvol par[1] = 0.000000; // ifield par[2] = 0.000000; // fieldm par[3] = 0.000000; // tmaxfd par[4] = 0.000000; // stemax par[5] = 0.000000; // deemax par[6] = 0.000000; // epsil par[7] = 0.000000; // stmin std::string stMedName = "mediumSensor"; stMedName.append( strId.str() ); TGeoMedium* pMed = new TGeoMedium( stMedName.c_str(), numed, pMat, par ); // Construct object shape // Shape: Box type: TGeoBBox // TGeo requires half-width of box side Double_t dx = siPlaneXSize( SensorId ) / 2.; Double_t dy = siPlaneYSize( SensorId ) / 2.; Double_t dz = siPlaneZSize( SensorId ) / 2.; TGeoShape *pBoxSensor = new TGeoBBox( "BoxSensor", dx, dy, dz ); std::cout << "Box for sensor: " << SensorId << " is: " << dx << "|" << dy << "|" << dz << '\n'; // Geometry navigation package requires following names for objects that have an ID name:ID std::string stVolName = "volume_SensorID:"; stVolName.append( strId.str() ); _planePath.insert( std::make_pair(SensorId, "/volume_World_1/"+stVolName+"_1") ); TGeoVolume* pvolumeSensor = new TGeoVolume( stVolName.c_str(), pBoxSensor, pMed ); pvolumeSensor->SetVisLeaves( kTRUE ); pvolumeWorld->AddNode(pvolumeSensor, 1/*(SensorId)*/, combi); //this line tells the pixel geometry manager to load the pixel geometry into the plane streamlog_out(DEBUG1) << " sensorID: " << SensorId << " " << stVolName << std::endl; std::string name = geoLibName(SensorId); if( name == "CAST" ) { _pixGeoMgr->addCastedPlane( SensorId, siPlaneXNpixels(SensorId), siPlaneYNpixels(SensorId), siPlaneXSize(SensorId), siPlaneYSize(SensorId), siPlaneZSize(SensorId), siPlaneRadLength(SensorId), stVolName); } else { _pixGeoMgr->addPlane( SensorId, name, stVolName); updatePlaneInfo(SensorId); } }
/** * Initialise ROOT geometry objects from GEAR objects * * @param geomName name of ROOT geometry object * @param dumpRoot dump automatically generated ROOT geometry file for further inspection */ void EUTelGeometryTelescopeGeoDescription::initializeTGeoDescription( std::string& geomName, bool dumpRoot = false ) { // #ifdef USE_TGEO // get access to ROOT's geometry manager if( _isGeoInitialized ) { streamlog_out( WARNING3 ) << "EUTelGeometryTelescopeGeoDescription: Geometry already initialized, using old initialization" << std::endl; return; } else { _geoManager = new TGeoManager("Telescope", "v0.1"); } if( !_geoManager ) { streamlog_out( ERROR3 ) << "Can't instantiate ROOT TGeoManager " << std::endl; return; } // Create top world volume containing telescope/DUT geometry // Create air mixture // see http://pdg.lbl.gov/2013/AtomicNuclearProperties/HTML_PAGES/104.html double air_density = 1.2e-3; // g/cm^3 double air_radlen = 36.62; // g/cm^2 TGeoMixture* pMatAir = new TGeoMixture("AIR",3,air_density); pMatAir->DefineElement(0, 14.007, 7., 0.755267 ); //Nitrogen pMatAir->DefineElement(1, 15.999, 8., 0.231781 ); //Oxygen pMatAir->DefineElement(2, 39.948, 18., 0.012827 ); //Argon pMatAir->DefineElement(3, 12.011, 6., 0.000124 ); //Carbon pMatAir->SetRadLen( air_radlen ); // Medium: medium_World_AIR TGeoMedium* pMedAir = new TGeoMedium("medium_World_AIR", 3, pMatAir ); // The World is the 10 x 10m x 10m box filled with air mixture Double_t dx,dy,dz; dx = 5000.000000; // [mm] dy = 5000.000000; // [mm] dz = 5000.000000; // [mm] TGeoShape *pBoxWorld = new TGeoBBox("Box_World", dx,dy,dz); // Volume: volume_World TGeoVolume* pvolumeWorld = new TGeoVolume("volume_World",pBoxWorld, pMedAir); pvolumeWorld->SetLineColor(4); pvolumeWorld->SetLineWidth(3); pvolumeWorld->SetVisLeaves(kTRUE); // Set top volume of geometry gGeoManager->SetTopVolume( pvolumeWorld ); // Iterate over registered GEAR objects and construct their TGeo representation const Double_t PI = 3.141592653589793; const Double_t DEG = 180./PI; double xc, yc, zc; // volume center position double alpha, beta, gamma; IntVec::const_iterator itrPlaneId; for ( itrPlaneId = _sensorIDVec.begin(); itrPlaneId != _sensorIDVec.end(); ++itrPlaneId ) { std::stringstream strId; strId << *itrPlaneId; // Get sensor center position xc = siPlaneXPosition( *itrPlaneId ); yc = siPlaneYPosition( *itrPlaneId ); zc = siPlaneZPosition( *itrPlaneId ); // Get sensor orientation alpha = siPlaneXRotation( *itrPlaneId ); // [rad] beta = siPlaneYRotation( *itrPlaneId ); // [rad] gamma = siPlaneZRotation( *itrPlaneId ); // [rad] // Spatial translations of the sensor center string stTranslationName = "matrixTranslationSensor"; stTranslationName.append( strId.str() ); TGeoTranslation* pMatrixTrans = new TGeoTranslation( stTranslationName.c_str(), xc, yc, zc ); //ALL clsses deriving from TGeoMatrix are not owned by the ROOT geometry manager, invoking RegisterYourself() transfers //ownership and thus ROOT will clean up pMatrixTrans->RegisterYourself(); // Spatial rotation around sensor center // TGeoRotation requires Euler angles in degrees string stRotationName = "matrixRotationSensorX"; stRotationName.append( strId.str() ); TGeoRotation* pMatrixRotX = new TGeoRotation( stRotationName.c_str(), 0., alpha*DEG, 0.); // around X axis stRotationName = "matrixRotationSensorY"; stRotationName.append( strId.str() ); TGeoRotation* pMatrixRotY = new TGeoRotation( stRotationName.c_str(), 90., beta*DEG, 0.); // around Y axis (combination of rotation around Z axis and new X axis) stRotationName = "matrixRotationSensorBackY"; stRotationName.append( strId.str() ); TGeoRotation* pMatrixRotY1 = new TGeoRotation( stRotationName.c_str(), -90., 0., 0.); // restoration of original orientation (valid in small angle approximataion ~< 15 deg) stRotationName = "matrixRotationSensorZ"; stRotationName.append( strId.str() ); TGeoRotation* pMatrixRotZ = new TGeoRotation( stRotationName.c_str(), 0. , 0., gamma*DEG); // around Z axis // Combined rotation in several steps TGeoRotation* pMatrixRot = new TGeoRotation( *pMatrixRotX ); pMatrixRot->MultiplyBy( pMatrixRotY ); pMatrixRot->MultiplyBy( pMatrixRotY1 ); pMatrixRot->MultiplyBy( pMatrixRotZ ); pMatrixRot->RegisterYourself(); pMatrixRotX->RegisterYourself(); pMatrixRotY->RegisterYourself(); pMatrixRotY1->RegisterYourself(); pMatrixRotZ->RegisterYourself(); // Combined translation and orientation TGeoCombiTrans* combi = new TGeoCombiTrans( *pMatrixTrans, *pMatrixRot ); combi->RegisterYourself(); // Construction of sensor objects // Construct object medium. Required for radiation length determination // assume SILICON, though all information except of radiation length is ignored double a = 28.085500; double z = 14.000000; double density = 2.330000; double radl = siPlaneMediumRadLen( *itrPlaneId ); double absl = 45.753206; string stMatName = "materialSensor"; stMatName.append( strId.str() ); TGeoMaterial* pMat = new TGeoMaterial( stMatName.c_str(), a, z, density, radl, absl ); pMat->SetIndex( 1 ); // Medium: medium_Sensor_SILICON int numed = 0; // medium number double par[8]; par[0] = 0.000000; // isvol par[1] = 0.000000; // ifield par[2] = 0.000000; // fieldm par[3] = 0.000000; // tmaxfd par[4] = 0.000000; // stemax par[5] = 0.000000; // deemax par[6] = 0.000000; // epsil par[7] = 0.000000; // stmin string stMedName = "mediumSensor"; stMedName.append( strId.str() ); TGeoMedium* pMed = new TGeoMedium( stMedName.c_str(), numed, pMat, par ); // Construct object shape // Shape: Box type: TGeoBBox // TGeo requires half-width of box side dx = siPlaneXSize( *itrPlaneId ) / 2.; dy = siPlaneYSize( *itrPlaneId ) / 2.; dz = siPlaneZSize( *itrPlaneId ) / 2.; TGeoShape *pBoxSensor = new TGeoBBox( "BoxSensor", dx, dy, dz ); // Volume: volume_Sensor1 // Geometry navigation package requires following names for objects that have an ID // name:ID string stVolName = "volume_SensorID:"; stVolName.append( strId.str() ); _planePath.insert( std::make_pair(*itrPlaneId, "/volume_World_1/"+stVolName+"_1") ); TGeoVolume* pvolumeSensor = new TGeoVolume( stVolName.c_str(), pBoxSensor, pMed ); pvolumeSensor->SetVisLeaves( kTRUE ); pvolumeWorld->AddNode(pvolumeSensor, 1/*(*itrPlaneId)*/, combi); //this line tells the pixel geometry manager to load the pixel geometry into the plane _pixGeoMgr->addPlane( *itrPlaneId, geoLibName( *itrPlaneId), stVolName); } // loop over sensorID _geoManager->CloseGeometry(); _isGeoInitialized = true; // Dump ROOT TGeo object into file if ( dumpRoot ) _geoManager->Export( geomName.c_str() ); // #endif //USE_TGEO return; }
void create_sfi_geo(const char* geoTag) { //fGlobalTrans->SetTranslation(0.0,0.0,0.0); // ------- Load media from media file ----------------------------------- FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader"); FairGeoInterface* geoFace = geoLoad->getGeoInterface(); TString geoPath = gSystem->Getenv("VMCWORKDIR"); TString medFile = geoPath + "/geometry/media_r3b.geo"; geoFace->setMediaFile(medFile); geoFace->readMedia(); gGeoMan = gGeoManager; // -------------------------------------------------------------------------- // ------- Geometry file name (output) ---------------------------------- TString geoFileName = geoPath + "/geometry/sfi_"; geoFileName = geoFileName + geoTag + ".geo.root"; // -------------------------------------------------------------------------- // ----------------- Get and create the required media ----------------- FairGeoMedia* geoMedia = geoFace->getMedia(); FairGeoBuilder* geoBuild = geoLoad->getGeoBuilder(); FairGeoMedium* mAir = geoMedia->getMedium("Air"); if ( ! mAir ) Fatal("Main", "FairMedium Air not found"); geoBuild->createMedium(mAir); TGeoMedium* pMed2 = gGeoMan->GetMedium("Air"); if ( ! pMed2 ) Fatal("Main", "Medium Air not found"); FairGeoMedium* mVac = geoMedia->getMedium("vacuum"); if ( ! mVac ) Fatal("Main", "FairMedium vacuum not found"); geoBuild->createMedium(mVac); TGeoMedium* pMed1 = gGeoMan->GetMedium("vacuum"); if ( ! pMed1 ) Fatal("Main", "Medium vacuum not found"); FairGeoMedium* mGfi = geoMedia->getMedium("plasticForGFI"); if ( ! mGfi ) Fatal("Main", "FairMedium plasticForGFI not found"); geoBuild->createMedium(mGfi); TGeoMedium* pMed35 = gGeoMan->GetMedium("plasticForGFI"); if ( ! pMed35 ) Fatal("Main", "Medium plasticForGFI not found"); FairGeoMedium* mAl = geoMedia->getMedium("aluminium"); if ( ! mAl ) Fatal("Main", "FairMedium aluminium not found"); geoBuild->createMedium(mAl); TGeoMedium* pMed21 = gGeoMan->GetMedium("aluminium"); if ( ! pMed21 ) Fatal("Main", "Medium aluminium not found"); // -------------------------------------------------------------------------- // -------------- Create geometry and top volume ------------------------- gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom"); gGeoMan->SetName("GFIgeom"); TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); // -------------------------------------------------------------------------- //LABPOS(GFI1,-73.274339,0.069976,513.649524) Float_t dx = -73.274339; //dE tracker, correction due to wrong angle Float_t dy = 0.069976; Float_t dz = 513.649524; TGeoRotation *pMatrix3 = new TGeoRotation(); //pMatrix3->RotateY(-16.7); TGeoCombiTrans* pMatrix2 = new TGeoCombiTrans("", dx,dy,dz,pMatrix3); //LABPOS(GFI2,-147.135037,0.069976,729.680342) dx = -147.135037; //dE tracker, correction due to wrong angle dy = 0.069976; dz = 729.680342; TGeoRotation *pMatrix5 = new TGeoRotation(); //pMatrix5->RotateY(-16.7); TGeoCombiTrans* pMatrix4 = new TGeoCombiTrans("", dx,dy,dz,pMatrix5); // World definition TGeoVolume* pWorld = gGeoManager->GetTopVolume(); pWorld->SetVisLeaves(kTRUE); // SHAPES, VOLUMES AND GEOMETRICAL HIERARCHY // Volume: GFILogWorld TGeoVolume* pGFILogWorld = new TGeoVolumeAssembly("GFILogWorld"); pGFILogWorld->SetVisLeaves(kTRUE); // Global positioning pWorld->AddNode( pGFILogWorld, 0, pMatrix2 ); Float_t detector_size = 5.120000; Float_t fiber_thickness = 0.020000; TGeoShape *pGFITube = new TGeoBBox("GFITube", fiber_thickness/2,detector_size/2,fiber_thickness/2); TGeoVolume* pGFILog = new TGeoVolume("SFILog",pGFITube, pMed35); TGeoShape *pGFITubeActive = new TGeoBBox("GFITubeActive", (fiber_thickness * .98)/2, detector_size/2-0.0001, (fiber_thickness * .98)/2); TGeoVolume* pGFILogActive = new TGeoVolume("SFI1Log",pGFITubeActive,pMed35); pGFILog->SetLineColor((Color_t) 1); pGFILog->SetVisLeaves(kTRUE); TGeoRotation *pMatrixTube = new TGeoRotation(); pMatrixTube->RotateZ(90); pGFILog -> AddNode(pGFILogActive, 0, new TGeoCombiTrans()); for(int fiber_id = 0; fiber_id < detector_size / fiber_thickness; fiber_id++) { pGFILogWorld->AddNode(pGFILog, fiber_id, new TGeoCombiTrans("", -detector_size / 2 + (fiber_id + .5) * fiber_thickness, 0, 0, new TGeoRotation() ) ); pGFILogWorld->AddNode(pGFILog, fiber_id + detector_size / fiber_thickness, new TGeoCombiTrans("", 0, -detector_size / 2 + (fiber_id + .5) * fiber_thickness, fiber_thickness, pMatrixTube ) ); } // Add the sensitive part // AddSensitiveVolume(pGFILog); // fNbOfSensitiveVol+=1; // --------------- Finish ----------------------------------------------- gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* geoFile = new TFile(geoFileName, "RECREATE"); top->Write(); geoFile->Close(); // -------------------------------------------------------------------------- }