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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final

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 >
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	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::HiggsDQM

(

const edm::ParameterSet &

ps

)

Definition at line 118 of file HiggsDQM.cc.

                                            {
  // 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;
}

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

~HiggsDQM()

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

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 224 of file HiggsDQM.cc.

                                                                     {
  // 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";
}

References pdwgLeptonRecoSkim_cfi::caloJetCollection, heavyionUCCDQM_cfi::caloMet, ewkDQM_cfi::caloMETCollection, MillePedeFileConverter_cfg::e, pdwgLeptonRecoSkim_cfi::electronCollection, muon::GlobalMuonPromptTight, mps_fire::i, reco::Muon::isGlobalMuon(), muon::isGoodMuon(), reco::Muon::isTrackerMuon(), edm::HandleBase::isValid(), dqmiolumiharvest::j, LogDebug, amptDefaultParameters_cff::mu, pdwgLeptonRecoSkim_cfi::muonCollection, reco::LeafCandidate::p(), RecoPFMET_cff::pfMet, B2GDQM_cfi::pfMETCollection, edm::Handle< T >::product(), DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), HLT_2018_cff::pt2, reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::LeafCandidate::pz(), reco::Muon::SegmentAndTrackArbitration, muon::segmentCompatibility(), mathSSE::sqrt(), findQualityFiles::v, and spclusmultinvestigator_cfi::vertexCollection.

bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 162 of file HiggsDQM.cc.

                                                                                           {
  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();
}

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

bookHistos()

void HiggsDQM::bookHistos ( DQMStore *  bei )

private

calcDeltaPhi()

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

private

Definition at line 105 of file HiggsDQM.cc.

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

References SiPixelRawToDigiRegional_cfi::deltaPhi.

Distance()

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

private

Definition at line 98 of file HiggsDQM.cc.

{ return deltaR(c1, c2); }

References alignmentValidation::c1, and PbPb_ZMuSkimMuonDPG_cff::deltaR.

DistancePhi()

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

private

Definition at line 100 of file HiggsDQM.cc.

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

References alignmentValidation::c1, SiPixelRawToDigiRegional_cfi::deltaPhi, and reco::Candidate::p4().

Member Data Documentation

h_caloMet

MonitorElement* HiggsDQM::h_caloMet

private

Definition at line 88 of file HiggsDQM.h.

h_caloMet_phi

MonitorElement* HiggsDQM::h_caloMet_phi

private

Definition at line 89 of file HiggsDQM.h.

h_dielemass

MonitorElement* HiggsDQM::h_dielemass

private

Definition at line 118 of file HiggsDQM.h.

h_dimumass_GMGM

MonitorElement* HiggsDQM::h_dimumass_GMGM

private

Definition at line 115 of file HiggsDQM.h.

h_dimumass_GMTM

MonitorElement* HiggsDQM::h_dimumass_GMTM

private

Definition at line 116 of file HiggsDQM.h.

h_dimumass_TMTM

MonitorElement* HiggsDQM::h_dimumass_TMTM

private

Definition at line 117 of file HiggsDQM.h.

h_eCombIso

MonitorElement* HiggsDQM::h_eCombIso

private

Definition at line 113 of file HiggsDQM.h.

h_eEta

MonitorElement* HiggsDQM::h_eEta

private

Definition at line 96 of file HiggsDQM.h.

h_eMultiplicity

MonitorElement* HiggsDQM::h_eMultiplicity

private

Definition at line 93 of file HiggsDQM.h.

h_ePhi

MonitorElement* HiggsDQM::h_ePhi

private

Definition at line 97 of file HiggsDQM.h.

h_ePt

MonitorElement* HiggsDQM::h_ePt

private

Definition at line 95 of file HiggsDQM.h.

h_jet2_et

MonitorElement* HiggsDQM::h_jet2_et

private

Definition at line 86 of file HiggsDQM.h.

h_jet_count

MonitorElement* HiggsDQM::h_jet_count

private

Definition at line 87 of file HiggsDQM.h.

h_jet_et

MonitorElement* HiggsDQM::h_jet_et

private

Definition at line 85 of file HiggsDQM.h.

h_lepcounts

MonitorElement* HiggsDQM::h_lepcounts

private

Definition at line 119 of file HiggsDQM.h.

h_mCombIso

MonitorElement* HiggsDQM::h_mCombIso

private

Definition at line 114 of file HiggsDQM.h.

h_mEta_GM

MonitorElement* HiggsDQM::h_mEta_GM

private

Definition at line 105 of file HiggsDQM.h.

h_mEta_GMPT

MonitorElement* HiggsDQM::h_mEta_GMPT

private

Definition at line 102 of file HiggsDQM.h.

h_mEta_GMTM

MonitorElement* HiggsDQM::h_mEta_GMTM

private

Definition at line 99 of file HiggsDQM.h.

h_mEta_STAM

MonitorElement* HiggsDQM::h_mEta_STAM

private

Definition at line 111 of file HiggsDQM.h.

h_mEta_TM

MonitorElement* HiggsDQM::h_mEta_TM

private

Definition at line 108 of file HiggsDQM.h.

h_mMultiplicity

MonitorElement* HiggsDQM::h_mMultiplicity

private

Definition at line 94 of file HiggsDQM.h.

h_mPhi_GM

MonitorElement* HiggsDQM::h_mPhi_GM

private

Definition at line 106 of file HiggsDQM.h.

h_mPhi_GMPT

MonitorElement* HiggsDQM::h_mPhi_GMPT

private

Definition at line 103 of file HiggsDQM.h.

