CgemLUTReader.cxx

Go to the documentation of this file.

// **************************************************************************/
// authors: L. Lavezzi (univ. of Torino & INFN, Italy) - R. Farinelli (INFN Ferrara, Italy)
//
 
 
//include system lib
#include <iostream>
#include <iomanip>
#include <string>
#include <cmath>
#include <csignal>
 
#include "TSystem.h"
 
#include "ReadCosmicRayData/CgemLUTReader.h"
 
using namespace std;
 
CgemLUTReader::CgemLUTReader() : tree(NULL) {
  CosmicRayLUTName = "";
}
 
CgemLUTReader::CgemLUTReader(TString lutname) : tree(NULL) {
  CosmicRayLUTName = lutname;
}
 
CgemLUTReader::~CgemLUTReader(){}
 
bool CgemLUTReader::ReadLUT() {
 
    const int nlayer = 3;                 // CHECK hardcoded 
    const int nsheet[nlayer] = {1, 2, 2}; // CHECK hardcoded
    const int nview = 2;                  // CHECK hardcoded 
    int nstrip[nlayer][nview] = {{856, 1173}, {630, 1077}, {832, 1395}}; // CHECK hardcoded
 
    // ----------------------------------------
    f = new TFile(CosmicRayLUTName, "READ");
    tree = (TTree*) f->Get("tree");
    tree->SetBranchAddress("ROC",&ROC);
    tree->SetBranchAddress("SW_FEB",&TIGER);
    tree->SetBranchAddress("FEB",&FEB);
    tree->SetBranchAddress("chip",&chip);
    tree->SetBranchAddress("channel",&channel);
    tree->SetBranchAddress("layer",&layer);
    tree->SetBranchAddress("side",&side);
    tree->SetBranchAddress("strip_x_graal",&strip_x_graal);
    tree->SetBranchAddress("strip_v_graal",&strip_v_graal);
    tree->SetBranchAddress("strip_x_boss",&strip_x_boss);
    tree->SetBranchAddress("strip_v_boss",&strip_v_boss);
    tree->SetBranchAddress("sheet",&sheet);
    tree->SetBranchAddress("calib_QDC_slope",&calib_QDC_slope);
    tree->SetBranchAddress("calib_QDC_const",&calib_QDC_const);
    tree->SetBranchAddress("calib_QDC_saturation",&calib_QDC_saturation);
    tree->SetBranchAddress("v_thr_T_mV",&v_thr_T_mV);
    tree->SetBranchAddress("v_thr_E_mV",&v_thr_E_mV);
    tree->SetBranchAddress("thr_T_fC",&thr_T_fC);
    tree->SetBranchAddress("thr_E_fC",&thr_E_fC);
    tree->SetBranchAddress("baseline_T_mV",&baseline_T_mV);
    tree->SetBranchAddress("baseline_E_mV",&baseline_E_mV);
    tree->SetBranchAddress("thr_fC",&thr_fC);
    tree->SetBranchAddress("noise_rate_Hz",&noise_rate_Hz);
    tree->SetBranchAddress("noise_meanQ_fC",&noise_meanQ_fC);
    tree->SetBranchAddress("signal_meanQ_fC",&signal_meanQ_fC);
    tree->SetBranchAddress("signal_maxQ_fC",&signal_maxQ_fC);
    tree->SetBranchAddress("signal_rate_Hz",&signal_rate_Hz);
    tree->SetBranchAddress("signal_startTime_ns",&signal_startTime_ns);
    tree->SetBranchAddress("signal_stopTime_ns",&signal_stopTime_ns);
    tree->SetBranchAddress("signal_sigmaTime_ns",&signal_sigmaTime_ns);
    tree->SetBranchAddress("signal_FEBstartTime_ns",&signal_FEBstartTime_ns);
    tree->SetBranchAddress("signal_FEBstopTime_ns",&signal_FEBstopTime_ns);
    tree->SetBranchAddress("signal_FEBsigmaTime_ns",&signal_FEBsigmaTime_ns);
    tree->SetBranchAddress("quality",&quality);
    // ----------------------------------------
    general = (TTree*) f->Get("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;
    general->SetBranchAddress("First_RUN", &First_RUN);
    general->SetBranchAddress("Last_RUN", &Last_RUN);
    general->SetBranchAddress("Integration_time", &Integration_time);
    general->SetBranchAddress("Energy_mode", &Energy_mode);
    general->SetBranchAddress("HV_L1_G3", &HV_L1_G3);
    general->SetBranchAddress("HV_L1_G2", &HV_L1_G2);
    general->SetBranchAddress("HV_L1_G1", &HV_L1_G1);
    general->SetBranchAddress("HV_L2_G3", &HV_L2_G3);
    general->SetBranchAddress("HV_L2_G2", &HV_L2_G2);
    general->SetBranchAddress("HV_L2_G1", &HV_L2_G1);
    general->SetBranchAddress("HV_L3_G3", &HV_L3_G3);
    general->SetBranchAddress("HV_L3_G2", &HV_L3_G2);
    general->SetBranchAddress("HV_L3_G1", &HV_L3_G1);
    general->SetBranchAddress("HV_L1_Dr", &HV_L1_Dr);
    general->SetBranchAddress("HV_L2_Dr", &HV_L2_Dr);
    general->SetBranchAddress("HV_L3_Dr", &HV_L3_Dr);
    general->SetBranchAddress("HV_L1_T1", &HV_L1_T1);
    general->SetBranchAddress("HV_L2_T1", &HV_L2_T1);
    general->SetBranchAddress("HV_L3_T1", &HV_L3_T1);
    general->SetBranchAddress("HV_L1_T2", &HV_L1_T2);
    general->SetBranchAddress("HV_L2_T2", &HV_L2_T2);
    general->SetBranchAddress("HV_L3_T2", &HV_L3_T2);
    general->SetBranchAddress("HV_L1_In", &HV_L1_In);
    general->SetBranchAddress("HV_L2_In", &HV_L2_In);
    general->SetBranchAddress("HV_L3_In", &HV_L3_In);
    general->GetEntry(0);
    // ----------------------------------------
 
    int nchannel_off=0;
      
    for(int ientry=0; ientry<tree->GetEntries(); ientry++) { 
      tree->GetEntry(ientry);
      
      std::pair<int, int> pairx(strip_x_boss, ientry);
      std::pair<int, int> pairv(strip_v_boss, ientry);
 
      if(strip_x_boss != -1) {
    if(layer == 0 && sheet == 0)      map_L1_S1_stripx_to_entry.insert(pairx);
        else if(layer == 1 && sheet == 0) map_L2_S1_stripx_to_entry.insert(pairx); // sheet 1 and 2 are swapped
    else if(layer == 1 && sheet == 1) map_L2_S2_stripx_to_entry.insert(pairx); // 1 now is 0 and viceversa
    else if(layer == 2 && sheet == 0) map_L3_S1_stripx_to_entry.insert(pairx);
    else if(layer == 2 && sheet == 1) map_L3_S2_stripx_to_entry.insert(pairx);
      }
      else if(strip_v_boss != -1) {
        if(layer == 0 && sheet == 0)      map_L1_S1_stripv_to_entry.insert(pairv);
    else if(layer == 1 && sheet == 0) map_L2_S1_stripv_to_entry.insert(pairv); // sheet 1 and 2 are swapped
    else if(layer == 1 && sheet == 1) map_L2_S2_stripv_to_entry.insert(pairv);  // 1 now is 0 and viceversa
    else if(layer == 2 && sheet == 0) map_L3_S1_stripv_to_entry.insert(pairv);
    else if(layer == 2 && sheet == 1) map_L3_S2_stripv_to_entry.insert(pairv);
      }
      else nchannel_off++;
 
      //      std::cout << "entry " << ientry << " channel " << channel << " layer " << layer << " sheet " << sheet << " ROC " << ROC << std::endl; 
    }
 
    std::cout <<"The number of channels not connected to strips is: " << nchannel_off << "(i.e. strip X&V = -1)" <<std::endl;
 
}
 
 
 
 
std::map<int, int> * CgemLUTReader::GetMap(int ilayer, int isheet, int iview) {
    
  if(iview==0) {
    if(ilayer == 0 && isheet == 0)       return &map_L1_S1_stripx_to_entry;
    else if(ilayer == 1 && isheet == 0)  return &map_L2_S1_stripx_to_entry; // unswapped - farinelli
    else if(ilayer == 1 && isheet == 1)  return &map_L2_S2_stripx_to_entry; // unswapped - farinelli
    else if(ilayer == 2 && isheet == 0)  return &map_L3_S1_stripx_to_entry;
    else if(ilayer == 2 && isheet == 1)  return &map_L3_S2_stripx_to_entry;
  }
  else if(iview==1) {
    if(ilayer == 0 && isheet == 0)       return &map_L1_S1_stripv_to_entry;
    else if(ilayer == 1 && isheet == 0)  return &map_L2_S1_stripv_to_entry; // unswapped - farinelli
    else if(ilayer == 1 && isheet == 1)  return &map_L2_S2_stripv_to_entry; // unswapped - farinelli
    else if(ilayer == 2 && isheet == 0)  return &map_L3_S1_stripv_to_entry;
    else if(ilayer == 2 && isheet == 1)  return &map_L3_S2_stripv_to_entry;
  }
  else return NULL;
}
 
bool CgemLUTReader::LoadStripInfo(int ilayer, int isheet, int iview, int istrip) {
  std::map<int, int> *mappa = GetMap(ilayer, isheet, iview);
  if(mappa == NULL) return false;
  std::map<int, int>::iterator it = mappa->find(istrip);
  if (it != mappa->end()) {
    int ientry = it->second;
    tree->GetEntry(ientry);
    return true;
  }
  return false;;
}
 
// --------------------------------
int CgemLUTReader::GetROC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return ROC;
}
 
