create_rootgeom_BD_run8.C

Go to the documentation of this file.

#include "TGeoCompositeShape.h"
#include "TGeoCone.h"
#include "TGeoManager.h"
#include "TGeoPara.h"
#include "TGeoTrd2.h"
#include "TGeoTube.h"
 
#include <iomanip>
#include <iostream>
 
R__ADD_INCLUDE_PATH($VMCWORKDIR)
#include "macro/run/bmnloadlibs.C"
//---------------------------
//---------------------------
 
TGeoManager* gGeoMan = NULL;
 
const Double_t strip_BD_Xpos = 0.0;
const Double_t strip_BD_Ypos = 0.0;
const Double_t strip_BD_Zpos = -2.5;   // cm
// const Double_t strip_BD_Zpos = -3.5; //cm
 
const Double_t tubeFe_Xpos = 0.0;
const Double_t tubeFe_Ypos = 0.0;
const Double_t tubeFe_Zpos = -8.45;   // cm
// const Double_t tubeFe_Zpos = -9.45; //cm
 
const Double_t tubeAl_Xpos = 0.0;
const Double_t tubeAl_Ypos = 0.0;
const Double_t tubeAl_Zpos = -14.3;   // cm
 
const Double_t roh_Xpos = 0.0;
const Double_t roh_Ypos = 0.0;
const Double_t roh_Zpos = -2.65;   // cm
 
const Double_t ZLength_BDActiveVolume = 15.;
const Double_t ZLength_tubeFe = 29.5;
const Double_t ZLength_tubeAl = 18.;
const Double_t ZLength_roh = 5.3;
 
const Double_t dphi1 = 360. / 40;
const Double_t rr = 4.95;
 
const Double_t MaxRadiusFe = 8.1;
const Double_t MinRadiusFe = 8.;
 
const Double_t MaxRadiusAl = 7.0;
const Double_t MinRadiusAl = 6.8;
 
const Double_t MaxRadiusRohacel = 8.0;
const Double_t MinRadiusRohacel = 7.0;
 
const Double_t Seg_tube_Min_dim = 10.7;
const Double_t Seg_tube_Max_dim = 50.0;
 
const Double_t Seg_tube_Max_phi = 270.0;
const Double_t Seg_tube_Min_phi = 90.0;
 
const Double_t Seg_tube_Thick = 1.0;
 
