create_scit_rootgeo_201505.C File Reference

Go to the source code of this file.

Functions
int	create_scit_rootgeo_201505 (Double_t distance=750)

Function Documentation

int create_scit_rootgeo_201505

(

Double_t

distance = 750

)

Definition at line 1 of file create_scit_rootgeo_201505.C.

References cos(), Double_t, dx, dz, fi, geobuild, geoFace, geoLoad, gGeoMan, m, Media, nmed, outfile, phi0, pi, sin(), top, and TString.

                                                       {
  // Forward tof geometry parameters
  //-----------------------------
  //-- macro created by Alicia S. based on step file
  //-- barrel-SciTil_18122012
  //-- edited by D. Steinschaden 
  //-- last update 04.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_201505.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);
  

//-------------------- creating Tile and SiPm and Readout Card --------------------

  //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", 1.425, 0.25, 1.425);// 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);
  
  TGeoBBox* SiPmShape= new TGeoBBox("SiPm_logicSiPm",0.215,0.1925,0.0725);
   TGeoVolume* SiPmVol= new TGeoVolume("SiPm",SiPmShape, gGeoMan->GetMedium("silicon"));
  SiPmVol->SetLineColor(kRed);

  TGeoBBox* CardShape= new TGeoBBox("SciTil_logicCard",2.,0.075,2.);
   TGeoVolume* CardVol= new TGeoVolume("Card",CardShape, gGeoMan->GetMedium("pcbmvd"));
   CardVol->SetLineColor(kGreen);
  

  // --------------creating module (Quad)-------------
  
  TGeoVolumeAssembly* QuadAssVol = new TGeoVolumeAssembly("SciTil_quad");// 4 tiles + 8 SiPm
          

  TGeoRotation* rot_SiPm = new TGeoRotation;
  TGeoRotation* rot_SiPm2 = new TGeoRotation;
  QuadAssVol->AddNode(TileVol,0,new TGeoCombiTrans(1.5, 0.,1.5 , new TGeoRotation (0)));
  QuadAssVol->AddNode(SiPmVol,0,new TGeoCombiTrans(2.2125, 0.,0.0025 , new TGeoRotation (0)));
  rot_SiPm.SetAngles(90.,90.,90.);
  QuadAssVol->AddNode(SiPmVol,1,new TGeoCombiTrans(0.0025., 0.,0.7875 , rot_SiPm) );

  QuadAssVol->AddNode(TileVol,1,new TGeoCombiTrans(1.5, 0.,-1.5, new TGeoRotation (0)) );
  QuadAssVol->AddNode(SiPmVol,2,new TGeoCombiTrans(0.0025, 0.,-2.2125 , rot_SiPm) );
  QuadAssVol->AddNode(SiPmVol,3,new TGeoCombiTrans(0.7875, 0.,-0.0025 , new TGeoRotation (0)));

  QuadAssVol->AddNode(TileVol,2,new TGeoCombiTrans(-1.5, 0.,-1.5, new TGeoRotation (0)));
  QuadAssVol->AddNode(SiPmVol,4,new TGeoCombiTrans(-2.2125, 0.,-0.0025 ,new TGeoRotation (0)));
  QuadAssVol->AddNode(SiPmVol,5,new TGeoCombiTrans(-0.0025, 0.,-0.7875 , rot_SiPm));

  QuadAssVol->AddNode(TileVol,3,new TGeoCombiTrans(-1.5, 0.,1.5 , new TGeoRotation (0)) );
  QuadAssVol->AddNode(SiPmVol,6,new TGeoCombiTrans(-0.0025, 0.,2.2125 , rot_SiPm));
  QuadAssVol->AddNode(SiPmVol,7,new TGeoCombiTrans(-0.7875, 0.,0.0025 , new TGeoRotation (0)));

  QuadAssVol->AddNode(CardVol,1,new TGeoCombiTrans(0.,0.325,0., new TGeoRotation (0)));

  
  //-----------------creating Super-Module----------------

  TGeoVolumeAssembly* SModAssVol = new TGeoVolumeAssembly("SciTil_mod");//90 quads 
  Double_t dx, dz;

  dx=0.;
  dz=0.;
  for(int j=0;j<30;j++){
    for(int l=0;l<3;l++){
      dz = -90+3+j*6;
      dx = -9+3+l*6; 

      SModAssVol->AddNode(QuadAssVol,j*3+l,new TGeoCombiTrans(dx, 0.,dz, new TGeoRotation (0)));

    }
  }
  
  // -----------------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(SModAssVol, m+1, new TGeoCombiTrans(dx_bbox, dy_bbox, 0, new TGeoRotation(rot_bbox)));   
          }
        
  
          top->AddNode(BarrelAssVol,0, new TGeoCombiTrans(0.,0.,32.,new TGeoRotation()));
          //top->AddNode(QuadAssVol,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;
  }