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(); // -------------------------------------------------------------------------- }
void Create_TOF_Geometry_v13_3b() { // Load the necessary FairRoot libraries gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C"); basiclibs(); gSystem->Load("libGeoBase"); gSystem->Load("libParBase"); gSystem->Load("libBase"); // Load needed material definition from media.geo file create_materials_from_media_file(); // Get the GeoManager for later usage gGeoMan = (TGeoManager*) gROOT->FindObject("FAIRGeom"); gGeoMan->SetVisLevel(5); // 2 = super modules gGeoMan->SetVisOption(0); // Create the top volume /* TGeoBBox* topbox= new TGeoBBox("", 1000., 1000., 1000.); TGeoVolume* top = new TGeoVolume("top", topbox, gGeoMan->GetMedium("air")); gGeoMan->SetTopVolume(top); */ TGeoVolume* top = new TGeoVolumeAssembly("TOP"); gGeoMan->SetTopVolume(top); TGeoVolume* tof = new TGeoVolumeAssembly(geoVersion); top->AddNode(tof, 1); for(Int_t counterType = 0; counterType < NumberOfDifferentCounterTypes; counterType++) { gCounter[counterType] = create_new_counter(counterType); } for(Int_t moduleType = 0; moduleType < NofModuleTypes; moduleType++) { gModules[moduleType] = create_new_tof_module(moduleType); gModules[moduleType]->SetVisContainers(1); } gPole = create_tof_pole(); position_side_tof_modules(1); // keep order !! position_inner_tof_modules(3); cout << "Outer Types "<<Outer_Module_Types[0][0]<<", "<<Outer_Module_Types[1][0] <<", col=1: "<<Outer_Module_Types[0][1]<<", "<<Outer_Module_Types[1][1] <<endl; cout << "Outer Number "<<Outer_Module_Number[0][0]<<", "<<Outer_Module_Number[1][0] <<", col=1: "<<Outer_Module_Number[0][1]<<", "<<Outer_Module_Number[1][1] <<endl; position_outer_tof_modules(4); position_tof_poles(0); position_tof_bars(0); gGeoMan->CloseGeometry(); gGeoMan->CheckOverlaps(0.001); gGeoMan->PrintOverlaps(); gGeoMan->Test(); TFile* outfile = new TFile(FileName,"RECREATE"); top->Write(); //gGeoMan->Write(); outfile->Close(); top->SetVisContainers(1); gGeoMan->SetVisLevel(5); top->Draw("ogl"); //top->Draw(); //gModules[0]->Draw("ogl"); // gModules[0]->Draw(""); gModules[0]->SetVisContainers(1); // gModules[1]->Draw(""); gModules[1]->SetVisContainers(1); //gModules[5]->Draw(""); // top->Raytrace(); }