HotlineDQM.cc

Go to the documentation of this file.

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Framework/interface/Frameworkfwd.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Common/interface/TriggerNames.h"
#include "DataFormats/HLTReco/interface/TriggerObject.h"
#include "DQMOffline/Trigger/interface/HotlineDQM.h"

HotlineDQM::HotlineDQM(const edm::ParameterSet& ps)
{

  edm::LogInfo("HotlineDQM") << "Constructor HotlineDQM::HotlineDQM " << std::endl;
  // Get parameters from configuration file
  theMuonCollection_ = consumes<reco::MuonCollection>(ps.getParameter<edm::InputTag>("muonCollection")); 
  thePfMETCollection_ = consumes<reco::PFMETCollection>(ps.getParameter<edm::InputTag>("pfMetCollection"));
  theMETCollection_ = consumes<reco::CaloMETCollection>(ps.getParameter<edm::InputTag>("caloMetCollection"));
  theCaloJetCollection_ = consumes<reco::CaloJetCollection>(ps.getParameter<edm::InputTag>("caloJetCollection"));
  thePhotonCollection_ = consumes<reco::PhotonCollection>(ps.getParameter<edm::InputTag>("photonCollection"));
  theTrigSummary_ = consumes<trigger::TriggerEvent>(ps.getParameter<edm::InputTag>("trigSummary"));
  triggerResults_ = consumes<edm::TriggerResults>(ps.getParameter<edm::InputTag>("triggerResults"));
  triggerPath_ = ps.getParameter<std::string>("triggerPath");
  triggerFilter_ = ps.getParameter<edm::InputTag>("triggerFilter");
  useMuons = ps.getParameter<bool>("useMuons");
  useMet = ps.getParameter<bool>("useMet");
  usePFMet = ps.getParameter<bool>("usePFMet");
  useHT = ps.getParameter<bool>("useHT");
  usePhotons = ps.getParameter<bool>("usePhotons");
}

HotlineDQM::~HotlineDQM()
{
   edm::LogInfo("HotlineDQM") << "Destructor HotlineDQM::~HotlineDQM " << std::endl;
}

void HotlineDQM::bookHistograms(DQMStore::IBooker & ibooker_, edm::Run const &, edm::EventSetup const &)
{
    edm::LogInfo("HotlineDQM") << "HotlineDQM::bookHistograms" << std::endl;

    ibooker_.cd();
    ibooker_.setCurrentFolder("HLT/Hotline/" + triggerPath_);

    //online quantities 
    h_MuPt = ibooker_.book1D("MuPt", "Muon Pt; GeV", 20, 0.0, 2000.0);
    h_PhotonPt = ibooker_.book1D("PhotonPt", "Photon Pt; GeV", 20, 0.0, 4000.0);
    h_HT = ibooker_.book1D("HT", "HT; GeV", 20, 0.0, 6000.0);
    h_MetPt = ibooker_.book1D("MetPt", "Calo MET; GeV", 20, 0.0, 2000);
    h_PFMetPt = ibooker_.book1D("PFMetPt", "PF MET; GeV", 20, 0.0, 2000);

    if(useMuons) h_OnlineMuPt = ibooker_.book1D("OnlineMuPt", "Online Muon Pt; GeV", 20, 0.0, 2000.0);
    if(usePhotons) h_OnlinePhotonPt = ibooker_.book1D("OnlinePhotonPt", "Online Photon Pt; GeV", 20, 0.0, 4000.0);
    if(useHT) h_OnlineHT = ibooker_.book1D("OnlineHT", "Online HT; GeV", 20, 0.0, 6000.0);
    if(useMet) h_OnlineMetPt = ibooker_.book1D("OnlineMetPt", "Online Calo MET; GeV", 20, 0.0, 2000);
    if(usePFMet) h_OnlinePFMetPt = ibooker_.book1D("OnlinePFMetPt", "Online PF MET; GeV", 20, 0.0, 2000);

    ibooker_.cd();
}

