MeasuredEcmsSvc Class Reference

#include <MeasuredEcmsSvc.h>

Inheritance diagram for MeasuredEcmsSvc:

Public Member Functions
	MeasuredEcmsSvc (const std::string &name, ISvcLocator *svcLoc)
virtual	~MeasuredEcmsSvc ()
virtual StatusCode	queryInterface (const InterfaceID &riid, void **ppvIF)
virtual StatusCode	initialize ()
virtual StatusCode	finalize ()
void	handle (const Incident &)
bool	isRunNoValid ()
double	getEcms ()
double	getBeamEnergy ()
double	getEcmsError ()
double	getAverageEcms ()
double	getAverageEcmsError ()
double	getSampleName ()
double	getPx ()
double	getPy ()
double	getPz ()
double	getBeamEnergyError ()
double	getEcms (int runNo)
double	getBeamEnergy (int runNo)
std::string	getRecord (int runNo, char *tab, char *col, char *min_col, char *max_col)
Public Member Functions inherited from IMeasuredEcmsSvc
virtual	~IMeasuredEcmsSvc ()
virtual bool	isRunNoValid ()=0
virtual double	getEcms ()=0
virtual double	getBeamEnergy ()=0
virtual double	getEcmsError ()=0
virtual double	getAverageEcms ()=0
virtual double	getAverageEcmsError ()=0
virtual double	getSampleName ()=0
virtual double	getPx ()=0
virtual double	getPy ()=0
virtual double	getPz ()=0
virtual double	getEcms (int runNo)=0
virtual std::string	getRecord (int runNo, char *tab, char *col, char *min_col, char *max_col)=0

Public Attributes
bool	m_readOneTime
int	m_runFrom
int	m_runTo
std::map< int, int >	m_mapBeamE
int	m_sele

Protected Member Functions
void	readDBTable (int runNo)
int	triggerDBTable (int runNo)
void	readDBTable (int runFrom, int runTo)
virtual void	readDBTable (int runNo)=0
virtual int	triggerDBTable (int runNo)=0

Friends
class	CnvFactory< MeasuredEcmsSvc >

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

Detailed Description

Class to read the data information from database, such as ECMS, ECMS Error ... ( beam_energy = ECMS / 2.0). Access to almost all the XYZ, 3773, and 4180 data.

Definition at line 24 of file MeasuredEcmsSvc.h.

Constructor & Destructor Documentation

MeasuredEcmsSvc()

MeasuredEcmsSvc::MeasuredEcmsSvc	(	const std::string &	name,
		ISvcLocator *	svcLoc
	)

Constructor

Definition at line 15 of file MeasuredEcmsSvc.cxx.

                                                                           : Service(name, svcLoc)
{
  // declare properties
  declareProperty("ReadOneTime", m_readOneTime=false);
  declareProperty("RunFrom", m_runFrom=8093);
  declareProperty("RunTo", m_runTo=9025);
  m_runflag = -1;
}

~MeasuredEcmsSvc()

MeasuredEcmsSvc::~MeasuredEcmsSvc ( )

virtual

Destructor

Definition at line 25 of file MeasuredEcmsSvc.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "destructor" << endreq;
}

Member Function Documentation

finalize()

StatusCode MeasuredEcmsSvc::finalize ( )

virtual

Finalize

Definition at line 74 of file MeasuredEcmsSvc.cxx.

{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "finalize begin here" << endreq;
 
  StatusCode status = Service::finalize();
  if (! status.isSuccess()) {
    log << MSG::ERROR << "can not finalize" << endreq;
    return status;
  }
 
  log << MSG::INFO << "finalized successfully" << endreq;
  return StatusCode::SUCCESS;
}

getAverageEcms()

double MeasuredEcmsSvc::getAverageEcms ( )

inlinevirtual

Return average Ecms

Implements IMeasuredEcmsSvc.

Definition at line 63 of file MeasuredEcmsSvc.h.

{ return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("Ave_Ecms") : 0.; }

getAverageEcmsError()

double MeasuredEcmsSvc::getAverageEcmsError ( )

inlinevirtual

Return average Ecms error

Implements IMeasuredEcmsSvc.

Definition at line 66 of file MeasuredEcmsSvc.h.

{ return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("Ave_Ecms_err") : 0.; }

getBeamEnergy() [1/2]

