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 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;
 
    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 SetPeak(double e, int n) { peak[n]=e; }
 
    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;
 
  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;
    // 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];
 
    std::string positionMode1;
    std::string positionMode2;
 
    // Shower energy correction
    TGraph2DErrors *dt;
    // Energy error
    TGraph2DErrors *dtErr;
};
 
#endif // EMC_REC_PARAMETER_A_H