CgemLUTReader.h

Go to the documentation of this file.

#ifndef CGEMLUTREADER
#define CGEMLUTREADER
 
#include <string>
#include "TFile.h"
#include "TString.h"
#include "TTree.h"
 
/**
#define MAXLAYER 3
#define MAXSHEET 2
#define MAXVIEW 2
#define MAXSTRIP 1395 // LAYER3
**/
 
class CgemLUTReader {
public:
     // Constructor
  CgemLUTReader();  
  CgemLUTReader(TString name);
  ~CgemLUTReader();
 
  bool ReadLUT();
  std::map<int, int> * GetMap(int ilayer, int isheet, int iview);
 
  bool LoadStripInfo(int ilayer, int isheet, int iview, int istrip);
 
  /** 
      Setup of each LAYER readout: strip -> TIGER -> FEB -> ROC
      > each ROC has up to 4 FEB
      > each FEB has 2 TIGER
      To select a strip use : layer, sheet, view and strip ID (in Cgemboss style) 
 
      More info at: https://docbes3.ihep.ac.cn/~cgem/index.php/Documentation
  **/
 
  // identifier of the ROC:
  // > LAYER1 has 4  ROC [0-3]   - [0,1]   west - [2.3]   east
  // > LAYER2 has 8  ROC [4-11]  - [4,7]   west - [8,11]  east
  // > LAYER3 HAS 10 ROC [12-21] - [12,16] west - [17,21] east 
  int GetROC(int ilayer, int isheet, int iview, int istrip);
 
  // identifier of the TIGER in a ROC
  // each ROC has up to 8 TIGER [0-8]
  int GetTIGER(int ilayer, int isheet, int iview, int istrip);
 
  // identifier of the FEB
  // > LAYER1 has 16 FEB [0-15]
  // > LAYER2 has 28 FEB [16-43]
  // > LAYER3 has 36 FEB [44,79]
  int GetFEB(int ilayer, int isheet, int iview, int istrip);
 
  // identifier of the TIGER chip in a FEB [0-1]
  int GetChip(int ilayer, int isheet, int iview, int istrip);
 
  // identifier of the TIGER channel [0-63]
  int GetChannel(int ilayer, int isheet, int iview, int istrip);
 
  // superfluous function, just for checking
  int GetLayer(int ilayer, int isheet, int iview, int istrip);
 
  // superfluous function, just for checking
  int GetSheet(int ilayer, int isheet, int iview, int istrip);
 
  // readout side: (fixed: +1->west/-1->east on 15 june 2020 lia lavezzi)
  // > +1 = west
  // > -1 = east
  int GetSide(int ilayer, int isheet, int iview, int istrip);
 
  // strip identifier in GRAAL
  int GetStrip_x_GRAAL(int ilayer, int isheet, int iview, int istrip);
  int GetStrip_v_GRAAL(int ilayer, int isheet, int iview, int istrip);
 
  // strip identifier in CgemBOSS (theone used as istrip!)
  int GetStrip_x_BOSS(int ilayer, int isheet, int iview, int istrip);
  int GetStrip_v_BOSS(int ilayer, int isheet, int iview, int istrip);
 
  // calibration line of sample-and-hold (E branch) 
  float GetCalib_QCD_slope(int ilayer, int isheet, int iview, int istrip);
  float GetCalib_QCD_const(int ilayer, int isheet, int iview, int istrip);
 
  // saturation value (QDC)
  float GetCalib_QCD_saturation(int ilayer, int isheet, int iview, int istrip);
 
  // threshold applied to the T branch (in mV)
  float GetV_thr_T_mV(int ilayer, int isheet, int iview, int istrip);
 
  // threshold applied to the E branch (in mV)
  float GetV_thr_E_mV(int ilayer, int isheet, int iview, int istrip);
 
  // threshold in the T branch (translated to fC)
  float Get_thr_T_fC(int ilayer, int isheet, int iview, int istrip);
 
  // threshold in the E branch (translated to fC)
  float Get_thr_E_fC(int ilayer, int isheet, int iview, int istrip);
 
  // baseline in the T branch (in mV) 
  float GetBaseline_T_mV(int ilayer, int isheet, int iview, int istrip);
 
  // baseline in the E branch (in mV)  
  float GetBaseline_E_mV(int ilayer, int isheet, int iview, int istrip);
 
