EgHLTOfflineClient.cc

Go to the documentation of this file.

#include "DQMOffline/Trigger/interface/EgHLTOfflineClient.h"
 
#include "FWCore/Framework/interface/Run.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
 
#include "FWCore/MessageLogger/interface/MessageLogger.h"
 
#include "DQMServices/Core/interface/DQMStore.h"
 
#include "DQMOffline/Trigger/interface/EgHLTTrigTools.h"
 
#include "HLTrigger/HLTcore/interface/HLTConfigProvider.h"
 
#include <boost/algorithm/string.hpp>
 
#include "TGraphAsymmErrors.h"
 
EgHLTOfflineClient::EgHLTOfflineClient(const edm::ParameterSet& iConfig) : isSetup_(false) {
  eleHLTFilterNames_ = iConfig.getParameter<std::vector<std::string> >("eleHLTFilterNames");
  eleHLTFilterNames2Leg_ = iConfig.getParameter<std::vector<std::string> >("eleHLTFilterNames2Leg");
  eleTightLooseTrigNames_ = iConfig.getParameter<std::vector<std::string> >("eleTightLooseTrigNames");
  phoHLTFilterNames_ = iConfig.getParameter<std::vector<std::string> >("phoHLTFilterNames");
  phoHLTFilterNames2Leg_ = iConfig.getParameter<std::vector<std::string> >("phoHLTFilterNames2Leg");
  phoTightLooseTrigNames_ = iConfig.getParameter<std::vector<std::string> >("phoTightLooseTrigNames");
 
  eleN1EffVars_ = iConfig.getParameter<std::vector<std::string> >("eleN1EffVars");
  eleSingleEffVars_ = iConfig.getParameter<std::vector<std::string> >("eleSingleEffVars");
  eleEffTags_ = iConfig.getParameter<std::vector<std::string> >("eleEffTags");
  eleTrigTPEffVsVars_ = iConfig.getParameter<std::vector<std::string> >("eleTrigTPEffVsVars");
  eleLooseTightTrigEffVsVars_ = iConfig.getParameter<std::vector<std::string> >("eleLooseTightTrigEffVsVars");
  eleHLTvOfflineVars_ = iConfig.getParameter<std::vector<std::string> >("eleHLTvOfflineVars");
 
  phoN1EffVars_ = iConfig.getParameter<std::vector<std::string> >("phoN1EffVars");
  phoSingleEffVars_ = iConfig.getParameter<std::vector<std::string> >("phoSingleEffVars");
  phoEffTags_ = iConfig.getParameter<std::vector<std::string> >("phoEffTags");
  phoTrigTPEffVsVars_ = iConfig.getParameter<std::vector<std::string> >("phoTrigTPEffVsVars");
  phoLooseTightTrigEffVsVars_ = iConfig.getParameter<std::vector<std::string> >("phoLooseTightTrigEffVsVars");
  phoHLTvOfflineVars_ = iConfig.getParameter<std::vector<std::string> >("phoHLTvOfflineVars");
 
  runClientEndLumiBlock_ = iConfig.getParameter<bool>("runClientEndLumiBlock");
  runClientEndRun_ = iConfig.getParameter<bool>("runClientEndRun");
  runClientEndJob_ = iConfig.getParameter<bool>("runClientEndJob");
 
  dirName_ = iConfig.getParameter<std::string>("DQMDirName");
 
  filterInactiveTriggers_ = iConfig.getParameter<bool>("filterInactiveTriggers");
  hltTag_ = iConfig.getParameter<std::string>("hltTag");
}
 
EgHLTOfflineClient::~EgHLTOfflineClient() = default;
 
void EgHLTOfflineClient::dqmEndJob(DQMStore::IBooker& ibooker_, DQMStore::IGetter& igetter_) {
  if (runClientEndJob_)
    runClient_(ibooker_, igetter_);
}
 
