SiPixelActionExecutor.cc

Go to the documentation of this file.

//#define printing false
//#define occupancyprinting false

#include "DQM/SiPixelMonitorClient/interface/SiPixelActionExecutor.h"
#include "DQM/SiPixelMonitorClient/interface/SiPixelUtility.h"
#include "DQM/SiPixelMonitorClient/interface/SiPixelInformationExtractor.h"
#include "DQM/SiPixelMonitorClient/interface/ANSIColors.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/FileInPath.h"
#include "DQMServices/Core/interface/DQMStore.h"
#include "DQMServices/Core/interface/MonitorElement.h"
#include "DQM/SiPixelCommon/interface/SiPixelFolderOrganizer.h"
#include "DataFormats/DetId/interface/DetId.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelName.h"
#include "DataFormats/SiPixelDetId/interface/PixelBarrelNameUpgrade.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapName.h"
#include "DataFormats/SiPixelDetId/interface/PixelEndcapNameUpgrade.h"
#include "DataFormats/SiPixelDetId/interface/PixelSubdetector.h"

#include <math.h>

#include <iostream>
using namespace std;
//=============================================================================================================
//
// -- Constructor
// 
SiPixelActionExecutor::SiPixelActionExecutor(bool offlineXMLfile, 
                                             bool Tier0Flag) : 
  offlineXMLfile_(offlineXMLfile), 
  Tier0Flag_(Tier0Flag) {
  edm::LogInfo("SiPixelActionExecutor") << 
    " Creating SiPixelActionExecutor " << "\n" ;
  configParser_ = 0;
  configWriter_ = 0;
  ndet_ = 0;
  //collationDone = false;
}
//=============================================================================================================
//
// --  Destructor
// 
SiPixelActionExecutor::~SiPixelActionExecutor() {
  edm::LogInfo("SiPixelActionExecutor") << 
    " Deleting SiPixelActionExecutor " << "\n" ;
  if (configParser_) delete configParser_;
  if (configWriter_) delete configWriter_;  
}
//=============================================================================================================
//
// -- Read Configuration File
//
void SiPixelActionExecutor::readConfiguration() {
  string localPath;
  if(offlineXMLfile_) localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
  else localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
  if (configParser_ == 0) {
    configParser_ = new SiPixelConfigParser();
    configParser_->getDocument(edm::FileInPath(localPath).fullPath());
  }
}
//=============================================================================================================
//
// -- Read Configuration File
//
bool SiPixelActionExecutor::readConfiguration(int& tkmap_freq, 
                                              int& sum_barrel_freq, 
                                              int& sum_endcap_freq, 
                          int& sum_grandbarrel_freq, 
                          int& sum_grandendcap_freq, 
                          int& message_limit_,
                          int& source_type_,
                          int& calib_type_) {
  //printing cout<<"Entering SiPixelActionExecutor::readConfiguration..."<<endl;
  string localPath;
  if(offlineXMLfile_) localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
  else localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
  if (configParser_ == 0) {
    configParser_ = new SiPixelConfigParser();
    configParser_->getDocument(edm::FileInPath(localPath).fullPath());
  }
    
  if (!configParser_->getFrequencyForTrackerMap(tkmap_freq)){
    cout << "SiPixelActionExecutor::readConfiguration: Failed to read TrackerMap configuration parameters!! ";
    return false;
  }
  if (!configParser_->getFrequencyForBarrelSummary(sum_barrel_freq)){
    edm::LogInfo("SiPixelActionExecutor") << "Failed to read Barrel Summary configuration parameters!! " << "\n" ;
    return false;
  }
  if (!configParser_->getFrequencyForEndcapSummary(sum_endcap_freq)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read Endcap Summary configuration parameters!! " << "\n" ;
    return false;
  }
  if (!configParser_->getFrequencyForGrandBarrelSummary(sum_grandbarrel_freq)){
    edm::LogInfo("SiPixelActionExecutor") << "Failed to read Grand Barrel Summary configuration parameters!! " << "\n" ;
    return false;
  }
  if (!configParser_->getFrequencyForGrandEndcapSummary(sum_grandendcap_freq)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read Grand Endcap Summary configuration parameters!! " << "\n" ;
    return false;
  }
  if (!configParser_->getMessageLimitForQTests(message_limit_)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read QTest Message Limit" << "\n" ;
    return false;
  }
  if (!configParser_->getSourceType(source_type_)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read Source Type" << "\n" ;
    return false;
  }
  if (!configParser_->getCalibType(calib_type_)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read Calib Type" << "\n" ;
    return false;
  }
  //printing cout<<"...leaving SiPixelActionExecutor::readConfiguration..."<<endl;
  return true;
}
//=============================================================================================================
bool SiPixelActionExecutor::readConfiguration(int& tkmap_freq, int& summary_freq) {
  //printing cout<<"Entering SiPixelActionExecutor::readConfiguration..."<<endl;
  string localPath;
  if(offlineXMLfile_) localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
  else localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
  if (configParser_ == 0) {
    configParser_ = new SiPixelConfigParser();
    configParser_->getDocument(edm::FileInPath(localPath).fullPath());
  }
    
  if (!configParser_->getFrequencyForTrackerMap(tkmap_freq)){
    cout << "SiPixelActionExecutor::readConfiguration: Failed to read TrackerMap configuration parameters!! ";
    return false;
  }
  if (!configParser_->getFrequencyForBarrelSummary(summary_freq)){
    edm::LogInfo("SiPixelActionExecutor") << "Failed to read Summary configuration parameters!! " << "\n" ;
    return false;
  }
  //printing cout<<"...leaving SiPixelActionExecutor::readConfiguration..."<<endl;
  return true;
}

//=============================================================================================================
void SiPixelActionExecutor::createSummary(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool isUpgrade) {
  //To be majorly overhauled and split into two, I guess.

  //cout<<"entering SiPixelActionExecutor::createSummary..."<<endl;
  string barrel_structure_name;
  vector<string> barrel_me_names;
  string localPath;
  if(offlineXMLfile_) localPath = string("DQM/SiPixelMonitorClient/test/sipixel_tier0_config.xml");
  else localPath = string("DQM/SiPixelMonitorClient/test/sipixel_monitorelement_config.xml");
//  cout<<"*********************ATTENTION! LOCALPATH= "<<localPath<<endl;
  if (configParser_ == 0) {
    configParser_ = new SiPixelConfigParser();
    configParser_->getDocument(edm::FileInPath(localPath).fullPath());
  }
  if (!configParser_->getMENamesForBarrelSummary(barrel_structure_name, barrel_me_names)){
    cout << "SiPixelActionExecutor::createSummary: Failed to read Barrel Summary configuration parameters!! ";
    return;
  }
  configParser_->getSourceType(source_type_); 

  iBooker.setCurrentFolder("Pixel/");
  iGetter.setCurrentFolder("Pixel/");
  fillSummary(iBooker, iGetter, barrel_structure_name, barrel_me_names, true, isUpgrade); // Barrel
  iBooker.setCurrentFolder("Pixel/");
  iGetter.setCurrentFolder("Pixel/");
  string endcap_structure_name;
  vector<string> endcap_me_names;
  if (!configParser_->getMENamesForEndcapSummary(endcap_structure_name, endcap_me_names)){
    edm::LogInfo("SiPixelActionExecutor")  << "Failed to read Endcap Summary configuration parameters!! " << "\n" ;
    return;
  }

//  cout << "--- Processing endcap" << endl;

  iBooker.setCurrentFolder("Pixel/");
  iGetter.setCurrentFolder("Pixel/");

  fillSummary(iBooker,iGetter, endcap_structure_name, endcap_me_names, false, isUpgrade); // Endcap
  iBooker.setCurrentFolder("Pixel/");
  iGetter.setCurrentFolder("Pixel/");

  if(source_type_==0||source_type_==5 || source_type_ == 20){//do this only if RawData source is present
    string federror_structure_name;
    vector<string> federror_me_names;
    if (!configParser_->getMENamesForFEDErrorSummary(federror_structure_name, federror_me_names)){
      cout << "SiPixelActionExecutor::createSummary: Failed to read FED Error Summary configuration parameters!! ";
      return;
    }
    iBooker.setCurrentFolder("Pixel/");
    iGetter.setCurrentFolder("Pixel/");

    fillFEDErrorSummary(iBooker, iGetter, federror_structure_name, federror_me_names);
    iBooker.setCurrentFolder("Pixel/");
    iGetter.setCurrentFolder("Pixel/");
  }
  if (configWriter_) delete configWriter_;
  configWriter_ = 0;
//  cout<<"leaving SiPixelActionExecutor::createSummary..."<<endl;
}

//=============================================================================================================
void SiPixelActionExecutor::bookDeviations(DQMStore::IBooker & iBooker, bool isUpgrade) {
  int nBPixModules;
  if (isUpgrade) {nBPixModules=1184;} else {nBPixModules=768;} 
  
  iBooker.cd();
  iBooker.setCurrentFolder("Pixel/Barrel");
  DEV_adc_Barrel = iBooker.book1D("DEV_adc_Barrel","Deviation from reference;Module;<adc_ref>-<adc>",nBPixModules,0.,nBPixModules);
  DEV_ndigis_Barrel = iBooker.book1D("DEV_ndigis_Barrel","Deviation from reference;Module;<ndigis_ref>-<ndigis>",nBPixModules,0.,nBPixModules);
  DEV_charge_Barrel = iBooker.book1D("DEV_charge_Barrel","Deviation from reference;Module;<charge_ref>-<charge>",nBPixModules,0.,nBPixModules);
  DEV_nclusters_Barrel = iBooker.book1D("DEV_nclusters_Barrel","Deviation from reference;Module;<nclusters_ref>-<nclusters>",nBPixModules,0.,nBPixModules);
  DEV_size_Barrel = iBooker.book1D("DEV_size_Barrel","Deviation from reference;Module;<size_ref>-<size>",nBPixModules,0.,nBPixModules);
  iBooker.cd();
  iBooker.setCurrentFolder("Pixel/Endcap");
  DEV_adc_Endcap = iBooker.book1D("DEV_adc_Endcap","Deviation from reference;Module;<adc_ref>-<adc>",672,0.,672.);
  DEV_ndigis_Endcap = iBooker.book1D("DEV_ndigis_Endcap","Deviation from reference;Module;<ndigis_ref>-<ndigis>",672,0.,672.);
  DEV_charge_Endcap = iBooker.book1D("DEV_charge_Endcap","Deviation from reference;Module;<charge_ref>-<charge>",672,0.,672.);
  DEV_nclusters_Endcap = iBooker.book1D("DEV_nclusters_Endcap","Deviation from reference;Module;<nclusters_ref>-<nclusters>",672,0.,672.);
  DEV_size_Endcap = iBooker.book1D("DEV_size_Endcap","Deviation from reference;Module;<size_ref>-<size>",672,0.,672.);  
  iBooker.cd();
}


void SiPixelActionExecutor::fillDeviations(DQMStore::IGetter& iGetter) {
  int n = 768;
  MonitorElement* me1; MonitorElement* me2; 
  MonitorElement* me3; MonitorElement* me4; 
  MonitorElement* me5; 
  TH1* ref1; TH1* ref2; 
  TH1* ref3; TH1* ref4; 
  TH1* ref5; 
  MonitorElement* dev1; MonitorElement* dev2; 
  MonitorElement* dev3; MonitorElement* dev4; 
  MonitorElement* dev5; 
  me1 = iGetter.get("Pixel/Barrel/SUMDIG_adc_Barrel");
  ref1 = me1->getRefTH1();
  dev1 = iGetter.get("Pixel/Barrel/DEV_adc_Barrel");
  me2 = iGetter.get("Pixel/Barrel/SUMDIG_ndigis_Barrel");
  ref2 = me2->getRefTH1();
  dev2 = iGetter.get("Pixel/Barrel/DEV_ndigis_Barrel");
  me3 = iGetter.get("Pixel/Barrel/SUMCLU_charge_Barrel");
  ref3 = me3->getRefTH1();
  dev3 = iGetter.get("Pixel/Barrel/DEV_charge_Barrel");
  me4 = iGetter.get("Pixel/Barrel/SUMCLU_nclusters_Barrel");
  ref4 = me4->getRefTH1();
  dev4 = iGetter.get("Pixel/Barrel/DEV_nclusters_Barrel");
  me5 = iGetter.get("Pixel/Barrel/SUMCLU_size_Barrel");
  ref5 = me5->getRefTH1();
  dev5 = iGetter.get("Pixel/Barrel/DEV_size_Barrel");
  for(int i=1; i!=n+1; i++){
    float ref_value; float new_value;
    // Barrel adc: 
    if(me1)if(ref1)if(dev1){
      new_value = me1->getBinContent(i);
      ref_value = ref1->GetBinContent(i); 
      dev1->setBinContent(i,ref_value-new_value);
    }
    //Barrel ndigis:
    if(me2)if(ref2)if(dev2){
      new_value = me2->getBinContent(i);
      ref_value = ref2->GetBinContent(i); 
      dev2->setBinContent(i,ref_value-new_value);
    }
    // Barrel cluster charge:
    if(me3)if(ref3)if(dev3){
      new_value = me3->getBinContent(i);
      ref_value = ref3->GetBinContent(i); 
      dev3->setBinContent(i,ref_value-new_value);
    }
    // Barrel nclusters:
    if(me4)if(ref4)if(dev4){
      new_value = me4->getBinContent(i);
      ref_value = ref4->GetBinContent(i); 
      dev4->setBinContent(i,ref_value-new_value);
    }
    // Barrel cluster size:
    if(me5)if(ref5)if(dev5){
      new_value = me5->getBinContent(i);
      ref_value = ref5->GetBinContent(i); 
      dev5->setBinContent(i,ref_value-new_value);
    }
  }

  int nn = 672;
  MonitorElement* me11; MonitorElement* me12; 
  MonitorElement* me13; MonitorElement* me14; 
  MonitorElement* me15; 
  TH1* ref11; TH1* ref12; 
  TH1* ref13; TH1* ref14; 
  TH1* ref15; 
  MonitorElement* dev11; MonitorElement* dev12; 
  MonitorElement* dev13; MonitorElement* dev14; 
  MonitorElement* dev15; 
  me11 = iGetter.get("Pixel/Endcap/SUMDIG_adc_Endcap");
  ref11 = me11->getRefTH1();
  dev11 = iGetter.get("Pixel/Endcap/DEV_adc_Endcap");
  me12 = iGetter.get("Pixel/Endcap/SUMDIG_ndigis_Endcap");
  ref12 = me12->getRefTH1();
  dev12 = iGetter.get("Pixel/Endcap/DEV_ndigis_Endcap");
  me13 = iGetter.get("Pixel/Endcap/SUMCLU_charge_Endcap");
  ref13 = me13->getRefTH1();
  dev13 = iGetter.get("Pixel/Endcap/DEV_charge_Endcap");
  me14 = iGetter.get("Pixel/Endcap/SUMCLU_nclusters_Endcap");
  ref14 = me14->getRefTH1();
  dev14 = iGetter.get("Pixel/Endcap/DEV_nclusters_Endcap");
  me15 = iGetter.get("Pixel/Endcap/SUMCLU_size_Endcap");
  ref15 = me15->getRefTH1();
  dev15 = iGetter.get("Pixel/Endcap/DEV_size_Endcap");
  for(int i=1; i!=nn+1; i++){
    float ref_value; float new_value;
    // Endcap adc: 
    if(me11)if(ref11)if(dev11){
      new_value = me11->getBinContent(i);
      ref_value = ref11->GetBinContent(i); 
      dev11->setBinContent(i,ref_value-new_value);
    }
    //Endcap ndigis:
    if(me12)if(ref12)if(dev12){
      new_value = me12->getBinContent(i);
      ref_value = ref12->GetBinContent(i); 
      dev12->setBinContent(i,ref_value-new_value);
    }
    // Endcap cluster charge:
    if(me13)if(ref13)if(dev13){
      new_value = me13->getBinContent(i);
      ref_value = ref13->GetBinContent(i); 
      dev13->setBinContent(i,ref_value-new_value);
    }
    // Endcap nclusters:
    if(me14)if(ref14)if(dev14){
      new_value = me14->getBinContent(i);
      ref_value = ref14->GetBinContent(i); 
      dev14->setBinContent(i,ref_value-new_value);
    }
    // Endcap cluster size:
    if(me15)if(ref15)if(dev15){
      new_value = me15->getBinContent(i);
      ref_value = ref15->GetBinContent(i); 
      dev15->setBinContent(i,ref_value-new_value);
    }
  }
}

