{
  // Open and draw things
  TGeoManager::Import("test.root");
  //TGeoManager::Import("test2.root");
  //gGeoManager->DefaultColors();
  gGeoManager->GetVolume("comet")->SetVisibility(false);
//  gGeoManager->GetVolume("ProdTgtSec")->SetVisibility();
//  gGeoManager->GetVolume("Tor1ToStopTgtSec")->SetVisibility();
//  gGeoManager->GetVolume("Torus1")->SetVisibility();
  gGeoManager->GetTopVolume()->Draw("ogl");
  gGeoManager->SetVisLevel(-1);
  //gGeoManager->SetVisOption(0);
  TGeoVolume* vol;
  //for(int i=0;i<5;i++){
  //        vol=gGeoManager->GetVolume(Form("StrawStation_%d",i));
  //        cout<<vol->GetName()<<endl;
  //        vol->SetVisContainers(false);
  //}
  //vol=gGeoManager->GetVolume("StrawTracker");
  //cout<<vol->GetName()<<endl;
  //vol->SetVisContainers();

  //double precision = 0.02; // cm
  //gGeoManager->CheckOverlaps(precision);
  //gGeoManager->PrintOverlaps();

  // Now setup the image
  TGLViewer* view=(TGLSAViewer*) gPad->GetViewer3D();
  if(!view) {
      cout<<"Unable to get view3D!"<<endl;
    return;
  }
  // Cliping
  TGLClipSet* clipping = view->GetClipSet();
  Double_t* coords[]={0,1,0,0};
  //TGLClip::EType clipType=0;
  //EClipType clipType=1;
  clipping->SetClipType(1 );
  //TGLPlaneSet_t planes(1);
  //planes[0]=new TGLPlane(0,1,0,0) ;
  //clipping->FillPlaneSet(planes);
  //clipping->SetClipState(clipType,coords );

  // Camera direction 
  TGLCamera* camera=&view->CurrentCamera();
  camera->SetCenterVecWarp(0,0,0);

}