  // measured effective threshold (in fC)
  float Get_thr_fC(int ilayer, int isheet, int iview, int istrip);
 
  // measured noise rate (Hz)
  float GetNoise_Rate_Hz(int ilayer, int isheet, int iview, int istrip);
 
  // measured noise level (fC)
  float GetNoise_MeanQ_fC(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal mean value (fC)
  float GetSignal_MeanQ_fC(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal maximum value (fC)
  float GetSignal_MaxQ_fC(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal rate (Hz)
  float GetSignal_Rate_Hz(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal rising edge (channel by channel) --> start time (ns)
  float GetSignal_StartTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal falling edge (channel by channel) --> stop time (ns)
  float GetSignal_StopTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal start time width (ns) (channel by channel)
  float GetSignal_SigmaTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal rising edge (FEB by FEB) --> start time (ns)
  float GetSignal_FEBStartTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal falling edge (FEB by FEB) --> stop time (ns)
  float GetSignal_FEBStopTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // measured signal start time width (ns) (FEB by FEB)
  float GetSignal_FEBSigmaTime_ns(int ilayer, int isheet, int iview, int istrip);
 
  // strip/channel quality
  // 0 = good
  // the following are with "minus" sign, but could be with "plus" sign: consider the abs value
  // -1 = channels not connected
  // -2 = low rate in signal time window (< 10 Hz)
  // -3 = max charge < 20 fC (only x view)
  // -4 = mean charge < 3 fC (only x view)
  // -5 = torn tail on L2 top
  // -6 = noise charge > 20 fC
  int GetQuality(int ilayer, int isheet, int iview, int istrip);
 
  bool PrintMap(int ilayer, int isheet, int iview);
     
  //private:
  TString CosmicRayLUTName;
     
  std::map<int, int> map_L1_S1_stripx_to_entry;
  std::map<int, int> map_L2_S1_stripx_to_entry;
  std::map<int, int> map_L2_S2_stripx_to_entry;
  std::map<int, int> map_L3_S1_stripx_to_entry;
  std::map<int, int> map_L3_S2_stripx_to_entry;
  std::map<int, int> map_L1_S1_stripv_to_entry;
  std::map<int, int> map_L2_S1_stripv_to_entry;
  std::map<int, int> map_L2_S2_stripv_to_entry;
  std::map<int, int> map_L3_S1_stripv_to_entry;
  std::map<int, int> map_L3_S2_stripv_to_entry;
 
 
  //public:
  TFile *f;
  TTree *tree;
  int ROC;
  int TIGER;
  // SW_FEB
  int FEB;
  int chip;
  int channel;
  int layer;
  int side;
  int strip_x_graal;
  int strip_v_graal;
  int strip_x_boss;
  int strip_v_boss;
  int sheet;
  float calib_QDC_slope;
  float calib_QDC_const;
  float calib_QDC_saturation;
  float v_thr_T_mV;
  float v_thr_E_mV;
  float thr_T_fC;
  float thr_E_fC;
  float baseline_T_mV;
  float baseline_E_mV;
  float thr_fC;
  float noise_rate_Hz;
  float noise_meanQ_fC;
  float signal_meanQ_fC;
  float signal_maxQ_fC;
  float signal_rate_Hz;
  float signal_startTime_ns;
  float signal_stopTime_ns;
  float signal_sigmaTime_ns;
  float signal_FEBstartTime_ns;
  float signal_FEBstopTime_ns;
  float signal_FEBsigmaTime_ns;
  int quality;
 
  // general
  TTree *general;
  int First_RUN;
  int Last_RUN;
  int Integration_time;
  int Energy_mode;
  int HV_L1_G3;
  int HV_L2_G3;
  int HV_L3_G3;
  int HV_L1_G2;
  int HV_L2_G2;
  int HV_L3_G2;
  int HV_L1_G1;
  int HV_L2_G1;
  int HV_L3_G1;
  int HV_L1_In;
  int HV_L2_In;
  int HV_L3_In;
  int HV_L1_T2;
  int HV_L2_T2;
  int HV_L3_T2;
  int HV_L1_T1;
  int HV_L2_T1;
  int HV_L3_T1;
  int HV_L1_Dr;
  int HV_L2_Dr;
  int HV_L3_Dr;
 
 
};
#endif