Updating geometry

alphacamm-geometry/ViewRESTGdml.C

@@ -0,0 +1,46 @@
+void ViewRESTGdml(string filename, float transparency = 50)
+{
+    // first import it and see if there is something wrong with the file
+    // we need TRestGDMLParser to visualize geometries that use external files
+
+    TRestGDMLParser *g = new TRestGDMLParser();
+    //TGeoManager::Import(filename.c_str());
+    //TGeoNode *node = gGeoManager->GetTopNode();
+    TGeoManager *gGeoManager = g->GetGeoManager(filename.c_str());
+    TGeoNode *node = gGeoManager->GetTopNode();
+    // check for overlaps
+    node->CheckOverlaps(0.0001);
+    // print overlaps
+    gGeoManager->PrintOverlaps();
+
+    for (auto i = 0; i < gGeoManager->GetListOfVolumes()->GetEntries(); i++)
+    {
+        TGeoVolume *geoVolume = gGeoManager->GetVolume(i); // https://root.cern/doc/v606/classTGeoVolume.html
+        if (!geoVolume)
+        {
+            continue;
+        }
+        TGeoMaterial *material = geoVolume->GetMaterial();
+        double density = material->GetDensity();
+        string materialName = material->GetName();
+
+        geoVolume->SetTransparency(transparency);
+        if (density < 0.10 /* g/cm3 */)
+        {
+            geoVolume->SetVisibility(kFALSE);
+        }
+    }
+
+    TEveManager::Create();
+
+    auto viewer = gEve->GetDefaultGLViewer();
+    viewer->SetClearColor(0);
+    viewer->SetCurrentCamera(TGLViewer::kCameraPerspXOZ);
+
+    TEveGeoTopNode *top_node = new TEveGeoTopNode(gGeoManager, node);
+    gEve->AddGlobalElement(top_node);
+
+    gEve->FullRedraw3D(kTRUE);
+
+    viewer->CurrentCamera().Reset();
+}