BesGeoSelector.cxx

Go to the documentation of this file.

/*
 *      2018/5/11   Long Peixun     IHEP
 *                  Geometry Selection for EventDisplay
 */
 
#include "BesGeoSelector.h"
#include "DetVerSvc/IDetVerSvc.h"
#include "GaudiKernel/Bootstrap.h"
#include "GaudiKernel/ISvcLocator.h"
#include <algorithm>
#include <iostream>
 
using namespace std;
 
BesGeoSelector::BesGeoSelector() {}
 
void BesGeoSelector::addGeomFile( const char* fname, int startno ) {
  cout << "BesGeoSelector: " << fname << " start phase " << startno << endl;
  fgeom.push_back( GeomFileInfo( fname, startno ) );
  sort( fgeom.begin(), fgeom.end() );
}
 
void BesGeoSelector::clear() { fgeom.clear(); }
 
// Given run number, select a suitable geometry file, and return file name.
const char* BesGeoSelector::getGeomFileNameFromRunNo( int run_no ) {
  auto dvs = Gaudi::svcLocator()->service<IDetVerSvc>( "DetVerSvc" );
  if ( !dvs )
  {
    cerr << "BesGeoSelector Error: DetVerSvc not found!" << endl;
    return "";
  }
  // dvs->initialize();
  int absno = ( run_no >= 0 ? run_no : -run_no );
 
  for ( vector<GeomFileInfo>::reverse_iterator i = fgeom.rbegin(); i != fgeom.rend(); ++i )
  {
    // if (absno >= i->startno)
    // {
    //     return i->filename.c_str();
    // }
 
    // In BOSS 7.0.3, use DecVerSvc package to choose geometry
    int phase = dvs->fromRun( absno );
    if ( phase >= i->startno ) { return i->filename.c_str(); }
  }
  cout << "BesGeoSelector Error: There is no suitable detector geometry!" << endl;
  return "";
}
 
// Given file name, find whether the list contains this file.
bool BesGeoSelector::hasGeomFile( const char* fname ) {
  string filename( fname );
  for ( vector<GeomFileInfo>::iterator i = fgeom.begin(); i != fgeom.end(); ++i )
  {
    if ( filename == i->filename ) return true;
  }
  return false;
}