void EgHLTOfflineClient::beginRun(const edm::Run& run, const edm::EventSetup& c) {
  if (!isSetup_) {
    if (filterInactiveTriggers_) {
      HLTConfigProvider hltConfig;
      bool changed = false;
      hltConfig.init(run, c, hltTag_, changed);
      std::vector<std::string> activeFilters;
      std::vector<std::string> activeEleFilters;
      std::vector<std::string> activeEle2LegFilters;
      std::vector<std::string> activePhoFilters;
      std::vector<std::string> activePho2LegFilters;
      egHLT::trigTools::getActiveFilters(
          hltConfig, activeFilters, activeEleFilters, activeEle2LegFilters, activePhoFilters, activePho2LegFilters);
 
      egHLT::trigTools::filterInactiveTriggers(eleHLTFilterNames_, activeEleFilters);
      egHLT::trigTools::filterInactiveTriggers(eleHLTFilterNames2Leg_, activeEle2LegFilters);
      egHLT::trigTools::filterInactiveTriggers(phoHLTFilterNames_, activePhoFilters);
      egHLT::trigTools::filterInactiveTightLooseTriggers(eleTightLooseTrigNames_, activeEleFilters);
      egHLT::trigTools::filterInactiveTightLooseTriggers(phoTightLooseTrigNames_, activePhoFilters);
    }
    isSetup_ = true;
  }
}
 
void EgHLTOfflineClient::dqmEndLuminosityBlock(DQMStore::IBooker& ibooker_,
                                               DQMStore::IGetter& igetter_,
                                               edm::LuminosityBlock const& iLumi,
                                               edm::EventSetup const& iSetup) {
  if (runClientEndLumiBlock_)
    runClient_(ibooker_, igetter_);
}
 
void EgHLTOfflineClient::runClient_(DQMStore::IBooker& ibooker, DQMStore::IGetter& igetter) {
  ibooker.setCurrentFolder(dirName_ + "/Client_Histos");
 
  std::vector<std::string> regions;
  regions.emplace_back("eb");
  regions.emplace_back("ee");
 
  for (auto const& eleHLTFilterName : eleHLTFilterNames_) {
    //std::cout<<"FilterName: "<<eleHLTFilterNames_[filterNr]<<std::endl;
    for (auto const& region : regions) {
      for (auto const& eleEffTag : eleEffTags_) {
        //----Morse----
        ibooker.setCurrentFolder(dirName_ + "/Client_Histos/" + eleHLTFilterName);
        //--------------
        createN1EffHists(
            eleHLTFilterName, eleHLTFilterName + "_gsfEle_" + eleEffTag, region, eleN1EffVars_, ibooker, igetter);
 
        createSingleEffHists(
            eleHLTFilterName, eleHLTFilterName + "_gsfEle_" + eleEffTag, region, eleSingleEffVars_, ibooker, igetter);
 
        createTrigTagProbeEffHistsNewAlgo(eleHLTFilterName, region, eleTrigTPEffVsVars_, "gsfEle", ibooker, igetter);
 
        createHLTvsOfflineHists(
            eleHLTFilterName, eleHLTFilterName + "_gsfEle_passFilter", region, eleHLTvOfflineVars_, ibooker, igetter);
      }
    }
  }
  for (auto& filterNr : eleHLTFilterNames2Leg_) {
    for (auto const& region : regions) {
      for (size_t effNr = 0; effNr < eleEffTags_.size(); effNr++) {
        std::string trigNameLeg1 = filterNr.substr(0, filterNr.find("::"));
 
        std::string trigNameLeg2 = filterNr.substr(filterNr.find("::") + 2);
 
        ibooker.setCurrentFolder(dirName_ + "/Client_Histos/" + trigNameLeg2);
        createTrigTagProbeEffHists2Leg(
            trigNameLeg1, trigNameLeg2, region, eleTrigTPEffVsVars_, "gsfEle", ibooker, igetter);
      }
    }
  }
 
  for (auto const& phoHLTFilterName : phoHLTFilterNames_) {
    for (auto const& region : regions) {
      for (auto const& phoEffTag : phoEffTags_) {
        //----Morse----
        ibooker.setCurrentFolder(dirName_ + "/Client_Histos/" + phoHLTFilterName);
        createN1EffHists(
            phoHLTFilterName, phoHLTFilterName + "_pho_" + phoEffTag, region, phoN1EffVars_, ibooker, igetter);
 
        createSingleEffHists(
            phoHLTFilterName, phoHLTFilterName + "_pho_" + phoEffTag, region, phoSingleEffVars_, ibooker, igetter);
 
        createTrigTagProbeEffHistsNewAlgo(phoHLTFilterName, region, phoTrigTPEffVsVars_, "pho", ibooker, igetter);
 
        createHLTvsOfflineHists(
            phoHLTFilterName, phoHLTFilterName + "_pho_passFilter", region, phoHLTvOfflineVars_, ibooker, igetter);
 
        //--------------
      }
    }
  }
 
  for (auto const& region : regions) {
    createLooseTightTrigEff(eleTightLooseTrigNames_, region, eleLooseTightTrigEffVsVars_, "gsfEle", ibooker, igetter);
 
    createLooseTightTrigEff(
        eleTightLooseTrigNames_, region, eleLooseTightTrigEffVsVars_, "gsfEle_trigCuts", ibooker, igetter);
 
    createLooseTightTrigEff(phoTightLooseTrigNames_, region, phoLooseTightTrigEffVsVars_, "pho", ibooker, igetter);
 
    createLooseTightTrigEff(
        phoTightLooseTrigNames_, region, phoLooseTightTrigEffVsVars_, "pho_trigCuts", ibooker, igetter);
  }
  //----Morse-----
  ibooker.setCurrentFolder(dirName_);
  //----------
}
 
