HiggsDQM Class Reference

#include <HiggsDQM.h>

Inheritance diagram for HiggsDQM:

Public Member Functions
	HiggsDQM (const edm::ParameterSet &ps)
	~HiggsDQM () override
Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistos (DQMStore *bei)
double	calcDeltaPhi (double phi1, double phi2)
double	Distance (const reco::Candidate &c1, const reco::Candidate &c2)
double	DistancePhi (const reco::Candidate &c1, const reco::Candidate &c2)

Private Attributes
MonitorElement *	h_caloMet
MonitorElement *	h_caloMet_phi
MonitorElement *	h_dielemass
MonitorElement *	h_dimumass_GMGM
MonitorElement *	h_dimumass_GMTM
MonitorElement *	h_dimumass_TMTM
MonitorElement *	h_eCombIso
MonitorElement *	h_eEta
MonitorElement *	h_eMultiplicity
MonitorElement *	h_ePhi
MonitorElement *	h_ePt
MonitorElement *	h_jet2_et
MonitorElement *	h_jet_count
MonitorElement *	h_jet_et
MonitorElement *	h_lepcounts
MonitorElement *	h_mCombIso
MonitorElement *	h_mEta_GM
MonitorElement *	h_mEta_GMPT
MonitorElement *	h_mEta_GMTM
MonitorElement *	h_mEta_STAM
MonitorElement *	h_mEta_TM
MonitorElement *	h_mMultiplicity
MonitorElement *	h_mPhi_GM
MonitorElement *	h_mPhi_GMPT
MonitorElement *	h_mPhi_GMTM
MonitorElement *	h_mPhi_STAM
MonitorElement *	h_mPhi_TM
MonitorElement *	h_mPt_GM
MonitorElement *	h_mPt_GMPT
MonitorElement *	h_mPt_GMTM
MonitorElement *	h_mPt_STAM
MonitorElement *	h_mPt_TM
MonitorElement *	h_pfMet
MonitorElement *	h_pfMet_phi
MonitorElement *	h_vertex_chi2
MonitorElement *	h_vertex_d0
MonitorElement *	h_vertex_number
MonitorElement *	h_vertex_numTrks
MonitorElement *	h_vertex_sumTrks
int	ievt
int	irun
int	leptonflavor
reco::CandidateCollection *	leptonscands_
unsigned long long	m_cacheID_
int	nElectron
int	nEvents_
int	nfourlept
int	nHiggs
int	nLepton
int	nLooseIsolEle
int	nLooseIsolMu
int	nMuon
int	nZEE
int	nZMuMu
float	pi
double	ptThrMu1_
double	ptThrMu2_
edm::InputTag	theCaloJetCollectionLabel_
edm::EDGetTokenT < reco::CaloJetCollection >	theCaloJetCollectionToken_
edm::EDGetTokenT < reco::CaloMETCollection >	theCaloMETCollectionToken_
std::string	theElecTriggerPathToPass
edm::EDGetTokenT < reco::GsfElectronCollection >	theElectronCollectionToken_
edm::EDGetTokenT < reco::MuonCollection >	theMuonCollectionToken_
std::string	theMuonTriggerPathToPass
edm::EDGetTokenT < reco::PFMETCollection >	thePfMETCollectionToken_
edm::EDGetTokenT < edm::TriggerResults >	theTriggerResultsCollection_
edm::InputTag	theTriggerResultsCollectionTag_
edm::EDGetTokenT < reco::VertexCollection >	vertexToken_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr < DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

Definition at line 42 of file HiggsDQM.h.

Constructor & Destructor Documentation

HiggsDQM::HiggsDQM

(

const edm::ParameterSet &

ps

)