//=============================================================================================================

void SiPixelActionExecutor::GetBladeSubdirs(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, vector<string>& blade_subdirs) {

    blade_subdirs.clear();
    vector<string> panels = iGetter.getSubdirs();
    vector<string> modules;
    for (vector<string>::const_iterator it = panels.begin(); it != panels.end(); it++) {
        iGetter.cd(*it);
        iBooker.cd(*it);
        modules = iGetter.getSubdirs();
        for (vector<string>::const_iterator m_it = modules.begin(); m_it != modules.end(); m_it++) {
            blade_subdirs.push_back(*m_it);
        }
    }
}


//=============================================================================================================

void SiPixelActionExecutor::fillSummary(DQMStore::IBooker& iBooker, DQMStore::IGetter & iGetter, string dir_name, vector<string>& me_names, bool isbarrel, bool isUpgrade)
{


  //cout<<"entering SiPixelActionExecutor::fillSummary..."<<endl;
  string currDir = iBooker.pwd();
  string prefix;
  if(source_type_==0) prefix="SUMRAW";
  else if (source_type_==1) prefix="SUMDIG";
  else if (source_type_==2) prefix="SUMCLU";
  else if (source_type_==3) prefix="SUMTRK";
  else if (source_type_==4) prefix="SUMHIT";
  else if (source_type_>=7 && source_type_<20) prefix="SUMCAL";
  else if (source_type_==20) prefix="SUMOFF";
  if (currDir.find(dir_name) != string::npos)  {
    vector<MonitorElement*> sum_mes;
    for (vector<string>::const_iterator iv = me_names.begin();
     iv != me_names.end(); iv++) {
      if(source_type_==5||source_type_==6){
    if((*iv)=="errorType"||(*iv)=="NErrors"||(*iv)=="fullType"||(*iv)=="chanNmbr"||
       (*iv)=="TBMType"||(*iv)=="EvtNbr"||(*iv)=="evtSize"||(*iv)=="linkId"||
       (*iv)=="ROCId"||(*iv)=="DCOLId"||(*iv)=="PXId"||(*iv)=="ROCNmbr"||
       (*iv)=="TBMMessage"||(*iv)=="Type36Hitmap") 
      prefix="SUMRAW";
    else if((*iv)=="ndigis"||(*iv)=="adc")
      prefix="SUMDIG";
    else if((*iv)=="nclusters"||(*iv)=="x"||(*iv)=="y"||(*iv)=="charge"||
        (*iv)=="size"||(*iv)=="sizeX"||(*iv)=="sizeY"||(*iv)=="minrow"||
        (*iv)=="maxrow"||(*iv)=="mincol"||(*iv)=="maxcol")
      prefix="SUMCLU";
    if(currDir.find("Track")!=string::npos) prefix="SUMTRK";
    else if((*iv)=="residualX"||(*iv)=="residualY")
      prefix="SUMTRK";
    else if((*iv)=="ClustX"||(*iv)=="ClustY"||(*iv)=="nRecHits"||(*iv)=="ErrorX"||(*iv)=="ErrorY")
      prefix="SUMHIT";
    else if((*iv)=="Gain1d"||(*iv)=="GainChi2NDF1d"||
        (*iv)=="GainChi2Prob1d"||(*iv)=="Pedestal1d"||
        (*iv)=="GainNPoints1d"||(*iv)=="GainHighPoint1d"||
        (*iv)=="GainLowPoint1d"||(*iv)=="GainEndPoint1d"||
        (*iv)=="GainFitResult2d"||(*iv)=="GainDynamicRange2d"||
        (*iv)=="GainSaturate2d"||
        (*iv)=="ScurveChi2ProbSummary"||(*iv)=="ScurveFitResultSummary"||
        (*iv)=="ScurveSigmasSummary"||(*iv)=="ScurveThresholdSummary"||
        (*iv)=="pixelAliveSummary"  || (*iv) == "SiPixelErrorsCalibDigis") 
      prefix="SUMCAL"; 
      }
      MonitorElement* temp; string tag;
      if((*iv).find("residual")!=string::npos){                           // track residuals
    tag = prefix + "_" + (*iv) + "_mean_" 
      + currDir.substr(currDir.find(dir_name));
    temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
    sum_mes.push_back(temp);
    tag = prefix + "_" + (*iv) + "_RMS_" 
      + currDir.substr(currDir.find(dir_name));
    temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
    sum_mes.push_back(temp);
      }else if(prefix == "SUMCAL"){                  // calibrations
    if((*iv)=="Gain1d" || (*iv)=="GainChi2NDF1d" || (*iv)=="GainChi2Prob1d" ||
       (*iv)=="GainNPoints1d" || (*iv)=="GainHighPoint1d" ||
       (*iv)=="GainLowPoint1d" || (*iv)=="GainEndPoint1d" || 
       (*iv)=="GainDynamicRange2d" || (*iv)=="GainSaturate2d" ||
       (*iv)=="Pedestal1d" ||
       (*iv)=="ScurveChi2ProbSummary" || (*iv)=="ScurveFitResultSummary" ||
       (*iv)=="ScurveSigmasSummary" || (*iv)=="ScurveThresholdSummary"){                    
      tag = prefix + "_" + (*iv) + "_mean_" 
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
      tag = prefix + "_" + (*iv) + "_RMS_" 
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
    }else if((*iv) == "SiPixelErrorsCalibDigis"){
      tag = prefix + "_" + (*iv) + "_NCalibErrors_"
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
    }else if((*iv)=="GainFitResult2d"){
      tag = prefix + "_" + (*iv) + "_NNegativeFits_"
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
    }else if((*iv)=="pixelAliveSummary"){
      tag = prefix + "_" + (*iv) + "_FracOfPerfectPix_"
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
      tag = prefix + "_" + (*iv) + "_mean_"
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
    }
      }else{
    tag = prefix + "_" + (*iv) + "_" + currDir.substr(currDir.find(dir_name));
    temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
    sum_mes.push_back(temp);
    if((*iv)=="ndigis"){
      tag = prefix + "_" + (*iv) + "FREQ_" 
        + currDir.substr(currDir.find(dir_name));
      temp = getSummaryME(iBooker, iGetter, tag, isUpgrade);
      sum_mes.push_back(temp);
    }
    if(prefix=="SUMDIG" && (*iv)=="adc"){
      tag = "ALLMODS_" + (*iv) + "COMB_" + currDir.substr(currDir.find(dir_name));
      temp = 0;
      string fullpathname = iBooker.pwd() + "/" + tag;
      temp = iGetter.get(fullpathname);
      if (temp) {
        temp->Reset();
      }else{
        temp = iBooker.book1D(tag.c_str(), tag.c_str(),128, 0., 256.);
      }
      sum_mes.push_back(temp);
    }
    if(prefix=="SUMCLU" && (*iv)=="charge"){
      tag = "ALLMODS_" + (*iv) + "COMB_" + currDir.substr(currDir.find(dir_name));
      temp = 0;
      string fullpathname = iBooker.pwd() + "/" + tag;
      temp = iGetter.get(fullpathname);
      if (temp) {
        temp->Reset();
      }else{
        temp = iBooker.book1D(tag.c_str(), tag.c_str(),100, 0., 200.); // To look to get the size automatically   
      }
      sum_mes.push_back(temp);
    }
      }
    }
    if (sum_mes.size() == 0) {
      edm::LogInfo("SiPixelActionExecutor") << " Summary MEs can not be created" << "\n" ;
      return;
    }
    vector<string> subdirs = iGetter.getSubdirs();
    //Blade
    if(dir_name.find("Blade_") == 0) GetBladeSubdirs(iBooker, iGetter, subdirs);
    
    int ndet = 0;
    for (vector<string>::const_iterator it = subdirs.begin();
     it != subdirs.end(); it++) {
      if (prefix!="SUMOFF" && (*it).find("Module_") == string::npos) continue;
      if (prefix=="SUMOFF" && (*it).find(isbarrel?"Layer_":"Disk_") == string::npos) continue;
      iBooker.cd(*it);
      iGetter.cd(*it);
      ndet++;

      vector<string> contents = iGetter.getMEs(); 
      
      for (vector<MonitorElement*>::const_iterator isum = sum_mes.begin();
       isum != sum_mes.end(); isum++) {
    for (vector<string>::const_iterator im = contents.begin();
         im != contents.end(); im++) {
      string sname = ((*isum)->getName());
      string tname = " ";
      tname = sname.substr(7,(sname.find("_",7)-6));
      if(sname.find("ALLMODS_adcCOMB_")!=string::npos) tname = "adc_";
      if(sname.find("ALLMODS_chargeCOMB_")!=string::npos) tname = "charge_";
      if(sname.find("_charge_")!=string::npos && sname.find("Track_")==string::npos) tname = "charge_";
      if(sname.find("_nclusters_")!=string::npos && sname.find("Track_")==string::npos) tname = "nclusters_";
      if(sname.find("_size_")!=string::npos && sname.find("Track_")==string::npos) tname = "size_";
      if(sname.find("_charge_OffTrack_")!=string::npos) tname = "charge_OffTrack_";
      if(sname.find("_nclusters_OffTrack_")!=string::npos) tname = "nclusters_OffTrack_";
      if(sname.find("_size_OffTrack_")!=string::npos) tname = "size_OffTrack_";
      if(sname.find("_sizeX_OffTrack_")!=string::npos) tname = "sizeX_OffTrack_";
      if(sname.find("_sizeY_OffTrack_")!=string::npos) tname = "sizeY_OffTrack_";
      if(sname.find("_charge_OnTrack_")!=string::npos) tname = "charge_OnTrack_";
      if(sname.find("_nclusters_OnTrack_")!=string::npos) tname = "nclusters_OnTrack_";
      if(sname.find("_size_OnTrack_")!=string::npos) tname = "size_OnTrack_";
      if(sname.find("_sizeX_OnTrack_")!=string::npos) tname = "sizeX_OnTrack_";
      if(sname.find("_sizeY_OnTrack_")!=string::npos) tname = "sizeY_OnTrack_";
      if(tname.find("FREQ")!=string::npos) tname = "ndigis_";
      if (((*im)).find(tname) == 0) {
        string fullpathname = iBooker.pwd() + "/" + (*im); 
        MonitorElement *  me = iGetter.get(fullpathname);
        
        if (me){ 
          if(sname.find("_charge")!=string::npos && sname.find("Track_")==string::npos && me->getName().find("Track_")!=string::npos) continue;
          if(sname.find("_nclusters_")!=string::npos && sname.find("Track_")==string::npos && me->getName().find("Track_")!=string::npos) continue;
          if(sname.find("_size")!=string::npos && sname.find("Track_")==string::npos && me->getName().find("Track_")!=string::npos) continue;
          // fill summary histos:
          if (sname.find("_RMS_")!=string::npos && 
          sname.find("GainDynamicRange2d")==string::npos && 
          sname.find("GainSaturate2d")==string::npos){
        (*isum)->Fill(ndet, me->getRMS());
          }else if (sname.find("GainDynamicRange2d")!=string::npos ||
            sname.find("GainSaturate2d")!=string::npos){
        float SumOfEntries=0.; float SumOfSquaredEntries=0.; int SumOfPixels=0;
        for(int cols=1; cols!=me->getNbinsX()+1; cols++) for(int rows=1; rows!=me->getNbinsY()+1; rows++){
          SumOfEntries+=me->getBinContent(cols,rows);
          SumOfSquaredEntries+=(me->getBinContent(cols,rows))*(me->getBinContent(cols,rows));
          SumOfPixels++;
          }

        float MeanInZ = SumOfEntries / float(SumOfPixels);
        float RMSInZ = sqrt(SumOfSquaredEntries/float(SumOfPixels));
        if(sname.find("_mean_")!=string::npos) (*isum)->Fill(ndet, MeanInZ);
        if(sname.find("_RMS_")!=string::npos) (*isum)->Fill(ndet, RMSInZ);
          }else if (sname.find("_FracOfPerfectPix_")!=string::npos){
        float nlast = me->getBinContent(me->getNbinsX());
        float nall = (me->getTH1F())->Integral(1,11);
        (*isum)->Fill(ndet, nlast/nall);
          }else if (sname.find("_NCalibErrors_")!=string::npos ||
            sname.find("FREQ_")!=string::npos){
        float nall = me->getEntries();
        (*isum)->Fill(ndet, nall);
          }else if (sname.find("GainFitResult2d")!=string::npos){
        int NegFitPixels=0;
        for(int cols=1; cols!=me->getNbinsX()+1; cols++) for(int rows=1; rows!=me->getNbinsY()+1; rows++){
          if(me->getBinContent(cols,rows)<0.) NegFitPixels++;
        }
        (*isum)->Fill(ndet, float(NegFitPixels));
          }else if (sname.find("ALLMODS_adcCOMB_")!=string::npos || 
                    (sname.find("ALLMODS_chargeCOMB_")!=string::npos && me->getName().find("Track_")==string::npos)){
        (*isum)->getTH1F()->Add(me->getTH1F());
          }else if (sname.find("_NErrors_")!=string::npos){
            string path1 = fullpathname;
        path1 = path1.replace(path1.find("NErrors"),7,"errorType");
        MonitorElement * me1 = iGetter.get(path1);
        bool notReset=true;
            if(me1){
              for(int jj=1; jj<16; jj++){
                if(me1->getBinContent(jj)>0.){
                  if(jj==6){ //errorType=30 (reset)
                    string path2 = path1;
            path2 = path2.replace(path2.find("errorType"),9,"TBMMessage");
                    MonitorElement * me2 = iGetter.get(path2);
                    if(me2) if(me2->getBinContent(6)>0. || me2->getBinContent(7)>0.) notReset=false;
              }
            }
          }
        }
        if(notReset) (*isum)->Fill(ndet, me1->getEntries());
          }else if ((sname.find("_charge_")!=string::npos && sname.find("Track_")==string::npos && 
                    me->getName().find("Track_")==string::npos) ||
            (sname.find("_charge_")!=string::npos && sname.find("_OnTrack_")!=string::npos && 
                    me->getName().find("_OnTrack_")!=string::npos) ||
            (sname.find("_charge_")!=string::npos && sname.find("_OffTrack_")!=string::npos && 
                    me->getName().find("_OffTrack_")!=string::npos) ||
            (sname.find("_nclusters_")!=string::npos && sname.find("Track_")==string::npos && 
                    me->getName().find("Track_")==string::npos) ||
            (sname.find("_nclusters_")!=string::npos && sname.find("_OnTrack_")!=string::npos && 
                    me->getName().find("_OnTrack_")!=string::npos) ||
            (sname.find("_nclusters_")!=string::npos && sname.find("_OffTrack_")!=string::npos && 
                    me->getName().find("_OffTrack_")!=string::npos) ||
            (sname.find("_size")!=string::npos && sname.find("Track_")==string::npos && 
                    me->getName().find("Track_")==string::npos) ||
            (sname.find("_size")!=string::npos && sname.find("_OnTrack_")!=string::npos && 
                    me->getName().find("_OnTrack_")!=string::npos) ||
            (sname.find("_size")!=string::npos && sname.find("_OffTrack_")!=string::npos && 
                    me->getName().find("_OffTrack_")!=string::npos)){
            (*isum)->Fill(ndet, me->getMean());
          }else if(sname.find("_charge_")==string::npos && sname.find("_nclusters_")==string::npos && sname.find("_size")==string::npos){
        (*isum)->Fill(ndet, me->getMean());
          }
          
          // set titles:
          if(prefix=="SUMOFF"){
        (*isum)->setAxisTitle(isbarrel?"Ladders":"Blades",1);
          }else if(sname.find("ALLMODS_adcCOMB_")!=string::npos){
        (*isum)->setAxisTitle("Digi charge [ADC]",1);
          }else if(sname.find("ALLMODS_chargeCOMB_")!=string::npos){
        (*isum)->setAxisTitle("Cluster charge [kilo electrons]",1);
          }else{
        (*isum)->setAxisTitle("Modules",1);
          }
          string title = " ";
          if (sname.find("_RMS_")!=string::npos){
        title = "RMS of " + sname.substr(7,(sname.find("_",7)-7)) + " per module"; 
          }else if (sname.find("_FracOfPerfectPix_")!=string::npos){
        title = "FracOfPerfectPix " + sname.substr(7,(sname.find("_",7)-7)) + " per module"; 
          }else if(sname.find("_NCalibErrors_")!=string::npos){
        title = "Number of CalibErrors " + sname.substr(7,(sname.find("_",7)-7)) + " per module"; 
          }else if(sname.find("_NNegativeFits_")!=string::npos){
        title = "Number of pixels with neg. fit result " + sname.substr(7,(sname.find("_",7)-7)) + " per module"; 
          }else if (sname.find("FREQ_")!=string::npos){
        title = "NEvents with digis per module"; 
          }else if (sname.find("ALLMODS_adcCOMB_")!=string::npos){
        title = "NDigis";
          }else if (sname.find("ALLMODS_chargeCOMB_")!=string::npos){
        title = "NClusters";
          }else if (sname.find("_NErrors_")!=string::npos){
            if(prefix=="SUMOFF" && isbarrel) title = "Total number of errors per Ladder";
        else if(prefix=="SUMOFF" && !isbarrel) title = "Total number of errors per Blade";
        else title = "Total number of errors per Module";
          }else{
        if(prefix=="SUMOFF") title = "Mean " + sname.substr(7,(sname.find("_",7)-7)) + (isbarrel?" per Ladder":" per Blade"); 
        else title = "Mean " + sname.substr(7,(sname.find("_",7)-7)) + " per Module"; 
          }
          (*isum)->setAxisTitle(title,2);
        }
        break;
      }
    }
      }
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
      if(dir_name.find("Blade") == 0) {
    iBooker.goUp(); // Going up a second time if we are processing the Blade
    iGetter.setCurrentFolder(iBooker.pwd());
      }
    } // end for it (subdirs)
  } else {  
    vector<string> subdirs = iGetter.getSubdirs();
    // printing cout << "#\t" << iBooker.pwd() << endl;
    if(isbarrel)
      {
            
    for (vector<string>::const_iterator it = subdirs.begin();
         it != subdirs.end(); it++) {
      //                 cout << "##\t" << iBooker.pwd() << "\t" << (*it) << endl;
      if((iBooker.pwd()).find("Endcap")!=string::npos ||
         (iBooker.pwd()).find("AdditionalPixelErrors")!=string::npos) {
        iBooker.goUp();
        iGetter.setCurrentFolder(iBooker.pwd());
      }
      iBooker.cd(*it);
      iGetter.cd(*it);
      if((*it).find("Endcap")!=string::npos ||
         (*it).find("AdditionalPixelErrors")!=string::npos) continue;
      fillSummary(iBooker, iGetter, dir_name, me_names, true, isUpgrade); // Barrel
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
    }
    string grandbarrel_structure_name;
    vector<string> grandbarrel_me_names;
    if (!configParser_->getMENamesForGrandBarrelSummary(grandbarrel_structure_name, grandbarrel_me_names)){
      cout << "SiPixelActionExecutor::createSummary: Failed to read Grand Barrel Summary configuration parameters!! ";
      return;
    }
    fillGrandBarrelSummaryHistos(iBooker, iGetter, grandbarrel_me_names, isUpgrade);
            
      }
    else // Endcap
      {
            
    for (vector<string>::const_iterator it = subdirs.begin();
         it != subdirs.end(); it++) {
      if((iBooker.pwd()).find("Barrel")!=string::npos ||
         (iBooker.pwd()).find("AdditionalPixelErrors")!=string::npos){
        iBooker.goUp();
        iGetter.setCurrentFolder(iBooker.pwd());
      }
      iBooker.cd(*it);
      iGetter.cd(*it);
      if ((*it).find("Barrel")!=string::npos ||
          (*it).find("AdditionalPixelErrors")!=string::npos) continue;
      fillSummary(iBooker, iGetter, dir_name, me_names, false, isUpgrade); // Endcap
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
    }
    string grandendcap_structure_name;
    vector<string> grandendcap_me_names;
    if (!configParser_->getMENamesForGrandEndcapSummary(grandendcap_structure_name, grandendcap_me_names)){
      cout << "SiPixelActionExecutor::createSummary: Failed to read Grand Endcap Summary configuration parameters!! ";
      return;
    }
    fillGrandEndcapSummaryHistos(iBooker, iGetter, grandendcap_me_names, isUpgrade);
            
            
      }
  }
//  cout<<"...leaving SiPixelActionExecutor::fillSummary!"<<endl;
    
}