void EgHLTOfflineClient::createHLTvsOfflineHists(const std::string& filterName,
                                                 const std::string& baseName,
                                                 const std::string& region,
                                                 const std::vector<std::string>& varNames,
                                                 DQMStore::IBooker& ibooker,
                                                 DQMStore::IGetter& igetter) {
  //need to do Energy manually to get SC Energy
  /*
  MonitorElement* numer = dbe_->get(dirName_+"/Source_Histos/"+filterName+"/"+baseName+"_HLTenergy"+"_"+region);
  MonitorElement* denom = dbe_->get(dirName_+"/Source_Histos/"+filterName+"/"+baseName+"_energy"+"_"+region);
 
  if(numer!=NULL && denom!=NULL){
    std::string effHistName(baseName+"_HLToverOfflineSC_energy_"+region);//std::cout<<"hltVSoffline:  "<<effHistName<<std::endl;
    std::string effHistTitle(effHistName);
    if(region=="eb" || region=="ee"){
      if(region=="eb") effHistTitle = "Barrel "+baseName+" HLToverOfflineSC Energy";
      if(region=="ee") effHistTitle = "Endcap "+baseName+" HLToverOfflineSC Energy";
      FillHLTvsOfflineHist(filterName,effHistName,effHistTitle,numer,denom);    
    }
  }//end Et
  */
 
  //now eta, phi automatically
  for (auto const& varName : varNames) {
    MonitorElement* numer =
        igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_HLT" + varName + "_" + region);
    MonitorElement* denom =
        igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_" + varName + "_" + region);
    if (numer != nullptr && denom != nullptr) {
      std::string effHistName(baseName + "_HLToverOffline_" + varName + "_" +
                              region);  //std::cout<<"hltVSoffline:  "<<effHistName<<std::endl;
      std::string effHistTitle(effHistName);
      if (region == "eb" || region == "ee") {
        if (region == "eb")
          effHistTitle = "Barrel " + baseName + " HLToverOffline " + varName;
        if (region == "ee")
          effHistTitle = "Endcap " + baseName + " HLToverOffline " + varName;
        FillHLTvsOfflineHist(filterName, effHistName, effHistTitle, numer, denom, ibooker, igetter);
      }
    }
  }  //end loop over varNames
}
 
EgHLTOfflineClient::MonitorElement* EgHLTOfflineClient::FillHLTvsOfflineHist(const std::string& filter,
                                                                             const std::string& name,
                                                                             const std::string& title,
                                                                             const MonitorElement* numer,
                                                                             const MonitorElement* denom,
                                                                             DQMStore::IBooker& ibooker,
                                                                             DQMStore::IGetter& igetter) {
  TH1F* num = numer->getTH1F();
  if (num->GetSumw2N() == 0)
    num->Sumw2();
 
  TH1F* den = denom->getTH1F();
  if (den->GetSumw2N() == 0)
    den->Sumw2();
 
  auto* h_eff = (TH1F*)num->Clone(name.c_str());
  h_eff->Divide(num, den, 1, 1, "B");
  h_eff->SetTitle(title.c_str());
  MonitorElement* eff = igetter.get(dirName_ + "/Client_Histos/" + filter + "/" + name);
  if (eff == nullptr) {
    eff = ibooker.book1D(name, h_eff);
  } else {  //I was having problems with collating the histograms, hence why I'm just resetting the histogram value
    *eff->getTH1F() = *h_eff;
    delete h_eff;
  }
  return eff;
}
 
