bak_ReadCosmicRayData-00-00-15/src/StripMapper.cxx

Go to the documentation of this file.

// **************************************************************************/
// authors: L. Lavezzi (univ. of Torino & INFN, Italy)
//
 
 
//include system lib
#include <fstream>
#include <iostream>
#include <iomanip>
#include <string>
#include <cmath>
#include <csignal>
 
#include "TSystem.h"
 
#include "ReadCosmicRayData/StripMapper.h"
 
//using namespace std;
 
StripMapper::StripMapper(TString dataset) {
  CosmicRayDataSetID = dataset;
}
 
StripMapper::~StripMapper(){}
 
bool StripMapper::FillMap() {
 
 
  if(CosmicRayDataSetID == "CR201909") {
 
    TString dir = gSystem->ExpandPathName("$READCOSMICRAYDATAROOT");
    dir += "/dat/CR201909/";
 
    // L1
    TString mapname = dir;
    ifstream input;
    std::string titleline;
 
    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
 
 
    for(int ilayer = 0; ilayer < nlayer; ilayer++) {
      for(int isheet = 0; isheet < nsheet[ilayer]; isheet++) {
    for(int iview = 0; iview < nview; iview++) {
 
      int layerid = ilayer + 1;
      int sheetid = isheet + 1;
 
      mapname = dir;
      mapname += "mapping_L";
      mapname += layerid;
      mapname += "_S";
      mapname += sheetid;
      if(iview == 0) mapname += "_x.csv";
      else mapname += "_v.csv";
      std::cout << "StripMapper::FillMap, reading " << mapname << std::endl;
 
      input.open(mapname);
      if(!input.good()) {
        std::cout << "StripMapper::FillMap, file not there!" << std::endl;
        continue;
      }
 
      std::getline(input, titleline);
 
      /**
         std::cout << mapname << std::endl;
         std::cout << "layer " << ilayer << std::endl;
         std::cout << "sheet " << isheet << std::endl;
         std::cout << "view  " << iview  << std::endl;
         std::cout << "titleline " << titleline << std::endl;
      **/
 
      for(int istrip = 0; istrip < nstrip[ilayer][iview]; istrip++) {
 
        // strip id / GEMROC / FEB / TIGER
        int stripid;
        input >> stripid
          >> mapping[istrip][ilayer][isheet][iview][0]  // gemroc
          >> mapping[istrip][ilayer][isheet][iview][1] // feb
          >> mapping[istrip][ilayer][isheet][iview][2]; // tiger
 
        //      std::cout << istrip << " <---> " << stripid << " " << mapping[istrip][ilayer][isheet][iview][0]  << " " << mapping[istrip][ilayer][isheet][iview][1] << " " << mapping[istrip][ilayer][isheet][iview][2] << std::endl;      
      }
 
      input.close();
    }
      }
    }
 
    return true;
  }
 
  return false;
 
 
}
 
 
int StripMapper::GetGEMROC(int strip, int type, int layer, int sheet) {
  return mapping[strip][layer][sheet][type][0];
}
 
int StripMapper::GetFEB(int strip, int type, int layer, int sheet) {
  return mapping[strip][layer][sheet][type][1];
}
 
int StripMapper::GetTIGER(int strip, int type, int layer, int sheet) {
  return mapping[strip][layer][sheet][type][2];
}