//=============================================================================================================
void SiPixelActionExecutor::fillFEDErrorSummary(DQMStore::IBooker& iBooker,
                        DQMStore::IGetter& iGetter,
                                                string dir_name,
                        vector<string>& me_names) {
  //printing cout<<"entering SiPixelActionExecutor::fillFEDErrorSummary..."<<endl;
  string currDir = iBooker.pwd();
  string prefix;
  if(source_type_==0) prefix="SUMRAW";
  else if(source_type_==20) prefix="SUMOFF";
  
  if (currDir.find(dir_name) != string::npos)  {
    vector<MonitorElement*> sum_mes;
    for (vector<string>::const_iterator iv = me_names.begin();
     iv != me_names.end(); iv++) {
      bool isBooked = false;
      vector<string> contents = iGetter.getMEs();
      for (vector<string>::const_iterator im = contents.begin(); im != contents.end(); im++)
        if ((*im).find(*iv) != string::npos) isBooked = true;
      if(source_type_==5||source_type_==6){
    if((*iv)=="errorType"||(*iv)=="NErrors"||(*iv)=="fullType"||(*iv)=="chanNmbr"||
       (*iv)=="TBMType"||(*iv)=="EvtNbr"||(*iv)=="evtSize"||(*iv)=="linkId"||
       (*iv)=="ROCId"||(*iv)=="DCOLId"||(*iv)=="PXId"||(*iv)=="ROCNmbr"||
       (*iv)=="TBMMessage"||(*iv)=="Type36Hitmap"||
       (*iv)=="FedChLErr"||(*iv)=="FedChNErr"||(*iv)=="FedETypeNErr") 
      prefix="SUMRAW";
      }
      if((*iv)=="errorType"||(*iv)=="NErrors"||(*iv)=="fullType"||(*iv)=="chanNmbr"||
     (*iv)=="TBMType"||(*iv)=="EvtNbr"||(*iv)=="evtSize"||(*iv)=="linkId"||
     (*iv)=="ROCId"||(*iv)=="DCOLId"||(*iv)=="PXId"||(*iv)=="ROCNmbr"||
     (*iv)=="TBMMessage"||(*iv)=="Type36Hitmap"){
        string tag = prefix + "_" + (*iv) + "_FEDErrors";
        MonitorElement* temp = getFEDSummaryME(iBooker, iGetter, tag);
        sum_mes.push_back(temp);
      }else if((*iv)=="FedChLErr"||(*iv)=="FedChNErr"||(*iv)=="FedETypeNErr"){
        string tag = prefix + "_" + (*iv);
    MonitorElement* temp;
    if((*iv)=="FedChLErr") {if (!isBooked) temp = iBooker.book2D("FedChLErr","Type of last error",40,-0.5,39.5,37,0.,37.);
      else{ 
        string fullpathname = iBooker.pwd() + "/" + (*iv);
        temp = iGetter.get(fullpathname);
        temp->Reset();}}  //If I don't reset this one, then I instead start adding error codes..
    if((*iv)=="FedChNErr") {if (!isBooked) temp = iBooker.book2D("FedChNErr","Total number of errors",40,-0.5,39.5,37,0.,37.);
      else{ 
        string fullpathname = iBooker.pwd() + "/" + (*iv);
        temp = iGetter.get(fullpathname);
        temp->Reset();}}  //If I don't reset this one, then I instead start adding error codes..
    if((*iv)=="FedETypeNErr"){
      if(!isBooked){
        temp = iBooker.book2D("FedETypeNErr","Number of each error type",40,-0.5,39.5,21,0.,21.);
        temp->setBinLabel(1,"ROC of 25",2);
        temp->setBinLabel(2,"Gap word",2);
        temp->setBinLabel(3,"Dummy word",2);
        temp->setBinLabel(4,"FIFO full",2);
        temp->setBinLabel(5,"Timeout",2);
        temp->setBinLabel(6,"Stack full",2);
        temp->setBinLabel(7,"Pre-cal issued",2);
        temp->setBinLabel(8,"Trigger clear or sync",2);
        temp->setBinLabel(9,"No token bit",2);
        temp->setBinLabel(10,"Overflow",2);
        temp->setBinLabel(11,"FSM error",2);
        temp->setBinLabel(12,"Invalid #ROCs",2);
        temp->setBinLabel(13,"Event number",2);
        temp->setBinLabel(14,"Slink header",2);
        temp->setBinLabel(15,"Slink trailer",2);
        temp->setBinLabel(16,"Event size",2);
        temp->setBinLabel(17,"Invalid channel#",2);
        temp->setBinLabel(18,"ROC value",2);
        temp->setBinLabel(19,"Dcol or pixel value",2);
        temp->setBinLabel(20,"Readout order",2);
        temp->setBinLabel(21,"CRC error",2);
      }
      else{
        string fullpathname = iBooker.pwd() + "/" + (*iv);
        temp = iGetter.get(fullpathname);
        temp->Reset();}  //If I don't reset this one, then I instead start adding error codes..
    }
    sum_mes.push_back(temp);
      }
    }
    if (sum_mes.size() == 0) {
      edm::LogInfo("SiPixelActionExecutor") << " Summary MEs can not be created" << "\n" ;
      return;
    }
    vector<string> subdirs = iGetter.getSubdirs();
    int ndet = 0;
    for (vector<string>::const_iterator it = subdirs.begin();
     it != subdirs.end(); it++) {
      if ( (*it).find("FED_") == string::npos) continue;
      iBooker.cd(*it);
      iGetter.cd(*it);
      string fedid = (*it).substr((*it).find("_")+1);
      std::istringstream isst;
      isst.str(fedid);
      isst>>ndet; ndet++;
      vector<string> contents = iGetter.getMEs(); 
            
      for (vector<MonitorElement*>::const_iterator isum = sum_mes.begin();
       isum != sum_mes.end(); isum++) {
    for (vector<string>::const_iterator im = contents.begin();
         im != contents.end(); im++) {
      if(((*im).find("FedChNErr")!=std::string::npos && (*isum)->getName().find("FedChNErr")!=std::string::npos) ||
         ((*im).find("FedChLErr")!=std::string::npos && (*isum)->getName().find("FedChLErr")!=std::string::npos) ||
         ((*im).find("FedETypeNErr")!=std::string::npos && (*isum)->getName().find("FedETypeNErr")!=std::string::npos)){
        string fullpathname = iBooker.pwd() + "/" + (*im); 
        MonitorElement *  me = iGetter.get(fullpathname);
        if(me){
          for(int i=0; i!=37; i++){
                if((*im).find("FedETypeNErr")!=std::string::npos && i<21) (*isum)->Fill(ndet-1,i,me->getBinContent(i+1));
                else (*isum)->Fill(ndet-1,i,me->getBinContent(i+1));
          }
        }
      }
      string sname = ((*isum)->getName());
      string tname = " ";
      tname = sname.substr(7,(sname.find("_",7)-6));
      if (((*im)).find(tname) == 0) {
        string fullpathname = iBooker.pwd() + "/" + (*im); 
        MonitorElement *  me = iGetter.get(fullpathname);
                        
        if (me){ 
          if(me->getMean()>0.){
            if (sname.find("_NErrors_")!=string::npos){
              string path1 = fullpathname;
          path1 = path1.replace(path1.find("NErrors"),7,"errorType");
          MonitorElement * me1 = iGetter.get(path1);
          bool notReset=true;
              if(me1){
                for(int jj=1; jj<16; jj++){
                  if(me1->getBinContent(jj)>0.){
                    if(jj==6){ //errorType=30 (reset)
                      string path2 = path1;
              path2 = path2.replace(path2.find("errorType"),9,"TBMMessage");
                      MonitorElement * me2 = iGetter.get(path2);
                      if(me2) if(me2->getBinContent(6)>0. || me2->getBinContent(7)>0.) notReset=false; 
                }
              }
            }
          }
          if(notReset) (*isum)->setBinContent(ndet, (*isum)->getBinContent(ndet) + me1->getEntries());
            }else (*isum)->setBinContent(ndet, (*isum)->getBinContent(ndet) + me->getEntries());
          }
          (*isum)->setAxisTitle("FED #",1);
          string title = " ";
          title = sname.substr(7,(sname.find("_",7)-7)) + " per FED"; 
          (*isum)->setAxisTitle(title,2);
        }
        break;
      }
    }
      }
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
    }
  } else {  
    vector<string> subdirs = iGetter.getSubdirs();
    for (vector<string>::const_iterator it = subdirs.begin();
     it != subdirs.end(); it++) {
      if((*it).find("Endcap")!=string::npos ||
     (*it).find("Barrel")!=string::npos) continue;
      iBooker.cd(*it);
      iGetter.cd(*it);
      fillFEDErrorSummary(iBooker,iGetter, dir_name, me_names);
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
    }
  }
  //printing cout<<"...leaving SiPixelActionExecutor::fillFEDErrorSummary!"<<endl;
}


