/home/bes3soft/bes3soft/Boss/7.0.2/dist/7.0.2/Reconstruction/EmcRec/EmcRec-01-02-57/EmcRec/EmcRecParameter.h

Go to the documentation of this file.

//
//  Parameters for Emc Reconstruction
//  
//  No Parameter is allowed to be hard coded into code!
//  
//  Created by Zhe Wang, May 31, 2004
//
#ifndef EMC_REC_PARAMETER_A_H
#define EMC_REC_PARAMETER_A_H
#include <pthread.h>
#include <iostream>
#include <vector>

using namespace std;

class TGraph2DErrors;

class EmcRecParameter
{
  private:
    // Constructors and destructors 
    EmcRecParameter();
    ~EmcRecParameter();

  public:
    //static member functions
    static EmcRecParameter& GetInstance();
    static bool Exist();
    static void Kill();
    static void lock(){
      if(pthread_mutex_lock(&m_pthread_lock) != 0){
        std::cerr << "LOCK MUTEX_LOCK @ GZFSSTREAM" << std::endl;
      }
    };
    static void unlock(){
      if(pthread_mutex_unlock(&m_pthread_lock) != 0){
        std::cerr << "UNLOCK MUTEX_LOCK @ GZFSSTREAM" << std::endl;
      }
    };



  private:
    //static data members
    static EmcRecParameter* fpInstance;
    static pthread_mutex_t m_pthread_lock;


  public:
    //access to each parameter
    double ElectronicsNoiseLevel() const;
    double EThresholdSeed() const;
    double EThresholdCluster() const;
    double LogPosOffset() const;

    double TimeMin() const;
    double TimeMax() const;
    double MethodMode() const;
    double PosCorr() const;
    double DataMode() const;
    int ElecSaturation() const;

    double MoliereRadius() const;
    double LateralProfile() const;

    double ECorr(int n) const;
    double SigE(int n) const;
    double SigTheta(int n) const;
    double SigPhi(int n) const;

    double HitNb(int n) const;
    double ElecBias(int n) const;
    double SmCut(int n) const;

    double Peak(int n) const;

    double EastLogThetaPara(int n, int m) const;
    double WestLogThetaPara(int n, int m) const;

    double EastLogPhiPara(int n, int m) const;
    double WestLogPhiPara(int n, int m) const;


    double EastDataLogThetaPara(int n, int m) const;
    double WestDataLogThetaPara(int n, int m) const;


    double EastLinThetaPara(int n, int m) const;
    double WestLinThetaPara(int n, int m) const;

    double EastLinPhiPara(int n, int m) const;
    double WestLinPhiPara(int n, int m) const;

    double BarrPosDataCor(int ntheta, int nphi) const;
    double WestPosDataCor(int ntheta, int nphi) const;
    double EastPosDataCor(int ntheta, int nphi) const;

    double BarrPosMCCor(int ntheta, int nphi) const;
    double WestPosMCCor(int ntheta, int nphi) const;
    double EastPosMCCor(int ntheta, int nphi) const;

    double BarrLogThetaPara(int n, int m) const;
    double BarrLogPhiPara(int n, int m) const;
 
    double BarrLoglinThetaPara(int n, int m) const;
    double BarrLoglinPhiPara(int n, int m) const;


    double BarrLinThetaPara(int n, int m) const;
    double BarrLinPhiPara(int n, int m) const;
 
    double BarrShLogThetaPara(int n, int m) const;
    double BarrShLogPhiPara(int n, int m) const;
 
    double BarrShLinThetaPara(int n, int m) const;
    double BarrShLinPhiPara(int n, int m) const;

    double BarrDataLogThetaPara(int n, int m) const;

    inline bool DigiCalib() const { return digiCalib; }
    inline void SetDigiCalib(bool digi) { digiCalib=digi; }

    inline void SetTimeMin(double min) { fTimeMin=min; }
    inline void SetTimeMax(double max) { fTimeMax=max; }

    inline void SetMethodMode(double en) { fMethodMode=en; }
    inline void SetPosCorr(double en) { fPosCorr=en; }
    inline void SetDataMode(double en) { fDataMode=en; }    
 
    inline void SetElecSaturation(int IO){ fElecSaturation=IO;}



    inline void SetPeak(double e, int n) { peak[n]=e; }

    inline void SetEastLogThetaPara(double p,int n,int m) { eastLogThetaPara[n][m]=p; }
    inline void SetWestLogThetaPara(double p,int n,int m) { westLogThetaPara[n][m]=p; }

    inline void SetEastLogPhiPara(double p,int n,int m) { eastLogPhiPara[n][m]=p; }
    inline void SetWestLogPhiPara(double p,int n,int m) { westLogPhiPara[n][m]=p; }

    inline void SetEastDataLogThetaPara(double p,int n,int m) { eastDataLogThetaPara[n][m]=p; }
    inline void SetWestDataLogThetaPara(double p,int n,int m) { westDataLogThetaPara[n][m]=p; }


    inline void SetEastLinThetaPara(double p,int n,int m) { eastLinThetaPara[n][m]=p; }
    inline void SetWestLinThetaPara(double p,int n,int m) { westLinThetaPara[n][m]=p; }