void HotlineDQM::analyze(edm::Event const& e, edm::EventSetup const& eSetup){
    edm::LogInfo("HotlineDQM") << "HotlineDQM::analyze" << std::endl;

    //-------------------------------
    //--- MET
    //-------------------------------
    edm::Handle<reco::PFMETCollection> pfMETCollection;
    e.getByToken(thePfMETCollection_, pfMETCollection);
    if ( !pfMETCollection.isValid() ){
        edm::LogError ("HotlineDQM") << "invalid collection: PFMET" << "\n";
        return;
    }
    edm::Handle<reco::CaloMETCollection> caloMETCollection;
    e.getByToken(theMETCollection_, caloMETCollection);
    if ( !caloMETCollection.isValid() ){
        edm::LogError ("HotlineDQM") << "invalid collection: CaloMET" << "\n";
        return;
    }

    //-------------------------------
    //--- Jets
    //-------------------------------
    edm::Handle<reco::CaloJetCollection> caloJetCollection;
    e.getByToken (theCaloJetCollection_,caloJetCollection);
    if ( !caloJetCollection.isValid() ){
        edm::LogError ("HotlineDQM") << "invalid collection: CaloJets" << "\n";
        return;
    }

    //-------------------------------
    //--- Muon
    //-------------------------------
    edm::Handle<reco::MuonCollection> MuonCollection;
    e.getByToken (theMuonCollection_, MuonCollection);
    if ( !MuonCollection.isValid() ){
        edm::LogError ("HotlineDQM") << "invalid collection: Muons " << "\n";
        return;
    }

    //-------------------------------
    //--- Photon 
    //-------------------------------
    edm::Handle<reco::PhotonCollection> PhotonCollection;
    e.getByToken (thePhotonCollection_, PhotonCollection);
    if ( !PhotonCollection.isValid() ){
        edm::LogError ("HotlineDQM") << "invalid collection: Photons " << "\n";
        return;
    }

    //-------------------------------
    //--- Trigger
    //-------------------------------
    edm::Handle<edm::TriggerResults> hltresults;
    e.getByToken(triggerResults_,hltresults);
    if(!hltresults.isValid()){
        edm::LogError ("HotlineDQM") << "invalid collection: TriggerResults" << "\n";
        return;
    }
    edm::Handle<trigger::TriggerEvent> triggerSummary;
    e.getByToken(theTrigSummary_, triggerSummary);
    if(!triggerSummary.isValid()) {
        edm::LogError ("HotlineDQM") << "invalid collection: TriggerSummary" << "\n";
        return;
    }

    bool hasFired = false;
    const edm::TriggerNames& trigNames = e.triggerNames(*hltresults);
    unsigned int numTriggers = trigNames.size();
    for( unsigned int hltIndex=0; hltIndex<numTriggers; ++hltIndex ){
        if (trigNames.triggerName(hltIndex).find(triggerPath_) != std::string::npos && hltresults->wasrun(hltIndex) && hltresults->accept(hltIndex)){
            hasFired = true;
        }
    }

    //get online objects
    float ptMuon=-1, ptPhoton=-1, met=-1, pfMet=-1, ht = 0;
    size_t filterIndex = triggerSummary->filterIndex( triggerFilter_ );
    trigger::TriggerObjectCollection triggerObjects = triggerSummary->getObjects();
    if( !(filterIndex >= triggerSummary->sizeFilters()) ){
        const trigger::Keys& keys = triggerSummary->filterKeys( filterIndex );
        for( size_t j = 0; j < keys.size(); ++j ){
            trigger::TriggerObject foundObject = triggerObjects[keys[j]];
            if(useMuons && fabs(foundObject.id()) == 13){ //muon
                if(foundObject.pt() > ptMuon) ptMuon = foundObject.pt();
            }
            else if(usePhotons && fabs(foundObject.id()) == 0){ //photon
                if(foundObject.pt() > ptPhoton) ptPhoton = foundObject.pt();
            }
            else if(useMet && fabs(foundObject.id()) == 0){ //MET
                met = foundObject.pt(); 
            }
            else if(usePFMet && fabs(foundObject.id()) == 0){ //PFMET
                pfMet = foundObject.pt();
            }
            else if(useHT && fabs(foundObject.id()) == 89){ //HT
                ht = foundObject.pt();
            }
        }
    }

    if(hasFired) {
        //fill appropriate online histogram
        if(useMuons) h_OnlineMuPt->Fill(ptMuon);
        if(usePhotons) h_OnlinePhotonPt->Fill(ptPhoton);
        if(useMet) h_OnlineMetPt->Fill(met);
        if(usePFMet) h_OnlinePFMetPt->Fill(pfMet);
        if(useHT) h_OnlineHT->Fill(ht);

        //fill muon pt histogram
        if(MuonCollection->size() > 0){
            float maxMuPt = -1.0;
            for(auto &mu : *MuonCollection){
                if(mu.pt() > maxMuPt) maxMuPt = mu.pt();
            }
            h_MuPt->Fill(maxMuPt);
        }

        //fill photon pt histogram
        if(PhotonCollection->size() > 0){
            float maxPhoPt = -1.0;
            for(auto &pho : *PhotonCollection){
                if(pho.pt() > maxPhoPt) maxPhoPt = pho.pt();
            }
            h_PhotonPt->Fill(maxPhoPt);
        }

        //fill HT histogram
        float caloHT = 0.0;
        for (reco::CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin(); i_calojet != caloJetCollection->end(); ++i_calojet){
            if (i_calojet->pt() < 40) continue;
            if (fabs(i_calojet->eta()) > 3.0) continue;
            caloHT += i_calojet->pt();
        }
        h_HT->Fill(caloHT);

        //fill CaloMET histogram
        h_MetPt->Fill(caloMETCollection->front().et()); 

        //fill PFMET histogram
        h_PFMetPt->Fill(pfMETCollection->front().et());
    }
}

void HotlineDQM::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("photonCollection", edm::InputTag("photons"));
  desc.add<edm::InputTag>("muonCollection", edm::InputTag("muons"));
  desc.add<edm::InputTag>("caloJetCollection", edm::InputTag("ak4CaloJets"));
  desc.add<edm::InputTag>("pfMetCollection", edm::InputTag("pfMet"));
  desc.add<edm::InputTag>("caloMetCollection", edm::InputTag("caloMet"));
  desc.add<edm::InputTag>("triggerResults",edm::InputTag("TriggerResults","","HLT"));
  desc.add<edm::InputTag>("trigSummary",edm::InputTag("hltTriggerSummaryAOD"));
  desc.add<std::string>("triggerPath","HLT_HT2000_v")->setComment("trigger path name");
  desc.add<edm::InputTag>("triggerFilter",edm::InputTag("hltHt2000","","HLT"))->setComment("name of the last filter in the path");
  desc.add<bool>("useMuons", false);
  desc.add<bool>("usePhotons", false);
  desc.add<bool>("useMet", false);
  desc.add<bool>("usePFMet", false);
  desc.add<bool>("useHT", false);
  descriptions.add("HotlineDQM",desc);
}

//define this as a plug-in
DEFINE_FWK_MODULE(HotlineDQM);