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::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, callgraph::m2, 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_FULL_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.