TriggerMatchMonitor Class Reference

#include <TriggerMatchMonitor.h>

Inheritance diagram for TriggerMatchMonitor:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
	TriggerMatchMonitor (const edm::ParameterSet &pSet)
	Constructor. More...
	~TriggerMatchMonitor () override
	Destructor. More...
Public Member Functions inherited from one::DQMEDAnalyzer< T >
	DQMEDAnalyzer ()=default
	DQMEDAnalyzer (DQMEDAnalyzer< T... > const &)=delete
	DQMEDAnalyzer (DQMEDAnalyzer< T... > &&)=delete
	~DQMEDAnalyzer () override=default

Private Attributes
edm::EDGetTokenT< reco::BeamSpot >	beamSpotToken_
MonitorElement *	h_passHLTPath1_eta_Tight
MonitorElement *	h_passHLTPath1_phi_Tight
MonitorElement *	h_passHLTPath1_pt_Tight
MonitorElement *	h_passHLTPath2_eta_Tight
MonitorElement *	h_passHLTPath2_phi_Tight
MonitorElement *	h_passHLTPath2_pt_Tight
MonitorElement *	h_totalHLTPath1_eta_Tight
MonitorElement *	h_totalHLTPath1_phi_Tight
MonitorElement *	h_totalHLTPath1_pt_Tight
MonitorElement *	h_totalHLTPath2_eta_Tight
MonitorElement *	h_totalHLTPath2_phi_Tight
MonitorElement *	h_totalHLTPath2_pt_Tight
std::vector< MonitorElement * >	matchHists
edm::ParameterSet	parameters
edm::EDGetTokenT< std::vector< reco::Vertex > >	primaryVerticesToken_
edm::EDGetTokenT< reco::BeamSpot >	theBeamSpotLabel_
std::string	theFolder
edm::EDGetTokenT< edm::View< reco::Muon > >	theMuonCollectionLabel_
edm::EDGetTokenT< edm::View< pat::Muon > >	thePATMuonCollectionLabel_
MuonServiceProxy *	theService
edm::EDGetTokenT< reco::VertexCollection >	theVertexLabel_
std::string	triggerHistName1_
std::string	triggerHistName2_
edm::EDGetTokenT< std::vector< pat::TriggerObjectStandAlone > >	triggerObjects_
std::string	triggerPathName1_
std::string	triggerPathName2_
double	triggerPtThresholdPath1_
double	triggerPtThresholdPath2_
edm::EDGetTokenT< edm::TriggerResults >	triggerResultsToken_

Detailed Description

Definition at line 39 of file TriggerMatchMonitor.h.

Constructor & Destructor Documentation

TriggerMatchMonitor::TriggerMatchMonitor

(

const edm::ParameterSet &

pSet

)

Constructor.

Definition at line 21 of file TriggerMatchMonitor.cc.

References LogTrace, metname, and MuonServiceProxy_cff::MuonServiceProxy.

                                                                    {
  LogTrace(metname)<<"[TriggerMatchMonitor] Parameters initialization";
  
  parameters = pSet;

  // the services
  theService = new MuonServiceProxy(parameters.getParameter<ParameterSet>("ServiceParameters"));
  
  beamSpotToken_ = consumes<reco::BeamSpot >(parameters.getUntrackedParameter<edm::InputTag>("offlineBeamSpot")),
  primaryVerticesToken_ = consumes<std::vector<reco::Vertex> >(parameters.getUntrackedParameter<edm::InputTag>("offlinePrimaryVertices")),
  theMuonCollectionLabel_ = consumes<edm::View<reco::Muon> >  (parameters.getParameter<edm::InputTag>("MuonCollection"));
  thePATMuonCollectionLabel_ = consumes<edm::View<pat::Muon> >  (parameters.getParameter<edm::InputTag>("patMuonCollection"));
  theVertexLabel_          = consumes<reco::VertexCollection>(parameters.getParameter<edm::InputTag>("VertexLabel"));
  theBeamSpotLabel_        = mayConsume<reco::BeamSpot>      (parameters.getParameter<edm::InputTag>("BeamSpotLabel"));
  triggerResultsToken_ = consumes<edm::TriggerResults>(parameters.getUntrackedParameter<edm::InputTag>("triggerResults"));
  triggerObjects_ = consumes<std::vector<pat::TriggerObjectStandAlone>>(parameters.getParameter<edm::InputTag>("triggerObjects"));

  triggerPathName1_ = parameters.getParameter<string>("triggerPathName1");
  triggerHistName1_ = parameters.getParameter<string>("triggerHistName1");
  triggerPtThresholdPath1_ = parameters.getParameter<double>("triggerPtThresholdPath1");
  triggerPathName2_ = parameters.getParameter<string>("triggerPathName2");
  triggerHistName2_ = parameters.getParameter<string>("triggerHistName2");
  triggerPtThresholdPath2_ = parameters.getParameter<double>("triggerPtThresholdPath2");
  theFolder = parameters.getParameter<string>("folder");
}

