{
  TGeoManager *geo = new TGeoManager("Test","Test Ellipsoid");
  geo->SetVisOption(0);
  TGeoMaterial *mat_fe = new TGeoMaterial("Fe",55.845,26,7.87);
  TGeoMedium *fe = new TGeoMedium("Fe",1,mat_fe);
  TGeoMaterial *mat_vac = new TGeoMaterial("Vacuum",0,0,0);
  TGeoMedium *vac= new TGeoMedium("Vacuum",1,mat_vac);
  
  TGeoVolume *root=geo->MakeBox("root",vac,1000.,1000.,1000.);

  TGeoSphere *tmp_sp= new TGeoSphere("sphere",0.,500.);
  TGeoScale  *tmp_scale=new TGeoScale("scale",0.2,0.4,1.);
  TGeoScaledShape *tmp_shape = new TGeoScaledShape("elle",tmp_sp,tmp_scale);
  TGeoVolume *ellip=new TGeoVolume("ellipsoid",tmp_shape,fe);
  root->AddNode(ellip,1,0);
  gGeoManager->SetTopVolume(root);
  gGeoManager->CloseGeometry();
  gGeoManager->GetTopVolume()->Draw();
}