#include <HiggsDQM.h>

Inheritance diagram for HiggsDQM:

Public Member Functions
	HiggsDQM (const edm::ParameterSet &ps)

virtual	~HiggsDQM ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

Public Member Functions inherited from edm::stream::EDAnalyzerBase
	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup)

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)

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

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)

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

Protected Member Functions inherited from edm::stream::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

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 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	nLumiSecs_

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 edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
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

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Detailed Description

Definition at line 45 of file HiggsDQM.h.

Constructor & Destructor Documentation

HiggsDQM::HiggsDQM ( const edm::ParameterSet & ps )

Definition at line 124 of file HiggsDQM.cc.

References edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_25ns14e33_v1_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");
 
   nLumiSecs_ = 0;
   nEvents_ = 0;
   pi = 3.14159265;
   // cout<<"...leaving  HiggsDQM::HiggsDQM. "<<endl;
 }

HiggsDQM::~HiggsDQM ( )

virtual

Definition at line 165 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
	)

protectedvirtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 261 of file HiggsDQM.cc.

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

                                                                      {
   // cout<<"[HiggsDQM::analyze()] "<<endl;
 
   //-------------------------------
   //--- Trigger Info
   //-------------------------------
   // Did it pass certain HLT path?
   bool passed_electron_HLT = true;
   bool passed_muon_HLT = true;
 
   //-------------------------------
   //--- Vertex Info
   //-------------------------------
   Handle<VertexCollection> vertexHandle;
   e.getByToken(vertexToken_, vertexHandle);
   if (vertexHandle.isValid()) {
     VertexCollection vertexCollection = *(vertexHandle.product());
     int vertex_number = vertexCollection.size();
     VertexCollection::const_iterator v = vertexCollection.begin();
     double vertex_chi2 = v->normalizedChi2();  // v->chi2();
     double vertex_d0 = sqrt(v->x() * v->x() + v->y() * v->y());
     // double vertex_ndof    = v->ndof();cout << "ndof="<<vertex_ndof<<endl;
     double vertex_numTrks = v->tracksSize();
     double vertex_sumTrks = 0.0;
     for (Vertex::trackRef_iterator vertex_curTrack = v->tracks_begin();
          vertex_curTrack != v->tracks_end(); vertex_curTrack++) {
       vertex_sumTrks += (*vertex_curTrack)->pt();
     }
     h_vertex_number->Fill(vertex_number);
     h_vertex_chi2->Fill(vertex_chi2);
     h_vertex_d0->Fill(vertex_d0);
     h_vertex_numTrks->Fill(vertex_numTrks);
     h_vertex_sumTrks->Fill(vertex_sumTrks);
   }
 
   //-------------------------------
   //--- Electrons
   //-------------------------------
   float nEle = 0;
   Handle<GsfElectronCollection> electronCollection;
   e.getByToken(theElectronCollectionToken_, electronCollection);
   if (electronCollection.isValid()) {
     int posEle = 0, negEle = 0;
     // If it passed electron HLT and the collection was found, find electrons
     // near Z mass
     if (passed_electron_HLT) {
       for (reco::GsfElectronCollection::const_iterator recoElectron =
                electronCollection->begin();
            recoElectron != electronCollection->end(); recoElectron++) {
         //      cout << "Electron with pt= " <<  recoElectron->pt() << " and
         // eta" << recoElectron->eta() << " p=" <<  recoElectron->p() << endl;
         h_ePt->Fill(recoElectron->pt());
         h_eEta->Fill(recoElectron->eta());
         h_ePhi->Fill(recoElectron->phi());
         if (recoElectron->charge() == 1) {
           posEle++;
         } else if (recoElectron->charge() == -1) {
           negEle++;
         }
         // Require electron to pass some basic cuts
         // if ( recoElectron->et() < 20 || fabs(recoElectron->eta())>2.5 )
         // continue;
         // Tighter electron cuts
         // if ( recoElectron->deltaPhiSuperClusterTrackAtVtx() > 0.58 ||
         //     recoElectron->deltaEtaSuperClusterTrackAtVtx() > 0.01 ||
         //     recoElectron->sigmaIetaIeta() > 0.027 ) continue;
       }  // end of loop over electrons
     }    // end if passed HLT
     nEle = posEle + negEle;
     if (nEle > 9.) nEle = 9.;
     h_eMultiplicity->Fill(nEle);
 
     // Z->ee:
     unsigned int eleCollectionSize = electronCollection->size();
     for (unsigned int i = 0; i < eleCollectionSize; i++) {
       const GsfElectron& ele = electronCollection->at(i);
       double pt = ele.pt();
       if (pt > ptThrMu1_) {
         for (unsigned int j = i + 1; j < eleCollectionSize; j++) {
           const GsfElectron& ele2 = electronCollection->at(j);
           double pt2 = ele2.pt();
           if (pt2 > ptThrMu2_) {
             const math::XYZTLorentzVector ZRecoEE(
                 ele.px() + ele2.px(), ele.py() + ele2.py(),
                 ele.pz() + ele2.pz(), ele.p() + ele2.p());
             h_dielemass->Fill(ZRecoEE.mass());
           }
         }
       }
     }
   }
 
   //-------------------------------
   //--- Muons
   //-------------------------------
   float nMu = 0;
   Handle<MuonCollection> muonCollection;
   e.getByToken(theMuonCollectionToken_, muonCollection);
   if (muonCollection.isValid()) {
     // Find the highest pt muons
     int posMu = 0, negMu = 0;
     TLorentzVector m1, m2;
     if (passed_muon_HLT) {
       for (reco::MuonCollection::const_iterator recoMuon =
                muonCollection->begin();
            recoMuon != muonCollection->end(); recoMuon++) {
         // cout << "Muon with pt= " <<  muIter->pt() << " and eta" <<
         // muIter->eta() << " p=" <<  muIter->p() << endl;
         if (recoMuon->isGlobalMuon() && recoMuon->isTrackerMuon()) {
           h_mPt_GMTM->Fill(recoMuon->pt());
           h_mEta_GMTM->Fill(recoMuon->eta());
           h_mPhi_GMTM->Fill(recoMuon->phi());
         } else if (recoMuon->isGlobalMuon() &&
                    (muon::isGoodMuon((*recoMuon),
                                      muon::GlobalMuonPromptTight))) {
           h_mPt_GMPT->Fill(recoMuon->pt());
           h_mEta_GMPT->Fill(recoMuon->eta());
           h_mPhi_GMPT->Fill(recoMuon->phi());
         } else if (recoMuon->isGlobalMuon()) {
           h_mPt_GM->Fill(recoMuon->pt());
           h_mEta_GM->Fill(recoMuon->eta());
           h_mPhi_GM->Fill(recoMuon->phi());
         } else if (recoMuon->isTrackerMuon() &&
                    (muon::segmentCompatibility(
                        (*recoMuon), reco::Muon::SegmentAndTrackArbitration))) {
           h_mPt_TM->Fill(recoMuon->pt());
           h_mEta_TM->Fill(recoMuon->eta());
           h_mPhi_TM->Fill(recoMuon->phi());
         } else if (recoMuon->isStandAloneMuon()) {
           h_mPt_STAM->Fill(recoMuon->pt());
           h_mEta_STAM->Fill(recoMuon->eta());
           h_mPhi_STAM->Fill(recoMuon->phi());
         }
         if (recoMuon->charge() == 1) {
           posMu++;
         } else if (recoMuon->charge() == -1) {
           negMu++;
         }
       }
       nMu = posMu + negMu;
       if (nMu > 9.) nMu = 9.;
       h_mMultiplicity->Fill(nMu);
     }
 
     // Z->mumu:
     unsigned int muonCollectionSize = muonCollection->size();
     for (unsigned int i = 0; i < muonCollectionSize; i++) {
       const Muon& mu = muonCollection->at(i);
       // if (!mu.isGlobalMuon()) continue;
       double pt = mu.pt();
       if (pt > ptThrMu1_) {
         for (unsigned int j = i + 1; j < muonCollectionSize; j++) {
           const Muon& mu2 = muonCollection->at(j);
           double pt2 = mu2.pt();
           if (pt2 > ptThrMu2_) {
             // Glb + Glb
             if (mu.isGlobalMuon() && mu2.isGlobalMuon()) {
               const math::XYZTLorentzVector ZRecoGMGM(
                   mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(),
                   mu.p() + mu2.p());
               h_dimumass_GMGM->Fill(ZRecoGMGM.mass());
             }
             // Glb + TM
             else if (mu.isGlobalMuon() && mu2.isTrackerMuon()) {
               const math::XYZTLorentzVector ZRecoGMTM(
                   mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(),
                   mu.p() + mu2.p());
               h_dimumass_GMTM->Fill(ZRecoGMTM.mass());
             }
             // TM + TM
             else if (mu.isTrackerMuon() && mu2.isTrackerMuon()) {
               const math::XYZTLorentzVector ZRecoTMTM(
                   mu.px() + mu2.px(), mu.py() + mu2.py(), mu.pz() + mu2.pz(),
                   mu.p() + mu2.p());
               h_dimumass_TMTM->Fill(ZRecoTMTM.mass());
             }
           }
         }
       }
     }
   }
 
   //-------------------------------
   //--- Jets
   //-------------------------------
   Handle<CaloJetCollection> caloJetCollection;
   e.getByToken(theCaloJetCollectionToken_, caloJetCollection);
   if (caloJetCollection.isValid()) {
     float jet_et = -8.0;
     //    float jet_eta   = -8.0; // UNUSED
     //    float jet_phi   = -8.0; // UNUSED
     int jet_count = 0;
     float jet2_et = -9.0;
     //    float jet2_eta  = -9.0; // UNUSED
     //    float jet2_phi  = -9.0; // UNUSED
     for (CaloJetCollection::const_iterator i_calojet =
              caloJetCollection->begin();
          i_calojet != caloJetCollection->end(); i_calojet++) {
       float jet_current_et = i_calojet->et();
       // if it overlaps with electron, it is not a jet
       // if ( electron_et>0.0 && fabs(i_calojet->eta()-electron_eta ) < 0.2 &&
       // calcDeltaPhi(i_calojet->phi(), electron_phi ) < 0.2) continue;
       // if ( electron2_et>0.0&& fabs(i_calojet->eta()-electron2_eta) < 0.2 &&
       // calcDeltaPhi(i_calojet->phi(), electron2_phi) < 0.2) continue;
       // if it has too low Et, throw away
       if (jet_current_et < 15) continue;
       jet_count++;
       if (jet_current_et > jet_et) {
         jet2_et = jet_et;  // 2nd highest jet get's et from current highest
                            //        jet2_eta = jet_eta; // UNUSED
                            //        jet2_phi = jet_phi; // UNUSED
         jet_et = i_calojet->et();  // current highest jet gets et from the new
                                    // highest
         //        jet_eta  = i_calojet->eta(); // UNUSED
         //        jet_phi  = i_calojet->phi(); // UNUSED
       } else if (jet_current_et > jet2_et) {
         jet2_et = i_calojet->et();
         //        jet2_eta = i_calojet->eta(); // UNUSED
         //        jet2_phi = i_calojet->phi(); // UNUSED
       }
     }
     if (jet_et > 0.0) {
       h_jet_et->Fill(jet_et);
       h_jet_count->Fill(jet_count);
     }
   }
 
   //-------------------------------
   //--- MET
   //-------------------------------
   Handle<CaloMETCollection> caloMETCollection;
   e.getByToken(theCaloMETCollectionToken_, caloMETCollection);
   if (caloMETCollection.isValid()) {
     float caloMet = caloMETCollection->begin()->et();
     float caloMet_phi = caloMETCollection->begin()->phi();
     h_caloMet->Fill(caloMet);
     h_caloMet_phi->Fill(caloMet_phi);
   }
   Handle<PFMETCollection> pfMETCollection;
   e.getByToken(thePfMETCollectionToken_, pfMETCollection);
   if (pfMETCollection.isValid()) {
     float pfMet = pfMETCollection->begin()->et();
     float pfMet_phi = pfMETCollection->begin()->phi();
     h_pfMet->Fill(pfMet);
     h_pfMet_phi->Fill(pfMet_phi);
   }
 
   //-------------------------------------
   //--- Events with more than 2 leptons:
   //-------------------------------------
   if (nMu + nEle > 2 && nMu + nEle < 10) {
     if (nMu == 0 && nEle == 3) h_lepcounts->Fill(0);
     if (nMu == 0 && nEle == 4) h_lepcounts->Fill(1);
     if (nMu == 0 && nEle == 5) h_lepcounts->Fill(2);
     if (nMu == 0 && nEle == 6) h_lepcounts->Fill(3);
     if (nMu == 0 && nEle == 7) h_lepcounts->Fill(4);
     if (nMu == 0 && nEle == 8) h_lepcounts->Fill(5);
     if (nMu == 0 && nEle == 9) h_lepcounts->Fill(6);
     if (nMu == 1 && nEle == 2) h_lepcounts->Fill(7);
     if (nMu == 1 && nEle == 3) h_lepcounts->Fill(8);
     if (nMu == 1 && nEle == 4) h_lepcounts->Fill(9);
     if (nMu == 1 && nEle == 5) h_lepcounts->Fill(10);
     if (nMu == 1 && nEle == 6) h_lepcounts->Fill(11);
     if (nMu == 1 && nEle == 7) h_lepcounts->Fill(12);
     if (nMu == 1 && nEle == 8) h_lepcounts->Fill(13);
     if (nMu == 2 && nEle == 1) h_lepcounts->Fill(14);
     if (nMu == 2 && nEle == 2) h_lepcounts->Fill(15);
     if (nMu == 2 && nEle == 3) h_lepcounts->Fill(16);
     if (nMu == 2 && nEle == 4) h_lepcounts->Fill(17);
     if (nMu == 2 && nEle == 5) h_lepcounts->Fill(18);
     if (nMu == 2 && nEle == 6) h_lepcounts->Fill(19);
     if (nMu == 2 && nEle == 7) h_lepcounts->Fill(20);
     if (nMu == 3 && nEle == 0) h_lepcounts->Fill(21);
     if (nMu == 3 && nEle == 1) h_lepcounts->Fill(22);
     if (nMu == 3 && nEle == 2) h_lepcounts->Fill(23);
     if (nMu == 3 && nEle == 3) h_lepcounts->Fill(24);
     if (nMu == 3 && nEle == 4) h_lepcounts->Fill(25);
     if (nMu == 3 && nEle == 5) h_lepcounts->Fill(26);
     if (nMu == 3 && nEle == 6) h_lepcounts->Fill(27);
     if (nMu == 4 && nEle == 0) h_lepcounts->Fill(28);
     if (nMu == 4 && nEle == 1) h_lepcounts->Fill(29);
     if (nMu == 4 && nEle == 2) h_lepcounts->Fill(30);
     if (nMu == 4 && nEle == 3) h_lepcounts->Fill(31);
     if (nMu == 4 && nEle == 4) h_lepcounts->Fill(32);
     if (nMu == 4 && nEle == 5) h_lepcounts->Fill(33);
     if (nMu == 5 && nEle == 0) h_lepcounts->Fill(34);
     if (nMu == 5 && nEle == 1) h_lepcounts->Fill(35);
     if (nMu == 5 && nEle == 2) h_lepcounts->Fill(36);
     if (nMu == 5 && nEle == 3) h_lepcounts->Fill(37);
     if (nMu == 5 && nEle == 4) h_lepcounts->Fill(38);
     if (nMu == 6 && nEle == 0) h_lepcounts->Fill(39);
     if (nMu == 6 && nEle == 1) h_lepcounts->Fill(40);
     if (nMu == 6 && nEle == 2) h_lepcounts->Fill(41);
     if (nMu == 6 && nEle == 3) h_lepcounts->Fill(42);
     if (nMu == 7 && nEle == 0) h_lepcounts->Fill(43);
     if (nMu == 7 && nEle == 1) h_lepcounts->Fill(44);
     if (nMu == 7 && nEle == 2) h_lepcounts->Fill(45);
     if (nMu == 8 && nEle == 0) h_lepcounts->Fill(46);
     if (nMu == 8 && nEle == 1) h_lepcounts->Fill(47);
     if (nMu == 9 && nEle == 0) h_lepcounts->Fill(48);
   }
   if ((nMu + nEle) >= 10)
     LogDebug("HiggsDQM") << "WARNING: " << nMu + nEle
                          << " leptons in this event: run=" << e.id().run()
                          << ", event=" << e.id().event() << "\n";
 }

