modifyRootFscGeometryFile.C

Go to the documentation of this file.

{
// Macro to create Fsc geometry file
// of the FSC detector for PANDA
// in root format
// with mechanical structure from CAD drawings
// fsc_cad2root_mech_support.root - obtained using CAD2RootConverter    
//03.10.2016 Dmitry Morozov

  // Fsc geometry parameters
  // in cm
  //-----------------------------
  Double_t     beamlineangle =  -TMath::DegToRad()*0.0; //in rad
  const Double_t     zpos          =  760.0;
  const Int_t        numcellsx     =  54;
  const Int_t        numcellsy     =  28;
  const Double_t     tyvekbundle   =  0.025;
  const Double_t     paperbundle   =  0.025;
  Double_t     crystalspace  =  tyvekbundle+paperbundle;
  const Double_t     Absorb_xsize  =  5.5;
  const Double_t     Absorb_ysize  =  5.5;
  const Double_t     Air_thickness =  0.0;
  const Double_t     Pb_thickness  =  0.0275;
  const Double_t     Sci_thickness =  0.15;
  const Double_t     Sci_xsize     =  5.5;
  const Double_t     Sci_ysize     =  5.5;
  const Int_t        EMCnumlayers  =  380;
  const Int_t        numfibers     =  6;
  const Double_t     fiber_r       =  0.06;
  const Double_t     fiber_hole_r  =  0.07;

  const Int_t        nClsSuperModuleX = 2;
  const Int_t        nClsSuperModuleY = 2;

  Double_t  cellxsize = TMath::Max(Absorb_xsize,Sci_xsize);
  Double_t  cellysize = TMath::Max(Absorb_ysize,Sci_ysize);
  Double_t  suprModXsize = 2.*cellxsize + 2.*crystalspace + Air_thickness;
  Double_t  suprModYsize = 2.*cellysize + 2.*crystalspace + Air_thickness;
  Double_t  xsize     = suprModXsize*numcellsx/nClsSuperModuleX;
  Double_t  ysize     = suprModYsize*numcellsy/nClsSuperModuleY;

  Double_t  cell_thickness = 2.0*Air_thickness+Pb_thickness+Sci_thickness;
  Double_t  thickness      = cell_thickness*EMCnumlayers + 2.*crystalspace;
  Double_t  holder_thickness = cell_thickness*EMCnumlayers;

  Double_t     hole_xsize    =  6.*cellxsize;
  Double_t     hole_ysize    =  6.*cellysize;
  //position of the hole means - coordinate of its center
  Double_t     hole_xpos     = 2.0*cellxsize; //Must be left to the beam line (check??)
  Double_t     hole_ypos     =  1.0*cellxsize;


  //--------------------------------------------------------------------
 //  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");

   // Load this example libraries
   gSystem->Load("libGeoBase");
   gSystem->Load("libParBase");
   gSystem->Load("libBase");
   gSystem->Load("libPndData");
   gSystem->Load("libPassive");
   TString outfile= "new_module5_fsc.root";
//   TString outfile= "emc_module5_fsc.root";
   TFile* fi = new TFile(outfile,"RECREATE");

   FairGeoLoader* geoLoad = new FairGeoLoader("TGeo","FairGeoLoader");
   FairGeoInterface *geoFace = geoLoad->getGeoInterface();
//   geoFace->setMediaFile("../../../geometry/media_pnd.geo");
   geoFace->setMediaFile("$VMCWORKDIR/geometry/media_pnd.geo");
//   geoFace->setMediaFile("/home/freez/fairsoft_may16p1/pandaroot/geometry/media_pnd.geo");
   geoFace->readMedia();
   geoFace->print();

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

   FairGeoMedium *MediumAir  = Media->getMedium("air");
   FairGeoMedium *MediumTyvek = Media->getMedium("tyvek");
   FairGeoMedium *MediumCarbon  = Media->getMedium("carbon");
   FairGeoMedium *MediumLead = Media->getMedium("lead");
   FairGeoMedium *MediumScint = Media->getMedium("FscScint");
   FairGeoMedium *MediumFiber = Media->getMedium("FscFiber");

   Int_t nmed=geobuild->createMedium(MediumAir);
   nmed=geobuild->createMedium(MediumTyvek);
   nmed=geobuild->createMedium(MediumCarbon);
   nmed=geobuild->createMedium(MediumLead);
   nmed=geobuild->createMedium(MediumScint);
   nmed=geobuild->createMedium(MediumFiber);

   TGeoManager* gGeoMan = (TGeoManager*)gROOT->FindObject("FAIRGeom");

   //building main container

   TGeoVolume *top = new TGeoVolumeAssembly("Emc5");

   gGeoMan->SetTopVolume(top);
//   TSrting name = "FscTopShape";
//
//   TGeoShape* TopShape = new TGeoBBox(name,xsize/2.,ysize/2.,thickness/2.);
//   TGeoVolume *TopVolume = new TGeoVolumeAssembly("FscBox");

 //  Building paper volume to place everything inside it
   TString name = "FscSuperModuleShape";
   TGeoShape* SuperModuleShape = new TGeoBBox(name,suprModXsize/2.,suprModYsize/2.,thickness/2.);
   name = "FscSuperModuleVolume";
   TString medium = "carbon";
   TGeoVolume* SuperModuleVolume = new TGeoVolume(name, SuperModuleShape, gGeoMan->GetMedium(medium));
   SuperModuleVolume->SetLineColor(2);

   //Building tyvek volume inside paper volume
   name = "FscTyvekShape";
   TGeoShape* TyvekShape = new TGeoBBox(name,suprModXsize/2. - paperbundle,suprModYsize/2. - paperbundle,thickness/2. - paperbundle);
   name = "FscTyvekVolume";
   medium = "tyvek";
   TGeoVolume* TyvekVolume = new TGeoVolume(name, TyvekShape, gGeoMan->GetMedium(medium));
   TyvekVolume->SetLineColor(3);
   SuperModuleVolume->AddNode(TyvekVolume,0);

   //Building module volume - 4 of them are inside tyvek volume
   name = "FscModuleShape";
   TGeoShape* ModuleShape = new TGeoBBox(name,cellxsize/2.,cellysize/2.,holder_thickness/2.);
   name = "FscModuleVolume";
   medium = "air";
   TGeoVolume* ModuleVolume = new TGeoVolume(name, ModuleShape, gGeoMan->GetMedium(medium));
   ModuleVolume->SetLineColor(2);
   TGeoCombiTrans* iCombTrans = NULL;
   TGeoRotation rot;
   Int_t iCopy = 0;
   Double_t x,y;

   for (int ix = 1; ix <= nClsSuperModuleX; ix++)
   {
       for (int iy = 1; iy <= nClsSuperModuleY; iy++)
       {
           x = -cellxsize + 2.*(cellxsize/(Double_t)nClsSuperModuleX)*(ix - 0.5);
           y = -cellxsize + 2.*(cellxsize/(Double_t)nClsSuperModuleY)*(iy - 0.5);
           cout<<"===Inside supermodule"<<endl;
           cout <<"x="<<x<<" y="<<y<<endl;
           cout <<"ix="<<ix<<" iy="<<iy<<endl;

           iCombTrans = new TGeoCombiTrans(x,y,0, new TGeoRotation(rot));
           name = "FscModuleVolume";
           name+=iCopy;
           iCombTrans->SetName(name);
           iCombTrans->RegisterYourself();
           TyvekVolume->AddNode(ModuleVolume,iCopy,iCombTrans);
           iCopy++;
       }
   }
   
   //Building and placing holder (for Scint+Lead layer) volumes inside module volume
   name = "FscLayerHolderShape";
   TGeoShape* LayerHolderShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, cell_thickness/2.);
   name = "FscLayerHolderVolume";
   medium = "Air";
 //  TGeoVolume* LayerHolderVolume = new TGeoVolume(name, LayerHolderShape, gGeoMan->GetMedium(medium));

   TGeoVolume* LayerHolderVolume = new TGeoVolumeAssembly(name);

 
      //Building and placing Lead + Scint layer inside holder Volume
      name = "FscLeadShape";
      TGeoShape* LeadShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, Pb_thickness/2.);
      name = "FscLeadVolume";
      medium = "lead";
      TGeoVolume* LeadVolume = new TGeoVolume(name, LeadShape, gGeoMan->GetMedium(medium));
      LeadVolume->SetLineColor(4);
      iCombTrans = new TGeoCombiTrans(0,0,-cell_thickness/2.+Pb_thickness/2., new TGeoRotation(rot));
      LayerHolderVolume->AddNode(LeadVolume,0,iCombTrans);

      name = "FscSciShape";
      TGeoShape* SciShape = new TGeoBBox(name, cellxsize/2. - crystalspace, cellysize/2. - crystalspace, Sci_thickness/2.);
      name = "FscSciVolume";
      medium = "FscScint";
      TGeoVolume* SciVolume = new TGeoVolume(name, SciShape, gGeoMan->GetMedium(medium));
      SciVolume->SetLineColor(5);
      iCombTrans = new TGeoCombiTrans(0,0,cell_thickness/2.-Sci_thickness/2., new TGeoRotation(rot));
      LayerHolderVolume->AddNode(SciVolume,0,iCombTrans);


   Double_t z;