TriggerMatchMonitor::~TriggerMatchMonitor ( )

override

Destructor.

Definition at line 46 of file TriggerMatchMonitor.cc.

                                          { 
  delete theService;
}

Member Function Documentation

void TriggerMatchMonitor::analyze ( const edm::Event &  iEvent, const edm::EventSetup &  iSetup  )

override

Mu50 test

Definition at line 93 of file TriggerMatchMonitor.cc.

References reco::BeamSpot::BeamWidthX(), reco::BeamSpot::BeamWidthY(), relativeConstraints::error, edm::Event::getByToken(), edm::HandleBase::isValid(), LogTrace, Max(), metname, extraflags_cff::muons, reco::BeamSpot::position(), position, reco::BeamSpot::sigmaZ(), TriggerAnalyzer::triggerObjects, and triggerResults.

                                                                                    {
  
  LogTrace(metname)<<"[TriggerMatchMonitor] Analyze the mu in different eta regions";
  theService->update(iSetup);
  
  edm::Handle<edm::View<reco::Muon> > muons; 
  iEvent.getByToken(theMuonCollectionLabel_,muons);

  edm::Handle<edm::View<pat::Muon> > PATmuons;
  iEvent.getByToken(thePATMuonCollectionLabel_,PATmuons);

  edm::Handle<std::vector<pat::TriggerObjectStandAlone> > triggerObjects;
  iEvent.getByToken(triggerObjects_, triggerObjects);

  Handle<edm::TriggerResults> triggerResults;
  iEvent.getByToken(triggerResultsToken_, triggerResults);

  reco::Vertex::Point posVtx;
  reco::Vertex::Error errVtx;
  Handle<std::vector<reco::Vertex> > recVtxs;
  iEvent.getByToken(primaryVerticesToken_,recVtxs);
  unsigned int theIndexOfThePrimaryVertex = 999.;
  for (unsigned int ind = 0; ind < recVtxs->size(); ++ind) {
     if ( (*recVtxs)[ind].isValid() && !((*recVtxs)[ind].isFake()) ) {
       theIndexOfThePrimaryVertex = ind;
       break;
     }
  }
  if (theIndexOfThePrimaryVertex<100) {
     posVtx = ((*recVtxs)[theIndexOfThePrimaryVertex]).position();
     errVtx = ((*recVtxs)[theIndexOfThePrimaryVertex]).error();
  }
  else {
     LogInfo("RecoMuonValidator") << "reco::PrimaryVertex not found, use BeamSpot position instead\n";
     Handle<reco::BeamSpot> recoBeamSpotHandle;
     iEvent.getByToken(beamSpotToken_,recoBeamSpotHandle);
     reco::BeamSpot bs = *recoBeamSpotHandle;
     posVtx = bs.position();
     errVtx(0,0) = bs.BeamWidthX();
     errVtx(1,1) = bs.BeamWidthY();
     errVtx(2,2) = bs.sigmaZ();
  }
  const reco::Vertex thePrimaryVertex(posVtx,errVtx);

  if(PATmuons.isValid()){//valid pat Muon
    for(const auto & patMuon : *PATmuons){//pat muon loop
      bool Isolated=patMuon.pfIsolationR04().sumChargedHadronPt + TMath::Max(0., patMuon.pfIsolationR04().sumNeutralHadronEt + patMuon.pfIsolationR04().sumPhotonEt - 0.5*patMuon.pfIsolationR04().sumPUPt)  / patMuon.pt() < 0.25;

      if(patMuon.isGlobalMuon() && Isolated && patMuon.isTightMuon(thePrimaryVertex)){//isolated tight muon

        TLorentzVector offlineMuon;
        offlineMuon.SetPtEtaPhiM(patMuon.pt(), patMuon.eta(), patMuon.phi(),0.0);

        const char *ptrmuPath1 = triggerPathName1_.c_str();        // 
        const char *ptrmuPath2 = triggerPathName2_.c_str();        // 
        if(patMuon.pt() > triggerPtThresholdPath1_){
        h_totalHLTPath1_eta_Tight->Fill(patMuon.eta());
        h_totalHLTPath1_pt_Tight->Fill(patMuon.pt());
        h_totalHLTPath1_phi_Tight->Fill(patMuon.phi());
        }
        if(patMuon.pt() > triggerPtThresholdPath2_){
        h_totalHLTPath2_eta_Tight->Fill(patMuon.eta());
        h_totalHLTPath2_pt_Tight->Fill(patMuon.pt());
        h_totalHLTPath2_phi_Tight->Fill(patMuon.phi());
        }
        if(patMuon.triggered(ptrmuPath1) && patMuon.hltObject() !=nullptr){
          TLorentzVector hltMuon;
          hltMuon.SetPtEtaPhiM(patMuon.hltObject()->pt(),patMuon.hltObject()->eta(),patMuon.hltObject()->phi(),0.0);
          double DelRrecoHLT=offlineMuon.DeltaR(hltMuon);

          matchHists[0]->Fill(DelRrecoHLT);   
          matchHists[1]->Fill(DelRrecoHLT);
          matchHists[2]->Fill(patMuon.hltObject()->pt()/patMuon.pt());
          if(DelRrecoHLT<0.2 && patMuon.pt() > triggerPtThresholdPath1_){
            h_passHLTPath1_eta_Tight->Fill(patMuon.eta());
            h_passHLTPath1_pt_Tight->Fill(patMuon.pt());
            h_passHLTPath1_phi_Tight->Fill(patMuon.phi());
          }
          if(patMuon.l1Object() !=nullptr){
            TLorentzVector L1Muon;
            L1Muon.SetPtEtaPhiM(patMuon.l1Object()->pt(),patMuon.l1Object()->eta(),patMuon.l1Object()->phi(),0.0);
            double DelRrecoL1=offlineMuon.DeltaR(L1Muon);
            matchHists[3]->Fill(DelRrecoL1);
            matchHists[4]->Fill(DelRrecoL1);
            matchHists[5]->Fill(patMuon.l1Object()->pt()/patMuon.pt());     
          }

        }

        if(patMuon.triggered(ptrmuPath2)){
          TLorentzVector hltMuon50;
          hltMuon50.SetPtEtaPhiM(patMuon.hltObject()->pt(),patMuon.hltObject()->eta(),patMuon.hltObject()->phi(),0.0);
          double DelRrecoHLT50=offlineMuon.DeltaR(hltMuon50);

          matchHists[6]->Fill(DelRrecoHLT50);
          matchHists[7]->Fill(DelRrecoHLT50);
          matchHists[8]->Fill(patMuon.hltObject()->pt()/patMuon.pt());
          if(DelRrecoHLT50<0.2 && patMuon.pt() > triggerPtThresholdPath2_ ){
            h_passHLTPath2_eta_Tight->Fill(patMuon.eta());
            h_passHLTPath2_pt_Tight->Fill(patMuon.pt());
            h_passHLTPath2_phi_Tight->Fill(patMuon.phi());
          }


          if(patMuon.l1Object() !=nullptr){
            TLorentzVector L1Muon50;
            L1Muon50.SetPtEtaPhiM(patMuon.l1Object()->pt(),patMuon.l1Object()->eta(),patMuon.l1Object()->phi(),0.0);
            double DelRrecoL150=offlineMuon.DeltaR(L1Muon50);
            matchHists[9]->Fill(DelRrecoL150);
            matchHists[10]->Fill(DelRrecoL150);
            matchHists[11]->Fill(patMuon.l1Object()->pt()/patMuon.pt());
          }
        }
      }//isolated tight muon
    } //pat muon loop
  } //valid pat muon
 
}