//=============================================================================================================
void SiPixelActionExecutor::fillGrandBarrelSummaryHistos(DQMStore::IBooker & iBooker,
                             DQMStore::IGetter & iGetter,
                                                         vector<string>& me_names,
                                                         bool isUpgrade) {
//  cout<<"Entering SiPixelActionExecutor::fillGrandBarrelSummaryHistos...:"<<me_names.size()<<endl;
  vector<MonitorElement*> gsum_mes;
  string currDir = iBooker.pwd();
  string path_name = iBooker.pwd();
  string dir_name =  path_name.substr(path_name.find_last_of("/")+1);
  if ((dir_name.find("DQMData") == 0) ||
      (dir_name.find("Pixel") == 0) ||
      (dir_name.find("AdditionalPixelErrors") == 0) ||
      (dir_name.find("Endcap") == 0) ||
      (dir_name.find("HalfCylinder") == 0) ||
      (dir_name.find("Disk") == 0) ||
      (dir_name.find("Blade") == 0) ||
      (dir_name.find("Panel") == 0) ) return;
  vector<string> subdirs = iGetter.getSubdirs();
  int nDirs = subdirs.size();
  int iDir =0;
  int nbin = 0;
  int nbin_i = 0; 
  int nbin_subdir = 0; 
  int cnt=0;
  bool first_subdir = true;
  for (vector<string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
    cnt++;
    iBooker.cd(*it);
    iGetter.cd(*it);
    vector<string> contents = iGetter.getMEs();
        
    iBooker.goUp();
    iGetter.setCurrentFolder(iBooker.pwd());
        
    string prefix;
    if(source_type_==0) prefix="SUMRAW";
    else if (source_type_==1) prefix="SUMDIG";
    else if (source_type_==2) prefix="SUMCLU";
    else if (source_type_==3) prefix="SUMTRK";
    else if (source_type_==4) prefix="SUMHIT";
    else if (source_type_>=7 && source_type_<20) prefix="SUMCAL";
    else if (source_type_==20) prefix="SUMOFF";
    
        
    for (vector<string>::const_iterator im = contents.begin();
     im != contents.end(); im++) {
      for (vector<string>::const_iterator iv = me_names.begin();
       iv != me_names.end(); iv++) {
    string var = "_" + (*iv) + "_";
    if ((*im).find(var) != string::npos) {
      if((var=="_charge_" || var=="_nclusters_" || var=="_size_" || var=="_sizeX_" || var=="_sizeY_") && 
         (*im).find("Track_")!=string::npos) continue;
      string full_path = (*it) + "/" +(*im);
      MonitorElement * me = iGetter.get(full_path.c_str());
      if (!me) continue; 
      if(source_type_==5||source_type_==6){
        if((*iv)=="errorType"||(*iv)=="NErrors"||(*iv)=="fullType"||(*iv)=="chanNmbr"||
           (*iv)=="TBMType"||(*iv)=="EvtNbr"||(*iv)=="evtSize"||(*iv)=="linkId"||
           (*iv)=="ROCId"||(*iv)=="DCOLId"||(*iv)=="PXId"||(*iv)=="ROCNmbr"||
           (*iv)=="TBMMessage"||(*iv)=="Type36Hitmap") 
          prefix="SUMRAW";
        else if((*iv)=="ndigis"||(*iv)=="adc" ||
            (*iv)=="ndigisFREQ" || (*iv)=="adcCOMB")
          prefix="SUMDIG";
        else if((*iv)=="nclusters"||(*iv)=="x"||(*iv)=="y"||(*iv)=="charge"||(*iv)=="chargeCOMB"||
            (*iv)=="size"||(*iv)=="sizeX"||(*iv)=="sizeY"||(*iv)=="minrow"||
            (*iv)=="maxrow"||(*iv)=="mincol"||(*iv)=="maxcol")
          prefix="SUMCLU";
        if(currDir.find("Track")!=string::npos) prefix="SUMTRK";
        else if((*iv)=="residualX_mean"||(*iv)=="residualY_mean"||
            (*iv)=="residualX_RMS"||(*iv)=="residualY_RMS")
          prefix="SUMTRK";
        else if((*iv)=="ClustX"||(*iv)=="ClustY"||(*iv)=="nRecHits"||(*iv)=="ErrorX"||(*iv)=="ErrorY")
          prefix="SUMHIT";
        else if((*iv)=="Gain1d_mean"||(*iv)=="GainChi2NDF1d_mean"||
            (*iv)=="GainChi2Prob1d_mean"||(*iv)=="Pedestal1d_mean"||
            (*iv)=="ScurveChi2ProbSummary_mean"||(*iv)=="ScurveFitResultSummary_mean"||
            (*iv)=="ScurveSigmasSummary_mean"||(*iv)=="ScurveThresholdSummary_mean"||
            (*iv)=="Gain1d_RMS"||(*iv)=="GainChi2NDF1d_RMS"||
            (*iv)=="GainChi2Prob1d_RMS"||(*iv)=="Pedestal1d_RMS"||
            (*iv)=="GainNPoints1d_mean" || (*iv)=="GainNPoints1d_RMS" ||
            (*iv)=="GainHighPoint1d_mean" || (*iv)=="GainHighPoint1d_RMS" ||
            (*iv)=="GainLowPoint1d_mean" || (*iv)=="GainLowPoint1d_RMS" ||
            (*iv)=="GainEndPoint1d_mean" || (*iv)=="GainEndPoint1d_RMS" ||
            (*iv)=="GainFitResult2d_mean" || (*iv)=="GainFitResult2d_RMS" ||
            (*iv)=="GainDynamicRange2d_mean" || (*iv)=="GainDynamicRange2d_RMS" ||
            (*iv)=="GainSaturate2d_mean" || (*iv)=="GainSaturate2d_RMS" ||
            (*iv)=="ScurveChi2ProbSummary_RMS"||(*iv)=="ScurveFitResultSummary_RMS"||
            (*iv)=="ScurveSigmasSummary_RMS"||(*iv)=="ScurveThresholdSummary_RMS"||
            (*iv)=="pixelAliveSummary_mean"||(*iv)=="pixelAliveSummary_FracOfPerfectPix" ||
            (*iv)=="SiPixelErrorsCalibDigis_NCalibErrors" )
          prefix="SUMCAL";
      } // end source_type if
      
      if (first_subdir && !isUpgrade){
        nbin = me->getTH1F()->GetNbinsX();        
        string me_name = prefix + "_" + (*iv) + "_" + dir_name;
        if((*iv)=="adcCOMB"||(*iv)=="chargeCOMB") me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
        else if(prefix=="SUMOFF" && dir_name=="Barrel") nbin=192;
        else if((*iv)=="adcCOMB") nbin=256;
        else if(dir_name=="Barrel") nbin=768;
        else if(prefix=="SUMOFF" && dir_name.find("Shell")!=string::npos) nbin=48;
        else if(dir_name.find("Shell")!=string::npos) nbin=192;
        else nbin=nbin*nDirs;

        getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
      } else if (first_subdir && isUpgrade){
        nbin = me->getTH1F()->GetNbinsX();        
        string me_name = prefix + "_" + (*iv) + "_" + dir_name;
        if((*iv)=="adcCOMB"||(*iv)=="chargeCOMB") me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
        else if(prefix=="SUMOFF" && dir_name=="Barrel") nbin=296;
        else if((*iv)=="adcCOMB") nbin=256;
        else if(dir_name=="Barrel") nbin=1184;
        else if(prefix=="SUMOFF" && dir_name.find("Shell")!=string::npos) nbin=74;
        else if(dir_name.find("Shell")!=string::npos) nbin=296;
        else nbin=nbin*nDirs;

        getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
      }

      for (vector<MonitorElement*>::const_iterator igm = gsum_mes.begin();
           igm != gsum_mes.end(); igm++) {
        if ((*igm)->getName().find(var) != string::npos) {
          if(prefix=="SUMOFF") (*igm)->setAxisTitle("Ladders",1);
          else if((*igm)->getName().find("adcCOMB_")!=string::npos) (*igm)->setAxisTitle("Digi charge [ADC]",1);
          else if((*igm)->getName().find("chargeCOMB_")!=string::npos) (*igm)->setAxisTitle("Cluster charge [kilo electrons]",1);
          else (*igm)->setAxisTitle("Modules",1);

          // Setting title

          string title="";
          if((*igm)->getName().find("NErrors_") != string::npos && prefix=="SUMOFF") title = "Total number of errors per Ladder";
          else if((*igm)->getName().find("NErrors_") != string::npos && prefix=="SUMRAW") title = "Total number of errors per Module";
          else if(prefix=="SUMOFF") title = "mean " + (*iv) + " per Ladder"; 
          else if((*igm)->getName().find("FREQ_") != string::npos && prefix!="SUMOFF") title = "NEvents with digis per Module"; 
          else if((*igm)->getName().find("FREQ_") != string::npos && prefix=="SUMOFF") title = "NEvents with digis per Ladder/Blade"; 
          else if((*igm)->getName().find("adcCOMB_") != string::npos) title = "NDigis";
          else if((*igm)->getName().find("chargeCOMB_") != string::npos) title = "NClusters";
          else title = "mean " + (*iv) + " per Module"; 
          (*igm)->setAxisTitle(title,2);
          
          // Setting binning
          if (!isUpgrade) {
          if((*igm)->getName().find("ALLMODS_adcCOMB_")!=string::npos){
        nbin_subdir=128;
          }else if((*igm)->getName().find("ALLMODS_chargeCOMB_")!=string::npos){
        nbin_subdir=100;
          }else if((*igm)->getName().find("Ladder") != string::npos){
        nbin_i=0; nbin_subdir=4;
          }else if((*igm)->getName().find("Layer") != string::npos){
        nbin_i=(cnt-1)*4; nbin_subdir=4;
          }else if((*igm)->getName().find("Shell") != string::npos){
        if(prefix!="SUMOFF"){
          if(iDir==0){ nbin_i=0; nbin_subdir=40; }
          else if(iDir==1){ nbin_i=40; nbin_subdir=64; }
          else if(iDir==2){ nbin_i=104; nbin_subdir=88; }
        }else{
          if(iDir==0){ nbin_i=0; nbin_subdir=10; }
          else if(iDir==1){ nbin_i=10; nbin_subdir=16; }
          else if(iDir==2){ nbin_i=26; nbin_subdir=22; }
        }
          }else if((*igm)->getName().find("Barrel") != string::npos){
        if(prefix!="SUMOFF"){
          if(iDir==0){ nbin_i=0; nbin_subdir=192; }
          else if(iDir==1){ nbin_i=192; nbin_subdir=192; }
          else if(iDir==2){ nbin_i=384; nbin_subdir=192; }
          else if(iDir==3){ nbin_i=576; nbin_subdir=192; }
        }else{
          if(iDir==0){ nbin_i=0; nbin_subdir=48; }
          else if(iDir==1){ nbin_i=48; nbin_subdir=48; }
          else if(iDir==2){ nbin_i=96; nbin_subdir=48; }
          else if(iDir==3){ nbin_i=144; nbin_subdir=48; }
        }
          }
          } else if (isUpgrade) {
            if((*igm)->getName().find("ALLMODS_adcCOMB_")!=string::npos){
          nbin_subdir=128;
            }else if((*igm)->getName().find("ALLMODS_chargeCOMB_")!=string::npos){
          nbin_subdir=100;
            }else if((*igm)->getName().find("Ladder") != string::npos){
          nbin_i=0; nbin_subdir=4;
            }else if((*igm)->getName().find("Layer") != string::npos){
          nbin_i=(cnt-1)*4; nbin_subdir=4;
            }else if((*igm)->getName().find("Shell") != string::npos){
          if(prefix!="SUMOFF"){
           if(iDir==0){ nbin_i=0; nbin_subdir=24; }//40(2*20)-->24(2*12)
            else if(iDir==1){ nbin_i=24; nbin_subdir=56; }//64(32*2)-->56(2*28)
            else if(iDir==2){ nbin_i=80; nbin_subdir=88; }//88(44*2)-->same88(44*2)
            else if(iDir==3){ nbin_i=168; nbin_subdir=128; }
          }else{
            if(iDir==0){ nbin_i=0; nbin_subdir=6; }//10-->6
            else if(iDir==1){ nbin_i=6; nbin_subdir=14; }//16-->14
            else if(iDir==2){ nbin_i=20; nbin_subdir=22; }//22-->same22
            else if(iDir==3){ nbin_i=42; nbin_subdir=32; }
          }
            }else if((*igm)->getName().find("Barrel") != string::npos){
          if(prefix!="SUMOFF"){
            if(iDir==0){ nbin_i=0; nbin_subdir=296; }//192=76 8/4-->296=1184/4
            else if(iDir==1){ nbin_i=296; nbin_subdir=296; }//296*2,*3,*4=1184
            else if(iDir==2){ nbin_i=592; nbin_subdir=296; }
            else if(iDir==3){ nbin_i=888; nbin_subdir=296; }
            else if(iDir==4){ nbin_i=1184; nbin_subdir=296; }
          }else{
            if(iDir==0){ nbin_i=0; nbin_subdir=74; }//48=192/4-->74=296/4
            else if(iDir==1){ nbin_i=74; nbin_subdir=74; }//74*2,...*4=296
            else if(iDir==2){ nbin_i=148; nbin_subdir=74; }
            else if(iDir==3){ nbin_i=222; nbin_subdir=74; }
            else if(iDir==4){ nbin_i=296; nbin_subdir=74; }
          }
            }
          }
          
          if((*igm)->getName().find("ndigisFREQ")==string::npos)
        { 
          if(((*igm)->getName().find("adcCOMB")!=string::npos && me->getName().find("adcCOMB")!=string::npos) 
             || ((*igm)->getName().find("chargeCOMB")!=string::npos && me->getName().find("chargeCOMB")!=string::npos))
            {
              (*igm)->getTH1F()->Add(me->getTH1F());
            }else if(((*igm)->getName().find("charge_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
                  me->getName().find("charge_")!=string::npos && me->getName().find("Track_")==string::npos) || 
                 ((*igm)->getName().find("nclusters_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
                  me->getName().find("nclusters_")!=string::npos && me->getName().find("Track_")==string::npos) || 
                 ((*igm)->getName().find("size_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
                  me->getName().find("size_")!=string::npos && me->getName().find("Track_")==string::npos) || 
                 ((*igm)->getName().find("charge_OffTrack_")!=string::npos && me->getName().find("charge_OffTrack_")!=string::npos) || 
                 ((*igm)->getName().find("nclusters_OffTrack_")!=string::npos && me->getName().find("nclusters_OffTrack_")!=string::npos) || 
                 ((*igm)->getName().find("size_OffTrack_")!=string::npos && me->getName().find("size_OffTrack_")!=string::npos) || 
                 ((*igm)->getName().find("charge_OnTrack_")!=string::npos && me->getName().find("charge_OnTrack_")!=string::npos) || 
                 ((*igm)->getName().find("nclusters_OnTrack_")!=string::npos && me->getName().find("nclusters_OnTrack_")!=string::npos) || 
                 ((*igm)->getName().find("size_OnTrack_")!=string::npos && me->getName().find("size_OnTrack_")!=string::npos) || 
                 ((*igm)->getName().find("charge_")==string::npos && (*igm)->getName().find("nclusters_")==string::npos && 
                  (*igm)->getName().find("size_")==string::npos)){
            for (int k = 1; k < nbin_subdir+1; k++) if(me->getBinContent(k) > 0) (*igm)->setBinContent(k+nbin_i, me->getBinContent(k));
          }
        }
          else if(me->getName().find("ndigisFREQ")!=string::npos)
        {
          for (int k = 1; k < nbin_subdir+1; k++) if(me->getBinContent(k) > 0) (*igm)->setBinContent(k+nbin_i, me->getBinContent(k));
        }
        } // end var in igm (gsum_mes)
      } // end igm loop
    } // end var in im (contents)
      } // end of iv loop
    } // end of im loop
    iDir++;
  first_subdir = false; // We are done processing the first directory, we don't add any new MEs in the future passes.   
  } // end of it loop (subdirs)
//  cout<<"...leaving SiPixelActionExecutor::fillGrandBarrelSummaryHistos!"<<endl;
}

//=============================================================================================================
void SiPixelActionExecutor::fillGrandEndcapSummaryHistos(DQMStore::IBooker& iBooker,
                             DQMStore::IGetter& iGetter,
                                                         vector<string>& me_names,
                                                         bool isUpgrade) {
  //printing cout<<"Entering SiPixelActionExecutor::fillGrandEndcapSummaryHistos..."<<endl;
  vector<MonitorElement*> gsum_mes;
  string currDir = iBooker.pwd();
  string path_name = iBooker.pwd();
  string dir_name =  path_name.substr(path_name.find_last_of("/")+1);
  if ((dir_name.find("DQMData") == 0) ||
      (dir_name.find("Pixel") == 0) ||
      (dir_name.find("AdditionalPixelErrors") == 0) ||
      (dir_name.find("Barrel") == 0) ||
      (dir_name.find("Shell") == 0) ||
      (dir_name.find("Layer") == 0) ||
      (dir_name.find("Ladder") == 0) ) return;
  vector<string> subdirs = iGetter.getSubdirs();
  int iDir =0;
  int nbin = 0;
  int nbin_i = 0; 
  int nbin_subdir = 0; 
  int cnt=0;
  bool first_subdir = true;  
  for (vector<string>::const_iterator it = subdirs.begin();
       it != subdirs.end(); it++) {
    cnt++;
    iBooker.cd(*it);
    iGetter.cd(*it);
    vector<string> contents = iGetter.getMEs();
    iBooker.goUp();
    iGetter.setCurrentFolder(iBooker.pwd());
        
    string prefix;
    if(source_type_==0) prefix="SUMRAW";
    else if (source_type_==1) prefix="SUMDIG";
    else if (source_type_==2) prefix="SUMCLU";
    else if (source_type_==3) prefix="SUMTRK";
    else if (source_type_==4) prefix="SUMHIT";
    else if (source_type_>=7 && source_type_<20) prefix="SUMCAL";
    else if (source_type_==20) prefix="SUMOFF";
        
    for (vector<string>::const_iterator im = contents.begin();
     im != contents.end(); im++) {
      for (vector<string>::const_iterator iv = me_names.begin();
       iv != me_names.end(); iv++) {
    string var = "_" + (*iv) + "_";
    if ((*im).find(var) != string::npos) {
      if((var=="_charge_" || var=="_nclusters_" || var=="_size_" || var=="_sizeX_" || var=="_sizeY_") && 
         (*im).find("Track_")!=string::npos) continue;
      string full_path = (*it) + "/" +(*im);
      MonitorElement * me = iGetter.get(full_path.c_str());
      if (!me) continue; 
      if(source_type_==5||source_type_==6){
        if((*iv)=="errorType"||(*iv)=="NErrors"||(*iv)=="fullType"||(*iv)=="chanNmbr"||
           (*iv)=="TBMType"||(*iv)=="EvtNbr"||(*iv)=="evtSize"||(*iv)=="linkId"||
           (*iv)=="ROCId"||(*iv)=="DCOLId"||(*iv)=="PXId"||(*iv)=="ROCNmbr"||
           (*iv)=="TBMMessage"||(*iv)=="Type36Hitmap") 
          prefix="SUMRAW";
        else if((*iv)=="ndigis"||(*iv)=="adc" ||
            (*iv)=="ndigisFREQ"||(*iv)=="adcCOMB")
          prefix="SUMDIG";
        else if((*iv)=="nclusters"||(*iv)=="x"||(*iv)=="y"||(*iv)=="charge"||(*iv)=="chargeCOMB"||
            (*iv)=="size"||(*iv)=="sizeX"||(*iv)=="sizeY"||(*iv)=="minrow"||
            (*iv)=="maxrow"||(*iv)=="mincol"||(*iv)=="maxcol")
          prefix="SUMCLU";
        if(currDir.find("Track")!=string::npos) prefix="SUMTRK";
        else if((*iv)=="residualX_mean"||(*iv)=="residualY_mean"||
            (*iv)=="residualX_RMS"||(*iv)=="residualY_RMS")
          prefix="SUMTRK";
        else if((*iv)=="ClustX"||(*iv)=="ClustY"||(*iv)=="nRecHits"||(*iv)=="ErrorX"||(*iv)=="ErrorY")
          prefix="SUMHIT";
        else if((*iv)=="Gain1d_mean"||(*iv)=="GainChi2NDF1d_mean"||
            (*iv)=="GainChi2Prob1d_mean"||(*iv)=="Pedestal1d_mean"||
            (*iv)=="ScurveChi2ProbSummary_mean"||(*iv)=="ScurveFitResultSummary_mean"||
            (*iv)=="ScurveSigmasSummary_mean"||(*iv)=="ScurveThresholdSummary_mean"||
            (*iv)=="Gain1d_RMS"||(*iv)=="GainChi2NDF1d_RMS"||
            (*iv)=="GainChi2Prob1d_RMS"||(*iv)=="Pedestal1d_RMS"||
            (*iv)=="GainNPoints1d_mean" || (*iv)=="GainNPoints1d_RMS" ||
            (*iv)=="GainHighPoint1d_mean" || (*iv)=="GainHighPoint1d_RMS" ||
            (*iv)=="GainLowPoint1d_mean" || (*iv)=="GainLowPoint1d_RMS" ||
            (*iv)=="GainEndPoint1d_mean" || (*iv)=="GainEndPoint1d_RMS" ||
            (*iv)=="GainFitResult2d_mean" || (*iv)=="GainFitResult2d_RMS" ||
            (*iv)=="GainDynamicRange2d_mean" || (*iv)=="GainDynamicRange2d_RMS" ||
            (*iv)=="GainSaturate2d_mean" || (*iv)=="GainSaturate2d_RMS" ||
            (*iv)=="ScurveChi2ProbSummary_RMS"||(*iv)=="ScurveFitResultSummary_RMS"||
            (*iv)=="ScurveSigmasSummary_RMS"||(*iv)=="ScurveThresholdSummary_RMS"||
            (*iv)=="pixelAliveSummary_mean"||(*iv)=="pixelAliveSummary_FracOfPerfectPix"|| 
            (*iv) == "SiPixelErrorsCalibDigis_NCalibErrors")
          prefix="SUMCAL"; 
      }

      if (first_subdir && !isUpgrade){
        nbin = me->getTH1F()->GetNbinsX();        
        string me_name = prefix + "_" + (*iv) + "_" + dir_name;
        if((*iv)=="adcCOMB"||(*iv)=="chargeCOMB") me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
        else if(prefix=="SUMOFF" && dir_name=="Endcap") nbin=96;
        else if(dir_name=="Endcap") nbin=672;
        else if(prefix=="SUMOFF" && dir_name.find("HalfCylinder")!=string::npos) nbin=24;
        else if(dir_name.find("HalfCylinder")!=string::npos) nbin=168;
        else if(prefix=="SUMOFF" && dir_name.find("Disk")!=string::npos) nbin=12;
        else if(dir_name.find("Disk")!=string::npos) nbin=84;
        else if(dir_name.find("Blade")!=string::npos) nbin=7;
        getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
      } else if(first_subdir && isUpgrade){
            nbin = me->getTH1F()->GetNbinsX();        
        string me_name = prefix + "_" + (*iv) + "_" + dir_name;
        if((*iv)=="adcCOMB"||(*iv)=="chargeCOMB") me_name = "ALLMODS_" + (*iv) + "_" + dir_name;
        else if(prefix=="SUMOFF" && dir_name=="Endcap") nbin=336;
        else if(dir_name=="Endcap") nbin=672;
        else if(prefix=="SUMOFF" && dir_name.find("HalfCylinder")!=string::npos) nbin=84;
        else if(dir_name.find("HalfCylinder")!=string::npos) nbin=168;
        else if(prefix=="SUMOFF" && dir_name.find("Disk")!=string::npos) nbin=28;
        else if(dir_name.find("Disk")!=string::npos) nbin=56;
        else if(dir_name.find("Blade")!=string::npos) nbin=2;
            getGrandSummaryME(iBooker, iGetter, nbin, me_name, gsum_mes);
          }

      for (vector<MonitorElement*>::const_iterator igm = gsum_mes.begin();
           igm != gsum_mes.end(); igm++) { 
        if ((*igm)->getName().find(var) != string::npos) {
          if(prefix=="SUMOFF") (*igm)->setAxisTitle("Blades",1);
          else if((*igm)->getName().find("adcCOMB_")!=string::npos) (*igm)->setAxisTitle("Digi charge [ADC]",1);
          else if((*igm)->getName().find("chargeCOMB_")!=string::npos) (*igm)->setAxisTitle("Cluster charge [kilo electrons]",1);
          else (*igm)->setAxisTitle("Modules",1);
          string title="";
          if((*igm)->getName().find("NErrors_") != string::npos && prefix=="SUMOFF") title = "Total number of errors per Blade";
          else if((*igm)->getName().find("NErrors_") != string::npos && prefix=="SUMRAW") title = "Total number of errors per Module";
          else if(prefix=="SUMOFF") title = "mean " + (*iv) + " per Blade"; 
          else if((*igm)->getName().find("FREQ_") != string::npos) title = "NEvents with digis per Module"; 
          else if((*igm)->getName().find("adcCOMB_")!=string::npos) title = "NDigis";
          else if((*igm)->getName().find("chargeCOMB_")!=string::npos) title = "NClusters";
          else title = "mean " + (*iv) + " per Module"; 
          (*igm)->setAxisTitle(title,2);
          nbin_i=0; 
          if (!isUpgrade) {
          if((*igm)->getName().find("ALLMODS_adcCOMB_")!=string::npos){
        nbin_subdir=128;
          }else if((*igm)->getName().find("ALLMODS_chargeCOMB_")!=string::npos){
        nbin_subdir=100;
          }else if((*igm)->getName().find("Panel_") != string::npos){
        nbin_subdir=7;
          }else if((*igm)->getName().find("Blade") != string::npos){
        if((*im).find("_1") != string::npos) nbin_subdir=4;
        if((*im).find("_2") != string::npos) {nbin_i=4; nbin_subdir=3;}
          }else if((*igm)->getName().find("Disk") != string::npos){
        nbin_i=((cnt-1)%12)*7; nbin_subdir=7;
          }else if((*igm)->getName().find("HalfCylinder") != string::npos){
        if(prefix!="SUMOFF"){
          nbin_subdir=84;
          if((*im).find("_2") != string::npos) nbin_i=84;
        }else{
          nbin_subdir=12;
          if((*im).find("_2") != string::npos) nbin_i=12;
        }
          }else if((*igm)->getName().find("Endcap") != string::npos){
        if(prefix!="SUMOFF"){
          nbin_subdir=168;
          if((*im).find("_mO") != string::npos) nbin_i=168;
          if((*im).find("_pI") != string::npos) nbin_i=336;
          if((*im).find("_pO") != string::npos) nbin_i=504;
        }else{
          nbin_subdir=24;
          if((*im).find("_mO") != string::npos) nbin_i=24;
          if((*im).find("_pI") != string::npos) nbin_i=48;
          if((*im).find("_pO") != string::npos) nbin_i=72;
        }
          }
              } else if (isUpgrade) {
                if((*igm)->getName().find("ALLMODS_adcCOMB_")!=string::npos){
          nbin_subdir=128;
            }else if((*igm)->getName().find("ALLMODS_chargeCOMB_")!=string::npos){
          nbin_subdir=100;
            }else if((*igm)->getName().find("Panel_") != string::npos){
          nbin_subdir=2;
            }else if((*igm)->getName().find("Blade") != string::npos){
          if((*im).find("_1") != string::npos) nbin_subdir=1;
          if((*im).find("_2") != string::npos) {nbin_i=1; nbin_subdir=1;}
            }else if((*igm)->getName().find("Disk") != string::npos){
          nbin_i=((cnt-1)%28)*2; nbin_subdir=2;
            }else if((*igm)->getName().find("HalfCylinder") != string::npos){
          if(prefix!="SUMOFF"){
            nbin_subdir=56;
            if((*im).find("_2") != string::npos) nbin_i=56;
                    if((*im).find("_3") != string::npos) nbin_i=112;
          }else{
            nbin_subdir=28;
            if((*im).find("_2") != string::npos) nbin_i=28;
                    if((*im).find("_3") != string::npos) nbin_i=56;
          }
            }else if((*igm)->getName().find("Endcap") != string::npos){
          if(prefix!="SUMOFF"){
            nbin_subdir=168;
            if((*im).find("_mO") != string::npos) nbin_i=168;
            if((*im).find("_pI") != string::npos) nbin_i=336;
            if((*im).find("_pO") != string::npos) nbin_i=504;
          }else{
            nbin_subdir=84;
            if((*im).find("_mO") != string::npos) nbin_i=84;
            if((*im).find("_pI") != string::npos) nbin_i=168;
            if((*im).find("_pO") != string::npos) nbin_i=252;
          }
            }
              }
                            
          if((*igm)->getName().find("ndigisFREQ")==string::npos){ 
        if(((*igm)->getName().find("adcCOMB")!=string::npos && me->getName().find("adcCOMB")!=string::npos) || ((*igm)->getName().find("chargeCOMB")!=string::npos && me->getName().find("chargeCOMB")!=string::npos)){
          (*igm)->getTH1F()->Add(me->getTH1F());
        }else if(((*igm)->getName().find("charge_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
               me->getName().find("charge_")!=string::npos && me->getName().find("Track_")==string::npos) || 
              ((*igm)->getName().find("nclusters_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
               me->getName().find("nclusters_")!=string::npos && me->getName().find("Track_")==string::npos) || 
              ((*igm)->getName().find("size_")!=string::npos && (*igm)->getName().find("Track_")==string::npos && 
               me->getName().find("size_")!=string::npos && me->getName().find("Track_")==string::npos) || 
              ((*igm)->getName().find("charge_OffTrack_")!=string::npos && me->getName().find("charge_OffTrack_")!=string::npos) || 
              ((*igm)->getName().find("nclusters_OffTrack_")!=string::npos && me->getName().find("nclusters_OffTrack_")!=string::npos) || 
              ((*igm)->getName().find("size_OffTrack_")!=string::npos && me->getName().find("size_OffTrack_")!=string::npos) || 
              ((*igm)->getName().find("charge_OnTrack_")!=string::npos && me->getName().find("charge_OnTrack_")!=string::npos) || 
              ((*igm)->getName().find("nclusters_OnTrack_")!=string::npos && me->getName().find("nclusters_OnTrack_")!=string::npos) || 
              ((*igm)->getName().find("size_OnTrack_")!=string::npos && me->getName().find("size_OnTrack_")!=string::npos) || 
              ((*igm)->getName().find("charge_")==string::npos && (*igm)->getName().find("nclusters_")==string::npos && 
               (*igm)->getName().find("size_")==string::npos)){
          for (int k = 1; k < nbin_subdir+1; k++) if(me->getBinContent(k) > 0) (*igm)->setBinContent(k+nbin_i, me->getBinContent(k));
        }
          }else if(me->getName().find("ndigisFREQ")!=string::npos){
        for (int k = 1; k < nbin_subdir+1; k++)  if(me->getBinContent(k) > 0) (*igm)->setBinContent(k+nbin_i, me->getBinContent(k));
          }
          //           }// for
                            
        }
      }
    }
      }
    }
        
    iDir++;
  first_subdir = false; // We are done processing the first directory, we don't add any new MEs in the future passes.   
  }  // end for it (subdirs)
}
//=============================================================================================================
//
// -- Get Summary ME
//
void SiPixelActionExecutor::getGrandSummaryME(DQMStore::IBooker& iBooker,
                          DQMStore::IGetter & iGetter,
                                              int nbin, 
                          string& me_name, 
                          vector<MonitorElement*> & mes) {
  //printing cout<<"Entering SiPixelActionExecutor::getGrandSummaryME for: "<<me_name<<endl;
  if((iBooker.pwd()).find("Pixel")==string::npos) return; // If one doesn't find pixel
  vector<string> contents = iGetter.getMEs();
    
  for (vector<string>::const_iterator it = contents.begin();
       it != contents.end(); it++) {
    //printing cout<<"in grand summary me: "<<me_name<<","<<(*it)<<endl;
    if ((*it).find(me_name) == 0) {
      string fullpathname = iBooker.pwd() + "/" + me_name;
      MonitorElement* me = iGetter.get(fullpathname);
            
      if (me) {
    me->Reset();
    mes.push_back(me);
    return;
      }
    }
  }

  MonitorElement* temp_me(0);
  if(me_name.find("ALLMODS_adcCOMB_")!=string::npos) temp_me = iBooker.book1D(me_name.c_str(),me_name.c_str(),128,0,256);
  else if(me_name.find("ALLMODS_chargeCOMB_")!=string::npos) temp_me = iBooker.book1D(me_name.c_str(),me_name.c_str(),100,0,200);
  else temp_me = iBooker.book1D(me_name.c_str(),me_name.c_str(),nbin,1.,nbin+1.);
  if (temp_me) mes.push_back(temp_me);
    
  //  if(temp_me) cout<<"finally found grand ME: "<<me_name<<endl;
}


//=============================================================================================================
//
// -- Get Summary ME
//
MonitorElement* SiPixelActionExecutor::getSummaryME(DQMStore::IBooker & iBooker,
                            DQMStore::IGetter & iGetter,
                                                    string me_name,
                                                    bool isUpgrade) {
  //printing cout<<"Entering SiPixelActionExecutor::getSummaryME for: "<<me_name<<endl;
  MonitorElement* me = 0;
  if((iBooker.pwd()).find("Pixel")==string::npos) return me;
  vector<string> contents = iGetter.getMEs();    
    
  for (vector<string>::const_iterator it = contents.begin();
       it != contents.end(); it++) {
    if ((*it).find(me_name) == 0) {
      string fullpathname = iBooker.pwd() + "/" + (*it); 
      me = iGetter.get(fullpathname);
      if (me) {
    me->Reset();
    return me;
      }
    }
  }
  contents.clear();
  if (!isUpgrade) {
  if(me_name.find("SUMOFF")==string::npos){
    if(me_name.find("Blade_")!=string::npos)me = iBooker.book1D(me_name.c_str(), me_name.c_str(),7,1.,8.);
    else me = iBooker.book1D(me_name.c_str(), me_name.c_str(),4,1.,5.);
  }else if(me_name.find("Layer_1")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),10,1.,11.);
  }else if(me_name.find("Layer_2")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),16,1.,17.);
  }else if(me_name.find("Layer_3")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),22,1.,23.);
  }else if(me_name.find("Disk_")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),12,1.,13.);
  }
  }//endifNOTUpgrade
  else if (isUpgrade) {
    if(me_name.find("SUMOFF")==string::npos){
      if(me_name.find("Blade_")!=string::npos)me = iBooker.book1D(me_name.c_str(), me_name.c_str(),2,1.,3.);
      else me = iBooker.book1D(me_name.c_str(), me_name.c_str(),1,1.,2.);
    }else if(me_name.find("Layer_1")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),6,1.,7.);
    }else if(me_name.find("Layer_2")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),14,1.,15.);
    }else if(me_name.find("Layer_3")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),22,1.,23.);
    }else if(me_name.find("Layer_4")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),32,1.,33.);
    }else if(me_name.find("Disk_")!=string::npos){ me = iBooker.book1D(me_name.c_str(), me_name.c_str(),28,1.,29.);
    }
  }//endifUpgrade
    
  return me;
}


