//
// $Id: d0StandaloneGeometry.cpp,v 1.2 2001/03/16 01:49:33 hobbs Exp $
//
// File: smt_get_geometry.cpp
// Purpose: 
// Created:  1-FEB-2001   John D. Hobbs
//
// $Revision: 1.2 $
//
//

// Include files
#include <string>
#include <vector>
#include "d0_geometry/base/d0StandaloneGeometry.hpp"

using namespace std;


d0StandaloneGeometry::d0StandaloneGeometry(const int argc, char **argv):
  _file(false), _ftype("EVPACK"), _fname("SiBaseGeometry")
{
  // Stash the command line arguments
  for( int i=0 ; i<argc ; i++ ) _argv.push_back(argv[i]);

  // Figure out which ones I can handle and flag 'em.  Also set up creation
  // information.  Don't actually do any work yet, though.
  int i=0;
  while( i<_argv.size() ) {
    if( _argv[i][0] == '-' ) {
      switch ( _argv[i][1] ) {
      case 'd':
	_used_args.push_back(i);
	_ftype="DSPACK";
	break;
      case 'f':
	_file=true;
      case 'r':    // This is the same for both -f and -r
	_used_args.push_back(i);
	if( argv[i][2] == '=' ) {
	  _fname=_argv[i].substr(3,_argv[i].size()-3);
	}
	else {
	  _fname=_argv[++i];
	  _used_args.push_back(i);
	}
        break;
      }
    }
    ++i;
  }

}


d0_Ref<BaseGeometry> d0StandaloneGeometry::get_geometry() const
{

  // If we've already gotten a geometry, use it
  if( !_d0det.is_null() ) return _d0det;

  // Get the geometry the first time around.
  if( !_file ) build_default();
  else read_geometry();

  // and hand it off...
  return _d0det;
}


void d0StandaloneGeometry::get_modified_args(vector<string> &partialArgv) const 
{
  int i=0,iskip=0;
  partialArgv.clear();
  // This loop relies on the entries in _used_args being sorted in increasing
  // order.  This is guaranteed by the parser in the constructor.
  while( i<_argv.size() ) {
    if( iskip>=_used_args.size() || _used_args[iskip] != i ) partialArgv.push_back(_argv[i]);
    else if( iskip<_used_args.size() ) ++iskip;
    ++i;
  }
}


#include "rcp/RCP.hpp"
#include "rcp/RCPManager.hpp"
using namespace edm;

void d0StandaloneGeometry::build_default() const
{  
  cout << "Building the geometry from RCP in " << _fname << endl;
  
  // Instantiate an RCP manager and read the silicon geometry RCP...
  RCPManager* rcpman = RCPManager::instance();
  RCP d0DefaultRCP = rcpman->extract("d0_geometry",_fname);
  
  // Instantiate a copy of the base geometry...
  _d0det = BaseGeometry::build_default(&d0DefaultRCP);

}

#include "d0stream/d0StreamFactory.hpp"
#include "d0stream/d0Stream.hpp"

void d0StandaloneGeometry::read_geometry() const 
{
  cout << "Reading the geometry from " << _fname << endl;
  _d0det = BaseGeometry::readGeometry(_fname.c_str(),_ftype.c_str());

}

ostream& operator <<(ostream& os, const d0StandaloneGeometry& me)
{
  os << "[-f=<geometry_file_name>] [-r=<geometry_rcp>] [-d]" << endl;
  return os;
}