void HiggsDQM::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 177 of file HiggsDQM.cc.

                                                                   {
   // cout<<"Entering HiggsDQM::beginLuminosityBlock: "<<endl;
 
   edm::LogInfo("HiggsDQM") << "[HiggsDQM]: Begin of LS transition";
 
   // cout<<"...leaving HiggsDQM::beginLuminosityBlock. "<<endl;
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 188 of file HiggsDQM.cc.

References DQMStore::IBooker::book1D(), DQMStore::IBooker::cd(), and 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 112 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 101 of file HiggsDQM.cc.

References deltaR().

                                                    {
   return deltaR(c1, c2);
 }

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

private

Definition at line 106 of file HiggsDQM.cc.

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

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

void HiggsDQM::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	c
	)

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 570 of file HiggsDQM.cc.

References edm::LuminosityBlockBase::luminosityBlock().

                                                                {
   //  cout<<"Entering HiggsDQM::endLuminosityBlock: "<<endl;
 
   edm::LogInfo("HiggsDQM") << "[HiggsDQM]: End of LS transition, performing "
                               "the DQM client operation";
 
   nLumiSecs_++;
   // cout << "nLumiSecs_: "<< nLumiSecs_ << endl;
 
   edm::LogInfo("HiggsDQM")
       << "====================================================== " << endl
       << " ===> Iteration # " << nLumiSecs_ << " " << lumiSeg.luminosityBlock()
       << endl
       << "====================================================== " << endl;
 
   //  cout<<"...leaving HiggsDQM::endLuminosityBlock. "<<endl;
 }

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