//=============================================================================================================
MonitorElement* SiPixelActionExecutor::getFEDSummaryME(DQMStore::IBooker & iBooker,
                               DQMStore::IGetter & iGetter,
                                                       string me_name) {
  //printing cout<<"Entering SiPixelActionExecutor::getFEDSummaryME..."<<endl;
  MonitorElement* me = 0;
  if((iBooker.pwd()).find("Pixel")==string::npos) return me;
  vector<string> contents = iGetter.getMEs();
    
  for (vector<string>::const_iterator it = contents.begin();
       it != contents.end(); it++) {
    if ((*it).find(me_name) == 0) {
      string fullpathname = iBooker.pwd() + "/" + (*it); 
            
      me = iGetter.get(fullpathname);
            
      if (me) {
    me->Reset();
    return me;
      }
    }
  }
  contents.clear();
  me = iBooker.book1D(me_name.c_str(), me_name.c_str(),40,-0.5,39.5);

  return me;
}

//=============================================================================================================
void SiPixelActionExecutor::bookOccupancyPlots(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool hiRes, bool isbarrel) // Polymorphism
{
  if(Tier0Flag_) return;
  vector<string> subdirs = iGetter.getSubdirs();
  for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++)
    {
      if(isbarrel && (*it).find("Barrel")==string::npos) continue;
      if(!isbarrel && (*it).find("Endcap")==string::npos) continue;
        
      if((*it).find("Module_")!=string::npos) continue;
      if((*it).find("Panel_")!=string::npos) continue;
      if((*it).find("Ladder_")!=string::npos) continue;
      if((*it).find("Blade_")!=string::npos) continue;
      if((*it).find("Layer_")!=string::npos) continue;
      if((*it).find("Disk_")!=string::npos) continue;
      iBooker.cd(*it);
      iGetter.cd(*it);
      bookOccupancyPlots(iBooker, iGetter, hiRes, isbarrel);
      if(!hiRes){
    //occupancyprinting cout<<"booking low res barrel occ plot now!"<<endl;
    OccupancyMap = iBooker.book2D((isbarrel?"barrelOccupancyMap":"endcapOccupancyMap"),"Barrel Digi Occupancy Map (4 pix per bin)",isbarrel?208:130,0.,isbarrel?416.:260.,80,0.,160.);
      }else{
    //occupancyprinting cout<<"booking high res barrel occ plot now!"<<endl;
    OccupancyMap = iBooker.book2D((isbarrel?"barrelOccupancyMap":"endcapOccupancyMap"),"Barrel Digi Occupancy Map (1 pix per bin)",isbarrel?416:260,0.,isbarrel?416.:260.,160,0.,160.);
      }
      OccupancyMap->setAxisTitle("Columns",1);
      OccupancyMap->setAxisTitle("Rows",2);
        
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
        
    }
    
    
    
}
//=============================================================================================================
void SiPixelActionExecutor::bookOccupancyPlots(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool hiRes) {
    
  if(Tier0Flag_) return;
  // Barrel
  iGetter.cd();
  iBooker.cd();
  iGetter.setCurrentFolder("Pixel");
  iBooker.setCurrentFolder("Pixel");
  this->bookOccupancyPlots(iBooker, iGetter, hiRes, true);
    
  // Endcap
  iGetter.cd();
  iBooker.cd();
  iGetter.setCurrentFolder("Pixel");
  iBooker.setCurrentFolder("Pixel");
  this->bookOccupancyPlots(iBooker, iGetter, hiRes, false);

}