void TriggerMatchMonitor::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

override

Definition at line 50 of file TriggerMatchMonitor.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), and DQMStore::IBooker::setCurrentFolder().

                                             {
  ibooker.cd();
  ibooker.setCurrentFolder(theFolder);

  // monitoring of trigger match parameter
   
  matchHists.push_back(ibooker.book1D("DelR_HLT_"+triggerHistName1_+"_v1", "DeltaR_(offline,HLT)_triggerPass("+triggerHistName1_+")", 500, 0.0, 0.5));
  matchHists.push_back(ibooker.book1D("DelR_HLT_"+triggerHistName1_+"_v2", "DeltaR_(offline,HLT)_triggerPass("+triggerHistName1_+")", 100, 0.5, 1.5));
  matchHists.push_back(ibooker.book1D("PtRatio_HLT_"+triggerHistName1_, "PtRatio_(HLTPt/OfflinePt)_"+triggerHistName1_, 200, -5., 5.0));
  
  matchHists.push_back(ibooker.book1D("DelR_L1_"+triggerHistName1_+"_v1", "DeltaR_(offline, L1)_triggerPass("+triggerHistName1_+")", 500, 0.0, 1.0));
  matchHists.push_back(ibooker.book1D("DelR_L1_"+triggerHistName1_+"_v2", "DeltaR_(offline, L1)_triggerPass("+triggerHistName1_+")", 500, 0.0, 2.0));  
  matchHists.push_back(ibooker.book1D("PtRatio_L1_"+triggerHistName1_, "PtRatio_(HLTPt/OfflinePt)_"+triggerHistName1_, 200, -5., 5.0));
   
  matchHists.push_back(ibooker.book1D("DelR_HLT_"+triggerHistName2_+"_v1", "DeltaR_(offline,HLT)_triggerPass("+triggerHistName2_+")", 500, 0.0, 0.5));
  matchHists.push_back(ibooker.book1D("DelR_HLT_"+triggerHistName2_+"_v2", "DeltaR_(offline,HLT)_triggerPass("+triggerHistName2_+")", 100, 0.5, 1.5));
  matchHists.push_back(ibooker.book1D("PtRatio_HLT_"+triggerHistName2_, "PtRatio_(HLTPt/OfflinePt)_"+triggerHistName2_, 200, -5., 5.0));

  matchHists.push_back(ibooker.book1D("DelR_L1_"+triggerHistName2_+"_v1", "DeltaR_(offline, L1)_triggerPass("+triggerHistName2_+")", 250, 0.0, 0.5));
  matchHists.push_back(ibooker.book1D("DelR_L1_"+triggerHistName2_+"_v2", "DeltaR_(offline, L1)_triggerPass("+triggerHistName2_+")", 100, 0.5, 1.5));
  matchHists.push_back(ibooker.book1D("PtRatio_L1_"+triggerHistName2_, "PtRatio_(HLTPt/OfflinePt)_"+triggerHistName2_, 200, -5., 5.0));

  ibooker.cd();
  ibooker.setCurrentFolder(theFolder+"/EfficiencyInput");
 
  h_passHLTPath1_eta_Tight          = ibooker.book1D("passHLT"+triggerHistName1_+"_eta_Tight", " HLT("+triggerHistName1_+") pass #eta", 8, -2.5, 2.5);
  h_passHLTPath1_pt_Tight          = ibooker.book1D("passHLT"+triggerHistName1_+"_pt_Tight", " HLT("+triggerHistName1_+") pass pt", 10, 20, 220);
  h_passHLTPath1_phi_Tight          = ibooker.book1D("passHLT"+triggerHistName1_+"_phi_Tight", " HLT("+triggerHistName1_+") pass phi", 8, -3.0, 3.0);
  h_totalHLTPath1_eta_Tight          = ibooker.book1D("totalHLT"+triggerHistName1_+"_eta_Tight", " HLT("+triggerHistName1_+") total #eta", 8, -2.5, 2.5);
  h_totalHLTPath1_pt_Tight          = ibooker.book1D("totalHLT"+triggerHistName1_+"_pt_Tight", " HLT("+triggerHistName1_+") total pt", 10, 20., 220);
  h_totalHLTPath1_phi_Tight          = ibooker.book1D("totalHLT"+triggerHistName1_+"_phi_Tight", " HLT("+triggerHistName1_+") total phi", 8, -3.0, 3.0);

  h_passHLTPath2_eta_Tight          = ibooker.book1D("passHLT"+triggerHistName2_+"_eta_Tight", " HLT("+triggerHistName2_+") pass #eta", 8, -2.5, 2.5);
  h_passHLTPath2_pt_Tight          = ibooker.book1D("passHLT"+triggerHistName2_+"_pt_Tight", " HLT("+triggerHistName2_+") pass pt", 10, 20., 220);
  h_passHLTPath2_phi_Tight          = ibooker.book1D("passHLT"+triggerHistName2_+"_phi_Tight", " HLT("+triggerHistName2_+") pass phi", 8, -3.0, 3.0);
  h_totalHLTPath2_eta_Tight          = ibooker.book1D("totalHLT"+triggerHistName2_+"_eta_Tight", " HLT("+triggerHistName2_+") total #eta", 8, -2.5, 2.5);
  h_totalHLTPath2_pt_Tight          = ibooker.book1D("totalHLT"+triggerHistName2_+"_pt_Tight", " HLT("+triggerHistName2_+") total pt", 10, 20, 220);
  h_totalHLTPath2_phi_Tight          = ibooker.book1D("totalHLT"+triggerHistName2_+"_phi_Tight", " HLT("+triggerHistName2_+") total phi", 8, -3.0, 3.0);


}

