MvdTBGeometry.C

Go to the documentation of this file.

/*
 * CreateLambdaGeometry.C
 *
 *  Created on: Mar 21, 2013
 *      Author: stockman
 */
MvdTBGeometry()
{
          //-----------------------------
        Double_t positionStage1 = 10.0;         //in cm
        Double_t positionStage2 = 16.0;         //in cm
        Double_t positionStage3 = 22.0;         //in cm
        Double_t positionStage4 = 28.0;         //in cm

          //--------------------------------------------------------------------

          
          TString outfile= "../../geometry/MvdTBGeometry.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 *CbmMediumSilicon  = Media->getMedium("silicon");
          
          Int_t nmed=geobuild->createMedium(CbmMediumSilicon);

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

          TGeoVolume *top = new TGeoVolumeAssembly("top");
                            
          gGeoMan->SetTopVolume(top);
          

        TGeoVolumeAssembly* TBSetup = new TGeoVolumeAssembly("MvdTestbeamSetup");
        TGeoVolumeAssembly* module = new TGeoVolumeAssembly("ToPix4Module");

        TGeoShape* ToPix4Shape = new TGeoBBox("ToPix4Shape", 0.2/2, 0.32/2, 0.0285/2);                  //todo: fix dimensions of ToPix4
        TGeoVolume* ToPix4Volume = new TGeoVolume("ToPix4Volume", ToPix4Shape, gGeoMan->GetMedium("silicon"));

        TGeoShape* EpiSensorShape = new TGeoBBox("EpiSensorShape", 0.2/2, 0.32/2, 0.0100/2);
        TGeoVolume* EpiSensorVolume = new TGeoVolume("EpiSensorVolume", EpiSensorShape, gGeoMan->GetMedium("silicon"));


        TGeoTranslation transToPix("ToPixTranslation", .0, .0, 0.0285/2 + 0.0100/2 + 0.002);

        module->AddNode(EpiSensorVolume, 0);
        module->AddNode(ToPix4Volume, 0, new TGeoHMatrix(transToPix));

        TBSetup->AddNode(module, 0, new TGeoTranslation(.0,.0,positionStage1));
        TBSetup->AddNode(module, 1, new TGeoTranslation(.0,.0,positionStage2));
        TBSetup->AddNode(module, 2, new TGeoTranslation(.0,.0,positionStage3));
        TBSetup->AddNode(module, 3, new TGeoTranslation(.0,.0,positionStage4));

        top->AddNode(TBSetup,0);

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