//   TGeoCombiTrans* iCombTrans = NULL;
//   TGeoRotation rot;
   for(int iz = 1; iz <= EMCnumlayers; iz++)
   {
           z = -holder_thickness/2. + 0.5 * cell_thickness *(2.*iz - 1.);
           iCombTrans = new TGeoCombiTrans(0,0,z, new TGeoRotation(rot));
           name = "FscLayerHolder";
           name += iz-1;
           ModuleVolume->AddNode(LayerHolderVolume,iz-1,iCombTrans);
   }


   //Building and placing Fibers inside FscModule Volume
   //Holes for fibers Bottom right corner is (ix = 1, iy = 1)
   name = "FscFibHoleSHape";
//   TGeoShape* FiberHoleShape = new TGeoTube(name, 0., fiber_hole_r, (thickness-crystalspace)/2.);
   TGeoShape* FiberHoleShape = new TGeoTube(name, 0., fiber_hole_r, (holder_thickness)/2.);
   name = "FscFibHoleVolume";
   medium = "air";
   TGeoVolume* FiberHoleVolume = new TGeoVolume(name, FiberHoleShape, gGeoMan->GetMedium(medium));
   for(int ix = 1; ix<= numfibers; ix++)
   {
           for(int iy = 1; iy<= numfibers; iy++)
           {
                   x = -Absorb_xsize/2. + (Absorb_xsize/(Double_t)numfibers)*(ix - 0.5);
                   y = -Absorb_ysize/2. + (Absorb_ysize/(Double_t)numfibers)*(iy - 0.5);
//                 iCombTrans = new TGeoCombiTrans(x,y,crystalspace, new TGeoRotation(rot));
                   iCombTrans = new TGeoCombiTrans(x,y,0, new TGeoRotation(rot));
                   name = "FscFibHoleVolume";
                   name+=iCopy;
                   iCombTrans->SetName(name);
                   iCombTrans->RegisterYourself();
//                 ModuleVolume->AddNode(FiberHoleVolume,iCopy,iCombTrans);
//                 ModuleVolume->AddNodeOverlap(FiberHoleVolume,iCopy,iCombTrans);
                   ModuleVolume->AddNodeOverlap(FiberHoleVolume,iCopy,iCombTrans);  
//                 TyvekVolume->AddNode(FiberHoleVolume,iCopy,iCombTrans);
                   iCopy++;
           }
   }
   //FIbers itself inside holes
   name = "FscFiberShape";