Member Data Documentation

edm::EDGetTokenT<reco::BeamSpot > TriggerMatchMonitor::beamSpotToken_

private

Definition at line 73 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath1_eta_Tight

private

Definition at line 78 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath1_phi_Tight

private

Definition at line 80 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath1_pt_Tight

private

Definition at line 79 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath2_eta_Tight

private

Definition at line 85 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath2_phi_Tight

private

Definition at line 87 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_passHLTPath2_pt_Tight

private

Definition at line 86 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath1_eta_Tight

private

Definition at line 81 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath1_phi_Tight

private

Definition at line 83 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath1_pt_Tight

private

Definition at line 82 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath2_eta_Tight

private

Definition at line 88 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath2_phi_Tight

private

Definition at line 90 of file TriggerMatchMonitor.h.

MonitorElement* TriggerMatchMonitor::h_totalHLTPath2_pt_Tight

private

Definition at line 89 of file TriggerMatchMonitor.h.

std::vector<MonitorElement*> TriggerMatchMonitor::matchHists

private

Definition at line 77 of file TriggerMatchMonitor.h.

edm::ParameterSet TriggerMatchMonitor::parameters

private

Definition at line 55 of file TriggerMatchMonitor.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

edm::EDGetTokenT<std::vector<reco::Vertex> > TriggerMatchMonitor::primaryVerticesToken_