Definition at line 118 of file HiggsDQM.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, and pi.

                                            {
  // cout<<"Entering  HiggsDQM::HiggsDQM: "<<endl;

  edm::LogInfo("HZZ4LeptonsDQM") << " Creating HZZ4LeptonsDQM "
                                 << "\n";

  typedef std::vector<edm::InputTag> vtag;
  // Get parameters from configuration file
  theElecTriggerPathToPass = ps.getParameter<string>("elecTriggerPathToPass");
  theMuonTriggerPathToPass = ps.getParameter<string>("muonTriggerPathToPass");
  theTriggerResultsCollectionTag_ = ps.getParameter<InputTag>("triggerResultsCollection");
  theCaloJetCollectionLabel_ = ps.getParameter<InputTag>("caloJetCollection");
  theTriggerResultsCollection_ = consumes<edm::TriggerResults>(theTriggerResultsCollectionTag_);
  theMuonCollectionToken_ = consumes<reco::MuonCollection>(ps.getParameter<InputTag>("muonCollection"));
  theElectronCollectionToken_ = consumes<reco::GsfElectronCollection>(ps.getParameter<InputTag>("electronCollection"));
  theCaloJetCollectionToken_ = consumes<reco::CaloJetCollection>(theCaloJetCollectionLabel_);
  theCaloMETCollectionToken_ = consumes<reco::CaloMETCollection>(ps.getParameter<InputTag>("caloMETCollection"));
  thePfMETCollectionToken_ = consumes<reco::PFMETCollection>(ps.getParameter<InputTag>("pfMETCollection"));
  vertexToken_ = consumes<reco::VertexCollection>(
      ps.getUntrackedParameter<InputTag>("vertexCollection", InputTag("offlinePrimaryVertices")));

  // cuts:
  ptThrMu1_ = ps.getUntrackedParameter<double>("PtThrMu1");
  ptThrMu2_ = ps.getUntrackedParameter<double>("PtThrMu2");

  nEvents_ = 0;
  pi = 3.14159265;
  // cout<<"...leaving  HiggsDQM::HiggsDQM. "<<endl;
}

HiggsDQM::~HiggsDQM ( )

override

Definition at line 150 of file HiggsDQM.cc.

                    {
  // cout<<"Entering HiggsDQM::~HiggsDQM: "<<endl;

  edm::LogInfo("HiggsDQM") << " Deleting HiggsDQM "
                           << "\n";

  // cout<<"...leaving HiggsDQM::~HiggsDQM. "<<endl;
}

Member Function Documentation

void HiggsDQM::analyze ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 224 of file HiggsDQM.cc.