double MeasuredEcmsSvc::getBeamEnergy ( )

inlinevirtual

Return beam energy

Implements IMeasuredEcmsSvc.

Definition at line 57 of file MeasuredEcmsSvc.h.

{ return (double)getEcms() / 2.00; }

getBeamEnergy() [2/2]

double MeasuredEcmsSvc::getBeamEnergy ( int  runNo )

inline

Return beam energy depend on runNo only. Here RunNo. can be set independ on event.

Definition at line 91 of file MeasuredEcmsSvc.h.

{ return (double)getEcms(runNo) / 2.00; }

getBeamEnergyError()

double MeasuredEcmsSvc::getBeamEnergyError ( )

inline

Return beam energy error

Definition at line 81 of file MeasuredEcmsSvc.h.

{ return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("beam_energy_err") : 0.; }

getEcms() [1/2]

double MeasuredEcmsSvc::getEcms ( )

virtual

Return Ecms ( Ecms = beam_energy * 2 )

Implements IMeasuredEcmsSvc.

Definition at line 130 of file MeasuredEcmsSvc.cxx.

{
    switch (m_tableOption) {
        case 1:
            return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("Ecms") : 0.00;
        break;
 
        case 2:
            return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("beam_energy") * 2.00 : 0.00;
        break;
 
        default:
            if (m_rowNo1 != 0) return (double)m_dbRecordVector[m_sele]->GetDouble("Ecms");
                        else if (m_rowNo2 != 0) return (double)m_dbRecordVector[m_sele]->GetDouble("beam_energy") * 2.00;
                        else return 0.00;
        break;
    }
}

Referenced by getBeamEnergy().

getEcms() [2/2]

double MeasuredEcmsSvc::getEcms ( int  runNo )

virtual

Return Ecms depend on runNo only. Here RunNo. can be set independ on event.

Implements IMeasuredEcmsSvc.

Definition at line 151 of file MeasuredEcmsSvc.cxx.

{
    m_runNo = std::abs(runNo);
 
    readDBTable(runNo);
 
    switch (m_tableOption) {
    case 1:
      return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("Ecms") : 0.00;
    break;
 
    case 2:
      return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("beam_energy") * 2.00 : 0.00;
    break;
 
    default:
      if (m_rowNo1 != 0) return (double)m_dbRecordVector[m_sele]->GetDouble("Ecms");
            else if (m_rowNo2 != 0) return (double)m_dbRecordVector[m_sele]->GetDouble("beam_energy") * 2.00;
            else return 0.00;
    break;
  }
}

getEcmsError()

double MeasuredEcmsSvc::getEcmsError ( )

inlinevirtual

Return error for Ecms

Implements IMeasuredEcmsSvc.

Definition at line 60 of file MeasuredEcmsSvc.h.

{ return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("Ecms_err") : 0.; }

getPx()

double MeasuredEcmsSvc::getPx ( )

inlinevirtual

Return x momentum

Implements IMeasuredEcmsSvc.

Definition at line 72 of file MeasuredEcmsSvc.h.

{ return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("px") : 0.011; }

getPy()

double MeasuredEcmsSvc::getPy ( )

inlinevirtual

Return y momentum

Implements IMeasuredEcmsSvc.

Definition at line 75 of file MeasuredEcmsSvc.h.

{ return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("py") : 0.; }

getPz()

double MeasuredEcmsSvc::getPz ( )

inlinevirtual

Return z momentum

Implements IMeasuredEcmsSvc.

Definition at line 78 of file MeasuredEcmsSvc.h.

