MeasuredEcmsSvc.h

Go to the documentation of this file.

/***********************************************************************
 * BES III Software                                                    *
 * Copyright(C) - BES III Collaboration                                *
 *                                                                     *
 * Author: The BESIII Collaboration                                    *
 * Contributor: Lianjin Wu                                             *
 *                                                                     *
 ***********************************************************************/
 
#pragma once
#ifndef Utilities_MeasuredEcmsSvc_H
#  define Utilities_MeasuredEcmsSvc_H
 
#  include "MeasuredEcmsSvc/IMeasuredEcmsSvc.h"
 
class IDatabaseSvc;
 
template <class TYPE> class CnvFactory;
 
/** 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.
 */
// class MeasuredEcmsSvc : public Service, virtual public IMeasuredEcmsSvc, virtual public
// IIncidentListener {
class MeasuredEcmsSvc : public extends<Service, IMeasuredEcmsSvc>,
                        virtual public IIncidentListener {
  friend class CnvFactory<MeasuredEcmsSvc>;
 
public:
  /** Constructor */
  MeasuredEcmsSvc( const std::string& name, ISvcLocator* svcLoc );
 
  /** Destructor */
  virtual ~MeasuredEcmsSvc();
 
  /** Query interface */
  // virtual StatusCode queryInterface(const InterfaceID& riid, void** ppvIF);
 
  /** Initialize */
  virtual StatusCode initialize();
 
  /** Finalize */
  virtual StatusCode finalize();
 
  /** Handle.
   * Process after initialization
   */
  void handle( const Incident& );
 
  /** Return true or false for runNo */
  bool isRunNoValid() { return m_isRunNoValid; }
 
  /** Return Ecms ( Ecms = beam_energy * 2 ) */
  double getEcms();
 
  /** Return beam energy */
  double getBeamEnergy() { return (double)getEcms() / 2.00; }
 
  /**   Return error for Ecms */
  double getEcmsError() {
    return ( m_rowNo1 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "Ecms_err" ) : 0.;
  }
 
  /** Return average Ecms */
  double getAverageEcms() {
    return ( m_rowNo1 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "Ave_Ecms" ) : 0.;
  }
 
  /** Return average Ecms error */
  double getAverageEcmsError() {
    return ( m_rowNo1 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "Ave_Ecms_err" )
                             : 0.;
  }
 
  /** Return sample name */
  double getSampleName() {
    return ( m_rowNo1 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "sample" ) : 0.;
  }
 
  /** Return x momentum */
  double getPx() {
    return ( m_rowNo2 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "px" ) : 0.011;
  }
 
  /** Return y momentum */
  double getPy() {
    return ( m_rowNo2 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "py" ) : 0.;
  }
 
  /** Return z momentum */
  double getPz() {
    return ( m_rowNo2 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "pz" ) : 0.;
  }
 
  /** Return beam energy error */
  double getBeamEnergyError() {
    return ( m_rowNo2 != 0 ) ? (double)m_dbRecordVector[m_sele]->GetDouble( "beam_energy_err" )
                             : 0.;
  }
 
  /** Return Ecms depend on runNo only.
   * Here RunNo. can be set independ on event.
   */
  double getEcms( int runNo );
 
  /** Return beam energy depend on runNo only.
   * Here RunNo. can be set independ on event.
   */
  double getBeamEnergy( int runNo ) { return (double)getEcms( runNo ) / 2.00; }
 
  /** Return the string depend on the parameters
   * @param runNo                           runNo. input
   * (independent)
   * @param tab                             table name in database,
   * such as MeasuredEcms2, RunParams
   * @param col                             the variable name which
   * you want to read, such as Ecms
   * @param min_col                     the down limits variable name,
   * such as Run_From
   * @param max_col           the upper limits variable name, such as Run_To
   */
  std::string getRecord( int runNo, char* tab, char* col, char* min_col, char* max_col );
 
  bool               m_readOneTime;
  int                m_runFrom;
  int                m_runTo;
  std::map<int, int> m_mapBeamE;
  int                m_sele;
 
protected:
  /** Read database table */
  void readDBTable( int runNo );
 
  /** trigger based on runNo
   * @option 1          MeasuredEcms2 table for XYZ (Charmonium data)
   * @option 2          RunParams   table for 3773, 4180 (Charm data)
   */
  int triggerDBTable( int runNo );
 
  void readDBTable( int runFrom, int runTo );
 
private:
  IDataProviderSvc*    m_eventSvc;
  IDatabaseSvc*        m_dbSvc;
  DatabaseRecordVector m_dbRecordVector;
 
  int  m_runNo;        /**< runNo. read from event */
  int  m_tableOption;  /**< triggerDBTable for (1) MeasuredEcms2 and (2) RunParams */
  bool m_isRunNoValid; /**< estimate whether the runNo valid */
  int  m_rowNo1;       /**< read MeasuredEcms2 */
  int  m_rowNo2;       /**< read RunParams */
 
  int m_runflag;
};
 
#endif