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