/afs/cern.ch/work/a/aaltunda/public/www/CMSSW_5_3_13_patch3/src/DQMOffline/Trigger/src/TopElectronHLTOfflineClient.cc

Go to the documentation of this file.

#include "DQMOffline/Trigger/interface/TopElectronHLTOfflineClient.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 "DQMServices/Core/interface/MonitorElement.h"


TopElectronHLTOfflineClient::TopElectronHLTOfflineClient(const edm::ParameterSet& iConfig) : dbe_(NULL)
{
        dbe_ = edm::Service<DQMStore>().operator->();
        
        if (!dbe_) 
        {
                edm::LogError("TopElectronHLTOfflineClient") << "unable to get DQMStore service, upshot is no client histograms will be made";
        }
        
        if(iConfig.getUntrackedParameter<bool>("DQMStore", false)) 
        {
                if (dbe_) dbe_->setVerbose(0);
        }
        
        dirName_=iConfig.getParameter<std::string>("DQMDirName");
        
        if (dbe_)
                dbe_->setCurrentFolder(dirName_);

        hltTag_ = iConfig.getParameter<std::string>("hltTag");
        
        electronIdNames_ = iConfig.getParameter<std::vector<std::string> >("electronIdNames");
        superTriggerNames_ = iConfig.getParameter<std::vector<std::string> >("superTriggerNames");
        electronTriggerNames_ = iConfig.getParameter<std::vector<std::string> >("electronTriggerNames");
        addExtraId_ = iConfig.getParameter<bool>("addExtraId");

        runClientEndLumiBlock_ = iConfig.getParameter<bool>("runClientEndLumiBlock");
        runClientEndRun_ = iConfig.getParameter<bool>("runClientEndRun");
        runClientEndJob_ = iConfig.getParameter<bool>("runClientEndJob");

}


TopElectronHLTOfflineClient::~TopElectronHLTOfflineClient()
{ 
}

void TopElectronHLTOfflineClient::beginJob()
{
        //compose the ME names we need
        
        // Eta regions
        std::vector<std::string> regions;
        regions.push_back("EB");
        regions.push_back("EE");
        
        // Electron IDs, including own extra ID
        std::vector<std::string> eleIdNames;
        for (size_t i = 0; i < electronIdNames_.size(); ++i)
        {
                eleIdNames.push_back(electronIdNames_[i]);
                if (addExtraId_)
                        eleIdNames.push_back(electronIdNames_[i]+"extraId");
        }
        
        std::vector<std::string> vars;
        vars.push_back("_et");
        vars.push_back("_eta");
        vars.push_back("_phi");
        vars.push_back("_isolEm");
        vars.push_back("_isolHad");
        vars.push_back("_minDeltaR");
        vars.push_back("_global_n30jets");
        vars.push_back("_global_sumEt");
        vars.push_back("_gsftrack_etaError");
        vars.push_back("_gsftrack_phiError");
        vars.push_back("_gsftrack_numberOfValidHits");
        vars.push_back("_gsftrack_dzPV");
        
        
        for (size_t i = 0; i < eleIdNames.size(); ++i)
                for (size_t j = 0; j < regions.size(); ++j)
                        for (size_t k = 0; k < vars.size(); ++k)
                                for (size_t l = 0; l < superTriggerNames_.size(); ++l)
                                {
                                        superMeNames_.push_back("ele_"+superTriggerNames_[l]+"_"+regions[j]+"_"+eleIdNames[i]+vars[k] );
                                        for (size_t m = 0; m < electronTriggerNames_.size(); ++m)
                                        {
                                                eleMeNames_.push_back("ele_"+superTriggerNames_[l]+"_"+electronTriggerNames_[m] +"_"+regions[j]+"_"+eleIdNames[i]+vars[k]);
                                        }
                                }

        

}

void TopElectronHLTOfflineClient::endJob() 
{
        if(runClientEndJob_)
                runClient_();
}

void TopElectronHLTOfflineClient::beginRun(const edm::Run& run, const edm::EventSetup& c)
{
}


void TopElectronHLTOfflineClient::endRun(const edm::Run& run, const edm::EventSetup& c)
{
        if(runClientEndRun_)
                runClient_();
}

//dummy analysis function
void TopElectronHLTOfflineClient::analyze(const edm::Event& iEvent,const edm::EventSetup& iSetup)
{
}

void TopElectronHLTOfflineClient::endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c)
{ 
        if(runClientEndLumiBlock_)
                runClient_();
}

void TopElectronHLTOfflineClient::runClient_()
{
        if (!dbe_) return; //we dont have the DQMStore so we cant do anything
        dbe_->setCurrentFolder(dirName_);
        
        size_t k = 0;
        for (size_t i = 0; i < superMeNames_.size(); ++i)
        {
                for (size_t j = 0; j < electronTriggerNames_.size(); ++j)
                {
                        if (k >= eleMeNames_.size())
                                continue;
                        createSingleEffHists(superMeNames_[i], eleMeNames_[k], eleMeNames_[k]+"_eff");
                        ++k;
                }
        }
        superTriggerNames_.size();
        electronTriggerNames_.size();


}

void TopElectronHLTOfflineClient::createSingleEffHists(const std::string& denomName, const std::string& nomName, const std::string& effName)
{ 
        MonitorElement* denom = dbe_->get(dirName_+"/"+denomName);
        
        MonitorElement* nom = dbe_->get(dirName_+"/"+nomName);
        if(nom!=NULL && denom!=NULL)
        {

                makeEffMonElemFromPassAndAll(effName, nom, denom);       
        }
}

        
MonitorElement* TopElectronHLTOfflineClient::makeEffMonElemFromPassAndAll(const std::string& name, const MonitorElement* pass, const MonitorElement* all)
{
        TH1F* passHist = pass->getTH1F();
        if(passHist->GetSumw2N()==0) 
                passHist->Sumw2();
        TH1F* allHist = all->getTH1F();
        if(allHist->GetSumw2N()==0)
                allHist->Sumw2();
        
        TH1F* effHist = (TH1F*) passHist->Clone(name.c_str());
        effHist->Divide(passHist,allHist,1,1,"B");

        MonitorElement* eff = dbe_->get(dirName_+"/"+name);
        if(eff==NULL)
        {
                eff= dbe_->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;
}