protectedvirtual

Reimplemented from edm::stream::EDAnalyzerBase.

Definition at line 591 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 100 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_caloMet_phi

private

Definition at line 101 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dielemass

private

Definition at line 131 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMGM

private

Definition at line 128 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_GMTM

private

Definition at line 129 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_dimumass_TMTM

private

Definition at line 130 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eCombIso

private

Definition at line 126 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eEta

private

Definition at line 109 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_eMultiplicity

private

Definition at line 106 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePhi

private

Definition at line 110 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_ePt

private

Definition at line 108 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet2_et

private

Definition at line 98 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_count

private

Definition at line 99 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_jet_et

private

Definition at line 97 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_lepcounts

private

Definition at line 132 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mCombIso

private

Definition at line 127 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GM

private

Definition at line 118 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMPT

private

Definition at line 115 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_GMTM

private

Definition at line 112 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_STAM

private

Definition at line 124 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mEta_TM

private

Definition at line 121 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mMultiplicity

private

Definition at line 107 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GM

private

Definition at line 119 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMPT

private

Definition at line 116 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_GMTM

private

Definition at line 113 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_STAM

private

Definition at line 125 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPhi_TM

private

Definition at line 122 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GM

private

Definition at line 117 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMPT

private

Definition at line 114 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_GMTM