void EgHLTOfflineClient::createN1EffHists(const std::string& filterName,
                                          const std::string& baseName,
                                          const std::string& region,
                                          const std::vector<std::string>& varNames,
                                          DQMStore::IBooker& ibooker,
                                          DQMStore::IGetter& igetter) {
  MonitorElement* numer =
      igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_allCuts_" + region);
 
  for (auto const& varName : varNames) {
    MonitorElement* denom =
        igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_n1_" + varName + "_" + region);
    if (numer != nullptr && denom != nullptr) {
      std::string effHistName(baseName + "_n1Eff_" + varName + "_" +
                              region);  //std::cout<<"N1:  "<<effHistName<<std::endl;
      //std::cout<<region<<"  ";
      //----Morse-----------
      std::string effHistTitle(effHistName);  //std::cout<<effHistTitle<<std::endl;
      if (region == "eb" || region == "ee") {
        if (region == "eb")
          effHistTitle = "Barrel " + baseName + " N1eff " + varName;
        if (region == "ee")
          effHistTitle = "Endcap " + baseName + " N1eff " + varName;
      }  //std::cout<<effHistTitle<<std::endl;
      makeEffMonElemFromPassAndAll(filterName, effHistName, effHistTitle, numer, denom, ibooker, igetter);
      //---------------------
    }
  }  //end loop over varNames
}
 
void EgHLTOfflineClient::createSingleEffHists(const std::string& filterName,
                                              const std::string& baseName,
                                              const std::string& region,
                                              const std::vector<std::string>& varNames,
                                              DQMStore::IBooker& ibooker,
                                              DQMStore::IGetter& igetter) {
  MonitorElement* denom = igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_noCuts_" + region);
 
  for (auto const& varName : varNames) {
    MonitorElement* numer =
        igetter.get(dirName_ + "/Source_Histos/" + filterName + "/" + baseName + "_single_" + varName + "_" + region);
    if (numer != nullptr && denom != nullptr) {
      std::string effHistName(baseName + "_singleEff_" + varName + "_" +
                              region);  //std::cout<<"Si:  "<<effHistName<<std::endl;
      //----Morse-----------
      std::string effHistTitle(effHistName);  //std::cout<<effHistTitle<<std::endl;
      if (region == "eb" || region == "ee") {
        if (region == "eb")
          effHistTitle = "Barrel " + baseName + " SingleEff " + varName;
        if (region == "ee")
          effHistTitle = "Endcap " + baseName + " SingleEff " + varName;
      }  //std::cout<<effHistTitle<<std::endl;
      makeEffMonElemFromPassAndAll(filterName, effHistName, effHistTitle, numer, denom, ibooker, igetter);
      //--------------------
    }
  }  //end loop over varNames
}
 