private

Definition at line 74 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<reco::BeamSpot> TriggerMatchMonitor::theBeamSpotLabel_

private

Definition at line 69 of file TriggerMatchMonitor.h.

std::string TriggerMatchMonitor::theFolder

private

Definition at line 92 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<edm::View<reco::Muon> > TriggerMatchMonitor::theMuonCollectionLabel_

private

Definition at line 66 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<edm::View<pat::Muon> > TriggerMatchMonitor::thePATMuonCollectionLabel_

private

Definition at line 67 of file TriggerMatchMonitor.h.

MuonServiceProxy* TriggerMatchMonitor::theService

private

Definition at line 54 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<reco::VertexCollection> TriggerMatchMonitor::theVertexLabel_

private

Definition at line 68 of file TriggerMatchMonitor.h.

std::string TriggerMatchMonitor::triggerHistName1_

private

Definition at line 59 of file TriggerMatchMonitor.h.

std::string TriggerMatchMonitor::triggerHistName2_

private

Definition at line 62 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<std::vector<pat::TriggerObjectStandAlone> > TriggerMatchMonitor::triggerObjects_

private

Definition at line 71 of file TriggerMatchMonitor.h.

std::string TriggerMatchMonitor::triggerPathName1_

private

Definition at line 58 of file TriggerMatchMonitor.h.

std::string TriggerMatchMonitor::triggerPathName2_

private

Definition at line 61 of file TriggerMatchMonitor.h.

double TriggerMatchMonitor::triggerPtThresholdPath1_

private

Definition at line 60 of file TriggerMatchMonitor.h.

double TriggerMatchMonitor::triggerPtThresholdPath2_

private

Definition at line 63 of file TriggerMatchMonitor.h.

edm::EDGetTokenT<edm::TriggerResults > TriggerMatchMonitor::triggerResultsToken_

private

Definition at line 70 of file TriggerMatchMonitor.h.