void Import_GDML_Export_ROOT()
{
	TGeoManager *gdml = new TGeoManager("gdml","My 3D Project");

	TGDMLParse parser;
	// Define your input GDML file HERE
	TGeoVolume* gdmlTop = parser.GDMLReadFile ("magnet_v10e.gdml");

	TGeoVolume* rootTop = new TGeoVolumeAssembly("geometryFromGDML");
	gGeoManager->SetTopVolume(rootTop);

	// Define your position HERE
	TGeoRotation* rot = new TGeoRotation ("rot", 0., 0., 0.);
	TGeoCombiTrans* posrot = new TGeoCombiTrans (0., 0., 0., rot);

	rootTop->AddNode (gdmlTop, 1, posrot);

	gGeoManager->CloseGeometry();
	// Define you output ROOT file HERE
	TFile* outfile = new TFile("magnet_v10e.root", "RECREATE");
	rootTop->Write();
	outfile->Close();
}