References CaloMET_cfi::caloMet, edm::EventID::event(), edm::Event::getByToken(), muon::GlobalMuonPromptTight, mps_fire::i, edm::EventBase::id(), reco::Muon::isGlobalMuon(), muon::isGoodMuon(), reco::Muon::isTrackerMuon(), edm::HandleBase::isValid(), dqmiolumiharvest::j, LogDebug, callgraph::m2, RPCpg::mu, reco::LeafCandidate::p(), edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), edm::EventID::run(), reco::Muon::SegmentAndTrackArbitration, muon::segmentCompatibility(), mathSSE::sqrt(), findQualityFiles::v, and GoodVertex_cfg::vertexCollection.

                                                                     {
  // cout<<"[HiggsDQM::analyze()] "<<endl;

  //-------------------------------
  //--- Trigger Info
  //-------------------------------
  // Did it pass certain HLT path?
  bool passed_electron_HLT = true;
  bool passed_muon_HLT = true;

  //-------------------------------
  //--- Vertex Info
  //-------------------------------
  Handle<VertexCollection> vertexHandle;
  e.getByToken(vertexToken_, vertexHandle);
  if (vertexHandle.isValid()) {
    VertexCollection vertexCollection = *(vertexHandle.product());
    int vertex_number = vertexCollection.size();
    VertexCollection::const_iterator v = vertexCollection.begin();
    double vertex_chi2 = v->normalizedChi2();  // v->chi2();
    double vertex_d0 = sqrt(v->x() * v->x() + v->y() * v->y());
    // double vertex_ndof    = v->ndof();cout << "ndof="<<vertex_ndof<<endl;
    double vertex_numTrks = v->tracksSize();
    double vertex_sumTrks = 0.0;
    for (Vertex::trackRef_iterator vertex_curTrack = v->tracks_begin(); vertex_curTrack != v->tracks_end();
         vertex_curTrack++) {
      vertex_sumTrks += (*vertex_curTrack)->pt();
    }
    h_vertex_number->Fill(vertex_number);
    h_vertex_chi2->Fill(vertex_chi2);
    h_vertex_d0->Fill(vertex_d0);
    h_vertex_numTrks->Fill(vertex_numTrks);
    h_vertex_sumTrks->Fill(vertex_sumTrks);
  }

  //-------------------------------
  //--- Electrons
  //-------------------------------
  float nEle = 0;
  Handle<GsfElectronCollection> electronCollection;
  e.getByToken(theElectronCollectionToken_, electronCollection);
  if (electronCollection.isValid()) {
    int posEle = 0, negEle = 0;
    // If it passed electron HLT and the collection was found, find electrons
    // near Z mass
    if (passed_electron_HLT) {
      for (reco::GsfElectronCollection::const_iterator recoElectron = electronCollection->begin();
           recoElectron != electronCollection->end();
           recoElectron++) {
        //      cout << "Electron with pt= " <<  recoElectron->pt() << " and
        // eta" << recoElectron->eta() << " p=" <<  recoElectron->p() << endl;
        h_ePt->Fill(recoElectron->pt());
        h_eEta->Fill(recoElectron->eta());
        h_ePhi->Fill(recoElectron->phi());
        if (recoElectron->charge() == 1) {
          posEle++;
        } else if (recoElectron->charge() == -1) {
          negEle++;
        }
        // Require electron to pass some basic cuts
        // if ( recoElectron->et() < 20 || fabs(recoElectron->eta())>2.5 )
        // continue;
        // Tighter electron cuts
        // if ( recoElectron->deltaPhiSuperClusterTrackAtVtx() > 0.58 ||
        //     recoElectron->deltaEtaSuperClusterTrackAtVtx() > 0.01 ||
        //     recoElectron->sigmaIetaIeta() > 0.027 ) continue;
      }  // end of loop over electrons
    }    // end if passed HLT
    nEle = posEle + negEle;
    if (nEle > 9.)
      nEle = 9.;
    h_eMultiplicity->Fill(nEle);

    // Z->ee:
    unsigned int eleCollectionSize = electronCollection->size();
    for (unsigned int i = 0; i < eleCollectionSize; i++) {
      const GsfElectron& ele = electronCollection->at(i);
      double pt = ele.pt();
      if (pt > ptThrMu1_) {
        for (unsigned int j = i + 1; j < eleCollectionSize; j++) {
          const GsfElectron& ele2 = electronCollection->at(j);
          double pt2 = ele2.pt();
          if (pt2 > ptThrMu2_) {
            const math::XYZTLorentzVector ZRecoEE(
                ele.px() + ele2.px(), ele.py() + ele2.py(), ele.pz() + ele2.pz(), ele.p() + ele2.p());
            h_dielemass->Fill(ZRecoEE.mass());
          }
        }
      }
    }
  }

  //-------------------------------
  //--- Muons
  //-------------------------------
  float nMu = 0;
  Handle<MuonCollection> muonCollection;
  e.getByToken(theMuonCollectionToken_, muonCollection);
  if (muonCollection.isValid()) {
    // Find the highest pt muons
    int posMu = 0, negMu = 0;
    TLorentzVector m1, m2;
    if (passed_muon_HLT) {
      for (reco::MuonCollection::const_iterator recoMuon = muonCollection->begin(); recoMuon != muonCollection->end();
           recoMuon++) {
        // cout << "Muon with pt= " <<  muIter->pt() << " and eta" <<
        // muIter->eta() << " p=" <<  muIter->p() << endl;
        if (recoMuon->isGlobalMuon() && recoMuon->isTrackerMuon()) {
          h_mPt_GMTM->Fill(recoMuon->pt());
          h_mEta_GMTM->Fill(recoMuon->eta());
          h_mPhi_GMTM->Fill(recoMuon->phi());
        } else if (recoMuon->isGlobalMuon() && (muon::isGoodMuon((*recoMuon), muon::GlobalMuonPromptTight))) {
          h_mPt_GMPT->Fill(recoMuon->pt());
          h_mEta_GMPT->Fill(recoMuon->eta());
          h_mPhi_GMPT->Fill(recoMuon->phi());
        } else if (recoMuon->isGlobalMuon()) {
          h_mPt_GM->Fill(recoMuon->pt());
          h_mEta_GM->Fill(recoMuon->eta());
          h_mPhi_GM->Fill(recoMuon->phi());
        } else if (recoMuon->isTrackerMuon() &&
                   (muon::segmentCompatibility((*recoMuon), reco::Muon::SegmentAndTrackArbitration))) {
          h_mPt_TM->Fill(recoMuon->pt());
          h_mEta_TM->Fill(recoMuon->eta());
          h_mPhi_TM->Fill(recoMuon->phi());
        } else if (recoMuon->isStandAloneMuon()) {
          h_mPt_STAM->Fill(recoMuon->pt());
          h_mEta_STAM->Fill(recoMuon->eta());
          h_mPhi_STAM->Fill(recoMuon->phi());
        }
        if (recoMuon->charge() == 1) {
          posMu++;
        } else if (recoMuon->charge() == -1) {
          negMu++;
        }
      }
      nMu = posMu + negMu;
      if (nMu > 9.)
        nMu = 9.;
      h_mMultiplicity->Fill(nMu);
    }

    // Z->mumu:
    unsigned int muonCollectionSize = muonCollection->size();
    for (unsigned int i = 0; i < muonCollectionSize; i++) {
      const Muon& mu = muonCollection->at(i);
      // if (!mu.isGlobalMuon()) continue;
      double pt = mu.pt();
      if (pt > ptThrMu1_) {
        for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
          const Muon& mu2 = muonCollection->at(j);
          double pt2 = mu2.pt();
          if (pt2 > ptThrMu2_) {
            // Glb + Glb
            if (mu.isGlobalMuon() && mu2.isGlobalMuon()) {
              const math::XYZTLorentzVector ZRecoGMGM(
                  mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(), mu.p() + mu2.p());
              h_dimumass_GMGM->Fill(ZRecoGMGM.mass());
            }
            // Glb + TM
            else if (mu.isGlobalMuon() && mu2.isTrackerMuon()) {
              const math::XYZTLorentzVector ZRecoGMTM(
                  mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(), mu.p() + mu2.p());
              h_dimumass_GMTM->Fill(ZRecoGMTM.mass());
            }
            // TM + TM
            else if (mu.isTrackerMuon() && mu2.isTrackerMuon()) {
              const math::XYZTLorentzVector ZRecoTMTM(
                  mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(), mu.p() + mu2.p());
              h_dimumass_TMTM->Fill(ZRecoTMTM.mass());
            }
          }
        }
      }
    }
  }

  //-------------------------------
  //--- Jets
  //-------------------------------
  Handle<CaloJetCollection> caloJetCollection;
  e.getByToken(theCaloJetCollectionToken_, caloJetCollection);
  if (caloJetCollection.isValid()) {
    float jet_et = -8.0;
    //    float jet_eta   = -8.0; // UNUSED
    //    float jet_phi   = -8.0; // UNUSED
    int jet_count = 0;
    float jet2_et = -9.0;
    //    float jet2_eta  = -9.0; // UNUSED
    //    float jet2_phi  = -9.0; // UNUSED
    for (CaloJetCollection::const_iterator i_calojet = caloJetCollection->begin();
         i_calojet != caloJetCollection->end();
         i_calojet++) {
      float jet_current_et = i_calojet->et();
      // if it overlaps with electron, it is not a jet
      // if ( electron_et>0.0 && fabs(i_calojet->eta()-electron_eta ) < 0.2 &&
      // calcDeltaPhi(i_calojet->phi(), electron_phi ) < 0.2) continue;
      // if ( electron2_et>0.0&& fabs(i_calojet->eta()-electron2_eta) < 0.2 &&
      // calcDeltaPhi(i_calojet->phi(), electron2_phi) < 0.2) continue;
      // if it has too low Et, throw away
      if (jet_current_et < 15)
        continue;
      jet_count++;
      if (jet_current_et > jet_et) {
        jet2_et = jet_et;          // 2nd highest jet get's et from current highest
                                   //        jet2_eta = jet_eta; // UNUSED
                                   //        jet2_phi = jet_phi; // UNUSED
        jet_et = i_calojet->et();  // current highest jet gets et from the new
                                   // highest
        //        jet_eta  = i_calojet->eta(); // UNUSED
        //        jet_phi  = i_calojet->phi(); // UNUSED
      } else if (jet_current_et > jet2_et) {
        jet2_et = i_calojet->et();
        //        jet2_eta = i_calojet->eta(); // UNUSED
        //        jet2_phi = i_calojet->phi(); // UNUSED
      }
    }
    if (jet_et > 0.0) {
      h_jet_et->Fill(jet_et);
      h_jet_count->Fill(jet_count);
    }
  }

  //-------------------------------
  //--- MET
  //-------------------------------
  Handle<CaloMETCollection> caloMETCollection;
  e.getByToken(theCaloMETCollectionToken_, caloMETCollection);
  if (caloMETCollection.isValid()) {
    float caloMet = caloMETCollection->begin()->et();
    float caloMet_phi = caloMETCollection->begin()->phi();
    h_caloMet->Fill(caloMet);
    h_caloMet_phi->Fill(caloMet_phi);
  }
  Handle<PFMETCollection> pfMETCollection;
  e.getByToken(thePfMETCollectionToken_, pfMETCollection);
  if (pfMETCollection.isValid()) {
    float pfMet = pfMETCollection->begin()->et();
    float pfMet_phi = pfMETCollection->begin()->phi();
    h_pfMet->Fill(pfMet);
    h_pfMet_phi->Fill(pfMet_phi);
  }

  //-------------------------------------
  //--- Events with more than 2 leptons:
  //-------------------------------------
  if (nMu + nEle > 2 && nMu + nEle < 10) {
    if (nMu == 0 && nEle == 3)
      h_lepcounts->Fill(0);
    if (nMu == 0 && nEle == 4)
      h_lepcounts->Fill(1);
    if (nMu == 0 && nEle == 5)
      h_lepcounts->Fill(2);
    if (nMu == 0 && nEle == 6)
      h_lepcounts->Fill(3);
    if (nMu == 0 && nEle == 7)
      h_lepcounts->Fill(4);
    if (nMu == 0 && nEle == 8)
      h_lepcounts->Fill(5);
    if (nMu == 0 && nEle == 9)
      h_lepcounts->Fill(6);
    if (nMu == 1 && nEle == 2)
      h_lepcounts->Fill(7);
    if (nMu == 1 && nEle == 3)
      h_lepcounts->Fill(8);
    if (nMu == 1 && nEle == 4)
      h_lepcounts->Fill(9);
    if (nMu == 1 && nEle == 5)
      h_lepcounts->Fill(10);
    if (nMu == 1 && nEle == 6)
      h_lepcounts->Fill(11);
    if (nMu == 1 && nEle == 7)
      h_lepcounts->Fill(12);
    if (nMu == 1 && nEle == 8)
      h_lepcounts->Fill(13);
    if (nMu == 2 && nEle == 1)
      h_lepcounts->Fill(14);
    if (nMu == 2 && nEle == 2)
      h_lepcounts->Fill(15);
    if (nMu == 2 && nEle == 3)
      h_lepcounts->Fill(16);
    if (nMu == 2 && nEle == 4)
      h_lepcounts->Fill(17);
    if (nMu == 2 && nEle == 5)
      h_lepcounts->Fill(18);
    if (nMu == 2 && nEle == 6)
      h_lepcounts->Fill(19);
    if (nMu == 2 && nEle == 7)
      h_lepcounts->Fill(20);
    if (nMu == 3 && nEle == 0)
      h_lepcounts->Fill(21);
    if (nMu == 3 && nEle == 1)
      h_lepcounts->Fill(22);
    if (nMu == 3 && nEle == 2)
      h_lepcounts->Fill(23);
    if (nMu == 3 && nEle == 3)
      h_lepcounts->Fill(24);
    if (nMu == 3 && nEle == 4)
      h_lepcounts->Fill(25);
    if (nMu == 3 && nEle == 5)
      h_lepcounts->Fill(26);
    if (nMu == 3 && nEle == 6)
      h_lepcounts->Fill(27);
    if (nMu == 4 && nEle == 0)
      h_lepcounts->Fill(28);
    if (nMu == 4 && nEle == 1)
      h_lepcounts->Fill(29);
    if (nMu == 4 && nEle == 2)
      h_lepcounts->Fill(30);
    if (nMu == 4 && nEle == 3)
      h_lepcounts->Fill(31);
    if (nMu == 4 && nEle == 4)
      h_lepcounts->Fill(32);
    if (nMu == 4 && nEle == 5)
      h_lepcounts->Fill(33);
    if (nMu == 5 && nEle == 0)
      h_lepcounts->Fill(34);
    if (nMu == 5 && nEle == 1)
      h_lepcounts->Fill(35);
    if (nMu == 5 && nEle == 2)
      h_lepcounts->Fill(36);
    if (nMu == 5 && nEle == 3)
      h_lepcounts->Fill(37);
    if (nMu == 5 && nEle == 4)
      h_lepcounts->Fill(38);
    if (nMu == 6 && nEle == 0)
      h_lepcounts->Fill(39);
    if (nMu == 6 && nEle == 1)
      h_lepcounts->Fill(40);
    if (nMu == 6 && nEle == 2)
      h_lepcounts->Fill(41);
    if (nMu == 6 && nEle == 3)
      h_lepcounts->Fill(42);
    if (nMu == 7 && nEle == 0)
      h_lepcounts->Fill(43);
    if (nMu == 7 && nEle == 1)
      h_lepcounts->Fill(44);
    if (nMu == 7 && nEle == 2)
      h_lepcounts->Fill(45);
    if (nMu == 8 && nEle == 0)
      h_lepcounts->Fill(46);
    if (nMu == 8 && nEle == 1)
      h_lepcounts->Fill(47);
    if (nMu == 9 && nEle == 0)
      h_lepcounts->Fill(48);
  }
  if ((nMu + nEle) >= 10)
    LogDebug("HiggsDQM") << "WARNING: " << nMu + nEle << " leptons in this event: run=" << e.id().run()
                         << ", event=" << e.id().event() << "\n";
}

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 162 of file HiggsDQM.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::NavigatorBase::cd(), and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                           {
  ibooker.setCurrentFolder("Physics/Higgs");

  h_vertex_number = ibooker.book1D("h_vertex_number", "Number of event vertices in collection", 10, -0.5, 9.5);
  h_vertex_chi2 = ibooker.book1D("h_vertex_chi2", "Event Vertex #chi^{2}/n.d.o.f.", 100, 0.0, 2.0);
  h_vertex_numTrks = ibooker.book1D("h_vertex_numTrks", "Event Vertex, number of tracks", 100, -0.5, 99.5);
  h_vertex_sumTrks = ibooker.book1D("h_vertex_sumTrks", "Event Vertex, sum of track pt", 100, 0.0, 100.0);
  h_vertex_d0 = ibooker.book1D("h_vertex_d0", "Event Vertex d0", 100, -10.0, 10.0);
  h_jet_et =
      ibooker.book1D("h_jet_et",
                     "Jet with highest E_{T} (from " + theCaloJetCollectionLabel_.label() + ");E_{T}(1^{st} jet) (GeV)",
                     20,
                     0.,
                     200.0);
  h_jet2_et = ibooker.book1D(
      "h_jet2_et",
      "Jet with 2^{nd} highest E_{T} (from " + theCaloJetCollectionLabel_.label() + ");E_{T}(2^{nd} jet) (GeV)",
      20,
      0.,
      200.0);
  h_jet_count = ibooker.book1D("h_jet_count",
                               "Number of " + theCaloJetCollectionLabel_.label() + " (E_{T} > 15 GeV);Number of Jets",
                               8,
                               -0.5,
                               7.5);
  h_caloMet = ibooker.book1D("h_caloMet", "Calo Missing E_{T}; GeV", 20, 0.0, 100);
  h_caloMet_phi = ibooker.book1D("h_caloMet_phi", "Calo Missing E_{T} #phi;#phi(MET)", 35, -3.5, 3.5);
  h_pfMet = ibooker.book1D("h_pfMet", "Pf Missing E_{T}; GeV", 20, 0.0, 100);
  h_pfMet_phi = ibooker.book1D("h_pfMet_phi", "Pf Missing E_{T} #phi;#phi(MET)", 35, -3.5, 3.5);
  h_eMultiplicity = ibooker.book1D("NElectrons", "# of electrons per event", 10, 0., 10.);
  h_mMultiplicity = ibooker.book1D("NMuons", "# of muons per event", 10, 0., 10.);
  h_ePt = ibooker.book1D("ElePt", "Pt of electrons", 50, 0., 100.);
  h_eEta = ibooker.book1D("EleEta", "Eta of electrons", 100, -5., 5.);
  h_ePhi = ibooker.book1D("ElePhi", "Phi of electrons", 100, -3.5, 3.5);
  h_mPt_GMTM = ibooker.book1D("MuonPt_GMTM", "Pt of global+tracker muons", 50, 0., 100.);
  h_mEta_GMTM = ibooker.book1D("MuonEta_GMTM", "Eta of global+tracker muons", 60, -3., 3.);
  h_mPhi_GMTM = ibooker.book1D("MuonPhi_GMTM", "Phi of global+tracker muons", 70, -3.5, 3.5);
  h_mPt_GMPT = ibooker.book1D("MuonPt_GMPT", "Pt of global prompt-tight muons", 50, 0., 100.);
  h_mEta_GMPT = ibooker.book1D("MuonEta_GMPT", "Eta of global prompt-tight muons", 60, -3., 3.);
  h_mPhi_GMPT = ibooker.book1D("MuonPhi_GMPT", "Phi of global prompt-tight muons", 70, -3.5, 3.5);
  h_mPt_GM = ibooker.book1D("MuonPt_GM", "Pt of global muons", 50, 0., 100.);
  h_mEta_GM = ibooker.book1D("MuonEta_GM", "Eta of global muons", 60, -3., 3.);
  h_mPhi_GM = ibooker.book1D("MuonPhi_GM", "Phi of global muons", 70, -3.5, 3.5);
  h_mPt_TM = ibooker.book1D("MuonPt_TM", "Pt of tracker muons", 50, 0., 100.);
  h_mEta_TM = ibooker.book1D("MuonEta_TM", "Eta of tracker muons", 60, -3., 3.);
  h_mPhi_TM = ibooker.book1D("MuonPhi_TM", "Phi of tracker muons", 70, -3.5, 3.5);
  h_mPt_STAM = ibooker.book1D("MuonPt_STAM", "Pt of STA muons", 50, 0., 100.);
  h_mEta_STAM = ibooker.book1D("MuonEta_STAM", "Eta of STA muons", 60, -3., 3.);
  h_mPhi_STAM = ibooker.book1D("MuonPhi_STAM", "Phi of STA muons", 70, -3.5, 3.5);
  h_eCombIso = ibooker.book1D("EleCombIso", "CombIso of electrons", 100, 0., 10.);
  h_mCombIso = ibooker.book1D("MuonCombIso", "CombIso of muons", 100, 0., 10.);
  h_dimumass_GMGM = ibooker.book1D("DimuMass_GMGM", "Invariant mass of GMGM pairs", 100, 0., 200.);
  h_dimumass_GMTM = ibooker.book1D("DimuMass_GMTM", "Invariant mass of GMTM pairs", 100, 0., 200.);
  h_dimumass_TMTM = ibooker.book1D("DimuMass_TMTM", "Invariant mass of TMTM pairs", 100, 0., 200.);
  h_dielemass = ibooker.book1D("DieleMass", "Invariant mass of EE pairs", 100, 0., 200.);
  h_lepcounts = ibooker.book1D("LeptonCounts", "LeptonCounts for multi lepton events", 49, 0., 49.);

  ibooker.cd();
}

