MeasuredEcmsSvc.cxx

Go to the documentation of this file.

//<[email protected]> 2015.11.28 Modified
#include "MeasuredEcmsSvc/MeasuredEcmsSvc.h"
using namespace std;
MeasuredEcmsSvc::MeasuredEcmsSvc(const std::string& name, ISvcLocator* svcLoc):Service(name, svcLoc){
    declareProperty("RunParamsMin", RunParamsMin=11414);
    declareProperty("RunParamsMax", RunParamsMax=23454);
    declareProperty("MeasuredEcmsMin", MeasuredEcmsMin=23463);
    declareProperty("MeasuredEcmsMax", MeasuredEcmsMax=38140);
}
 
MeasuredEcmsSvc::~MeasuredEcmsSvc(){}
 
StatusCode MeasuredEcmsSvc::initialize(){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc initialize begin here"<<endreq;
  StatusCode status=Service::initialize();
  if(!status.isSuccess()){
    log<<MSG::ERROR<<"MeasuredEcmsSvc Can not initialize"<<endreq;
    return status;
  }
  IIncidentSvc* incsvc;
  status=service("IncidentSvc", incsvc);
  int priority = 100;
  if(status.isSuccess()) incsvc->addListener(this, "NewRun", priority);
  status=serviceLocator()->service("DatabaseSvc", m_dbSvc, true);
  if(!status.isSuccess()){
    log<<MSG::ERROR<<"MeasuredEcmsSvc Unable to find DatabaseSvc"<<endreq;
    return status;
  }
  status=serviceLocator()->service("EventDataSvc", m_eventSvc, true);
  if(!status.isSuccess()){
    log<<MSG::ERROR<<"MeasuredEcmsSvc Unable to find EventDataSvc "<<endreq;
    return status;
  }
  log<<MSG::INFO<<"MeasuredEcmsSvc initialized successfully"<<endreq;
  return StatusCode::SUCCESS;
}
 
StatusCode MeasuredEcmsSvc::finalize(){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc finalize begin here"<<endreq;
  StatusCode status= Service::finalize();
  if(!status.isSuccess()){
    log<<MSG::ERROR<<"MeasuredEcmsSvc Can not finalize"<<endreq;
    return status;
  }
  log<<MSG::INFO<<"MeasuredEcmsSvc finalized successfully"<<endreq;
  return StatusCode::SUCCESS;
}
 
StatusCode MeasuredEcmsSvc::queryInterface(const InterfaceID& riid, void** ppvIF){
  if(IMeasuredEcmsSvc::interfaceID().versionMatch(riid)) *ppvIF = dynamic_cast<IMeasuredEcmsSvc*>(this);
  else return Service::queryInterface(riid, ppvIF);
  addRef();
  return StatusCode::SUCCESS;
}
 
void MeasuredEcmsSvc::handle(const Incident& inc){
  MsgStream log(messageService(), name());
  log<<MSG::DEBUG<<"handle: "<<inc.type()<<endreq;
    m_callRead=false; m_recallRead=true;
  if(inc.type()=="NewRun"){
    log<<MSG::INFO<<"RunNo is new, call the readMeasuredEcms again"<<endreq;
    m_callRead=true;
  }
}
 