private

Definition at line 111 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_STAM

private

Definition at line 123 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_mPt_TM

private

Definition at line 120 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet

private

Definition at line 102 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_pfMet_phi

private

Definition at line 103 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_chi2

private

Definition at line 93 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_d0

private

Definition at line 94 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_number

private

Definition at line 92 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_numTrks

private

Definition at line 95 of file HiggsDQM.h.

MonitorElement* HiggsDQM::h_vertex_sumTrks

private

Definition at line 96 of file HiggsDQM.h.

int HiggsDQM::ievt

private

Definition at line 70 of file HiggsDQM.h.

int HiggsDQM::irun

private

Definition at line 70 of file HiggsDQM.h.

int HiggsDQM::leptonflavor

private

Definition at line 72 of file HiggsDQM.h.

reco::CandidateCollection* HiggsDQM::leptonscands_

private

Definition at line 71 of file HiggsDQM.h.

unsigned long long HiggsDQM::m_cacheID_

private

Definition at line 68 of file HiggsDQM.h.

int HiggsDQM::nElectron

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nEvents_

private

Definition at line 70 of file HiggsDQM.h.

int HiggsDQM::nfourlept

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nHiggs

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nLepton

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolEle

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nLooseIsolMu

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nLumiSecs_

private

Definition at line 69 of file HiggsDQM.h.