{ return (m_rowNo2 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("pz") : 0.; }

getRecord()

std::string MeasuredEcmsSvc::getRecord ( int  runNo, char *  tab, char *  col, char *  min_col, char *  max_col  )

virtual

Return the string depend on the parameters

Parameters

runNo	runNo. input (independent)
tab	table name in database, such as MeasuredEcms2, RunParams
col	the variable name which you want to read, such as Ecms
min_col	the down limits variable name, such as Run_From
max_col	the upper limits variable name, such as Run_To

Implements IMeasuredEcmsSvc.

Definition at line 331 of file MeasuredEcmsSvc.cxx.

{
    std::string strRecord;
    m_isRunNoValid = false;
 
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "MeasuredEcmsSvc begin read database " << tab << "-" << col << endreq;
 
  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[500];
  snprintf(stmt, 1024, "select %s from %s where %s <= %d and %s >= %d", col, tab, min_col, runNum, max_col, runNum);
 
  DatabaseRecordVector res;
  int row_no = m_dbSvc->query("offlinedb", stmt, res);
 
  if(row_no == 0) {
        log << MSG::ERROR << "RowNo can not be 0 for " << tab << ", failed to read DB for " << runNo << endreq;
    }
  else {
    log << MSG::INFO << "MeasuredEcmsSvc first read the runNo, read the " << tab << " database " << runNo << endreq;
    strRecord = res[m_sele]->GetString(col);
        m_isRunNoValid = true;
  } 
 
    return strRecord;
}

getSampleName()

double MeasuredEcmsSvc::getSampleName ( )

inlinevirtual

Return sample name

Implements IMeasuredEcmsSvc.

Definition at line 69 of file MeasuredEcmsSvc.h.

{ return (m_rowNo1 != 0) ? (double)m_dbRecordVector[m_sele]->GetDouble("sample") : 0.; }

handle()

void MeasuredEcmsSvc::handle

(

const Incident &

inc

)

Handle. Process after initialization

Definition at line 101 of file MeasuredEcmsSvc.cxx.

{
  MsgStream log(messageService(), name());
  log << MSG::DEBUG << "handle: " << inc.type() << endreq;
 
    SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc, "/Event/EventHeader");
    int runNo = eventHeader->runNumber();
    m_runNo = std::abs(runNo);
 
  if (inc.type() == "NewRun") {
    log << MSG::INFO << "RunNo is new, call the readMeasuredEcms again" << endreq;
        if(!m_readOneTime) {
           m_sele=0;
       readDBTable(runNo);
        } else {
           SmartDataPtr<Event::EventHeader> eventHeader(m_eventSvc,"/Event/EventHeader");
           int run = eventHeader->runNumber();
           if(run<0) run=-run;
           if(m_mapBeamE[run]>=0) {
              m_isRunNoValid = true;
              m_sele=m_mapBeamE[run];
           } else {
              std::cout<<"MeasuredEcmsSvc, could not get Ecm infor in handle new run"<<std::endl;
           }
        }
  }
}

initialize()

StatusCode MeasuredEcmsSvc::initialize ( )

virtual

Initialize

Definition at line 32 of file MeasuredEcmsSvc.cxx.

{
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "initialize begin here" << endreq;
 
  StatusCode status = Service::initialize();
  if (! status.isSuccess()) {
    log << MSG::ERROR << "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<< "unable to find DatabaseSvc" << endreq;
    return status;
  }
 
  status = serviceLocator()->service("EventDataSvc", m_eventSvc, true);
  if (! status.isSuccess()) {
    log << MSG::ERROR << "unable to find EventDataSvc " << endreq;
    return status;
  }
 
  if(m_readOneTime){
    if(m_runFrom>=8093){
        readDBTable(m_runFrom, m_runTo);
    }
    else
      std::cout<<"MeasuredEcmsSvc, invalid RunFrom, RunFrom should be >=8093"<<std::endl;
  }
 
  m_sele = 0;
  log << MSG::INFO << "initialized successfully" << endreq;
  return StatusCode::SUCCESS;
}

isRunNoValid()

bool MeasuredEcmsSvc::isRunNoValid ( )

inlinevirtual

Return true or false for runNo

Implements IMeasuredEcmsSvc.

Definition at line 51 of file MeasuredEcmsSvc.h.

{ return m_isRunNoValid; }

queryInterface()

StatusCode MeasuredEcmsSvc::queryInterface ( const InterfaceID &  riid, void **  ppvIF  )

virtual

Query interface

Definition at line 90 of file MeasuredEcmsSvc.cxx.

{
  if (IMeasuredEcmsSvc::interfaceID().versionMatch(riid)) *ppvIF = dynamic_cast<IMeasuredEcmsSvc*>(this);
  else return Service::queryInterface(riid, ppvIF);
 
  addRef();
 
  return StatusCode::SUCCESS;
}

readDBTable() [1/2]

void MeasuredEcmsSvc::readDBTable ( int  runFrom, int  runTo  )

protected

Definition at line 247 of file MeasuredEcmsSvc.cxx.

