EmcRecGeoSvc.cxx

Go to the documentation of this file.

#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/PropertyMgr.h"
#include "GaudiKernel/StatusCode.h"
 
#include "EmcRecGeoSvc.h"
#include "EmcRecGeoSvc/EmcRecROOTGeo.h"
 
DECLARE_COMPONENT( EmcRecGeoSvc )
 
// static SvcFactory<EmcRecGeoSvc> s_factory;
// const ISvcFactory& EmcRecGeoSvcFactory = s_factory;
 
EmcRecGeoSvc::EmcRecGeoSvc( const std::string& name, ISvcLocator* svcloc )
    : base_class( name, svcloc ) {
  fGdml = true;
 
  // Declare the properties
  declareProperty( "Gdml", fGdml );
 
  if ( fGdml ) { fROOTGeo = new EmcRecROOTGeo(); }
}
 
EmcRecGeoSvc::~EmcRecGeoSvc() {
  if ( fGdml )
  {
    if ( fROOTGeo ) delete fROOTGeo;
  }
}
 
/*StatusCode EmcRecGeoSvc::queryInterface (const InterfaceID& riid, void** ppvInterface ){
 
  if ( IID_IEmcRecGeoSvc.versionMatch(riid) ) {
    *ppvInterface = static_cast<IEmcRecGeoSvc*> (this);
  } else {
    return Service::queryInterface(riid, ppvInterface) ;
  }
  return StatusCode::SUCCESS;
}
*/
StatusCode EmcRecGeoSvc::initialize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << ": Start of run initialisation" << endmsg;
 
  StatusCode sc = Service::initialize();
  if ( sc.isFailure() ) return sc;
 
  // initialize geometry from gdml
  if ( fGdml ) { fROOTGeo->InitFromXML(); }
  return StatusCode::SUCCESS;
}
 
StatusCode EmcRecGeoSvc::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << ": End of Run" << endmsg;
  return StatusCode::SUCCESS;
}
 
EmcRecCrystal EmcRecGeoSvc::GetCrystal( const Identifier& id ) const {
  if ( fGdml ) { return fROOTGeo->GetCrystal( id ); }
  else
  {
    EmcRecCrystal cry;
 
    if ( EmcID::is_barrel( id ) ) { cry = fBarrel.GetCrystal( id ); }
    else { cry = fEndCap.GetCrystal( id ); }
 
    return cry;
  }
}
 
HepPoint3D EmcRecGeoSvc::GetCrystalPoint( const Identifier& id, const int i ) const {
  EmcRecCrystal cry;
  if ( fGdml ) { cry = fROOTGeo->GetCrystal( id ); }
  else
  {
    if ( EmcID::is_barrel( id ) ) { cry = fBarrel.GetCrystal( id ); }
    else { cry = fEndCap.GetCrystal( id ); }
  }
  return cry.Get( i ) / cm;
}
 
HepPoint3D EmcRecGeoSvc::GetCCenter( const Identifier& id ) const {
  if ( fGdml ) { return fROOTGeo->GetCCenter( id ) / cm; }
  else
  {
    HepPoint3D center;
 
    if ( EmcID::is_barrel( id ) ) { center = fBarrel.GetCCenter( id ); }
    else { center = fEndCap.GetCCenter( id ); }
 
    return center / cm;
  }
}
 
HepPoint3D EmcRecGeoSvc::GetCFrontCenter( const Identifier& id ) const {
  if ( fGdml ) { return fROOTGeo->GetCFrontCenter( id ) / cm; }
  else
  {
    HepPoint3D frontCenter;
 
    if ( EmcID::is_barrel( id ) ) { frontCenter = fBarrel.GetCFrontCenter( id ); }
    else { frontCenter = fEndCap.GetCFrontCenter( id ); }
 
    return frontCenter / cm;
  }
}
 
double EmcRecGeoSvc::GetBarrelR() const { return fBarrel.GetBarrelR() / cm; }
 
double EmcRecGeoSvc::GetBarrelOffset1() const { return fBarrel.GetBarrelOffset1() / cm; }
 
double EmcRecGeoSvc::GetBarrelOffset2() const { return fBarrel.GetBarrelOffset2() / cm; }
 
double EmcRecGeoSvc::GetBarrelh1() const { return fBarrel.GetBarrelh1() / cm; }
 
double EmcRecGeoSvc::GetBarrelh2() const { return fBarrel.GetBarrelh2() / cm; }
 
double EmcRecGeoSvc::GetBarrelh3() const { return fBarrel.GetBarrelh3() / cm; }
 
double EmcRecGeoSvc::GetBarrelL() const { return fBarrel.GetBarrelL() / cm; }
 
int EmcRecGeoSvc::GetBarrelNPhiMax() const { return fBarrel.GetBarrelNPhiMax(); }
 
int EmcRecGeoSvc::GetBarrelNThetaMax() const { return fBarrel.GetBarrelNThetaMax(); }