void HiggsDQM::bookHistos ( DQMStore *  bei )

private

double HiggsDQM::calcDeltaPhi ( double  phi1, double  phi2  )

private

Definition at line 105 of file HiggsDQM.cc.

References srCondWrite_cfg::deltaPhi.

                                                      {
  double deltaPhi = phi1 - phi2;
  if (deltaPhi < 0)
    deltaPhi = -deltaPhi;
  if (deltaPhi > 3.1415926) {
    deltaPhi = 2 * 3.1415926 - deltaPhi;
  }
  return deltaPhi;
}

double HiggsDQM::Distance ( const reco::Candidate &  c1, const reco::Candidate &  c2  )

private

Definition at line 98 of file HiggsDQM.cc.

References HLT_FULL_cff::deltaR.

{ return deltaR(c1, c2); }

double HiggsDQM::DistancePhi ( const reco::Candidate &  c1, const reco::Candidate &  c2  )

private

Definition at line 100 of file HiggsDQM.cc.

References srCondWrite_cfg::deltaPhi, and reco::Candidate::p4().

                                                                             {
  return deltaPhi(c1.p4().phi(), c2.p4().phi());
}

Member Data Documentation

MonitorElement* HiggsDQM::h_caloMet

private

Definition at line 88 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_caloMet_phi

