APVGainHelpers.cc

Go to the documentation of this file.

#include "CalibTracker/SiStripChannelGain/interface/APVGainHelpers.h"
#include "DataFormats/DetId/interface/DetId.h"

#include "DataFormats/SiStripDetId/interface/StripSubdetector.h"

int APVGain::subdetectorId(uint32_t det_id) {
    return DetId(det_id).subdetId();
};


int APVGain::subdetectorId(const std::string& tag) {
    std::string d = tag.substr(0,3);
    if ( d=="TIB" ) return 3;
    if ( d=="TID" ) return 4;
    if ( d=="TOB" ) return 5;
    if ( d=="TEC" ) return 6;
    return 0;
};

int APVGain::subdetectorSide(uint32_t det_id, const TrackerTopology* topo) {
    return topo->side( det_id );
}


int APVGain::subdetectorSide(const std::string& tag) {
    std::size_t m = tag.find("minus");
    std::size_t p = tag.find("plus");
    if (m!=std::string::npos) return 1;
    if (p!=std::string::npos) return 2;
    return 0;
}


int APVGain::subdetectorPlane(uint32_t det_id, const TrackerTopology* topo) {
  if( topo ) {
    if( APVGain::subdetectorId(det_id)==StripSubdetector::TIB )      return topo->tibLayer(det_id);
    else if( APVGain::subdetectorId(det_id)==StripSubdetector::TID ) return (2*topo->tidSide(det_id)-3)*topo->tidWheel(det_id);
    else if( APVGain::subdetectorId(det_id)==StripSubdetector::TOB ) return topo->tobLayer(det_id);
    else if( APVGain::subdetectorId(det_id)==StripSubdetector::TEC ) return (2*topo->tecSide(det_id)-3)*topo->tecWheel(det_id);
  }
  return 0;
};

int APVGain::subdetectorPlane(const std::string& tag) {
   std::size_t p = (tag.find("layer")!=std::string::npos)? tag.find("layer") : tag.find("wheel");
   if (p!=std::string::npos) {
       std::size_t start = tag.find("_",p+1) + 1;
       std::size_t stop  = tag.find('_',start);
       std::string plane = tag.substr(start,stop-start);
       return atoi( plane.c_str());
   }
   return 0;
};


std::vector<MonitorElement*> APVGain::FetchMonitor(std::vector<APVGain::APVmon> histos, uint32_t det_id, 
                           const TrackerTopology* topo) {
                       
  std::vector<MonitorElement*> found = std::vector<MonitorElement*>();
  int sId    = APVGain::subdetectorId(det_id);
  int sPlane = APVGain::subdetectorPlane(det_id, topo);
  int sSide  = APVGain::subdetectorSide(det_id, topo);
  auto it = histos.begin();

  LogDebug("APVGainHelpers")<<"sId: "<<sId<<" sPlane: "<<sPlane<<" sSide: "<<sSide<<std::endl; 

  while (it!=histos.end()) {
    std::string tag = (*it).getMonitor()->getName();
    int subdetectorId    = (*it).getSubdetectorId();
    int subdetectorSide  = (*it).getSubdetectorSide();
    int subdetectorPlane = (*it).getSubdetectorPlane();
    
    bool match = (subdetectorId==0  || subdetectorId==sId) && (subdetectorPlane==0 || subdetectorPlane==sPlane) && (subdetectorSide==0  || subdetectorSide==sSide);

    if (match) {
      found.emplace_back((*it).getMonitor());
      LogDebug("APVGainHelpers")<<det_id<<" found: "<< tag << std::endl;
      (*it).printAll();
    }
    it++;
  }
  return found;
}

std::vector<unsigned int> APVGain::FetchIndices(std::map<unsigned int,APVloc> theMap, uint32_t det_id, const TrackerTopology* topo){

  std::vector<unsigned int> found_indices = std::vector<unsigned int>();

  int sId    = APVGain::subdetectorId(det_id);
  int sPlane = APVGain::subdetectorPlane(det_id, topo);
  int sSide  = APVGain::subdetectorSide(det_id, topo);
 
  for(auto &element : theMap){

    int subdetectorId    = element.second.m_subdetectorId;
    int subdetectorSide  = element.second.m_subdetectorSide;
    int subdetectorPlane = element.second.m_subdetectorPlane;

    bool match = (subdetectorId==0  || subdetectorId==sId) && (subdetectorPlane==0 || subdetectorPlane==sPlane) && (subdetectorSide==0  || subdetectorSide==sSide);

      if (match ){
      found_indices.push_back(element.first);
    }
  }
  return found_indices;
}

std::vector<std::pair<std::string,std::string>> 
APVGain::monHnames(std::vector<std::string> VH, bool allPlanes, const char* tag) {

  std::vector<std::pair<std::string,std::string>> out;

  // total number of measurement layers/wheels in the Strips Tracker
  // 4(TIB) + 6(TOB) + 3(TID+) + 3(TID-) + 9(TEC+) + 9(TEC-)
  constexpr int countOfPlanes = 34;  

  int re = (allPlanes)? countOfPlanes + VH.size() : VH.size();
  out.reserve( re );

  std::string Tag = tag;
  if (Tag.length())  Tag = "__" + Tag;

  std::string h_tag = "";
  std::string h_tit = "";

  if (allPlanes) {
    // Names of monitoring histogram for TIB layers
    constexpr int TIBlayers = 4;  //number of TIB layers.
    for(int i=1; i<=TIBlayers; i++) {
      h_tag = "TIB_layer_" + std::to_string(i) + Tag;
      h_tit = h_tag; std::replace(h_tit.begin(),h_tit.end(),'_',' ');
      out.push_back(std::pair<std::string,std::string>(h_tag,h_tit));
    }
    // Names of monitoring histogram for TOB layers
    constexpr int TOBlayers = 6;  //number of TOB layers
    for(int i=1; i<=TOBlayers; i++) {
      h_tag = "TOB_layer_" + std::to_string(i) + Tag;
      h_tit = h_tag; std::replace(h_tit.begin(),h_tit.end(),'_',' ');
      out.push_back(std::pair<std::string,std::string>(h_tag,h_tit));
    }
    // Names of monitoring histogram for TID wheels
    constexpr int TIDwheels = 3;  //number of TID wheels
    for(int i=-TIDwheels; i<=TIDwheels; i++) {
      if (i==0) continue;
      if (i<0)  h_tag = "TIDminus_wheel_" + std::to_string(i) + Tag;
      else      h_tag = "TIDplus_wheel_" + std::to_string(i) + Tag;
      h_tit = h_tag; std::replace(h_tit.begin(),h_tit.end(),'_',' ');
      out.push_back(std::pair<std::string,std::string>(h_tag,h_tit));
    }
    // Names of monitoring histogram for TEC wheels
    constexpr int TECwheels = 9;  //number of TEC wheels
    for(int i=-TECwheels; i<=TECwheels; i++) {
      if (i==0) continue;
      if (i<0) h_tag = "TECminus_wheel_" + std::to_string(i) + Tag;
      else     h_tag = "TECplus_wheel_" + std::to_string(i) + Tag;
      h_tit = h_tag; std::replace(h_tit.begin(),h_tit.end(),'_',' ');
      out.push_back(std::pair<std::string,std::string>(h_tag,h_tit));
    }
  }
  
  for(unsigned int i=0; i<VH.size();i++) {
    h_tag = VH[i] + Tag;
    h_tit = h_tag; std::replace(h_tit.begin(),h_tit.end(),'_',' ');
    out.push_back(std::pair<std::string,std::string>(h_tag,h_tit));
  }
  
  return out;
}