er  dev
Import_GDML_Export_ROOT.c
1 void Import_GDML_Export_ROOT()
2 {
3  TString vmcDir = TString(gSystem->Getenv("VMCWORKDIR"));
4 
5  TGeoManager *gdml = new TGeoManager("gdml", "FAIRGeom");
6 
7  TGDMLParse parser;
8  // Define your input GDML file HERE
9  TGeoVolume* gdmlTop = parser.GDMLReadFile ("gadast.gdml");
10  TGeoVolume* rootTop = new TGeoVolumeAssembly("TOP");
11 
12  gGeoManager->SetTopVolume(gdmlTop);
13 
14 
15  // Define your position HERE
16  TGeoRotation* rot = new TGeoRotation ("rot", 0., 0., 0.);
17  TGeoCombiTrans* posrot = new TGeoCombiTrans (0., 0., 0., rot);
18 
19  rootTop->AddNode (gdmlTop, 1, posrot);
20 
21  gGeoManager->CloseGeometry();
22 
23  // Define you output ROOT file HERE
24  TFile* outfile = new TFile(vmcDir + "/geometry/gadast.root", "RECREATE");
25  rootTop->Write();
26  outfile->Close();
27 
28 }
29