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((uint32_t)det_id);
    int sPlane = APVGain::subdetectorPlane((uint32_t)det_id, topo);
    int sSide  = APVGain::subdetectorSide((uint32_t)det_id, topo);
    std::vector<APVGain::APVmon>::iterator it= histos.begin();
    
    while (it!=histos.end()) {
        //std::string tag = (*it)->getName();
        int subdetectorId    = (*it).subdetectorId;
        int subdetectorSide  = (*it).subdetectorSide;
        int subdetectorPlane = (*it).subdetectorPlane;

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

        if (match) {
            found.push_back((*it).monitor);
        }
        it++;
    }
    return found;
}



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;
    int re = (allPlanes)? 34 + 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
        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
        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
        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
        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;
}