//   TGeoShape* FiberShape = new TGeoTube(name, 0., fiber_r, (thickness-crystalspace)/2.);
   TGeoShape* FiberShape = new TGeoTube(name, 0., fiber_r, (holder_thickness)/2.);
   name = "FscFiberVolume";
   medium = "FscFiber";
   TGeoVolume* FiberVolume = new TGeoVolume(name, FiberShape, gGeoMan->GetMedium(medium));
   FiberVolume->SetLineColor(5);
   iCombTrans = new TGeoCombiTrans(0,0,0, new TGeoRotation(rot));
   FiberHoleVolume->AddNode(FiberVolume,0,iCombTrans);


   iCopy =0;
   Int_t CopyCode=0;
   // placing SuperModules inside top FSC
   // Bottom right corner is (ix = 1, iy = 1)
   for(int ix = 1; ix<= numcellsx/nClsSuperModuleX; ix++)
   {
           for(int iy = 1; iy<= numcellsy/nClsSuperModuleY; iy++)
           {
                   x = 0.5 * suprModXsize  * (2. * ix - numcellsx/nClsSuperModuleX - 1.);
                   y = 0.5 * suprModYsize  * (2. * iy - numcellsy/nClsSuperModuleY - 1.);
                   cout<<"======Supermodules"<<endl;
                   cout <<"x="<<x<<" y="<<y<<endl;
                   cout <<"ix="<<ix<<" iy="<<iy<<endl;
                   if(TMath::Abs(x-hole_xpos)<hole_xsize/2. && TMath::Abs(y-hole_ypos)<hole_ysize/2.) continue;
                   iCombTrans = new TGeoCombiTrans(x,y,0., new TGeoRotation(rot));
                   name = "FscSuperModuleVolume";
                   name+=iCopy;
                   CopyCode = 100*ix+iy;

                   iCombTrans->SetName(name);
                   iCombTrans->RegisterYourself();
//                 top->AddNode(ModuleVolume,iCopy,iCombTrans);
                   top->AddNode(SuperModuleVolume,CopyCode,iCombTrans);
                   iCopy++;

           }
   }


 
   //Adding mechanical support structure ==================================================================
