MeasuredEcmsSvc Class Reference

#include <MeasuredEcmsSvc.h>

Inheritance diagram for MeasuredEcmsSvc:

Public Member Functions
	MeasuredEcmsSvc (const std::string &name, ISvcLocator *svcLoc)
virtual	~MeasuredEcmsSvc ()
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 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)

Friends
class	CnvFactory< MeasuredEcmsSvc >

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 25 of file MeasuredEcmsSvc.h.

Constructor & Destructor Documentation

MeasuredEcmsSvc()

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

Constructor

Definition at line 16 of file MeasuredEcmsSvc.cxx.

                                                                           : base_class(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 26 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 72 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 ( )

inline

Return average Ecms

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 ( )

inline

Return average Ecms error

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 ( )

inline

Return beam energy

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

(

)

Return Ecms ( Ecms = beam_energy * 2 )

Definition at line 127 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

)

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

Definition at line 148 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 ( )

inline

Return error for Ecms

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 ( )

inline

Return x momentum

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 ( )

inline

Return y momentum

Definition at line 75 of file MeasuredEcmsSvc.h.

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

getPz()

double MeasuredEcmsSvc::getPz ( )

inline

Return z momentum

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
	)

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

Definition at line 327 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 ( )

inline

Return sample name

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 99 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

Query interface Initialize

Definition at line 33 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 ( )

inline

Return true or false for runNo

Definition at line 51 of file MeasuredEcmsSvc.h.

{ return m_isRunNoValid; }

readDBTable() [1/2]

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

protected

Definition at line 243 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) || (runFrom >= 70522 && runFrom <= 73929)) {
      if ((runTo >= 11414 && runTo <= 23454) || (runTo >= 43716 && runTo <= 47066) || (runFrom >= 70522 && runFrom <= 73929)) {
         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 )

protected

Read database table

Definition at line 172 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 )

protected

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

Definition at line 303 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)|| (runNum >= 70522 && runNum <= 73929)) {
        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(), and initialize().

m_runFrom

int MeasuredEcmsSvc::m_runFrom

Definition at line 104 of file MeasuredEcmsSvc.h.

Referenced by initialize(), and readDBTable().

m_runTo

int MeasuredEcmsSvc::m_runTo

Definition at line 105 of file MeasuredEcmsSvc.h.

Referenced by initialize(), 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