void SiPixelActionExecutor::createOccupancy(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter) {
  //std::cout<<"entering SiPixelActionExecutor::createOccupancy..."<<std::endl;
  if(Tier0Flag_) return;
  iBooker.cd();
  iGetter.cd();
  fillOccupancy(iBooker,iGetter, true);
  iBooker.cd();
  iGetter.cd();
  fillOccupancy(iBooker,iGetter, false);
  iBooker.cd();
  iGetter.cd();

  //std::cout<<"leaving SiPixelActionExecutor::createOccupancy..."<<std::endl;
}

//=============================================================================================================

void SiPixelActionExecutor::fillOccupancy(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool isbarrel)
{
  //occupancyprinting cout<<"entering SiPixelActionExecutor::fillOccupancy..."<<std::endl;
  if(Tier0Flag_) return;
  string currDir = iBooker.pwd();
  string dname = currDir.substr(currDir.find_last_of("/")+1);
    
  if(dname.find("Layer_")!=string::npos || dname.find("Disk_")!=string::npos){ 
    vector<string> meVec = iGetter.getMEs();
    for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
      string full_path = currDir + "/" + (*it);
      if(full_path.find("hitmap_siPixelDigis")!=string::npos){ // If we have the hitmap ME
    MonitorElement * me = iGetter.get(full_path);
    if (!me) continue;
    string path = full_path;
    while (path.find_last_of("/") != 5) // Stop before Pixel/
      {
        path = path.substr(0,path.find_last_of("/"));
        //                          cout << "\t" << path << endl;
        OccupancyMap = iGetter.get(path + "/" + (isbarrel?"barrel":"endcap") + "OccupancyMap");
                    
        if(OccupancyMap){ 
        for(int i=1; i!=me->getNbinsX()+1; i++) for(int j=1; j!=me->getNbinsY()+1; j++){
          float previous = OccupancyMap->getBinContent(i,j);
          OccupancyMap->setBinContent(i,j,previous + me->getBinContent(i,j));
        }
          OccupancyMap->getTH2F()->SetEntries(OccupancyMap->getTH2F()->Integral());
        }       
                    
      }
      }
            
    }
  } else {  
    vector<string> subdirs = iGetter.getSubdirs();
    for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
      iGetter.cd(*it);
      iBooker.cd(*it);
      if(*it != "Pixel" && ((isbarrel && (*it).find("Barrel")==string::npos) || (!isbarrel && (*it).find("Endcap")==string::npos))) continue;
      fillOccupancy(iBooker,iGetter, isbarrel);
      iBooker.goUp();
      iGetter.setCurrentFolder(iBooker.pwd());
    }
  }
    
  //occupancyprinting cout<<"leaving SiPixelActionExecutor::fillOccupancy..."<<std::endl;
    
}

//=============================================================================================================

void SiPixelActionExecutor::normaliseAvDigiOcc(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter){
  //occupancyprinting cout<<"entering SiPixelActionExecutor::normaliseAvDigiOcc..."<<std::endl;
  
  iGetter.cd();

  MonitorElement * roccupancyPlot = iGetter.get("Pixel/averageDigiOccupancy");
  
  float totalDigisBPIX = 0.;
  float totalDigisFPIX = 0.;
  for (int i = 1; i !=41; i++){
    if (i < 33) totalDigisBPIX += roccupancyPlot->getBinContent(i);
    else totalDigisFPIX += roccupancyPlot->getBinContent(i);
  } 
  float averageBPIXOcc = totalDigisBPIX/32.;
  float averageFPIXOcc = totalDigisFPIX/8.;
  for (int i = 1; i !=41; i++){
    if (i < 33) roccupancyPlot->setBinContent(i,roccupancyPlot->getBinContent(i)/averageBPIXOcc);
    else roccupancyPlot->setBinContent(i,roccupancyPlot->getBinContent(i)/averageFPIXOcc);
  }
  
  iGetter.setCurrentFolder(iBooker.pwd());
  
}

//=============================================================================================================

void SiPixelActionExecutor::normaliseAvDigiOccVsLumi(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, int lumisec){

  iGetter.cd();

  MonitorElement* avgfedDigiOccvsLumi = iGetter.get("Pixel/avgfedDigiOccvsLumi");

  float totalDigisBPIX = 0.; 
  float totalDigisFPIX = 0.;
  for (int i = 1; i !=41; i++){
    if (i < 33) totalDigisBPIX += avgfedDigiOccvsLumi->getBinContent(lumisec,i);
    else        totalDigisFPIX += avgfedDigiOccvsLumi->getBinContent(lumisec,i);
  }  
  float averageBPIXOcc = totalDigisBPIX/32.;
  float averageFPIXOcc = totalDigisFPIX/8.;
  for (int i = 1; i !=41; i++){
    if (i < 33) avgfedDigiOccvsLumi->setBinContent(lumisec,i,avgfedDigiOccvsLumi->getBinContent(lumisec,i)/averageBPIXOcc);
    else        avgfedDigiOccvsLumi->setBinContent(lumisec,i,avgfedDigiOccvsLumi->getBinContent(lumisec,i)/averageFPIXOcc);
  }

  iGetter.setCurrentFolder(iBooker.pwd());

}

//=============================================================================================================