void EgHLTOfflineClient::createTrigTagProbeEffHists(const std::string& filterName,
                                                    const std::string& region,
                                                    const std::vector<std::string>& vsVarNames,
                                                    const std::string& objName,
                                                    DQMStore::IBooker& ibooker,
                                                    DQMStore::IGetter& igetter) {
  for (auto const& vsVarName : vsVarNames) {
    std::string allName(dirName_ + "/Source_Histos/" + filterName + "/" + filterName + "_trigTagProbe_" + objName +
                        "_all_" + vsVarName + "_" + region);
    MonitorElement* all = igetter.get(allName);
    if (all == nullptr) {
      continue;
    }
    std::string passName(dirName_ + "/Source_Histos/" + filterName + "/" + filterName + "_trigTagProbe_" + objName +
                         "_pass_" + vsVarName + "_" + region);
    MonitorElement* pass = igetter.get(passName);
    if (pass == nullptr) {
      continue;
    }
    //----Morse-----
    std::string effHistTitle(filterName + "_trigTagProbeEff_" + objName + "_vs_" + vsVarName + "_" +
                             region);  //std::cout<<effHistTitle<<std::endl;
    if (region == "eb" || region == "ee") {
      if (region == "eb")
        effHistTitle = "Barrel " + filterName + "_" + objName + " TrigTagProbeEff vs " + vsVarName;
      if (region == "ee")
        effHistTitle = "Endcap " + filterName + "_" + objName + " TrigTagProbeEff vs " + vsVarName;
    }
    //------------
    makeEffMonElemFromPassAndAll(filterName,
                                 filterName + "_trigTagProbeEff_" + objName + "_vs_" + vsVarName + "_" + region,
                                 effHistTitle,
                                 pass,
                                 all,
                                 ibooker,
                                 igetter);
 
  }  //end loop over vsVarNames
}
 
void EgHLTOfflineClient::createTrigTagProbeEffHistsNewAlgo(const std::string& filterName,
                                                           const std::string& region,
                                                           const std::vector<std::string>& vsVarNames,
                                                           const std::string& objName,
                                                           DQMStore::IBooker& ibooker,
                                                           DQMStore::IGetter& igetter) {
  for (auto const& vsVarName : vsVarNames) {
    /* 
       std::string allName(dirName_+"/Source_Histos/"+filterName+"/"+filterName+"_trigTagProbe_"+objName+"_all_"+vsVarNames[varNr]+"_"+region);
       MonitorElement* all = dbe_->get(allName); 
       if(all==NULL){
       //edm::LogInfo("EgHLTOfflineClient") <<" couldnt get hist "<<allName;
       continue;
       }*/
    std::string passName(dirName_ + "/Source_Histos/" + filterName + "/" + filterName + "_trigTagProbe_" + objName +
                         "_passNotTag_" + vsVarName + "_" + region);
    MonitorElement* passNotTag = igetter.get(passName);
    if (passNotTag == nullptr) {
      //edm::LogInfo("EgHLTOfflineClient") <<" couldnt get hist "<<passName;
      continue;
    }
    std::string passTagTagName(dirName_ + "/Source_Histos/" + filterName + "/" + filterName + "_trigTagProbe_" +
                               objName + "_passTagTag_" + vsVarName + "_" + region);
    MonitorElement* passTagTag = igetter.get(passTagTagName);
    if (passTagTag == nullptr) {
      //edm::LogInfo("EgHLTOfflineClient") <<" couldnt get hist "<<passTagTagName;
      continue;
    }
    std::string failName(dirName_ + "/Source_Histos/" + filterName + "/" + filterName + "_trigTagProbe_" + objName +
                         "_fail_" + vsVarName + "_" + region);
    MonitorElement* fail = igetter.get(failName);
    if (fail == nullptr) {
      //edm::LogInfo("EgHLTOfflineClient") <<" couldnt get hist "<<failName;
      continue;
    }
    //----Morse-----
    std::string effHistTitle(filterName + "_trigTagProbeEff_" + objName + "_vs_" + vsVarName + "_" +
                             region);  //std::cout<<effHistTitle<<std::endl;
    if (region == "eb" || region == "ee") {
      if (region == "eb")
        effHistTitle = "Barrel " + filterName + "_" + objName + " TrigTagProbeEff vs " + vsVarName;
      if (region == "ee")
        effHistTitle = "Endcap " + filterName + "_" + objName + " TrigTagProbeEff vs " + vsVarName;
    }  //std::cout<<effHistTitle<<std::endl;
    //------------
    makeEffMonElemFromPassAndFailAndTagTag(
        filterName,
        filterName + "_trigTagProbeEff_" + objName + "_vs_" + vsVarName + "_" + region,
        effHistTitle,
        passNotTag,
        fail,
        passTagTag,
        ibooker,
        igetter);
  }  //end loop over vsVarNames
}
 