{
   MsgStream log(msgSvc(), name());
   log << MSG::INFO << "begin read database MeasuredEcmsSvc from " << m_runFrom << " to " << m_runTo << endreq;
   m_dbRecordVector.clear();
 
   m_rowNo1 = 0;
   m_rowNo2 = 0;
   m_isRunNoValid = false;
   char stmt[400];
 
   if ((runFrom >= 23463 && runFrom <= 38140) || (runFrom >= 47543 && runFrom <= 51498)) {
      if ((runTo >= 23463 && runTo <= 38140) || (runTo >= 47543 && runTo <= 51498)) {
         log << MSG::INFO << "read table MeasuredEcms2 for " << runFrom <<":" << runTo << endreq;
     snprintf(stmt, 1024, "select sample, RunFrom, RunTo, Ecms, Ecms_err, Ave_Ecms, Ave_Ecms_err from MeasuredEcms2 where RunFrom >= %d and RunTo <= %d", runFrom, runTo);
         m_tableOption = 1;
      }
   } 
   else if ((runFrom >= 11414 && runFrom <= 23454) || (runFrom >= 43716 && runFrom <= 47066)) {
      if ((runTo >= 11414 && runTo <= 23454) || (runTo >= 43716 && runTo <= 47066)) {
         log << MSG::INFO << "read table RunParams for " << runFrom <<":" << runTo << endreq;
     snprintf(stmt, 1024, "select run_number, beam_energy, px, py, pz, beam_energy_err from RunParams where run_number >= %d and run_number <= %d", runFrom, runTo);
         m_tableOption = 2;
      }
   }
   else {
         m_tableOption = -1;
   }   
 
   if (m_tableOption == 1 || m_tableOption == 2) {
      int rowNo = m_dbSvc->query("offlinedb", stmt, m_dbRecordVector);
      if (rowNo == 0) { 
     log << MSG::ERROR << "RowNo can not be 0, failed to read DB for " << runFrom <<":" << runTo << endreq; 
         cout << "RowNo can not be 0, failed to read DB for " << runFrom <<":" << runTo << endl;
      } else {
     if (m_tableOption == 1) m_rowNo1 = 1;
     if (m_tableOption == 2) m_rowNo2 = 1;
     log << MSG::INFO << "initialize the DatabaseRecordVector OK for " << runFrom <<":" << runTo << endreq;
         cout << "MeasuredEcmsSvc: initialize the DatabaseRecordVector OK for " << runFrom <<":" << runTo << ", Nrecord = " << rowNo << "." << endl;
     m_isRunNoValid = true;
      }
      if (m_rowNo2==1) {
         for ( int i=0; i<rowNo; i++ ) {
            int runnumber = m_dbRecordVector[i]->GetInt("run_number");
            m_mapBeamE[runnumber] = i;
         }
      }
      if (m_rowNo1==1) {
         for ( int i=0; i<rowNo; i++ ) {
            int run0 = m_dbRecordVector[i]->GetInt("RunFrom");
            int run1 = m_dbRecordVector[i]->GetInt("RunTo");
            for ( int j = run0; j<=run1; j++) {
               m_mapBeamE[j] = i;
            }
         }
      }
   }
}

readDBTable() [2/2]

void MeasuredEcmsSvc::readDBTable ( int  runNo )

protectedvirtual

Read database table

Implements IMeasuredEcmsSvc.

Definition at line 175 of file MeasuredEcmsSvc.cxx.