void SiPixelActionExecutor::bookEfficiency(DQMStore::IBooker & iBooker, bool isUpgrade){
  // Barrel
  iBooker.cd();
  iBooker.setCurrentFolder("Pixel/Barrel");
  if (!isUpgrade) {
  if(Tier0Flag_){
    HitEfficiency_L1 = iBooker.book2D("HitEfficiency_L1","Hit Efficiency in Barrel_Layer1;Module;Ladder",9,-4.5,4.5,21,-10.5,10.5);
    HitEfficiency_L2 = iBooker.book2D("HitEfficiency_L2","Hit Efficiency in Barrel_Layer2;Module;Ladder",9,-4.5,4.5,33,-16.5,16.5);
    HitEfficiency_L3 = iBooker.book2D("HitEfficiency_L3","Hit Efficiency in Barrel_Layer3;Module;Ladder",9,-4.5,4.5,45,-22.5,22.5);
  }else{
    HitEfficiency_L1 = iBooker.book2D("HitEfficiency_L1","Hit Efficiency in Barrel_Layer1;Module;Ladder",9,-4.5,4.5,21,-10.5,10.5);
    HitEfficiency_L2 = iBooker.book2D("HitEfficiency_L2","Hit Efficiency in Barrel_Layer2;Module;Ladder",9,-4.5,4.5,33,-16.5,16.5);
    HitEfficiency_L3 = iBooker.book2D("HitEfficiency_L3","Hit Efficiency in Barrel_Layer3;Module;Ladder",9,-4.5,4.5,45,-22.5,22.5); 
  }
  }//endifNOTUpgrade
  else if (isUpgrade) {
      if(Tier0Flag_){
      HitEfficiency_L1 = iBooker.book2D("HitEfficiency_L1","Hit Efficiency in Barrel_Layer1;z-side;Ladder",2,-1.,1.,12,-6.,6.);
      HitEfficiency_L2 = iBooker.book2D("HitEfficiency_L2","Hit Efficiency in Barrel_Layer2;z-side;Ladder",2,-1.,1.,28,-14.,14.);
      HitEfficiency_L3 = iBooker.book2D("HitEfficiency_L3","Hit Efficiency in Barrel_Layer3;z-side;Ladder",2,-1.,1.,44,-22.,22.);
      HitEfficiency_L4 = iBooker.book2D("HitEfficiency_L4","Hit Efficiency in Barrel_Layer4;z-side;Ladder",2,-1.,1.,64,-32.,32.);
    }else{
      HitEfficiency_L1 = iBooker.book2D("HitEfficiency_L1","Hit Efficiency in Barrel_Layer1;Module;Ladder",8,-4.,4.,12,-6.,6.);
      HitEfficiency_L2 = iBooker.book2D("HitEfficiency_L2","Hit Efficiency in Barrel_Layer2;Module;Ladder",8,-4.,4.,28,-14.,14.);
      HitEfficiency_L3 = iBooker.book2D("HitEfficiency_L3","Hit Efficiency in Barrel_Layer3;Module;Ladder",8,-4.,4.,44,-22.,22.);
      HitEfficiency_L4 = iBooker.book2D("HitEfficiency_L4","Hit Efficiency in Barrel_Layer4;Module;Ladder",8,-4.,4.,64,-32.,32.);
    }
  }//endifUpgrade
  // Endcap
  iBooker.cd();
  iBooker.setCurrentFolder("Pixel/Endcap");
  if (!isUpgrade) {
    if(Tier0Flag_){
      HitEfficiency_Dp1 = iBooker.book2D("HitEfficiency_Dp1","Hit Efficiency in Endcap_Disk_p1;Blade;Panel",26,-13.,13.,2,0.5,2.5);
      HitEfficiency_Dp2 = iBooker.book2D("HitEfficiency_Dp2","Hit Efficiency in Endcap_Disk_p2;Blade;Panel",26,-13.,13.,2,0.5,2.5);
      HitEfficiency_Dm1 = iBooker.book2D("HitEfficiency_Dm1","Hit Efficiency in Endcap_Disk_m1;Blade;Panel",26,-13.,13.,2,0.5,2.5);
      HitEfficiency_Dm2 = iBooker.book2D("HitEfficiency_Dm2","Hit Efficiency in Endcap_Disk_m2;;Blade;Panel",26,-13.,13.,2,0.5,2.5);   
    }else{
      HitEfficiency_Dp1 = iBooker.book2D("HitEfficiency_Dp1","Hit Efficiency in Endcap_Disk_p1;Blades;Modules",24,-12.,12.,7,1.,8.);
      HitEfficiency_Dp2 = iBooker.book2D("HitEfficiency_Dp2","Hit Efficiency in Endcap_Disk_p2;Blades;Modules",24,-12.,12.,7,1.,8.);
      HitEfficiency_Dm1 = iBooker.book2D("HitEfficiency_Dm1","Hit Efficiency in Endcap_Disk_m1;Blades;Modules",24,-12.,12.,7,1.,8.);
      HitEfficiency_Dm2 = iBooker.book2D("HitEfficiency_Dm2","Hit Efficiency in Endcap_Disk_m2;Blades;Modules",24,-12.,12.,7,1.,8.);
    }
  } else if (isUpgrade) {
    if(Tier0Flag_){
      HitEfficiency_Dp1 = iBooker.book2D("HitEfficiency_Dp1","Hit Efficiency in Endcap_Disk_p1;Blades;",28,-17.,11.,1,0.,1.);
      HitEfficiency_Dp2 = iBooker.book2D("HitEfficiency_Dp2","Hit Efficiency in Endcap_Disk_p2;Blades;",28,-17.,11.,1,0.,1.);
      HitEfficiency_Dp3 = iBooker.book2D("HitEfficiency_Dp3","Hit Efficiency in Endcap_Disk_p3;Blades;",28,-17.,11.,1,0.,1.);
      HitEfficiency_Dm1 = iBooker.book2D("HitEfficiency_Dm1","Hit Efficiency in Endcap_Disk_m1;Blades;",28,-17.,11.,1,0.,1.);
      HitEfficiency_Dm2 = iBooker.book2D("HitEfficiency_Dm2","Hit Efficiency in Endcap_Disk_m2;Blades;",28,-17.,11.,1,0.,1.);
      HitEfficiency_Dm3 = iBooker.book2D("HitEfficiency_Dm3","Hit Efficiency in Endcap_Disk_m3;Blades;",28,-17.,11.,1,0.,1.);
    }else{
      HitEfficiency_Dp1 = iBooker.book2D("HitEfficiency_Dp1","Hit Efficiency in Endcap_Disk_p1;Blades;Modules",28,-17.,11.,2,1.,3.);
      HitEfficiency_Dp2 = iBooker.book2D("HitEfficiency_Dp2","Hit Efficiency in Endcap_Disk_p2;Blades;Modules",28,-17.,11.,2,1.,3.);
      HitEfficiency_Dp3 = iBooker.book2D("HitEfficiency_Dp3","Hit Efficiency in Endcap_Disk_p3;Blades;Modules",28,-17.,11.,2,1.,3.);
      HitEfficiency_Dm1 = iBooker.book2D("HitEfficiency_Dm1","Hit Efficiency in Endcap_Disk_m1;Blades;Modules",28,-17.,11.,2,1.,3.);
      HitEfficiency_Dm2 = iBooker.book2D("HitEfficiency_Dm2","Hit Efficiency in Endcap_Disk_m2;Blades;Modules",28,-17.,11.,2,1.,3.);
      HitEfficiency_Dm3 = iBooker.book2D("HitEfficiency_Dm3","Hit Efficiency in Endcap_Disk_m3;Blades;Modules",28,-17.,11.,2,1.,3.);
    }
  }//endif(isUpgrade)
  iBooker.cd();
  iBooker.cd("Pixel/");
  string bins[] = {"Layer1","Layer2","Layer3","Disk1+","Disk2+","Disk1-","Disk2-"};
  HitEfficiencySummary = iBooker.book1D("HitEfficiencySummary","Hit efficiency per sub detector",7,0,7);
  HitEfficiencySummary->setAxisTitle("Sub Detector", 1);
  HitEfficiencySummary->setAxisTitle("Efficiency (%)",2);
  for (int i = 1; i < 8; i++){
    HitEfficiencySummary->setBinLabel(i,bins[i-1]);
  }
}

//=============================================================================================================

void SiPixelActionExecutor::createEfficiency(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool isUpgrade){
  //std::cout<<"entering SiPixelActionExecutor::createEfficiency..."<<std::endl;
  iGetter.cd();
  iBooker.cd();
  fillEfficiency(iBooker, iGetter, true, isUpgrade); // Barrel
  iGetter.cd();
  iBooker.cd();
  fillEfficiency(iBooker, iGetter, false, isUpgrade); // Endcap
  iGetter.cd();
  iBooker.cd();
  //std::cout<<"leaving SiPixelActionExecutor::createEfficiency..."<<std::endl;
}

//=============================================================================================================

int SiPixelActionExecutor::getLadder(const std::string & dname_){
  int biny_= 0;
  string lad=dname_.substr(dname_.find("Ladder_") + 7, 2);
  if(dname_.find(lad)!=string::npos){ biny_ = atoi(lad.c_str());}
  return biny_;
}

//=============================================================================================================

int SiPixelActionExecutor::getBlade(const std::string & dname_){
  int binx_= 0;
  string blad=dname_.substr(dname_.find("Blade_") + 6, 2);
  if(dname_.find(blad)!=string::npos){ binx_ = atoi(blad.c_str());}
  return binx_;
}

//=============================================================================================================

