lmd_geohandler_benchmark.C

Go to the documentation of this file.

// compile with
// g++ $(root-config --cflags --glibs) -lGeom -I${PANDAROOT}/geobase -L${PANDAROOT}/buildPanda/lib -lGeoBase create_Lmd_HV_MAPS.C -o /tmp/create_geometry -I${PANDAROOT}/lmd/LmdMC
// and run with
// /tmp/create_geometry

#include<TROOT.h>
#include<FairGeoLoader.h>
#include<FairGeoInterface.h>
#include<FairGeoBuilder.h>
#include<FairGeoPcon.h>
#include<FairGeoMedia.h>
#include<TGeoCompositeShape.h>
#include<TGeoMatrix.h>
#include<TGeoVolume.h>
#include<TGeoTube.h>
#include<TGeoArb8.h>
#include<TGeoPara.h>
#include<TSystem.h>
#include<TFile.h>
#include<TGeoPcon.h>
#include<TGeoManager.h>
#include<TGeoTorus.h>
#include<TGeoSphere.h>
#include<TMath.h>
#include<TRandom.h>
#include<PndLmdDim.h>
#include<TGeoPhysicalNode.h>
#include<TGeoCone.h>
#include<PndGeoHandling.h>

#include<iostream>
#include<cmath>
#include<sstream>

//using namespace LmdDim;

// to check the position and dimensions
// one may enable a beam pipe dummy to be displayed
const bool show_beam_pipe_dummy = true;

const bool include_box = false;

void create_HV_MAPS(bool misalign = false) {

//        TString outfile= "geo/HV_MAPS-Design.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 *FairMediumSilicon = Media->getMedium("silicon");

          Int_t nmed=geobuild->createMedium(FairMediumSilicon);

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

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

          gGeoMan->SetTopVolume(top);

          PndLmdDim& lmddim = PndLmdDim::Get_instance();

          // for testing purposes
          lmddim.Generate_rootgeom(*top, misalign);

          //gGeoMan->LocalToMaster();
          LocalToMasterShortId

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

#include<TApplication.h>

int main() {
        TApplication myapp("myapp", 0, 0);
        create_HV_MAPS();
        myapp.Run();
        return 0;
}