void Geo() {
  gSystem->Load("libGeom");
 
  TGeoManager* m_Geometry = new TGeoManager("Strange", "behaviour");

  TGeoMaterial* Al = new TGeoMaterial("Al", 26.98,13, 2.7);
  TGeoMedium* Medium = new TGeoMedium("Medium", 1, Al);

  // Level 0
  TGeoTranslation T0(0, 0, 0);
  TGeoRotation R0("rLevel0", 90, 0, 90, 90, 0, 0);
  TGeoCombiTrans* C0 = new TGeoCombiTrans(T0, R0);
  //TGeoBBox* Shape0 = new TGeoBBox(5, 5, 5);
  //TGeoVolume* Volume0 = new TGeoVolume("vLevel0", Shape0, 0);
  TGeoVolume* Volume0 = m_Geometry->MakeBox("vLevel0", Medium, 5, 5, 5);
  m_Geometry->SetTopVolume(Volume0);
  m_Geometry->SetTopVisible(kTRUE);

  // Level 1
  TGeoTranslation T1(0, 0, 0);
  TGeoRotation R1("rLevel1", 90, 0, 90, 90, 0, 0);
  TGeoCombiTrans* C1 = new TGeoCombiTrans(T1, R1);
  //TGeoBBox* Shape1 = new TGeoBBox(4, 4, 4);
  //TGeoVolume* Volume1 = new TGeoVolume("vLevel1", Shape1, 0);
  TGeoVolume* Volume1 = m_Geometry->MakeBox("vLevel1", Medium, 4, 4, 4);
  Volume1->SetVisibility(kTRUE);
  Volume1->SetLineColor(kGreen);
  Volume0->AddNode(Volume1, 1, C1);

  // Level 2
  TGeoTranslation T2(0, 0, 0);
  TGeoRotation R2("rLevel2", 90, 0, 90, 90, 0, 0);
  TGeoCombiTrans* C2 = new TGeoCombiTrans(T2, R2);
  //TGeoBBox* Shape2 = new TGeoBBox(3, 3, 3);
  //TGeoVolume* Volume2 = new TGeoVolume("vLevel2", Shape2, 0);
  TGeoVolume* Volume2 = m_Geometry->MakeBox("vLevel2", Medium, 3, 3, 3);
  Volume2->SetVisibility(kTRUE);
  Volume2->SetLineColor(kBlue);
  Volume1->AddNode(Volume2, 2, C2);

  // Level 3
  TGeoTranslation T3(0, 0, 0);
  TGeoRotation R3("rLevel3", 90, 0, 90, 90, 0, 0);
  TGeoCombiTrans* C3 = new TGeoCombiTrans(T3, R3);
  //TGeoBBox* Shape3 = new TGeoBBox(2, 2, 2);
  //TGeoVolume* Volume3 = new TGeoVolume("vLevel3", Shape3, 0);
  TGeoVolume* Volume3 = m_Geometry->MakeBox("vLevel3", Medium, 2, 2, 2);
  Volume3->SetVisibility(kTRUE);
  Volume3->SetLineColor(kRed);
  Volume2->AddNode(Volume3, 3, C3);

  m_Geometry->CloseGeometry();
  //m_Geometry->SetVisLevel(1000);
  //m_Geometry->SetVisDensity(0.0);
  //m_Geometry->SetMaxVisNodes(40000);
  
  if (m_Geometry->GetTopVolume() != 0) m_Geometry->GetTopVolume()->Draw("ogle");
  
  cout<<"Volume below top volume: "<<m_Geometry->GetTopVolume()->GetNode(0)->GetName()<<endl;
}