{
    if (runNo == m_runflag) return;
    m_runflag = runNo;
 
  MsgStream log(msgSvc(), name());
  log << MSG::INFO << "begin read database MeasuredEcmsSvc" << endreq;
 
    m_dbRecordVector.clear();
 
    m_rowNo1 = 0;
    m_rowNo2 = 0;
 
    m_isRunNoValid = false;
    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;
 
    triggerDBTable(runNum);
 
  char stmt[400];
 
    char stmt1[400];
    char stmt2[400];
 
    switch(m_tableOption) {
        case 1:
            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);
        break;
 
        case 2:
            snprintf(stmt, 1024, "select beam_energy, px, py, pz, beam_energy_err from RunParams where run_number = %d", runNum);
        break;
 
        default:
            snprintf(stmt1, 1024, "select sample, RunFrom, RunTo, Ecms, Ecms_err, Ave_Ecms, Ave_Ecms_err from MeasuredEcms2 where RunFrom <= %d and RunTo >= %d", runNum, runNum);
            snprintf(stmt2, 1024, "select beam_energy, px, py, pz, beam_energy_err from RunParams where run_number = %d", runNum);
        break;
    }
 
    if (m_tableOption == 1 || m_tableOption == 2) {
        int rowNo = m_dbSvc->query("offlinedb", stmt, m_dbRecordVector);
 
        if (rowNo == 0) { 
            log << MSG::ERROR << "RowNo can not be 0, failed to read DB for " << runNo << endreq; 
        }
        else {
            if (m_tableOption == 1) m_rowNo1 = 1;
            if (m_tableOption == 2) m_rowNo2 = 1;
            log << MSG::INFO << "initialize the DatabaseRecordVector OK for " << runNo << endreq;
            m_isRunNoValid = true;
        }
    }
    else {
        m_rowNo1 = m_dbSvc->query("offlinedb", stmt1, m_dbRecordVector);
        if (m_rowNo1 == 0) {
            m_dbRecordVector.clear();
            m_rowNo2 = m_dbSvc->query("offlinedb", stmt2, m_dbRecordVector);
        }
 
        if (m_rowNo1 == 0 && m_rowNo2 == 0) {
            log << MSG::ERROR << "RowNo can not be 0, failed to read DB for " << runNo << endreq;
        }
        else {
            log << MSG::INFO << "initialize the DatabaseRecordVector OK for " << runNo << endreq;
            m_isRunNoValid = true;
        }
    }
}

Referenced by getEcms(), handle(), and initialize().

triggerDBTable()

int MeasuredEcmsSvc::triggerDBTable ( int  runNo )

protectedvirtual

trigger based on runNo @option 1 MeasuredEcms2 table for XYZ (Charmonium data) @option 2 RunParams table for 3773, 4180 (Charm data)

Implements IMeasuredEcmsSvc.

Definition at line 307 of file MeasuredEcmsSvc.cxx.

{
    MsgStream log(msgSvc(), name());
    log << MSG::INFO << "triggerDBTable based on runNo" << endl;
 
    int runNum = std::abs(runNo);
 
  if ((runNum >= 23463 && runNum <= 38140) || (runNum >= 47543 && runNum <= 51498)) {
        log << MSG::INFO << "read table MeasuredEcms2 for " << runNo << endreq;
        m_tableOption = 1;
    }
 
    else if ((runNum >= 11414 && runNum <= 23454) || (runNum >= 43716 && runNum <= 47066)) {
        log << MSG::INFO << "read table RunParams for " << runNo << endreq;
        m_tableOption = 2;
    }
 
    else m_tableOption = -1;
 
    return m_tableOption;
}

Referenced by readDBTable().

Friends And Related Function Documentation

CnvFactory< MeasuredEcmsSvc >

friend class CnvFactory< MeasuredEcmsSvc >

friend

Definition at line 1 of file MeasuredEcmsSvc.h.

Member Data Documentation

m_mapBeamE

std::map<int,int> MeasuredEcmsSvc::m_mapBeamE

Definition at line 106 of file MeasuredEcmsSvc.h.

Referenced by handle(), and readDBTable().

m_readOneTime

bool MeasuredEcmsSvc::m_readOneTime

Definition at line 103 of file MeasuredEcmsSvc.h.

Referenced by handle(), initialize(), and MeasuredEcmsSvc().

m_runFrom

int MeasuredEcmsSvc::m_runFrom

Definition at line 104 of file MeasuredEcmsSvc.h.

Referenced by initialize(), MeasuredEcmsSvc(), and readDBTable().

m_runTo

int MeasuredEcmsSvc::m_runTo

Definition at line 105 of file MeasuredEcmsSvc.h.

Referenced by initialize(), MeasuredEcmsSvc(), and readDBTable().

m_sele

int MeasuredEcmsSvc::m_sele

Definition at line 107 of file MeasuredEcmsSvc.h.

Referenced by getAverageEcms(), getAverageEcmsError(), getBeamEnergyError(), getEcms(), getEcmsError(), getPx(), getPy(), getPz(), getRecord(), getSampleName(), handle(), and initialize().

---

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

• MeasuredEcmsSvc.h
• MeasuredEcmsSvc.cxx