DedxSimSvc.cxx

Go to the documentation of this file.

#include "DedxSimSvc/DedxSimSvc.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/StatusCode.h"
#include "GaudiKernel/SvcFactory.h"
#include "GaudiKernel/MsgStream.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/DataSvc.h"
 
#include "GaudiKernel/IIncidentSvc.h"
#include "GaudiKernel/Incident.h"
#include "GaudiKernel/IIncidentListener.h"
 
#include "Identifier/Identifier.h"
#include "CLHEP/Vector/ThreeVector.h"
#include "CalibData/Dedx/DedxCalibData.h"
#include "CalibData/CalibModel.h"
#include <math.h>
#include <vector>
#include <iostream>
#include <fstream>
#include <string>
#include "EventModel/EventHeader.h"
#include "TH1F.h"
#include "CalibData/Dedx/DedxSimData.h"
#include "CalibDataSvc/CalibDataSvc.h"
 
using namespace std;
using CLHEP::Hep3Vector;
 
 
DedxSimSvc::DedxSimSvc( const string& name, ISvcLocator* svcloc) :
    Service (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(messageService(), name());
    log << MSG::INFO << name() << "DedxSimSvc::initialize()" << endreq;
 
    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 " << endreq;
      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(messageService(), name());
    log << MSG::INFO << name() << "DedxSimSvc::finalize()" << endreq;
    return StatusCode::SUCCESS;
}
 
void DedxSimSvc::handle(const Incident& inc){
    MsgStream log( messageService(), name() );
    log << MSG::DEBUG << "handle: " << inc.type() << endreq;
 
    if ( inc.type() == "NewRun" ){
        log << MSG::DEBUG << "New Run" << endreq;
        SmartDataPtr<Event::EventHeader> evt(m_eventSvc,"/Event/EventHeader");
        int runNo;
        if( evt ){
            runNo = evt -> runNumber();
            log << MSG::DEBUG <<"The runNumber of current event is  "<<runNo<<endreq;
        }
        else{
            log << MSG::ERROR << "ERROR accessing Event" <<endreq;
        }
        if(runNo<m_runfrom || runNo>m_runto){
            update_param_svc();
        }
    }
}
 
 
void
DedxSimSvc::update_param_svc() {
    MsgStream log(messageService(), name());
    IDataProviderSvc* pCalibDataSvc;
    StatusCode sc = service("CalibDataSvc", pCalibDataSvc, true);
    if ( !sc.isSuccess() ) {
        log << MSG::ERROR 
            << "Could not get IDataProviderSvc interface of CalibXmlCnvSvc" 
            << endreq;
        return ;
    } else {
        log << MSG::DEBUG 
            << "Retrieved IDataProviderSvc interface of CalibXmlCnvSvc" 
            << endreq;
    }
 
    std::string fullPath = "/Calib/DedxSim";
    CalibDataSvc* m_calibDataSvc;
    m_calibDataSvc = dynamic_cast<CalibDataSvc*>(pCalibDataSvc);
    SmartDataPtr<CalibData::DedxSimData> pDedxSimData(m_calibDataSvc, fullPath);
    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" << endreq;
}