    inline void SetEastLinPhiPara(double p,int n,int m) { eastLinPhiPara[n][m]=p; }
    inline void SetWestLinPhiPara(double p,int n,int m) { westLinPhiPara[n][m]=p; }


    inline void SetBarrLogThetaPara(double p,int n,int m) {barrLogThetaPara[n][m]=p; }
    inline void SetBarrLogPhiPara(double p,int n,int m) { barrLogPhiPara[n][m]=p; }
 
    inline void SetBarrLoglinThetaPara(double p,int n,int m) {barrLoglinThetaPara[n][m]=p; }
    inline void SetBarrLoglinPhiPara(double p,int n,int m) { barrLoglinPhiPara[n][m]=p; }

    inline void SetBarrLinThetaPara(double p,int n,int m) { barrLinThetaPara[n][m]=p; }
    inline void SetBarrLinPhiPara(double p,int n,int m) { barrLinPhiPara[n][m]=p; }
    
    inline void SetBarrShLogThetaPara(double p,int n,int m) { barrShLogThetaPara[n][m]=p; }
    inline void SetBarrShLogPhiPara(double p,int n,int m) { barrShLogPhiPara[n][m]=p; }
 
    inline void SetBarrShLinThetaPara(double p,int n,int m) { barrShLinThetaPara[n][m]=p; }
    inline void SetBarrShLinPhiPara(double p,int n,int m) { barrShLinPhiPara[n][m]=p; }

    inline void SetBarrPosDataCor(double p, int nphi,int ntheta) {barrPosDataCorPara[ntheta][nphi]=p; }
    inline void SetWestPosDataCor(double p, int nphi,int ntheta) {westPosDataCorPara[ntheta][nphi]=p; }
    inline void SetEastPosDataCor(double p, int nphi,int ntheta) {eastPosDataCorPara[ntheta][nphi]=p; }

    inline void SetBarrPosMCCor(double p, int nphi,int ntheta) {barrPosMCCorPara[ntheta][nphi]=p; }
    inline void SetWestPosMCCor(double p, int nphi,int ntheta) {westPosMCCorPara[ntheta][nphi]=p; }
    inline void SetEastPosMCCor(double p, int nphi,int ntheta) {eastPosMCCorPara[ntheta][nphi]=p; }

    inline void SetBarrDataLogThetaPara(double p,int n,int m) {barrDataLogThetaPara[n][m]=p; }

    inline std::string PositionMode1() const { return positionMode1; }
    inline std::string PositionMode2() const { return positionMode2; }
    void SetPositionMode(std::vector<std::string>& mode);

    double ECorrMC(double eg, double theid) const;
    double ErrMC(double eg, double theid) const;
    double E25min(int n) const;
    double E25max(int n) const;

    double e25min[28];
    double e25max[28];
  private:
    //each parameter
    // ElectronicsNoiseLevel
    double fElectronicsNoiseLevel;
    // Energy threshold for seed search
    double fEThresholdSeed;
    // Energy threshold for cluster search
    double fEThresholdCluster;
    // LogPosOffset
    double fLogPosOffset;

    // Time window left
    double fTimeMin;
    // Time window right
    double fTimeMax;
    //Position correction method
    double fMethodMode;

   //MC reconstruction for electronics saturation
   // 0: using Bhabha calibration constants
   // 1: using Emaxdata from data
   int fElecSaturation;

    //Position correction or not
    double fPosCorr;
    double fDataMode;
    // Moliere radius
    double fMoliereRadius;
    // Lateral profile
    double fLateralProfile;

    // Energy correction
    double eCorr[4];
    // Energy error
    double sigE[3];
    // Theta error
    double sigTheta[2];
    // Phi error
    double sigPhi[2];

    // Hit number
    double hitNb[3];
    // Correction of electronics bias
    double elecBias[5];
    // Cluster splitting cut with second moment
    double smCut[4];

    // Digi calibration
    bool digiCalib;

    // Shower energy correction VS theta
    double peak[56];

    double barrLogThetaPara[66][5];
    double barrLogPhiPara[66][5];

    double eastLogThetaPara[18][5];
    double westLogThetaPara[18][5];

    double eastLogPhiPara[3][5];
    double westLogPhiPara[3][5];


    double eastDataLogThetaPara[6][5];
    double westDataLogThetaPara[6][5];


    double eastLinThetaPara[6][5];
    double westLinThetaPara[6][5];

    double eastLinPhiPara[1][5];
    double westLinPhiPara[1][5];

    double barrLoglinThetaPara[22][5];
    double barrLoglinPhiPara[1][5];
 
    double barrLinThetaPara[66][5];
    double barrLinPhiPara[3][5];
 
    double barrShLogThetaPara[66][5];
    double barrShLogPhiPara[3][5];
 
    double barrShLinThetaPara[66][5];
    double barrShLinPhiPara[3][5];

    double barrDataLogThetaPara[22][5];

    double barrPosDataCorPara[44][120];
    double westPosDataCorPara[6][100];
    double eastPosDataCorPara[6][100];

    double barrPosMCCorPara[44][120];
    double westPosMCCorPara[6][100];
    double eastPosMCCorPara[6][100];

    std::string positionMode1;
    std::string positionMode2;

    // Shower energy correction
    TGraph2DErrors *dt;
    // Energy error
    TGraph2DErrors *dtErr;
};

#endif // EMC_REC_PARAMETER_A_H

---