h_mPhi_GMTM

MonitorElement* HiggsDQM::h_mPhi_GMTM

private

Definition at line 100 of file HiggsDQM.h.

h_mPhi_STAM

MonitorElement* HiggsDQM::h_mPhi_STAM

private

Definition at line 112 of file HiggsDQM.h.

h_mPhi_TM

MonitorElement* HiggsDQM::h_mPhi_TM

private

Definition at line 109 of file HiggsDQM.h.

h_mPt_GM

MonitorElement* HiggsDQM::h_mPt_GM

private

Definition at line 104 of file HiggsDQM.h.

h_mPt_GMPT

MonitorElement* HiggsDQM::h_mPt_GMPT

private

Definition at line 101 of file HiggsDQM.h.

h_mPt_GMTM

MonitorElement* HiggsDQM::h_mPt_GMTM

private

Definition at line 98 of file HiggsDQM.h.

h_mPt_STAM

MonitorElement* HiggsDQM::h_mPt_STAM

private

Definition at line 110 of file HiggsDQM.h.

h_mPt_TM

MonitorElement* HiggsDQM::h_mPt_TM

private

Definition at line 107 of file HiggsDQM.h.

h_pfMet

MonitorElement* HiggsDQM::h_pfMet

private

Definition at line 90 of file HiggsDQM.h.

h_pfMet_phi

MonitorElement* HiggsDQM::h_pfMet_phi

private

Definition at line 91 of file HiggsDQM.h.

h_vertex_chi2

MonitorElement* HiggsDQM::h_vertex_chi2

private

Definition at line 81 of file HiggsDQM.h.

h_vertex_d0

MonitorElement* HiggsDQM::h_vertex_d0

private

Definition at line 82 of file HiggsDQM.h.

h_vertex_number

MonitorElement* HiggsDQM::h_vertex_number

private

Definition at line 80 of file HiggsDQM.h.

h_vertex_numTrks

MonitorElement* HiggsDQM::h_vertex_numTrks

private

Definition at line 83 of file HiggsDQM.h.

h_vertex_sumTrks

MonitorElement* HiggsDQM::h_vertex_sumTrks

private

Definition at line 84 of file HiggsDQM.h.

ievt

int HiggsDQM::ievt

private

Definition at line 59 of file HiggsDQM.h.

irun

int HiggsDQM::irun

private

Definition at line 59 of file HiggsDQM.h.

leptonflavor

int HiggsDQM::leptonflavor

private

Definition at line 61 of file HiggsDQM.h.

leptonscands_

reco::CandidateCollection* HiggsDQM::leptonscands_

private

Definition at line 60 of file HiggsDQM.h.

m_cacheID_

unsigned long long HiggsDQM::m_cacheID_

private

Definition at line 58 of file HiggsDQM.h.

nElectron

int HiggsDQM::nElectron

private

Definition at line 92 of file HiggsDQM.h.

nEvents_

int HiggsDQM::nEvents_

private

Definition at line 59 of file HiggsDQM.h.

nfourlept

int HiggsDQM::nfourlept

private

Definition at line 92 of file HiggsDQM.h.

nHiggs

int HiggsDQM::nHiggs

private

Definition at line 92 of file HiggsDQM.h.

nLepton

int HiggsDQM::nLepton

private

Definition at line 92 of file HiggsDQM.h.

nLooseIsolEle

int HiggsDQM::nLooseIsolEle

private

Definition at line 92 of file HiggsDQM.h.

nLooseIsolMu

int HiggsDQM::nLooseIsolMu

private

Definition at line 92 of file HiggsDQM.h.

nMuon

int HiggsDQM::nMuon

private

Definition at line 92 of file HiggsDQM.h.

nZEE

int HiggsDQM::nZEE

private

Definition at line 92 of file HiggsDQM.h.

nZMuMu

int HiggsDQM::nZMuMu

private

Definition at line 92 of file HiggsDQM.h.

pi

float HiggsDQM::pi

private

Definition at line 62 of file HiggsDQM.h.

ptThrMu1_

double HiggsDQM::ptThrMu1_

private

Definition at line 76 of file HiggsDQM.h.

ptThrMu2_

double HiggsDQM::ptThrMu2_

private

Definition at line 77 of file HiggsDQM.h.

theCaloJetCollectionLabel_

edm::InputTag HiggsDQM::theCaloJetCollectionLabel_

private

Definition at line 68 of file HiggsDQM.h.

theCaloJetCollectionToken_

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

private

Definition at line 73 of file HiggsDQM.h.

theCaloMETCollectionToken_

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

private

Definition at line 74 of file HiggsDQM.h.

theElecTriggerPathToPass

std::string HiggsDQM::theElecTriggerPathToPass

private

Definition at line 65 of file HiggsDQM.h.

theElectronCollectionToken_

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

private

Definition at line 72 of file HiggsDQM.h.

theMuonCollectionToken_

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

private

Definition at line 71 of file HiggsDQM.h.

theMuonTriggerPathToPass

std::string HiggsDQM::theMuonTriggerPathToPass

private

Definition at line 66 of file HiggsDQM.h.

thePfMETCollectionToken_

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

private

Definition at line 75 of file HiggsDQM.h.

theTriggerResultsCollection_

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

private

Definition at line 69 of file HiggsDQM.h.

theTriggerResultsCollectionTag_

edm::InputTag HiggsDQM::theTriggerResultsCollectionTag_

private

Definition at line 67 of file HiggsDQM.h.

vertexToken_

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

private

Definition at line 70 of file HiggsDQM.h.