//   TString FileName = "/home/freez/fairsoft_may16p1/pandaroot/macro/emc/dedicated/fsc_cad2root_mech_support_addnodeoverlap.root";
//   TString FileName = "/home/freez/fairsoft_may16p1/pandaroot/macro/emc/dedicated/fsc_cad2root_mech_support.root";
//   TString FileName = "/home/freez/fairsoft_may16p1/pandaroot/macro/emc/dedicated/fsc_cad2root_mech_support_reduced_overlaps.root";
//   TString FileName = "/home/freez/fairsoft_may16p1/pandaroot/macro/emc/dedicated/fsc_cad2root_mech_support_reduced_overlaps2.root";
//   TString FileName = "/home/freez/fairsoft_may16p1/pandaroot/macro/emc/dedicated/fsc_cad2root_mech_support_reduced_overlaps3.root";
   TString FileName = "fsc_cad2root_mech_support.root";
   TFile *fb = new TFile(FileName);
   TGeoVolume *MechFsc=(TGeoVolume *)fb->Get("topNode");
   fb->Close();
   MechFsc->Print();
   MechFsc->PrintNodes();

   TGeoNode *n=MechFsc->GetNode(0);
//   gGeoMan->AddVolume(MechFsc);
   
   n->Print();

   TGeoVolume* MFscVol = n->GetVolume();
    
   MFscVol->Print();

   MFscVol->PrintNodes();
//==========================================================================================================

   
   
   TGeoRotation rotMechFsc;
   rotMechFsc.RotateX(-90.);
   rotMechFsc.RotateY(180.);
   //DAM version
//   top->AddNode(MFscVol,0, new TGeoCombiTrans(11.0, 5.5, -38.9, new TGeoRotation(rotMechFsc)));
   //Slava Moiseev version
   top->AddNode(MFscVol,0, new TGeoCombiTrans(11.0, 4.8, -38.9, new TGeoRotation(rotMechFsc)));
//   top->AddNode(MFscVol,0, new TGeoCombiTrans(0, -18.55, 16.6, new TGeoRotation(rotMechFsc)));
   
   
   
   gGeoMan->CloseGeometry();
   fi->cd();
   top->Write();
   fi->Close();
  //gGeoMan->Export(outfile);//
//   gGeoManager->SetVisLevel(5);
   top->Draw("ogl");//

}