private

Definition at line 89 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dielemass

private

Definition at line 118 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMGM

private

Definition at line 115 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMTM

private

Definition at line 116 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_TMTM

private

Definition at line 117 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eCombIso

private

Definition at line 113 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eEta

private

Definition at line 96 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eMultiplicity

private

Definition at line 93 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePhi

private

Definition at line 97 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePt

private

Definition at line 95 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet2_et

private

Definition at line 86 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_count

private

Definition at line 87 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_et

private

Definition at line 85 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_lepcounts

private

Definition at line 119 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mCombIso

private

Definition at line 114 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GM

private

Definition at line 105 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMPT

private

Definition at line 102 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMTM

private

Definition at line 99 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_STAM

private

Definition at line 111 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_TM

private

Definition at line 108 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mMultiplicity

private

Definition at line 94 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GM

private

Definition at line 106 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMPT

private

Definition at line 103 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMTM

private

Definition at line 100 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_STAM

private

Definition at line 112 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_TM

private

Definition at line 109 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GM

private

Definition at line 104 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMPT

private

Definition at line 101 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMTM

private

Definition at line 98 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_STAM

private

Definition at line 110 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_TM

private

Definition at line 107 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet

private

Definition at line 90 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet_phi

private

Definition at line 91 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_chi2