void EgHLTOfflineClient::createTrigTagProbeEffHists2Leg(const std::string& filterNameLeg1,
                                                        const std::string& filterNameLeg2,
                                                        const std::string& region,
                                                        const std::vector<std::string>& vsVarNames,
                                                        const std::string& objName,
                                                        DQMStore::IBooker& ibooker,
                                                        DQMStore::IGetter& igetter) {
  for (auto const& vsVarName : vsVarNames) {
    std::string allName(dirName_ + "/Source_Histos/" + filterNameLeg2 + "/" + filterNameLeg2 + "_trigTagProbe_" +
                        objName + "_all_" + vsVarName + "_" + region);
    MonitorElement* all = igetter.get(allName);
    if (all == nullptr) {
      edm::LogInfo("EgHLTOfflineClient") << " couldnt get hist " << allName;
      continue;
    }
 
    std::string Leg2NotLeg1SourceName(dirName_ + "/Source_Histos/" + filterNameLeg2 + "/" + filterNameLeg2 +
                                      "_trigTagProbe_" + objName + "_passLeg2failLeg1_" + vsVarName + "_" + region);
    MonitorElement* Leg2NotLeg1Source = igetter.get(Leg2NotLeg1SourceName);
    if (Leg2NotLeg1Source == nullptr) {
      edm::LogInfo("EgHLTOfflineClient") << " couldnt get hist " << Leg2NotLeg1SourceName;
      continue;
    }
 
    std::string Leg1EffName(dirName_ + "/Client_Histos/" + filterNameLeg1 + "/" + filterNameLeg1 + "_trigTagProbeEff_" +
                            objName + "_vs_" + vsVarName + "_" + region);
    MonitorElement* Leg1Eff = igetter.get(Leg1EffName);
    if (Leg1Eff == nullptr) {
      edm::LogInfo("EgHLTOfflineClient") << " couldnt get hist " << Leg1EffName;
      continue;
    }
 
    std::string effHistTitle(filterNameLeg2 + "_trigTagProbeEff2Leg_" + objName + "_vs_" + vsVarName + "_" +
                             region);  //std::cout<<effHistTitle<<std::endl;
    if (region == "eb" || region == "ee") {
      if (region == "eb")
        effHistTitle = "Barrel " + filterNameLeg2 + "_" + objName + " TrigTagProbeEff2Leg vs " + vsVarName;
      if (region == "ee")
        effHistTitle = "Endcap " + filterNameLeg2 + "_" + objName + " TrigTagProbeEff2Leg vs " + vsVarName;
    }  //std::cout<<effHistTitle<<std::endl;
    makeEffMonElem2Leg(filterNameLeg2,
                       filterNameLeg2 + "_trigTagProbeEff2Leg_" + objName + "_vs_" + vsVarName + "_" + region,
                       effHistTitle,
                       Leg1Eff,
                       Leg2NotLeg1Source,
                       all,
                       ibooker,
                       igetter);
  }  //end loop over vsVarNames
}
 