int CgemLUTReader::GetTIGER(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return TIGER;
}
 
int CgemLUTReader::GetFEB(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return FEB;
}
 
int CgemLUTReader::GetChip(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return chip;
}
 
int CgemLUTReader::GetChannel(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return channel;
}
 
int CgemLUTReader::GetLayer(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return layer;
}
 
int CgemLUTReader::GetSheet(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return sheet;
}
 
int CgemLUTReader::GetSide(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return side;
}
 
int CgemLUTReader::GetStrip_x_GRAAL(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return strip_x_graal;
}
 
int CgemLUTReader::GetStrip_v_GRAAL(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return strip_v_graal;
}
 
int CgemLUTReader::GetStrip_x_BOSS(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return strip_x_boss;
}
 
int CgemLUTReader::GetStrip_v_BOSS(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return strip_v_boss;
}
 
float CgemLUTReader::GetCalib_QCD_slope(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return calib_QDC_slope;
}
 
float CgemLUTReader::GetCalib_QCD_const(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return calib_QDC_const;
}
 
float CgemLUTReader::GetCalib_QCD_saturation(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return calib_QDC_saturation;
}
 
float CgemLUTReader::GetV_thr_T_mV(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return v_thr_T_mV;
}
 
float CgemLUTReader::GetV_thr_E_mV(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return v_thr_E_mV;
}
 
