create_scit_rootgeo_pSuperModule.C

Go to the documentation of this file.

int create_scit_rootgeo_pSuperModule(Double_t distance = 750){
  // Forward tof geometry parameters
  //-----------------------------
  //-- macro created by Alicia S. based on step file
  //-- barrel-SciTil_18122012
  //-- edited by D. Steinschaden 
  //-- last update 08.2015


 //--------------------------------------------------------------------
  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  
  // Load this libraries
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libPndData");
  gSystem->Load("libPassive");
  
  TString outfile= "../../../geometry/SciTil_pSuperModule.root";
  //TString outfile= "../../../geometry/SciTil_test.root";
  TFile* fi = new TFile(outfile,"RECREATE");  
  
  FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
  FairGeoInterface *geoFace = geoLoad->getGeoInterface();
  geoFace->setMediaFile("../../../geometry/media_pnd.geo");
  geoFace->readMedia();
  geoFace->print();

  FairGeoMedia *Media =  geoFace->getMedia();
  FairGeoBuilder *geobuild=geoLoad->getGeoBuilder();

  FairGeoMedium *MediumSci  = Media->getMedium("polyvinyltoluene"); 
  FairGeoMedium *MediumSiPm = Media->getMedium("silicon"); 
  FairGeoMedium *MediumCard = Media->getMedium("pcbmvd");

  const Double_t pi = 3.1415926535;
  Int_t nmed=geobuild->createMedium(MediumSci); // ----why is the return value saved?
   nmed=geobuild->createMedium(MediumSiPm);
  nmed=geobuild->createMedium(MediumCard);

   TGeoManager* gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom");
   TGeoVolume *top = new TGeoVolumeAssembly("SciT");
   
                    
  gGeoMan->SetTopVolume(top);
  

  //Volumes containing SENSOR in the name will be taken into account in the MC hit production
  //check PndSciT.h and PndSciT.cxx and the listofsensitives 

  TGeoBBox* TileShape;
  TileShape = new TGeoBBox("SciTil_logicTile", 9.0, 0.25, 90.0);// dimension is always given from the mid to the end e.g. only the half length

  TGeoVolume* TileVol;
  TileVol = new TGeoVolume("SciTil_SENSOR",TileShape, gGeoMan->GetMedium("polyvinyltoluene"));
  TileVol->SetLineColor(kBlue);


  // -----------------creating Barrel --------------------------------------------
  TGeoVolumeAssembly* BarrelAssVol = new TGeoVolumeAssembly("SciTil_barrel");
  
  //                                       ------------Values of Dirc Detector------------
  Double_t radius        =  50.55      ;  // 50. radius in middle of the barbox (x and y)
  Double_t hthick        =  1.7/2.;       // 1.7/2. half thickness of the bars
  Double_t barnum        =  6;            // 6 number of bars per barbox
  Int_t bbnum         =  16;      //16. total number of sides = barboxes
  Double_t bbGap         =  1.5;          //1.5 gap btw the neighboring barboxes (at the middle height)
  Double_t pipehAngle    =  3.6;          //3.6 [degrees] half of the angular space needed for the target pipe
   
  
  Double_t phi0 = (180.-2.*pipehAngle)/bbnum + pipehAngle;
  Double_t dphi = (180.-2.*pipehAngle)/bbnum*2.;  

//taking angles from dirc geometry 
          Double_t dx_bbox, dy_bbox, dz_bbox, phi_curr;   
             
          for(Int_t m = 0; m < bbnum; m ++){
            phi_curr = (90. - phi0 - dphi*m)/180.*pi;   
            if(m > bbnum/2-1){ phi_curr = (90. - phi0 - dphi*m - 2.*pipehAngle)/180.*pi; }
            dx_bbox = radius * cos(phi_curr);
            dy_bbox = radius * sin(phi_curr);
            
            TGeoRotation rot_bbox;   
            rot_bbox.RotateZ( -phi0 - m*dphi - (TMath::Floor(2.*m/bbnum))*(2.*pipehAngle));   
            BarrelAssVol->AddNode(TileVol, m, new TGeoCombiTrans(dx_bbox, dy_bbox, 0, new TGeoRotation(rot_bbox)));   
          }
        
  
          top->AddNode(BarrelAssVol,0, new TGeoCombiTrans(0.,0.,32.,new TGeoRotation()));
          //top->AddNode(SModAssVol,0, new TGeoCombiTrans(0.,0.,32.,new TGeoRotation()));
          //top->AddNode(LineAssVol,0, new TGeoCombiTrans(0.,0.,32.,new TGeoRotation()));

  gGeoMan->SetVisLevel(5);
          gGeoMan->CloseGeometry();
          top->Write();
          fi->Close();
          // gGeoManager->Export(outfile);
          top->Draw("ogl");
          
  return 0;
  }