createLmdGeometry.C File Reference

#include <iostream>
#include "detectors/lmd/LmdMC/PndLmdGeometryFactory.h"
#include "FairGeoLoader.h"
#include "FairGeoInterface.h"
#include "FairGeoBuilder.h"
#include "FairGeoMedia.h"
#include "TFile.h"
#include "TROOT.h"
#include "TGeoManager.h"
#include "TGeoVolume.h"
#include <boost/property_tree/json_parser.hpp>

Go to the source code of this file.

Functions
void	createLmdGeometry (std::string output_filename="Luminosity-Detector.root", std::string geo_params_config_file_url="../lmd-geo-params.json")

Function Documentation

void createLmdGeometry	(	std::string	output_filename = "Luminosity-Detector.root",
		std::string	geo_params_config_file_url = "../lmd-geo-params.json"
	)

Definition at line 19 of file createLmdGeometry.C.

References file, PndLmdGeometryFactory::generateLmdGeometry(), geobuild, geoFace, geoLoad, gGeoMan, PndLmdGeometryFactory::init(), Media, pandaroot_dir(), and top.

                                                                                 {
        // load parameters
        ptree geo_property_tree;
        read_json(geo_params_config_file_url, geo_property_tree);

        FairGeoLoader* geoLoad = new FairGeoLoader("TGeo", "FairGeoLoader");
        FairGeoInterface *geoFace = geoLoad->getGeoInterface();
        std::string pandaroot_dir(getenv("VMCWORKDIR"));
        geoFace->setMediaFile((pandaroot_dir + "/geometry/media_pnd.geo").c_str());
        geoFace->readMedia();
        FairGeoMedia *Media = geoFace->getMedia();
        FairGeoBuilder *geobuild = geoLoad->getGeoBuilder();

        TFile file(output_filename.c_str(), "RECREATE");

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

        PndLmdGeometryFactory lmd_geo_factory(geo_property_tree);
        lmd_geo_factory.init(geoLoad);
        TGeoVolumeAssembly *top = lmd_geo_factory.generateLmdGeometry();

        gGeoMan->CheckOverlaps(0.0001); // 1mu
        gGeoMan->PrintOverlaps();

        top->Write();
        file.Close();
        
        gGeoMan->SetVisLevel(20);
        top->Draw("ogle");
        
}