float CgemLUTReader::Get_thr_T_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return thr_T_fC;
}
 
float CgemLUTReader::Get_thr_E_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return thr_E_fC;
}
 
float CgemLUTReader::GetBaseline_T_mV(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return baseline_T_mV;
}
 
float CgemLUTReader::GetBaseline_E_mV(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return baseline_E_mV;
}
 
float CgemLUTReader::Get_thr_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return thr_fC;
}
 
float CgemLUTReader::GetNoise_Rate_Hz(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return noise_rate_Hz;
}
 
float CgemLUTReader::GetNoise_MeanQ_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return noise_meanQ_fC;
}
 
float CgemLUTReader::GetSignal_MeanQ_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_meanQ_fC;
}
 
float CgemLUTReader::GetSignal_MaxQ_fC(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_maxQ_fC;
}
 
float CgemLUTReader::GetSignal_Rate_Hz(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_rate_Hz;
}
 
float CgemLUTReader::GetSignal_StartTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_startTime_ns;
}
 
float CgemLUTReader::GetSignal_StopTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_stopTime_ns;
}
 
float CgemLUTReader::GetSignal_SigmaTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_sigmaTime_ns;
}
 
float CgemLUTReader::GetSignal_FEBStartTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_FEBstartTime_ns;
}
 
float CgemLUTReader::GetSignal_FEBStopTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_FEBstopTime_ns;
}
 
float CgemLUTReader::GetSignal_FEBSigmaTime_ns(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return 0.;
  return signal_FEBsigmaTime_ns;
}
 
int CgemLUTReader::GetQuality(int ilayer, int isheet, int iview, int istrip) {
  bool doesexist = LoadStripInfo(ilayer, isheet, iview, istrip);
  if(doesexist == false) return -1;
  return quality;
}
// ------------------------------------------------
bool CgemLUTReader::PrintMap(int ilayer, int isheet, int iview) {
 
  std::map<int, int> *mappa = GetMap(ilayer, isheet, iview);
  if(mappa == NULL) return false;
 
  int nstrip = mappa->size();
  std::cout << "MAPPING STRIPS of view " << iview << " on sheet " << isheet << " of layer " << ilayer << std::endl;
  std::cout << "N. OF STRIP " << nstrip << std::endl;
  
  std::map<int, int>::iterator it;
  for(int istrip=0; istrip < nstrip; istrip++) {
    it = mappa->find(istrip);
    int ientry = it->second;
    tree->GetEntry(ientry);
    std::cout << "strip " << istrip << " channel " << channel << " ROC " << ROC << " TIGER " << TIGER << " FEB " << FEB << std::endl;
  }
 
  return true;
}