void EgHLTOfflineClient::createLooseTightTrigEff(const std::vector<std::string>& tightLooseTrigNames,
                                                 const std::string& region,
                                                 const std::vector<std::string>& vsVarNames,
                                                 const std::string& objName,
                                                 DQMStore::IBooker& ibooker,
                                                 DQMStore::IGetter& igetter) {
  for (auto const& vsVarName : vsVarNames) {
    for (auto const& tightLooseTrigName : tightLooseTrigNames) {
      std::vector<std::string> splitString;
      boost::split(splitString, tightLooseTrigName, boost::is_any_of(std::string(":")));
      if (splitString.size() != 2) {
        continue;  //format incorrect
      }
 
      const std::string& tightTrig = splitString[0];
      const std::string& looseTrig = splitString[1];
      MonitorElement* fail = igetter.get(dirName_ + "/Source_Histos/" + tightTrig + "_" + looseTrig + "_" + objName +
                                         "_failTrig_" + vsVarName + "_" + region);
      if (fail == nullptr) {
        continue;
      }
 
      MonitorElement* pass = igetter.get(dirName_ + "/Source_Histos/" + tightTrig + "_" + looseTrig + "_" + objName +
                                         "_passTrig_" + vsVarName + "_" + region);
      if (pass == nullptr) {
        continue;
      }
 
      const std::string newHistName(tightTrig + "_trigEffTo_" + looseTrig + "_" + objName + "_vs_" + vsVarName + "_" +
                                    region);
      //----Morse-----
      std::string effHistTitle(newHistName);  //std::cout<<effHistTitle<<std::endl;
      if (region == "eb" || region == "ee") {
        if (region == "eb")
          effHistTitle = "Barrel " + tightTrig + "_TrigEffTo_" + looseTrig + "_" + objName + " vs " + vsVarName;
        if (region == "ee")
          effHistTitle = "Endcap " + tightTrig + "_TrigEffTo_" + looseTrig + "_" + objName + " vs " + vsVarName;
      }
      //------------
      makeEffMonElemFromPassAndFail("LooseTight", newHistName, effHistTitle, pass, fail, ibooker, igetter);
 
    }  //end loop over trigger pairs
  }    //end loop over vsVarNames
}
//-----Morse-------
EgHLTOfflineClient::MonitorElement* EgHLTOfflineClient::makeEffMonElemFromPassAndAll(const std::string& filterName,
                                                                                     const std::string& name,
                                                                                     const std::string& title,
                                                                                     const MonitorElement* pass,
                                                                                     const MonitorElement* all,
                                                                                     DQMStore::IBooker& ibooker,
                                                                                     DQMStore::IGetter& igetter) {
  TH1F* passHist = pass->getTH1F();
  if (passHist->GetSumw2N() == 0)
    passHist->Sumw2();
  TH1F* allHist = all->getTH1F();
  if (allHist->GetSumw2N() == 0)
    allHist->Sumw2();
  auto* effHist = (TH1F*)passHist->Clone(name.c_str());
  effHist->Divide(passHist, allHist, 1, 1, "B");
  //----Morse---------
  effHist->SetTitle(title.c_str());
  //------------------
  MonitorElement* eff = igetter.get(dirName_ + "/Client_Histos/" + filterName + "/" + name);
  if (eff == nullptr) {
    eff = ibooker.book1D(name, effHist);
  } else {  //I was having problems with collating the histograms, hence why I'm just resetting the histogram value
    *eff->getTH1F() = *effHist;
    delete effHist;
  }
  return eff;
}
 
EgHLTOfflineClient::MonitorElement* EgHLTOfflineClient::makeEffMonElemFromPassAndFailAndTagTag(
    const std::string& filter,
    const std::string& name,
    const std::string& title,
    const MonitorElement* passNotTag,
    const MonitorElement* fail,
    const MonitorElement* tagtag,
    DQMStore::IBooker& ibooker,
    DQMStore::IGetter& igetter) {
  TH1F* passHist = passNotTag->getTH1F();
  if (passHist->GetSumw2N() == 0)
    passHist->Sumw2();
  TH1F* failHist = fail->getTH1F();
  if (failHist->GetSumw2N() == 0)
    failHist->Sumw2();
  TH1F* tagtagHist = tagtag->getTH1F();
  if (tagtagHist->GetSumw2N() == 0)
    tagtagHist->Sumw2();
  auto* numer = (TH1F*)passHist->Clone(name.c_str());
  if (numer->GetSumw2N() == 0)
    numer->Sumw2();
  numer->Add(tagtagHist, passHist, 2, 1);
  auto* denom = (TH1F*)passHist->Clone(name.c_str());
  if (denom->GetSumw2N() == 0)
    denom->Sumw2();
  denom->Add(tagtagHist, passHist, 2, 1);
  denom->Add(failHist, 1);
  auto* effHist = (TH1F*)passHist->Clone(name.c_str());
  //TGraphAsymmErrors *effHist = new TGraphAsymmErrors(numer,denom,"cl=0.683 b(1,1) mode");
  effHist->Divide(numer, denom, 1, 1, "B");
  //effHist->Divide(numer,denom,"cl=0.683 b(1,1) mode");
  //----Morse---------
  effHist->SetTitle(title.c_str());
  //------------------
  MonitorElement* eff = igetter.get(dirName_ + "/Client_Histos/" + filter + "/" + name);
  if (eff == nullptr) {
    eff = ibooker.book1D(name, effHist);
  } else {  //I was having problems with collating the histograms, hence why I'm just resetting the histogram value
    *eff->getTH1F() = *effHist;
    //*eff->getTGraphAsymmErrors()=*effHist;
    delete effHist;
  }
  return eff;
}
 
