DedxCurSvc Class Reference

#include <DedxCurSvc.h>

Inheritance diagram for DedxCurSvc:

Public Member Functions
	DedxCurSvc (const std::string &name, ISvcLocator *svcloc)
	~DedxCurSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvUnknown)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &)
void	getDedxCurveInfo ()
const double	getCurve (int i)
const double	getSigma (int i)
const int	getCurveSize ()
const int	getSigmaSize ()
virtual const double	getSigma (int i)=0
virtual const double	getCurve (int i)=0
virtual const int	getCurveSize ()=0
virtual const int	getSigmaSize ()=0

Additional Inherited Members
Static Public Member Functions inherited from IDedxCurSvc
static const InterfaceID &	interfaceID ()

Detailed Description

Definition at line 16 of file DedxCurSvc.h.

Constructor & Destructor Documentation

DedxCurSvc()

DedxCurSvc::DedxCurSvc	(	const std::string &	name,
		ISvcLocator *	svcloc
	)

Definition at line 32 of file DedxCurSvc.cxx.

                                                                   : Service(name, svcloc) {
declareProperty("BossRelease",m_bossRelease="default");
declareProperty("BossVer",m_sftver="default");
declareProperty("ParVer",m_calParVer="default");
declareProperty("DbStatus",m_dbStatus="OK");
declareProperty("Type",m_type="default");
}

~DedxCurSvc()

DedxCurSvc::~DedxCurSvc

(

)

Definition at line 98 of file DedxCurSvc.cxx.

                       {
}

Member Function Documentation

finalize()

StatusCode DedxCurSvc::finalize ( )

virtual

Definition at line 91 of file DedxCurSvc.cxx.

                                  {
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << name() << ": End of Run" << endreq;
  return StatusCode::SUCCESS;
}

getCurve()

const double DedxCurSvc::getCurve ( int  i )

inlinevirtual

Implements IDedxCurSvc.

Definition at line 31 of file DedxCurSvc.h.

{return m_curve[i];}

getCurveSize()

const int DedxCurSvc::getCurveSize ( )

inlinevirtual

Implements IDedxCurSvc.

Definition at line 33 of file DedxCurSvc.h.

{return m_curve_size;}

getDedxCurveInfo()

void DedxCurSvc::getDedxCurveInfo

(

)

Definition at line 117 of file DedxCurSvc.cxx.

                                 {
  MsgStream log(messageService(), name());
  SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc,"/Event/EventHeader");
  int  run = eventHeader->runNumber();
  if(run<0){
    log << MSG::INFO << "This data is the MC sample with the Run Number: " << run << endreq;
    if(m_type=="Sim"){
      run=-run;
      log << MSG::INFO << "Reverse the sign of Run Number" << endreq;
    }
  }
  //  unsigned long *lengths;
  MYSQL_RES *res_set;
  MYSQL_ROW row;
  char stmt1[400];
 
  std::string sftver;
  
  if(m_sftver == "default" && m_bossRelease == "default")
  {
    log << MSG::FATAL << " BossRelease and BossVer for DedxCurve not set!" << endreq;
    exit(1);
  }
 
  if( (run<m_runFromMax || run>m_runToMin) ) 
  { 
    if(m_sftver == "default")
    {
      sprintf(stmt1,"select RunFrom, RunTo, SftVer, ParVer from CalVtxLumVer where BossRelease = '%s' and RunFrom <= %d and RunTo >= %d and DataType = 'DedxCurve'", m_bossRelease.c_str(), run, run);
      std::cout<<stmt1<<std::endl;
      DatabaseRecordVector rest;
      int row_no = m_dbsvc->query("offlinedb",stmt1,rest);
      if(row_no==0) {
        log << MSG::FATAL << "can not find result for DedxCur with: " << stmt1 << endreq;
        exit(1);
      }
      if(row_no>1) {
        log << MSG::FATAL << "find more than 1 results for DedxCur with: " << stmt1 << endreq;
        exit(1);
      }
      DatabaseRecord* recordst = rest[0];
      sftver = recordst->GetString("SftVer");
      m_calParVer = recordst->GetString("ParVer");
      m_runFromMax = atoi((recordst->GetString("RunFrom")).c_str());
      m_runToMin = atoi((recordst->GetString("RunTo")).c_str());
      log << MSG::INFO << "get from CalVtxLumVer, m_runFromMax: "<< m_runFromMax << " m_runToMin: "<<m_runToMin<< endreq;
    }
    else
      sftver = m_sftver;
  
    if(m_calParVer!="default")
      sprintf(stmt1,"select RunFrom, RunTo, DedxCurvePar,DedxSigmaPar from DedxCurvePar where SftVer = '%s' and RunFrom <= %d and RunTo >= %d and DedxCurveParVer = %s and Status='%s'",sftver.c_str(),run,run,m_calParVer.c_str(), m_dbStatus.c_str());
    else
      sprintf(stmt1,"select RunFrom, RunTo, DedxCurvePar,DedxSigmaPar from DedxCurvePar where SftVer = '%s' and RunFrom <= %d and RunTo >= %d and Status='%s' order by DedxCurveParVer desc",sftver.c_str(),run,run, m_dbStatus.c_str());
    std::cout<<stmt1<<std::endl; 
    
    DatabaseRecordVector res;
    int row_no  = m_dbsvc->query("offlinedb", stmt1, res);
    if(row_no==0){
     log << MSG::FATAL << "can not find result for DedxCur with: "<< stmt1 << endreq;
     exit(1);
    }
    if(row_no>1){
      log << MSG::FATAL << "find more than 1 results for DedxCur with: "<< stmt1 << endreq;
      exit(1);
    }
    
    DatabaseRecord* records = res[0];
    int runFrom = atoi((records->GetString("RunFrom")).c_str());
    int runTo = atoi((records->GetString("RunTo")).c_str());
    if(m_sftver=="default")
    {
      m_runFromMax = runFrom>m_runFromMax? runFrom:m_runFromMax;
      m_runToMin = runTo<m_runToMin? runTo:m_runToMin;
    }
    else
    {
      m_runFromMax = runFrom;
      m_runToMin = runTo;
    }
    log << MSG::INFO << "get from DedxCurvePar runFrom: "<< runFrom << " runTo: "<<runTo<< endreq;
    log << MSG::INFO << "m_runFromMax: "<< m_runFromMax << " m_runToMin: "<<m_runToMin<< endreq;
 
    TBuffer *buf1 = new TBufferFile(TBuffer::kRead);
    buf1->SetBuffer((*records)["DedxCurvePar"],327680,kFALSE);
    TBuffer *buf2 = new TBufferFile(TBuffer::kRead);
    buf2->SetBuffer((*records)["DedxSigmaPar"],327680,kFALSE);
    
    
    TTree* curvetree = new TTree(); 
    curvetree->Streamer(*buf1); 
    
    TTree* sigmatree = new TTree(); 
    sigmatree->Streamer(*buf2);
    
    double curve[50];
    double sigma[50];
    int SigmaSize;
    int CurveSize;
    curvetree -> SetBranchAddress("curve", curve);
    curvetree -> SetBranchAddress("CurveSize", &CurveSize);
    sigmatree -> SetBranchAddress("sigma", sigma);
    sigmatree -> SetBranchAddress("SigmaSize", &SigmaSize);
    
    Int_t nentries_curve = (Int_t)curvetree ->GetEntries();
    Int_t nentries_sigma = (Int_t)sigmatree ->GetEntries();
    curvetree -> GetEntry(0);
    if(CurveSize>50){
      log << MSG::ERROR<<"CurveSize larger than designed number"<<endreq;
      exit(1);}  
    for(int i=0; i<CurveSize;i++){
        m_curve[i]=curve[i];
       }
    
    sigmatree-> GetEntry(0);
    if(SigmaSize>50){
      log << MSG::ERROR<<"SigmaSize larger than designed number"<<endreq;
      exit(1);}
    for(int i=0; i<SigmaSize;i++){
        m_sigma[i]=sigma[i];}
      m_curve_size=CurveSize;
      m_sigma_size=SigmaSize;
     delete curvetree;
     delete sigmatree;
  }
   return;
}