private

Definition at line 81 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_d0

private

Definition at line 82 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_number

private

Definition at line 80 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_numTrks

private

Definition at line 83 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_sumTrks

private

Definition at line 84 of file HiggsDQM.h.

int HiggsDQM::ievt

private

Definition at line 59 of file HiggsDQM.h.

int HiggsDQM::irun

private

Definition at line 59 of file HiggsDQM.h.

int HiggsDQM::leptonflavor

private

Definition at line 61 of file HiggsDQM.h.

reco::CandidateCollection* HiggsDQM::leptonscands_

private

Definition at line 60 of file HiggsDQM.h.

unsigned long long HiggsDQM::m_cacheID_

private

Definition at line 58 of file HiggsDQM.h.

int HiggsDQM::nElectron

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nEvents_

private

Definition at line 59 of file HiggsDQM.h.

int HiggsDQM::nfourlept

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nHiggs

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nLepton

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolEle

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolMu

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nMuon

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nZEE

private

Definition at line 92 of file HiggsDQM.h.

int HiggsDQM::nZMuMu

private

Definition at line 92 of file HiggsDQM.h.

float HiggsDQM::pi

private

Definition at line 62 of file HiggsDQM.h.

double HiggsDQM::ptThrMu1_

private

Definition at line 76 of file HiggsDQM.h.

