//
// $Id: D0MaterialGeometry_t.cpp,v 1.12 2000/12/05 16:52:16 hobbs Exp $
//
// File: test_material.cpp
// Purpose: 
// Created: 08-JUN-1998  Harry Melanson
//
// $Revision: 1.12 $
//
//

// Include files

#include "d0_geometry/base/BaseGeometry.hpp"
#include "d0_geometry/material/MaterialGeometry.hpp"

#include <cft_util/CFTDetector.hpp>
#include <silicon_geometry/base/SiGeometerDecoupler.hpp>
#include <cft_geometry/base/CftGeometerDecoupler.hpp>
#include <cps_geometry/base/CpsGeometerDecoupler.hpp>
#include <fps_geometry/base/FPSGeometerDecoupler.hpp>
#include <calorimeter_geometry/base/CalGeometerDecoupler.hpp>
#include <muon_geometry/base/MuoGeometerDecoupler.hpp>

#include <silicon_geometry/material/SiMaterialGeometer.hpp>
#include <cft_geometry/material/CftMaterialGeometer.hpp>

#include <rcp/RCPManager.hpp>                          // RCP stuff
#include <rcp/RCP.hpp>
#include <stdexcept>
#include <iostream>

using namespace std;
using edm::RCPManager;
using edm::RCP;

int main(int argc, char **argv)
//
// Purpose: Test MaterialGeometry
//
{
  d0om_init("materialG_testing");

  //  Base Geometry

  RCPManager* rcpman = RCPManager::instance();
  d0_Ref<BaseGeometry> base;

  try {
    // Build the default base geometry
    RCP DefaultRCP = rcpman->extract("d0_geometry","BaseGeometry");
    base = BaseGeometry::build_default(&DefaultRCP);

    // Update each of the geometers in turn.  Some Material geometries 
    // require a valid base geometry as part of initialization...
    SiGeometerDecoupler smt;
    d0_Ref<const SiBaseGeometry> smtgeom = base->get_silicon();
    smt.refresh(smtgeom);
 
    CftGeometerDecoupler cft;
    d0_Ref<const CftBaseGeometry> cftgeom = base->get_cft();
    cft.refresh(cftgeom);
 
    CpsGeometerDecoupler cps;
    d0_Ref<const CpsBaseGeometry> cpsgeom = base->get_cps();
    cps.refresh(cpsgeom);
 
    FPSGeometerDecoupler fps;
    d0_Ref<const FPSBaseGeometry> fpsgeom = base->get_fps();
    fps.refresh(fpsgeom);

    CalGeometerDecoupler cal;
    d0_Ref<const CalBaseGeometry> calgeom = base->get_calorimeter();
    cal.refresh(calgeom);
 
    MuoGeometerDecoupler muo;
    d0_Ref<const MuoBaseGeometry> muogeom  = base->get_muon();
    muo.refresh(muogeom);
 
  }
  catch(runtime_error mistake) {
    cerr << "Caught error: " << mistake.what() << endl;
    return(1);
  }
  catch(const edm::XRCP& x) {
    cerr << "Caught error: " << x.logMessage() << endl;
    return(1);
  }

  //  Material Geometry

  MaterialGeometry* material_geom = new MaterialGeometry();
  try {
    RCP mainrcp = rcpman->extract("d0_geometry","MaterialGeometry");
    material_geom->build(&mainrcp);
  }
  catch(runtime_error mistake) {    cerr << "Caught error: " << mistake.what() << endl;
    return(1);
  }
  catch(const edm::XRCP& x) {
    cerr << "Caught error: " << x.logMessage() << endl;
    return(1);
  }

  cout << "Get Material Geometer" << endl;
  CftMaterialGeometer* cft_mat = CftMaterialGeometer::get_instance();
  cout << "CFT Material Geometry" << endl;
  cft_mat->print_geometry();
  return 0;

}