Referenced by handle().

getSigma()

const double DedxCurSvc::getSigma ( int  i )

inlinevirtual

Implements IDedxCurSvc.

Definition at line 32 of file DedxCurSvc.h.

{return m_sigma[i];}

getSigmaSize()

const int DedxCurSvc::getSigmaSize ( )

inlinevirtual

Implements IDedxCurSvc.

Definition at line 34 of file DedxCurSvc.h.

{return m_sigma_size;}

handle()

void DedxCurSvc::handle

(

const Incident &

inc

)

Definition at line 103 of file DedxCurSvc.cxx.

                                          {
  MsgStream log( messageService(), name() );
  log << MSG::DEBUG << "handle: " << inc.type() << endreq;
 
  SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc,"/Event/EventHeader");
  int  run = eventHeader->runNumber();
 
  if ( inc.type() == "NewRun" ){
    log << MSG::DEBUG << "New Run:" << run << endreq;
    getDedxCurveInfo();
  }
}

initialize()

StatusCode DedxCurSvc::initialize ( )

virtual

Definition at line 51 of file DedxCurSvc.cxx.

                                    {
  
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << name() << ": Start of run initialisation" << 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("DatabaseSvc",m_dbsvc,true);
  if (sc .isFailure() ) {
    log << MSG::ERROR << "Unable to find DatabaseSvc " << endreq;
    return sc;
  }
 
 
  sc = serviceLocator()->service("EventDataSvc", m_eventSvc, true);
  if (sc .isFailure() ) {
    log << MSG::ERROR << "Unable to find EventDataSvc " << endreq;
    return sc;
  }
  
  StatusCode scc;
  log<<MSG::INFO << "setProperties()" << endreq;
  scc = setProperties();
 
  m_runFromMax=0;
  m_runToMin=0;
 
  return sc;
}  

queryInterface()

StatusCode DedxCurSvc::queryInterface ( const InterfaceID &  riid, void **  ppvUnknown  )

virtual

Definition at line 41 of file DedxCurSvc.cxx.

                                                                                   {
  
  if ( IID_IDedxCurSvc.versionMatch(riid) ) { 
    *ppvInterface = static_cast<IDedxCurSvc*> (this); 
  } 
  else return Service::queryInterface(riid, ppvInterface) ; 
  return StatusCode::SUCCESS;
}

---

The documentation for this class was generated from the following files:

• DedxCurSvc.h
• DedxCurSvc.cxx