double HiggsDQM::ptThrMu2_

private

Definition at line 77 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theCaloJetCollectionLabel_

private

Definition at line 68 of file HiggsDQM.h.

edm::EDGetTokenT<reco::CaloJetCollection> HiggsDQM::theCaloJetCollectionToken_

private

Definition at line 73 of file HiggsDQM.h.

edm::EDGetTokenT<reco::CaloMETCollection> HiggsDQM::theCaloMETCollectionToken_

private

Definition at line 74 of file HiggsDQM.h.

std::string HiggsDQM::theElecTriggerPathToPass

private

Definition at line 65 of file HiggsDQM.h.

edm::EDGetTokenT<reco::GsfElectronCollection> HiggsDQM::theElectronCollectionToken_

private

Definition at line 72 of file HiggsDQM.h.

edm::EDGetTokenT<reco::MuonCollection> HiggsDQM::theMuonCollectionToken_

private

Definition at line 71 of file HiggsDQM.h.

std::string HiggsDQM::theMuonTriggerPathToPass

private

Definition at line 66 of file HiggsDQM.h.

edm::EDGetTokenT<reco::PFMETCollection> HiggsDQM::thePfMETCollectionToken_

private

Definition at line 75 of file HiggsDQM.h.

edm::EDGetTokenT<edm::TriggerResults> HiggsDQM::theTriggerResultsCollection_

private

Definition at line 69 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theTriggerResultsCollectionTag_

private

Definition at line 67 of file HiggsDQM.h.

edm::EDGetTokenT<reco::VertexCollection> HiggsDQM::vertexToken_

private

Definition at line 70 of file HiggsDQM.h.