CgemCalibFunSvc.h

Go to the documentation of this file.

#ifndef CGEMCALIBFUNSVC_H_
#define CGEMCALIBFUNSVC_H_
 
#include "GaudiKernel/IInterface.h"
#include "GaudiKernel/Kernel.h"
#include "GaudiKernel/Service.h"
#include "GaudiKernel/IIncidentListener.h"
 
#include "GaudiKernel/IDataProviderSvc.h"
#include "GaudiKernel/SmartDataPtr.h"
#include "GaudiKernel/DataSvc.h"
#include "CalibData/CalibModel.h"
/* #include "CalibData/Cgem/CgemCalibData.h" */
#include "CalibDataSvc/ICalibRootSvc.h"
 
#include "CgemGeomSvc/ICgemGeomSvc.h"
#include "CgemCalibFunSvc/ICgemCalibFunSvc.h"
#include "ReadCosmicRayData/CgemLUTReader.h"
 
#include <string>
#include <vector>
 
class CgemCalibFunSvc: public Service, virtual public ICgemCalibFunSvc, virtual public IIncidentListener{
public:
     CgemCalibFunSvc( const std::string& name, ISvcLocator* svcloc );
     ~CgemCalibFunSvc();
 
     virtual StatusCode queryInterface(const InterfaceID& riid, void** ppvUnknown);
     virtual StatusCode initialize();
     virtual StatusCode finalize();
 
     // Incident handler
     void handle(const Incident&);
 
     /* Access by users */
 
     /* xvFlag: 0 -> x, 1 -> v;  readoutMode: 0 -> charge centroid, 1 -> micro-TPC */
     double getSigma(int layer, int xvFlag, int readoutMode, double angle, double Q, double T) const;
 
     /* xvFlag: 0 -> x, 1 -> v; */
     double getTimeRising(int layer, int xvFlag, int sheet, int stripID, double Q=100., double z=0.) const;
 
     /* xvFlag: 0 -> x, 1 -> v; */
     double getTimeFalling(int layer, int xvFlag, int sheet, int stripID, double Q=100., double z=0.) const;
 
     double getTimeWalk(int layer, int xvFlag, int sheet, int stripID, double Q) const;
     double getTimeWalk(double Q, double threshold) const;
private:
     bool initCalibConst();
 
     int                m_run;
     bool               m_getRunNo;
     ICgemGeomSvc*      m_pCgemGeomSvc;
     IDataProviderSvc*  m_pCalDataSvc;
     ICalibRootSvc*     m_pRootSvc;
 
     static const int NLAYER = 3;
     static const int NXV = 2;
     static const int NSHEET = 2;
     static const int NTHRESHODS = 20;
 
     bool m_printInfo;
     std::string m_timeFitFile;
     std::string m_timeFitFilePath;
     double m_tRising[NLAYER][NXV][NSHEET];
     double m_tFalling[NLAYER][NXV][NSHEET];
 
     std::string m_timeWalkFileName;
     std::vector<double> m_tw;
     std::vector<double> m_charge;
     std::vector<double> m_threshold;
 
     int m_readLUT;
     std::string m_fileLUT;
     CgemLUTReader* m_cgemLUTReader;
     std::vector<double> m_valThreshold;
     std::vector<int> m_thIdFrom;
     std::vector<int> m_thIdEnd;
};
 
#endif /* CGEMCALIBFUNSVC_H_ */