void create_rootgeom_BD_run8()
{
 
    // -------   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.geo";
    geoFace->setMediaFile(medFile);
    geoFace->readMedia();
    gGeoMan = gGeoManager;
    // --------------------------------------------------------------------------
 
    // -------   Geometry file name (output)   ----------------------------------
    const TString geoDetectorName = "BD";
    const TString geoDetectorVersion = "run8_v1";
    TString geoFileName = geoPath + "/geometry/" + geoDetectorName + "_" + geoDetectorVersion + ".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* pMedAir = gGeoMan->GetMedium("air");
    if (!pMedAir)
        Fatal("Main", "Medium air not found");
 
    FairGeoMedium* mBC408 = geoMedia->getMedium("BC408");
    if (!mBC408)
        Fatal("Main", "FairMedium BC408 not found");
    geoBuild->createMedium(mBC408);
    TGeoMedium* pMedBC408 = gGeoMan->GetMedium("BC408");
    if (!pMedBC408)
        Fatal("Main", "Medium BC408 not found");
 
    FairGeoMedium* miron = geoMedia->getMedium("iron");
    if (!miron)
        Fatal("Main", "FairMedium iron   not found");
    geoBuild->createMedium(miron);
    TGeoMedium* pMediron = gGeoMan->GetMedium("iron");
    if (!pMediron)
        Fatal("Main", "Medium iron  not found");
 
    FairGeoMedium* maluminium = geoMedia->getMedium("aluminium");
    if (!maluminium)
        Fatal("Main", "FairMedium aluminium   not found");
    geoBuild->createMedium(maluminium);
    TGeoMedium* pMedaluminium = gGeoMan->GetMedium("aluminium");
    if (!pMedaluminium)
        Fatal("Main", "Medium aluminium  not found");
 
    FairGeoMedium* mpolypropylene = geoMedia->getMedium("polypropylene");
    if (!mpolypropylene)
        Fatal("Main", "FairMedium polypropylene   not found");
    geoBuild->createMedium(mpolypropylene);
    TGeoMedium* pMedpolypropylene = gGeoMan->GetMedium("polypropylene");
    if (!pMedpolypropylene)
        Fatal("Main", "Medium polypropylene  not found");
 
    FairGeoMedium* mlead = geoMedia->getMedium("lead");
    if (!mlead)
        Fatal("Main", "FairMedium lead   not found");
    geoBuild->createMedium(mlead);
    TGeoMedium* pMedlead = gGeoMan->GetMedium("lead");
    if (!pMedlead)
        Fatal("Main", "Medium lead  not found");
 
    FairGeoMedium* mVac = geoMedia->getMedium("vacuum");
    if (!mVac)
        Fatal("Main", "FairMedium vacuum not found");
    geoBuild->createMedium(mVac);
    TGeoMedium* pMedVac = gGeoMan->GetMedium("vacuum");
    if (!pMedVac)
        Fatal("Main", "Medium vacuum not found");
 
    // --------------------------------------------------------------------------
 
    // --------------   Create geometry and top volume  -------------------------
    gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom");
    gGeoMan->SetName(geoDetectorName + "_geom");
    TGeoVolume* top = new TGeoVolumeAssembly("TOP");
    top->SetMedium(pMedAir);
    gGeoMan->SetTopVolume(top);
 
    // --------------------------------------------------------------------------
 
    TGeoTranslation* BD_Position = new TGeoTranslation("BD_Position", 0, 0, 0);
 
    // Define TOP Geometry
    TGeoVolume* BD = new TGeoVolumeAssembly("BD");
    BD->SetMedium(pMedVac);
 
    TGeoTranslation* strips_BD1Pos_trans =
        new TGeoTranslation("strips_BD1Pos_trans", strip_BD_Xpos, strip_BD_Ypos, strip_BD_Zpos);
    TGeoTranslation* tubefePos_trans = new TGeoTranslation("tubefePos_trans", tubeFe_Xpos, tubeFe_Ypos, tubeFe_Zpos);
    TGeoTranslation* tubealPos_trans = new TGeoTranslation("tubefePos_trans", tubeAl_Xpos, tubeAl_Ypos, tubeAl_Zpos);
    TGeoTranslation* Postubepb_trans = new TGeoTranslation("Postubepb_trans", 0.0, 0.0, -4);
    TGeoTranslation* Seg_tube_trans = new TGeoTranslation("Seg_tube_trans", 0, 0, 4.6);
 
    TGeoBBox* strip = new TGeoBBox("srtip", 0.7 / 2., 0.7 / 2., 15. / 2.);
    TGeoTube* tubefeContainerS = new TGeoTube("tubefeContainerS", MinRadiusFe / 2, MaxRadiusFe / 2, ZLength_tubeFe / 2);
    TGeoTube* tubealContainerS = new TGeoTube("tubealContainerS", MinRadiusAl / 2, MaxRadiusAl / 2, ZLength_tubeAl / 2);
    TGeoTube* tubepbContainerS = new TGeoTube("tubepbContainerS", 7.6 / 2, 8 / 2, 18.8 / 2);
    TGeoTubeSeg* Seg_tube = new TGeoTubeSeg("Seg_tube", Seg_tube_Min_dim / 2, Seg_tube_Max_dim / 2, Seg_tube_Thick / 2,
                                            Seg_tube_Min_phi, Seg_tube_Max_phi);
 
    Double_t phi = 90 * TMath::DegToRad(), delta_phi = 360. / 40., dphi = delta_phi * TMath::DegToRad();
 
    TGeoCombiTrans* combi[41];
 
    for (Int_t i = 1; i < 41; i++) {
        combi[i] = new TGeoCombiTrans(rr * TMath::Cos(dphi * (i + 0.5)), rr * TMath::Sin(dphi * (i + 0.5)),
                                      strip_BD_Zpos, new TGeoRotation("rot", dphi1 * (i + 0.5), 0., 0.));
    }
 
    TGeoVolume* BD1ActiveVolumeV1 = new TGeoVolume("BD1ActiveVolumeV1", strip);
    BD1ActiveVolumeV1->SetMedium(pMedBC408);
    BD1ActiveVolumeV1->SetLineColor(kYellow - 10);
 
    TGeoVolume* tubefeMat = new TGeoVolume("tubefeMat", tubefeContainerS);
    tubefeMat->SetMedium(pMediron);
    tubefeMat->SetLineColor(kOrange - 1);
 
    TGeoVolume* tubealMat = new TGeoVolume("tubealMat", tubealContainerS);
    tubealMat->SetMedium(pMedaluminium);
    tubealMat->SetLineColor(kCyan - 10);
 
    TGeoVolume* tubepb_prot = new TGeoVolume("tubepb_prot", tubepbContainerS);
    tubepb_prot->SetMedium(pMedlead);
    tubepb_prot->SetLineColor(kRed);
 
    TGeoVolume* seg_tube_prot = new TGeoVolume("seg_tube_prot", Seg_tube);
    seg_tube_prot->SetMedium(pMedlead);
    seg_tube_prot->SetLineColor(kRed);
 
    for (Int_t i = 1; i < 41; i++) {
        BD->AddNode(BD1ActiveVolumeV1, i, combi[i]);
    }
 
    BD->AddNode(tubefeMat, 41, tubefePos_trans);
    BD->AddNode(tubealMat, 42, tubealPos_trans);
    BD->AddNode(tubepb_prot, 43, Postubepb_trans);
    BD->AddNode(seg_tube_prot, 44, Seg_tube_trans);
 
    // Adding volumes to the TOP Volume
    top->AddNode(BD, 1, BD_Position);
 
    top->SetVisContainers(kTRUE);
 
    // ---------------   Finish   -----------------------------------------------
    gGeoMan->CloseGeometry();
    gGeoMan->CheckOverlaps(0.001);
    gGeoMan->PrintOverlaps();
    gGeoMan->Test();
 
    TFile* geoFile = new TFile(geoFileName, "RECREATE");
    top->Write();
    geoFile->Close();
    top->Draw("ogl");
    // top->Draw("");
}