DedxSimSvc.cxx

Go to the documentation of this file.

#include "CLHEP/Vector/ThreeVector.h"
#include "GaudiKernel/DataSvc.h"
#include "GaudiKernel/IIncidentListener.h"
#include "GaudiKernel/IIncidentSvc.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/Incident.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/StatusCode.h"
#include "TH1F.h"
#include <fstream>
#include <iostream>
#include <math.h>
#include <string>
#include <vector>
 
#include "CalibData/CalibModel.h"
#include "CalibData/Dedx/DedxCalibData.h"
#include "CalibData/Dedx/DedxSimData.h"
#include "EventModel/EventHeader.h"
#include "Identifier/Identifier.h"
 
#include "DedxSimSvc.h"
 
using namespace std;
using CLHEP::Hep3Vector;
 
DECLARE_COMPONENT( DedxSimSvc )
 
DedxSimSvc::DedxSimSvc( const std::string& name, ISvcLocator* svcloc )
    : base_class( name, svcloc ) {}
 
DedxSimSvc::~DedxSimSvc() {}
/*
StatusCode DedxSimSvc::queryInterface(const InterfaceID& riid, void** ppvInterface){
    if( IID_IDedxSimSvc.versionMatch(riid) ){
        *ppvInterface = static_cast<IDedxSimSvc*> (this);
    } else{
        return Service::queryInterface(riid, ppvInterface);
    }
    return StatusCode::SUCCESS;
}
*/
 
StatusCode DedxSimSvc::initialize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << "DedxSimSvc::initialize()" << endmsg;
 
  StatusCode sc = Service::initialize();
  if ( sc.isFailure() ) return sc;
 
  IIncidentSvc* incsvc;
  sc           = service( "IncidentSvc", incsvc );
  int priority = 100;
  if ( sc.isSuccess() ) { incsvc->addListener( this, "NewRun", priority ); }
  sc = serviceLocator()->service( "EventDataSvc", m_eventSvc, true );
  if ( sc.isFailure() )
  {
    log << MSG::ERROR << "Unable to find EventDataSvc " << endmsg;
    return sc;
  }
 
  m_runfrom      = -999999;
  m_runto        = -999999;
  m_version      = -1;
  m_numDedxHists = 0;
  m_numBg        = 0;
  m_numTheta     = 0;
  m_dedx_hists   = nullptr;
  m_bgRange      = nullptr;
  return StatusCode::SUCCESS;
}
 
StatusCode DedxSimSvc::finalize() {
  MsgStream log( msgSvc(), name() );
  log << MSG::INFO << name() << "DedxSimSvc::finalize()" << endmsg;
  return StatusCode::SUCCESS;
}
 
void DedxSimSvc::handle( const Incident& inc ) {
  MsgStream log( msgSvc(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endmsg;
 
  if ( inc.type() == "NewRun" )
  {
    log << MSG::DEBUG << "New Run" << endmsg;
    SmartDataPtr<Event::EventHeader> evt( m_eventSvc, "/Event/EventHeader" );
    int runNo=0;
    if ( evt )
    {
      runNo = evt->runNumber();
      log << MSG::DEBUG << "The runNumber of current event is  " << runNo << endmsg;
    }
    else { log << MSG::ERROR << "ERROR accessing Event" << endmsg; }
    if ( runNo < m_runfrom || runNo > m_runto ) { update_param_svc(); }
  }
}
 
void DedxSimSvc::update_param_svc() {
  MsgStream         log( msgSvc(), name() );
  IDataProviderSvc* pCalibDataSvc;
  StatusCode        sc = Gaudi::svcLocator()->service( "CalibDataSvc", pCalibDataSvc, true );
  if ( !sc.isSuccess() )
  {
    log << MSG::ERROR << "Could not get IDataProviderSvc interface of CalibXmlCnvSvc"
        << endmsg;
    return;
  }
  else
  { log << MSG::DEBUG << "Retrieved IDataProviderSvc interface of CalibXmlCnvSvc" << endmsg; }
 
  SmartDataPtr<CalibData::DedxSimData> pDedxSimData( pCalibDataSvc, "/Calib/DedxSim" );
  m_runfrom      = pDedxSimData->getrunfrm();
  m_runto        = pDedxSimData->getrunto();
  m_version      = pDedxSimData->getVersion();
  m_numDedxHists = pDedxSimData->gethistNo();
  m_numBg        = pDedxSimData->getRangeNo();
  if ( m_version == 0 ) m_numTheta = 10;
  else m_numTheta = pDedxSimData->getThetaNo();
  m_dedx_hists = pDedxSimData->getHist();
  m_bgRange    = pDedxSimData->getRange();
  log << MSG::DEBUG << "DedxSimSvc::update_param_svc() finish" << endmsg;
}