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
virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)
void	endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) final
void	endLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &, edm::EventSetup const &) final
void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndRuns () const final
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
	~ProducerBase () noexcept(false) override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

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
void	dqmEndRun (edm::Run const &run, edm::EventSetup const &eSetup) override
Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

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::one::EDProducerBase
typedef EDProducerBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_

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_2018_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 pdwgLeptonRecoSkim_cfi::caloJetCollection, heavyionUCCDQM_cfi::caloMet, ewkDQM_cfi::caloMETCollection, pdwgLeptonRecoSkim_cfi::electronCollection, 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, 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(), edm::EventID::run(), reco::Muon::SegmentAndTrackArbitration, muon::segmentCompatibility(), mathSSE::sqrt(), findQualityFiles::v, and spclusmultinvestigator_cfi::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::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::cd(), and dqm::dqmstoreimpl::DQMStore::IBooker::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 SiPixelRawToDigiRegional_cfi::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 PbPb_ZMuSkimMuonDPG_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 SiPixelRawToDigiRegional_cfi::deltaPhi, and reco::Candidate::p4().

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

void HiggsDQM::dqmEndRun ( edm::Run const &  run, edm::EventSetup const &  eSetup  )

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 576 of file HiggsDQM.cc.

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

  // edm::LogVerbatim ("HiggsDQM") <<"[HiggsDQM]: End of Run, saving  DQM output
  // ";
  // int iRun = run.run();

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

Member Data Documentation

MonitorElement* HiggsDQM::h_caloMet

private

Definition at line 89 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_caloMet_phi

private

Definition at line 90 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dielemass

private

Definition at line 119 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMGM

private

Definition at line 116 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMTM

private

Definition at line 117 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_TMTM

private

Definition at line 118 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eCombIso

private

Definition at line 114 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eEta

private

Definition at line 97 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eMultiplicity

private

Definition at line 94 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePhi

private

Definition at line 98 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePt

private

Definition at line 96 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet2_et

private

Definition at line 87 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_count

private

Definition at line 88 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_et

private

Definition at line 86 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_lepcounts

private

Definition at line 120 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mCombIso

private

Definition at line 115 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GM

private

Definition at line 106 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMPT

private

Definition at line 103 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMTM

private

Definition at line 100 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_STAM

private

Definition at line 112 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_TM

private

Definition at line 109 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mMultiplicity

private

Definition at line 95 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GM

private

Definition at line 107 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMPT

private

Definition at line 104 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMTM

private

Definition at line 101 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_STAM

private

Definition at line 113 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_TM

private

Definition at line 110 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GM

private

Definition at line 105 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMPT

private

Definition at line 102 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMTM

private

Definition at line 99 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_STAM

private

Definition at line 111 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_TM

private

Definition at line 108 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet

private

Definition at line 91 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet_phi

private

Definition at line 92 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_chi2

private

Definition at line 82 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_d0

private

Definition at line 83 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_number

private

Definition at line 81 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_numTrks

private

Definition at line 84 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_sumTrks

private

Definition at line 85 of file HiggsDQM.h.

int HiggsDQM::ievt

private

Definition at line 60 of file HiggsDQM.h.

int HiggsDQM::irun

private

Definition at line 60 of file HiggsDQM.h.

int HiggsDQM::leptonflavor

private

Definition at line 62 of file HiggsDQM.h.

reco::CandidateCollection* HiggsDQM::leptonscands_

private

Definition at line 61 of file HiggsDQM.h.

unsigned long long HiggsDQM::m_cacheID_

private

Definition at line 59 of file HiggsDQM.h.

int HiggsDQM::nElectron

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nEvents_

private

Definition at line 60 of file HiggsDQM.h.

int HiggsDQM::nfourlept

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nHiggs

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nLepton

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolEle

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolMu

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nMuon

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nZEE

private

Definition at line 93 of file HiggsDQM.h.

int HiggsDQM::nZMuMu

private

Definition at line 93 of file HiggsDQM.h.

float HiggsDQM::pi

private

Definition at line 63 of file HiggsDQM.h.

double HiggsDQM::ptThrMu1_

private

Definition at line 77 of file HiggsDQM.h.

double HiggsDQM::ptThrMu2_

private

Definition at line 78 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theCaloJetCollectionLabel_

private

Definition at line 69 of file HiggsDQM.h.

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

private

Definition at line 74 of file HiggsDQM.h.

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

private

Definition at line 75 of file HiggsDQM.h.

std::string HiggsDQM::theElecTriggerPathToPass

private

Definition at line 66 of file HiggsDQM.h.

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

private

Definition at line 73 of file HiggsDQM.h.

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

private

Definition at line 72 of file HiggsDQM.h.

std::string HiggsDQM::theMuonTriggerPathToPass

private

Definition at line 67 of file HiggsDQM.h.

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

private

Definition at line 76 of file HiggsDQM.h.

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

private

Definition at line 70 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theTriggerResultsCollectionTag_

private

Definition at line 68 of file HiggsDQM.h.

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

private

Definition at line 71 of file HiggsDQM.h.