int HiggsDQM::nMuon

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nZEE

private

Definition at line 104 of file HiggsDQM.h.

int HiggsDQM::nZMuMu

private

Definition at line 104 of file HiggsDQM.h.

float HiggsDQM::pi

private

Definition at line 73 of file HiggsDQM.h.

double HiggsDQM::ptThrMu1_

private

Definition at line 88 of file HiggsDQM.h.

double HiggsDQM::ptThrMu2_

private

Definition at line 89 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theCaloJetCollectionLabel_

private

Definition at line 79 of file HiggsDQM.h.

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

private

Definition at line 85 of file HiggsDQM.h.

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

private

Definition at line 86 of file HiggsDQM.h.

std::string HiggsDQM::theElecTriggerPathToPass

private

Definition at line 76 of file HiggsDQM.h.

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

private

Definition at line 84 of file HiggsDQM.h.

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

private

Definition at line 83 of file HiggsDQM.h.

std::string HiggsDQM::theMuonTriggerPathToPass

private

Definition at line 77 of file HiggsDQM.h.

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

private

Definition at line 87 of file HiggsDQM.h.

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

private

Definition at line 80 of file HiggsDQM.h.

edm::InputTag HiggsDQM::theTriggerResultsCollectionTag_

private

Definition at line 78 of file HiggsDQM.h.

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

private

Definition at line 82 of file HiggsDQM.h.

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation