CgemLUTWriter.h

Go to the documentation of this file.

#ifndef CGEMLUTWRITER
#define CGEMLUTWRITER
 
#include "CgemGeomSvc.h"
#include "RealizationSvc/RealizationSvc.h"
#include "RealizationSvc/IRealizationSvc.h"
#include <vector>
#include <string>
#include "TFile.h"
#include "TTree.h"
#include "TString.h"
#include "Identifier/Identifier.h"
#include "TGraph.h"
#include "TCanvas.h"
#include <iostream>
#include <fstream>
 
#define N_FEB 80
#define N_CHIP 2
#define N_CHANNEL 64
 
class IDataProviderSvc;
 
class CgemLUTWriter : public Algorithm{
 public:
  // Constructor
  CgemLUTWriter(const std::string& name, ISvcLocator* pSvcLocator);
  ~CgemLUTWriter();
  // Executable
  StatusCode initialize();
  StatusCode execute   ();
  StatusCode finalize  ();
  
 private:
  //Variables Code
  int FEB_i,chip_i, ch_i;
  bool PrintCanvas;
  TString DrawOption;
  string OutputPdf;
  string OutputRoot;
  string InputFileData;
  string InputFileMapping;
  string InputFileTiger;
  string InputFileHvSetting;
  ifstream mapping_file, tiger_file, hv_file;
  TCanvas *canvas;
  bool DisableMeasures;
  bool DisableLowEntries;
  int  MinEntries;
  bool is_failing;
  bool Debug;
  //Variables LUT
  struct LUT_digi{
    int   layer;
    int   gemroc;
    int   sw_feb;
    int   side;
    int   sheet;
    int   HW_FEB_ID;
    float qdc_slope;
    float qdc_constant;
    float qdc_saturation;
    int   strip_x;
    int   strip_v;
    float noise_rate;
    float noise_charge;
    float thr_effective;
    float thr_width;
    float signal_rate;
    float signal_charge_mean;
    float signal_charge_saturation;
    float signal_saturation_ratio;    
    float baseline_e;
    float baseline_t;
    float vth1;
    float vth2;
    float delta_vth1_baseline;
    float delta_vth2_baseline;
    float thr_t_fC;
    float thr_e_fC;
  };
  LUT_digi LUT_digi_coll[N_FEB][N_CHIP][N_CHANNEL];
  TString ofile_name;
  TTree *otree, *ootree;
  TFile *ofile;
  int t_roc,t_sw_feb,t_feb, t_chip, t_ch,t_layer,t_side,t_sheet,t_strip_x_boss, t_strip_v_boss,t_strip_x_graal,t_strip_v_graal;
  int t_quality;
  float t_noise_Hz, t_noise_Q, t_thr, t_thr_wid, t_sig_meanQ, t_sig_maxQ, t_sig_Hz;
  float t_rate_sat_max;
  float t_qdc_slope,t_qdc_constant,t_qdc_saturation;
  float t_baseline_e,t_baseline_t,t_vth1,t_vth2,t_delta_vth1_baseline,t_delta_vth2_baseline,t_thr_t_fC,t_thr_e_fC;
  float t_time_start,t_time_stop,t_time_sigma,t_time_startFEB,t_time_stopFEB,t_time_sigmaFEB;
  int t_First_RUN,t_Last_RUN,t_Integration_time,t_Energy_mode,t_HV_L1_G3,t_HV_L1_G2,t_HV_L1_G1,t_HV_L2_G3,t_HV_L2_G2,t_HV_L2_G1,t_HV_L3_G3,t_HV_L3_G2,t_HV_L3_G1,t_HV_L1_Dr,t_HV_L2_Dr,t_HV_L3_Dr,t_HV_L1_T1,t_HV_L2_T1,t_HV_L3_T1,t_HV_L1_T2,t_HV_L2_T2,t_HV_L3_T2,t_HV_L1_In,t_HV_L2_In,t_HV_L3_In;
  //Function
  TString Get_Cut(int);
  void Read_HV_Setting();
  void Read_Mapping();
  void Read_Tiger();
  void Measure_Noise();
  void Measure_Threshold();
  void Measure_Signal();
  void Define_TTree();
  void Fill_TTree();
};
#endif