void SiPixelActionExecutor::fillEfficiency(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter, bool isbarrel, bool isUpgrade){
  //cout<<"entering SiPixelActionExecutor::fillEfficiency..."<<std::endl;
  string currDir = iBooker.pwd();
  string dname = currDir.substr(currDir.find_last_of("/")+1);
  //cout<<"currDir= "<<currDir<< " , dname= "<<dname<<std::endl;
  
  if(Tier0Flag_){ // Offline    
    if(isbarrel && dname.find("Ladder_")!=string::npos){ 
      if (!isUpgrade) {
      vector<string> meVec = iGetter.getMEs();
      for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
        string full_path = currDir + "/" + (*it);

    if(full_path.find("missingMod_")!=string::npos){ // If we have missing hits ME
      
      //Get the MEs that contain missing and valid hits
      MonitorElement * missing = iGetter.get(full_path);
      if (!missing) continue;
      string new_path = full_path.replace(full_path.find("missing"),7,"valid");
      MonitorElement * valid = iGetter.get(new_path);
      if (!valid) continue;
      int binx = 0; 
      int biny = 0;
      //get the ladder number
       biny = getLadder(dname);//Current
      if(currDir.find("Shell_mO")!=string::npos || currDir.find("Shell_pO")!=string::npos){ 
        biny=-biny;
      }
      const int nMod = 4; 
      for(int i=1; i<nMod+1;i++){
        float hitEfficiency = -1.0;
        float missingHits=0;
        float validHits=0;
        binx=i; //Module
        if(currDir.find("Shell_m")!=string::npos) binx=-binx;

        missingHits=missing->getBinContent(i);
        validHits=valid->getBinContent(i);
        if(validHits + missingHits > 0.) hitEfficiency = validHits / (validHits + missingHits);

        if(currDir.find("Layer_1")!=string::npos){
          HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
          if(HitEfficiency_L1) HitEfficiency_L1->setBinContent(HitEfficiency_L1->getTH2F()->FindBin(binx, biny),(float)hitEfficiency);
        }
        else if(currDir.find("Layer_2")!=string::npos){
          HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
          if(HitEfficiency_L2) HitEfficiency_L2->setBinContent(HitEfficiency_L2->getTH2F()->FindBin(binx, biny),(float)hitEfficiency);
        }
        else if(currDir.find("Layer_3")!=string::npos){
          HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");      
          if(HitEfficiency_L3) HitEfficiency_L3->setBinContent(HitEfficiency_L3->getTH2F()->FindBin(binx, biny),(float)hitEfficiency);
        }     
        
      }
      
    }
      }
      }//endifNOTUpgradeInBPix
      else if (isUpgrade) {
        vector<string> meVec = iGetter.getMEs();
        for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
          string full_path = currDir + "/" + (*it);
          if(full_path.find("missing_")!=string::npos){ // If we have missing hits ME
        MonitorElement * me = iGetter.get(full_path);
        if (!me) continue;
        float missingHits = me->getEntries();
        string new_path = full_path.replace(full_path.find("missing"),7,"valid");
        me = iGetter.get(new_path);
        if (!me) continue;
        float validHits = me->getEntries();
        float hitEfficiency = -1.;
        if(validHits + missingHits > 0.) hitEfficiency = validHits / (validHits + missingHits);
        int binx = 0; int biny = 0;
        biny = getLadder(dname);
        if(currDir.find("Shell_mO")!=string::npos || currDir.find("Shell_pO")!=string::npos){
          biny=-biny;
        }
        if(currDir.find("Shell_m")!=string::npos){ binx = 1;}else{ binx = 2;} //x-axis: z-side
        if(currDir.find("Layer_1")!=string::npos){
          HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
          if(HitEfficiency_L1) HitEfficiency_L1->setBinContent(binx, biny,(float)hitEfficiency);
        }else if(currDir.find("Layer_2")!=string::npos){
          HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
          if(HitEfficiency_L2) HitEfficiency_L2->setBinContent(binx, biny,(float)hitEfficiency);
        }else if(currDir.find("Layer_3")!=string::npos){
          HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");
          if(HitEfficiency_L3) HitEfficiency_L3->setBinContent(binx, biny,(float)hitEfficiency);
        }else if(currDir.find("Layer_4")!=string::npos){
          HitEfficiency_L4 = iGetter.get("Pixel/Barrel/HitEfficiency_L4");
          if(HitEfficiency_L4) HitEfficiency_L4->setBinContent(binx, biny,(float)hitEfficiency);
        }
          }
        }
      }//endifUpgradeInBPix
    }else if(!isbarrel && dname.find("Blade_")!=string::npos && !isUpgrade){ 
      vector<string> meVec = iGetter.getMEs();
      for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
        string full_path = currDir + "/" + (*it);
        if(full_path.find("missing_")!=string::npos){ // If we have missing hits ME
      MonitorElement * missing = iGetter.get(full_path);
      if (!missing) continue;
      //float missingHits = missing->getEntries();
      string new_path = full_path.replace(full_path.find("missing"),7,"valid");
      MonitorElement * valid = iGetter.get(new_path);
      if (!valid) continue;
      //float validHits = valid->getEntries();
      int binx = 0; int biny = 0;
      binx = getBlade(dname);
      if(currDir.find("HalfCylinder_mI")!=string::npos || currDir.find("HalfCylinder_pI")!=string::npos){ binx = binx + 14;}
      else{  binx = 13 -binx;}
      const int nPanel=2;
      for (int i=1; i<nPanel+1; i++){
        float hitEfficiency = -1.;
        float missingHits=0;
        float validHits=0;
        biny=i;
        missingHits=missing->getBinContent(i);
        validHits=valid->getBinContent(i);
        if(validHits + missingHits > 0.) hitEfficiency = validHits / (validHits + missingHits);
        if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
          HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
          if(HitEfficiency_Dm1) HitEfficiency_Dm1->setBinContent(binx, biny, (float)hitEfficiency);
        }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
          HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
          if(HitEfficiency_Dm2) HitEfficiency_Dm2->setBinContent(binx, biny, (float)hitEfficiency);
        }else if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
          HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
          if(HitEfficiency_Dp1) HitEfficiency_Dp1->setBinContent(binx, biny, (float)hitEfficiency);
        }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
          HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
          if(HitEfficiency_Dp2) HitEfficiency_Dp2->setBinContent(binx, biny, (float)hitEfficiency);
        }
      }//EndOfFor
    }
      } 
    }else if(!isbarrel && dname.find("Blade_")!=string::npos && isUpgrade){ 
      vector<string> meVec = iGetter.getMEs();
      for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
        string full_path = currDir + "/" + (*it);
        if(full_path.find("missing_")!=string::npos){ // If we have missing hits ME
      MonitorElement * me = iGetter.get(full_path);
      if (!me) continue;
      float missingHits = me->getEntries();
      string new_path = full_path.replace(full_path.find("missing"),7,"valid");
      me = iGetter.get(new_path);
      if (!me) continue;
      float validHits = me->getEntries();
      float hitEfficiency = -1.;
      if(validHits + missingHits > 0.) hitEfficiency = validHits / (validHits + missingHits);
      int binx = 0; int biny = 1;
      binx = getBlade(dname);
      if(currDir.find("HalfCylinder_mI")!=string::npos || currDir.find("HalfCylinder_pI")!=string::npos){ binx = binx + 12;}
      else{ 
        if(binx==1) binx = 17;
        else if(binx==2) binx = 16;
        else if(binx==3) binx = 15;
        else if(binx==4) binx = 14;
        else if(binx==5) binx = 13;
        else if(binx==6) binx = 12;
        else if(binx==7) binx = 11;
        else if(binx==8) binx = 10;
        else if(binx==9) binx = 9;
        else if(binx==10) binx = 8;
        else if(binx==11) binx = 7;
        else if(binx==12) binx = 6;
            else if(binx==13) binx = 5;
        else if(binx==14) binx = 4;
        else if(binx==15) binx = 3;
        else if(binx==16) binx = 2;
        else if(binx==17) binx = 1;
      }
      if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
        if(HitEfficiency_Dm1) HitEfficiency_Dm1->setBinContent(binx, biny, (float)hitEfficiency);
      }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
        if(HitEfficiency_Dm2) HitEfficiency_Dm2->setBinContent(binx, biny, (float)hitEfficiency);
      }else if(currDir.find("Disk_3")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm3");
        if(HitEfficiency_Dm3) HitEfficiency_Dm3->setBinContent(binx, biny, (float)hitEfficiency);
      }else if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
        if(HitEfficiency_Dp1) HitEfficiency_Dp1->setBinContent(binx, biny, (float)hitEfficiency);
      }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
        if(HitEfficiency_Dp2) HitEfficiency_Dp2->setBinContent(binx, biny, (float)hitEfficiency);
          }else if(currDir.find("Disk_3")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp3");
        if(HitEfficiency_Dp3) HitEfficiency_Dp3->setBinContent(binx, biny, (float)hitEfficiency);
          }
      //std::cout<<"EFFI: "<<currDir<<" , x: "<<binx<<" , y: "<<biny<<std::endl;
    }
      } 
    }else{  
      vector<string> subdirs = iGetter.getSubdirs();
      for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
        iBooker.cd(*it);
    iGetter.cd(*it);
        if(*it != "Pixel" && ((isbarrel && (*it).find("Barrel")==string::npos) || (!isbarrel && (*it).find("Endcap")==string::npos))) continue;
        fillEfficiency(iBooker, iGetter, isbarrel, isUpgrade);
        iBooker.goUp();
    iGetter.setCurrentFolder(iBooker.pwd());
      }
    }
  }else{ // Online
    if(dname.find("Module_")!=string::npos){ 
      vector<string> meVec = iGetter.getMEs();
      for (vector<string>::const_iterator it = meVec.begin(); it != meVec.end(); it++) {
        string full_path = currDir + "/" + (*it);
        if(full_path.find("missing_")!=string::npos){ // If we have missing hits ME
      MonitorElement * me = iGetter.get(full_path);
      if (!me) continue;
      float missingHits = me->getEntries();
      string new_path = full_path.replace(full_path.find("missing"),7,"valid");
      me = iGetter.get(new_path);
      if (!me) continue;
      float validHits = me->getEntries();
      float hitEfficiency = -1.;
      if(validHits + missingHits > 0.) hitEfficiency = validHits / (validHits + missingHits);
      int binx = 0; int biny = 0;
      if(isbarrel){
        if(currDir.find("Shell_m")!=string::npos){
          if(currDir.find("Module_4")!=string::npos){ binx = 1;}else if(currDir.find("Module_3")!=string::npos){ binx = 2;}
          if(currDir.find("Module_2")!=string::npos){ binx = 3;}else if(currDir.find("Module_1")!=string::npos){ binx = 4;}
        }else if(currDir.find("Shell_p")!=string::npos){
          if(currDir.find("Module_1")!=string::npos){ binx = 5;}else if(currDir.find("Module_2")!=string::npos){ binx = 6;}
          if(currDir.find("Module_3")!=string::npos){ binx = 7;}else if(currDir.find("Module_4")!=string::npos){ binx = 8;}
        }
        if (!isUpgrade) {
        if(currDir.find("01")!=string::npos){ biny = 1;}else if(currDir.find("02")!=string::npos){ biny = 2;}
        else if(currDir.find("03")!=string::npos){ biny = 3;}else if(currDir.find("04")!=string::npos){ biny = 4;}
        else if(currDir.find("05")!=string::npos){ biny = 5;}else if(currDir.find("06")!=string::npos){ biny = 6;}
        else if(currDir.find("07")!=string::npos){ biny = 7;}else if(currDir.find("08")!=string::npos){ biny = 8;}
        else if(currDir.find("09")!=string::npos){ biny = 9;}else if(currDir.find("10")!=string::npos){ biny = 10;}
        else if(currDir.find("11")!=string::npos){ biny = 11;}else if(currDir.find("12")!=string::npos){ biny = 12;}
        else if(currDir.find("13")!=string::npos){ biny = 13;}else if(currDir.find("14")!=string::npos){ biny = 14;}
        else if(currDir.find("15")!=string::npos){ biny = 15;}else if(currDir.find("16")!=string::npos){ biny = 16;}
        else if(currDir.find("17")!=string::npos){ biny = 17;}else if(currDir.find("18")!=string::npos){ biny = 18;}
        else if(currDir.find("19")!=string::npos){ biny = 19;}else if(currDir.find("20")!=string::npos){ biny = 20;}
        else if(currDir.find("21")!=string::npos){ biny = 21;}else if(currDir.find("22")!=string::npos){ biny = 22;}
        if(currDir.find("Shell_mO")!=string::npos || currDir.find("Shell_pO")!=string::npos){
          if(currDir.find("Layer_1")!=string::npos){ biny = biny + 10;}
          else if(currDir.find("Layer_2")!=string::npos){ biny = biny + 16;}
          else if(currDir.find("Layer_3")!=string::npos){ biny = biny + 22;}
        }
        }
        else if (isUpgrade) {
          if(currDir.find("01")!=string::npos){ biny = 1;}else if(currDir.find("02")!=string::npos){ biny = 2;}
          else if(currDir.find("03")!=string::npos){ biny = 3;}else if(currDir.find("04")!=string::npos){ biny = 4;}
          else if(currDir.find("05")!=string::npos){ biny = 5;}else if(currDir.find("06")!=string::npos){ biny = 6;}
          else if(currDir.find("07")!=string::npos){ biny = 7;}else if(currDir.find("08")!=string::npos){ biny = 8;}
          else if(currDir.find("09")!=string::npos){ biny = 9;}else if(currDir.find("10")!=string::npos){ biny = 10;}
          else if(currDir.find("11")!=string::npos){ biny = 11;}else if(currDir.find("12")!=string::npos){ biny = 12;}
          else if(currDir.find("13")!=string::npos){ biny = 13;}else if(currDir.find("14")!=string::npos){ biny = 14;}
          else if(currDir.find("15")!=string::npos){ biny = 15;}else if(currDir.find("16")!=string::npos){ biny = 16;}
          else if(currDir.find("17")!=string::npos){ biny = 17;}else if(currDir.find("18")!=string::npos){ biny = 18;}
          else if(currDir.find("19")!=string::npos){ biny = 19;}else if(currDir.find("20")!=string::npos){ biny = 20;}
          else if(currDir.find("21")!=string::npos){ biny = 21;}else if(currDir.find("22")!=string::npos){ biny = 22;}
          else if(currDir.find("23")!=string::npos){ biny = 23;}else if(currDir.find("24")!=string::npos){ biny = 24;}
          else if(currDir.find("25")!=string::npos){ biny = 25;}else if(currDir.find("25")!=string::npos){ biny = 25;}
          else if(currDir.find("26")!=string::npos){ biny = 26;}else if(currDir.find("27")!=string::npos){ biny = 27;}
          else if(currDir.find("28")!=string::npos){ biny = 28;}else if(currDir.find("29")!=string::npos){ biny = 29;}
          else if(currDir.find("30")!=string::npos){ biny = 30;}else if(currDir.find("31")!=string::npos){ biny = 31;}
          else if(currDir.find("32")!=string::npos){ biny = 32;}
          if(currDir.find("Shell_mO")!=string::npos || currDir.find("Shell_pO")!=string::npos){
            if(currDir.find("Layer_1")!=string::npos){ biny = biny + 6;}
            else if(currDir.find("Layer_2")!=string::npos){ biny = biny + 14;}
            else if(currDir.find("Layer_3")!=string::npos){ biny = biny + 22;}
            else if(currDir.find("Layer_4")!=string::npos){ biny = biny + 32;}
          }
        }
      }else{ //endcap
        if (!isUpgrade) {
        if(currDir.find("01")!=string::npos){ binx = 1;}else if(currDir.find("02")!=string::npos){ binx = 2;}
        else if(currDir.find("03")!=string::npos){ binx = 3;}else if(currDir.find("04")!=string::npos){ binx = 4;}
        else if(currDir.find("05")!=string::npos){ binx = 5;}else if(currDir.find("06")!=string::npos){ binx = 6;}
        else if(currDir.find("07")!=string::npos){ binx = 7;}else if(currDir.find("08")!=string::npos){ binx = 8;}
        else if(currDir.find("09")!=string::npos){ binx = 9;}else if(currDir.find("10")!=string::npos){ binx = 10;}
        else if(currDir.find("11")!=string::npos){ binx = 11;}else if(currDir.find("12")!=string::npos){ binx = 12;}
        if(currDir.find("HalfCylinder_mO")!=string::npos || currDir.find("HalfCylinder_pO")!=string::npos){ binx = binx + 12;}
        if(currDir.find("Panel_1/Module_1")!=string::npos){ biny = 1;}else if(currDir.find("Panel_2/Module_1")!=string::npos){ biny = 2;}
        else if(currDir.find("Panel_1/Module_2")!=string::npos){ biny = 3;}else if(currDir.find("Panel_2/Module_2")!=string::npos){ biny = 4;}
        else if(currDir.find("Panel_1/Module_3")!=string::npos){ biny = 5;}else if(currDir.find("Panel_2/Module_3")!=string::npos){ biny = 6;}
        else if(currDir.find("Panel_1/Module_4")!=string::npos){ biny = 7;}
        } else if (isUpgrade) {
          if(currDir.find("01")!=string::npos){ binx = 1;}else if(currDir.find("02")!=string::npos){ binx = 2;}
          else if(currDir.find("03")!=string::npos){ binx = 3;}else if(currDir.find("04")!=string::npos){ binx = 4;}
          else if(currDir.find("05")!=string::npos){ binx = 5;}else if(currDir.find("06")!=string::npos){ binx = 6;}
          else if(currDir.find("07")!=string::npos){ binx = 7;}else if(currDir.find("08")!=string::npos){ binx = 8;}
          else if(currDir.find("09")!=string::npos){ binx = 9;}else if(currDir.find("10")!=string::npos){ binx = 10;}
          else if(currDir.find("11")!=string::npos){ binx = 11;}else if(currDir.find("12")!=string::npos){ binx = 12;}
          else if(currDir.find("13")!=string::npos){ binx = 13;}else if(currDir.find("14")!=string::npos){ binx = 14;}
          else if(currDir.find("15")!=string::npos){ binx = 15;}else if(currDir.find("16")!=string::npos){ binx = 16;}
          else if(currDir.find("17")!=string::npos){ binx = 17;}
          if(currDir.find("HalfCylinder_mO")!=string::npos || currDir.find("HalfCylinder_pO")!=string::npos){ binx = binx + 17;}
          if(currDir.find("Panel_1/Module_1")!=string::npos){ biny = 1;}else if(currDir.find("Panel_2/Module_1")!=string::npos){ biny = 2;}
        }//endif(isUpgrade)
      }
      
      if(currDir.find("Layer_1")!=string::npos){
        HitEfficiency_L1 = iGetter.get("Pixel/Barrel/HitEfficiency_L1");
        if(HitEfficiency_L1) HitEfficiency_L1->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Layer_2")!=string::npos){
        HitEfficiency_L2 = iGetter.get("Pixel/Barrel/HitEfficiency_L2");
        if(HitEfficiency_L2) HitEfficiency_L2->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Layer_3")!=string::npos){
        HitEfficiency_L3 = iGetter.get("Pixel/Barrel/HitEfficiency_L3");
        if(HitEfficiency_L3) HitEfficiency_L3->setBinContent(binx, biny,(float)hitEfficiency);
      }else if( isUpgrade && (currDir.find("Layer_4")!=string::npos) ){
        HitEfficiency_L4 = iGetter.get("Pixel/Barrel/HitEfficiency_L4");
        if(HitEfficiency_L4) HitEfficiency_L4->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm1");
        if(HitEfficiency_Dm1) HitEfficiency_Dm1->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm2");
        if(HitEfficiency_Dm2) HitEfficiency_Dm2->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_3")!=string::npos && currDir.find("HalfCylinder_m")!=string::npos){
        HitEfficiency_Dm3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dm3");
        if(HitEfficiency_Dm3) HitEfficiency_Dm3->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_1")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp1 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp1");
        if(HitEfficiency_Dp1) HitEfficiency_Dp1->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_2")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp2 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp2");
        if(HitEfficiency_Dp2) HitEfficiency_Dp2->setBinContent(binx, biny,(float)hitEfficiency);
      }else if(currDir.find("Disk_3")!=string::npos && currDir.find("HalfCylinder_p")!=string::npos){
        HitEfficiency_Dp3 = iGetter.get("Pixel/Endcap/HitEfficiency_Dp3");
        if(HitEfficiency_Dp3) HitEfficiency_Dp3->setBinContent(binx, biny,(float)hitEfficiency);
          }
        }
      }
    }else{  
      //cout<<"finding subdirs now"<<std::endl;
      vector<string> subdirs = iGetter.getSubdirs();
      for (vector<string>::const_iterator it = subdirs.begin(); it != subdirs.end(); it++) {
        iBooker.cd(*it);
    iGetter.cd(*it);
        if(*it != "Pixel" && ((isbarrel && (*it).find("Barrel")==string::npos) || (!isbarrel && (*it).find("Endcap")==string::npos))) continue;
        fillEfficiency(iBooker, iGetter, isbarrel, isUpgrade);
        iBooker.goUp();
    iGetter.setCurrentFolder(iBooker.pwd());
      }
    }
  } // end online/offline
  //cout<<"leaving SiPixelActionExecutor::fillEfficiency..."<<std::endl;
}

//=============================================================================================================

void SiPixelActionExecutor::fillEfficiencySummary(DQMStore::IBooker & iBooker, DQMStore::IGetter & iGetter){
  //cout<<"entering SiPixelActionExecutor::fillEfficiencySummary..."<<std::endl;
  //First we get the summary plot"
  if (!Tier0Flag_) return;
  HitEfficiencySummary = iGetter.get("Pixel/HitEfficiencySummary");
  //Now we will loop over the hit efficiency plots and fill it"
  string hitEfficiencyPostfix[] = {"L1","L2","L3","Dp1","Dp2","Dm1","Dm2"};
  std::vector<std::vector<float> > ignoreXBins = {{-4,2}, {4,4,-1,-3,3,-4,-3,-2,-1,-4,-3,-2,-1,1,-4},    {1,-4,1},{},{},{},{}};
  std::vector<std::vector<float> > ignoreYBins = {{-9,-3},{1,16,1,-13,-13,-5,-5,-5,-5,-6,-6,-6,-6,-8,-8},{3,14,6},{},{},{},{}};
  
  for (int i = 0; i < 7; i++){
    string subdetName = "Endcap/";
    if (i < 3) subdetName = "Barrel/";
    char meName [50];
    sprintf(meName,"Pixel/%sHitEfficiency_%s",subdetName.c_str(),hitEfficiencyPostfix[i].c_str());
    MonitorElement * tempHitEffMap = iGetter.get(meName);
    float totalEff = 0.;
    int totalBins = 0;
    TH1 * hitEffMap = tempHitEffMap->getTH1();
    for (int xBin = 1; xBin < tempHitEffMap->getNbinsX() + 1; xBin++){
      if (fabs(hitEffMap->GetXaxis()->GetBinCenter(xBin)) < 1.) continue;
      for (int yBin = 1; yBin < tempHitEffMap->getNbinsY() + 1; yBin++){
    if (fabs(hitEffMap->GetYaxis()->GetBinCenter(yBin)) < 0.5) continue;
    bool ignoreBin = false;
    for (unsigned int j = 0; j < ignoreXBins[i].size(); j++){
      if (hitEffMap->GetXaxis()->GetBinCenter(xBin) == ignoreXBins[i][j] && hitEffMap->GetYaxis()->GetBinCenter(yBin) == ignoreYBins[i][j]){
        ignoreBin = true;
        break;
      }
    }
    if (ignoreBin) continue;
    if (!(tempHitEffMap->getBinContent(xBin,yBin) < 0.)) totalEff += tempHitEffMap->getBinContent(xBin,yBin);
    totalBins++;
      }
    }
    float overalEff = 0.;
    //    std::cout << i << " " << totalEff << " " << totalBins << std::endl;
    if (totalBins > 0) overalEff = totalEff/float(totalBins);
    HitEfficiencySummary->setBinContent(i+1,overalEff);
  }
  
}