void MeasuredEcmsSvc::readMeasuredEcms(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc begin read database MeasuredEcmsSvc"<<endreq;
  m_sampleName=0; m_ecms=0; m_ecmsErr=0;
  m_aveEcms=0; m_aveEcmsErr=0; m_runFrom=0;
  m_runTo=0; m_runNoValid=false;
    m_px=0; m_py=0; m_pz=0;
  int runNum=std::abs(runNo);
  if(runNo<0) log<<MSG::INFO<<"MeasuredEcmsSvc for MC events"<<endreq;
  else log<<MSG::INFO<<"MeasuredEcmsSvc for DT events"<<endreq;
  char stmt[400];
  snprintf(stmt, 1024, "select sample, RunFrom, RunTo, Ecms, Ecms_err, Ave_Ecms, Ave_Ecms_err from MeasuredEcms2 where RunFrom <= %d and RunTo >= %d", runNum, runNum);
  DatabaseRecordVector res;
  int rowNo=m_dbSvc->query("offlinedb", stmt, res);
  if(rowNo==0){ log<<MSG::ERROR<<"RowNo can not be 0 for MeasuredEcmsSvc, failed to read DB"<<endreq; }
  else{
    DatabaseRecord* records=res[0];
    log<<MSG::INFO<<"MeasuredEcmsSvc first read the runNo, read the MeasuredEcms database"<<runNo<<endreq;
    m_runFrom=records->GetDouble("RunFrom");
    m_runTo=records->GetDouble("RunTo");
    m_sampleName=records->GetDouble("sample");
    m_ecms=records->GetDouble("Ecms");
    m_ecmsErr=records->GetDouble("Ecms_err");
    m_aveEcms=records->GetDouble("Ave_Ecms");
    m_aveEcmsErr=records->GetDouble("Ave_Ecms_err");
    m_runNoValid=true;
  }
}
 
void MeasuredEcmsSvc::readRunParams(int runNo){
    MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc begin read database RunParams"<<endreq;
    m_sampleName=0; m_ecms=0; m_ecmsErr=0;
  m_aveEcms=0; m_aveEcmsErr=0; m_runFrom=0;
  m_runTo=0; m_runNoValid=false;
  m_px=0; m_py=0; m_pz=0;
    int runNum=std::abs(runNo);
  if(runNo<0) log<<MSG::INFO<<"MeasuredEcmsSvc for MC events"<<endreq;
  else log<<MSG::INFO<<"MeasuredEcmsSvc for DT events"<<endreq;
    char stmt1[400];
  snprintf(stmt1, 1024, "select beam_energy, px, py, pz from RunParams where run_number = %d", runNum);
  DatabaseRecordVector res;
  int row_no = m_dbSvc->query("offlinedb", stmt1, res);
  if(row_no==0){log<<MSG::ERROR<<"RowNo can not be 0 for RunParams, failed to read DB"<<endreq;}
  else{
        DatabaseRecord* records = res[0];
        log<<MSG::INFO<<"MeasuredEcmsSvc first read the runNo, read the  RunParams database"<<runNo<<endreq;
    double bE=0;
    bE=records->GetDouble("beam_energy");
    m_px=records->GetDouble("px");
    m_py=records->GetDouble("py");
    m_pz=records->GetDouble("pz");
    m_ecms=bE*2.00;
        m_runNoValid=true;
    }
}
 
void MeasuredEcmsSvc::readDB(int runNo){
    if(std::abs(runNo)>=RunParamsMin&&std::abs(runNo)<=RunParamsMax) readRunParams(runNo);
    if(std::abs(runNo)>=MeasuredEcmsMin&&std::abs(runNo)<=MeasuredEcmsMax) readMeasuredEcms(runNo);
}
 
int MeasuredEcmsSvc::getSampleName(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the sample name"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_sampleName;
}
 
int MeasuredEcmsSvc::getRunFrom(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the RunFrom"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_runFrom;
}
 
int MeasuredEcmsSvc::getRunTo(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the RunTo"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_runTo;
}
 
double MeasuredEcmsSvc::getEcms(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the ecms"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_ecms;
}
 
double MeasuredEcmsSvc::getEcmsErr(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the error of ecms"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_ecmsErr;
}
 
double MeasuredEcmsSvc::getAveEcms(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the average ecms"<<endreq;
    if(m_recallRead&&m_callRead){
    readDB(runNo);
    m_recallRead=false;
  }
  return m_aveEcms;
}
 
double MeasuredEcmsSvc::getAveEcmsErr(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc read the error of the average ecms"<<endreq;
  if(m_recallRead&&m_callRead){
        readDB(runNo);
        m_recallRead=false;
    }
  return m_aveEcmsErr;
}
 
bool MeasuredEcmsSvc::isRunNoValid(int runNo){
  MsgStream log(msgSvc(), name());
  log<<MSG::INFO<<"MeasuredEcmsSvc judge the valid of runNo"<<endreq;
  if(m_callRead){
        readDB(runNo);
        m_recallRead=false;
    }
  return m_runNoValid;
}