EgHLTOfflineClient::MonitorElement* EgHLTOfflineClient::makeEffMonElem2Leg(const std::string& filter,
                                                                           const std::string& name,
                                                                           const std::string& title,
                                                                           const MonitorElement* Leg1Eff,
                                                                           const MonitorElement* Leg2NotLeg1Source,
                                                                           const MonitorElement* all,
                                                                           DQMStore::IBooker& ibooker,
                                                                           DQMStore::IGetter& igetter) {
  TH1F* allHist = all->getTH1F();
  if (allHist->GetSumw2() == nullptr)
    allHist->Sumw2();
  TH1F* Leg2NotLeg1SourceHist = Leg2NotLeg1Source->getTH1F();
  if (Leg2NotLeg1SourceHist->GetSumw2() == nullptr)
    Leg2NotLeg1SourceHist->Sumw2();
 
  auto* effHistLeg2NotLeg1 = (TH1F*)allHist->Clone(name.c_str());
  if (effHistLeg2NotLeg1->GetSumw2() == nullptr)
    effHistLeg2NotLeg1->Sumw2();
  effHistLeg2NotLeg1->Divide(Leg2NotLeg1SourceHist, allHist, 1, 1, "B");
 
  TH1F* Leg1EffHist = Leg1Eff->getTH1F();
  if (Leg1EffHist->GetSumw2() == nullptr)
    Leg1EffHist->Sumw2();
 
  auto* effHistTerm1 = (TH1F*)allHist->Clone(name.c_str());
  if (effHistTerm1->GetSumw2() == nullptr)
    effHistTerm1->Sumw2();
  effHistTerm1->Multiply(Leg1EffHist, Leg1EffHist, 1, 1, "B");
 
  auto* effHistTerm2 = (TH1F*)allHist->Clone(name.c_str());
  if (effHistTerm2->GetSumw2() == nullptr)
    effHistTerm2->Sumw2();
  effHistTerm2->Multiply(Leg1EffHist, effHistLeg2NotLeg1, 1, 1, "B");
  effHistTerm2->Scale(2);
 
  auto* effHist = (TH1F*)allHist->Clone(name.c_str());
  if (effHist->GetSumw2() == nullptr)
    effHist->Sumw2();
  effHist->Add(effHistTerm1, effHistTerm2, 1, 1);
  effHist->SetTitle(title.c_str());
 
  MonitorElement* eff = igetter.get(dirName_ + "/Client_Histos/" + filter + "/" + name);
  if (eff == nullptr) {
    eff = ibooker.book1D(name, effHist);
  } else {  //I was having problems with collating the histograms, hence why I'm just resetting the histogram value
    *eff->getTH1F() = *effHist;
    delete effHist;
  }
  return eff;
}
 
//-----Morse-------
EgHLTOfflineClient::MonitorElement* EgHLTOfflineClient::makeEffMonElemFromPassAndFail(const std::string& filterName,
                                                                                      const std::string& name,
                                                                                      const std::string& title,
                                                                                      const MonitorElement* pass,
                                                                                      const MonitorElement* fail,
                                                                                      DQMStore::IBooker& ibooker,
                                                                                      DQMStore::IGetter& igetter) {
  TH1F* failHist = fail->getTH1F();
  if (failHist->GetSumw2N() == 0)
    failHist->Sumw2();
  TH1F* passHist = pass->getTH1F();
  if (passHist->GetSumw2N() == 0)
    passHist->Sumw2();
 
  auto* effHist = (TH1F*)passHist->Clone(name.c_str());
  effHist->Add(failHist);
  effHist->Divide(passHist, effHist, 1, 1, "B");
  //----Morse---------
  effHist->SetTitle(title.c_str());
  //------------------
  MonitorElement* eff = igetter.get(dirName_ + "/Client_Histos/" + filterName + "/" + name);
  if (eff == nullptr) {
    eff = ibooker.book1D(name, effHist);
  } else {  //I was having problems with collating the histograms, hence why I'm just reseting the histogram value
    *eff->getTH1F() = *effHist;
    delete effHist;
  }
  return eff;
}