B2GDQM Class Reference

#include <B2GDQM.h>

Inheritance diagram for B2GDQM:

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

Protected Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override
virtual void	analyzeAllHad (edm::Event const &e, edm::EventSetup const &eSetup)
virtual void	analyzeJets (edm::Event const &e, edm::EventSetup const &eSetup)
virtual void	analyzeSemiE (edm::Event const &e, edm::EventSetup const &eSetup)
virtual void	analyzeSemiMu (edm::Event const &e, edm::EventSetup const &eSetup)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistograms (DQMStore::IBooker &bei, edm::Run const &, edm::EventSetup const &) override

Private Attributes
MonitorElement *	allHad_mass0
MonitorElement *	allHad_mass1
MonitorElement *	allHad_mttbar
MonitorElement *	allHad_phi0
MonitorElement *	allHad_phi1
MonitorElement *	allHad_pt0
MonitorElement *	allHad_pt1
MonitorElement *	allHad_wMass0
MonitorElement *	allHad_wMass1
MonitorElement *	allHad_y0
MonitorElement *	allHad_y1
double	allHadDeltaPhiCut_
double	allHadPtCut_
double	allHadRapidityCut_
std::vector< MonitorElement * >	boostedJet_massDrop
std::vector< MonitorElement * >	boostedJet_subjetM
std::vector< MonitorElement * >	boostedJet_subjetN
std::vector< MonitorElement * >	boostedJet_subjetPhi
std::vector< MonitorElement * >	boostedJet_subjetPt
std::vector< MonitorElement * >	boostedJet_subjetY
std::vector< MonitorElement * >	boostedJet_wMass
std::shared_ptr< StringCutObjectSelector< reco::GsfElectron > >	elecSelect_
edm::InputTag	electronLabel_
edm::EDGetTokenT< edm::View< reco::GsfElectron > >	electronToken_
HLTConfigProvider	hltConfigProvider_
int	ievt
int	irun
bool	isValidHltConfig_
std::vector< edm::InputTag >	jetLabels_
std::vector< double >	jetPtMins_
std::vector< edm::EDGetTokenT< edm::View< reco::Jet > > >	jetTokens_
edm::InputTag	muonLabel_
std::shared_ptr< StringCutObjectSelector< reco::Muon > >	muonSelect_
edm::EDGetTokenT< edm::View< reco::Muon > >	muonToken_
int	nEvents_
int	nLumiSecs_
std::vector< MonitorElement * >	pfJet_cemf
std::vector< MonitorElement * >	pfJet_chef
std::vector< MonitorElement * >	pfJet_m
std::vector< MonitorElement * >	pfJet_nemf
std::vector< MonitorElement * >	pfJet_nhef
std::vector< MonitorElement * >	pfJet_phi
std::vector< MonitorElement * >	pfJet_pt
std::vector< MonitorElement * >	pfJet_y
std::string	PFJetCorService_
MonitorElement *	pfMet_phi
MonitorElement *	pfMet_pt
edm::InputTag	PFMETLabel_
edm::EDGetTokenT< std::vector< reco::PFMET > >	PFMETToken_
edm::InputTag	sdjetLabel_
edm::EDGetTokenT< edm::View< reco::BasicJet > >	sdjetToken_
double	semiE_dphiHadCut_
double	semiE_dRMin_
MonitorElement *	semiE_eDRMin
MonitorElement *	semiE_eEta
MonitorElement *	semiE_ePhi
MonitorElement *	semiE_ePt
MonitorElement *	semiE_ePtRel
MonitorElement *	semiE_hadJetDR
MonitorElement *	semiE_hadJetMass
MonitorElement *	semiE_hadJetPhi
MonitorElement *	semiE_hadJetPt
double	semiE_HadJetPtCut_
MonitorElement *	semiE_hadJetWMass
MonitorElement *	semiE_hadJetY
double	semiE_LepJetPtCut_
MonitorElement *	semiE_mttbar
double	semiE_ptRel_
double	semiMu_dphiHadCut_
double	semiMu_dRMin_
MonitorElement *	semiMu_hadJetDR
MonitorElement *	semiMu_hadJetMass
MonitorElement *	semiMu_hadJetPhi
MonitorElement *	semiMu_hadJetPt
double	semiMu_HadJetPtCut_
MonitorElement *	semiMu_hadJetWMass
MonitorElement *	semiMu_hadJetY
double	semiMu_LepJetPtCut_
MonitorElement *	semiMu_mttbar
MonitorElement *	semiMu_muDRMin
MonitorElement *	semiMu_muEta
MonitorElement *	semiMu_muPhi
MonitorElement *	semiMu_muPt
MonitorElement *	semiMu_muPtRel
double	semiMu_ptRel_
edm::InputTag	theTriggerResultsCollection
edm::Handle< edm::TriggerResults >	triggerResults_
edm::EDGetTokenT< edm::TriggerResults >	triggerToken_

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

Detailed Description

Definition at line 78 of file B2GDQM.h.

Constructor & Destructor Documentation

B2GDQM()

B2GDQM::B2GDQM

(

const edm::ParameterSet &

ps

)

Definition at line 77 of file B2GDQM.cc.

                                        {
  edm::LogInfo("B2GDQM") << " Starting B2GDQM "
                         << "\n";
 
  typedef std::vector<edm::InputTag> vtag;
 
  // Get parameters from configuration file
  // Trigger
  theTriggerResultsCollection = ps.getParameter<InputTag>("triggerResultsCollection");
  triggerToken_ = consumes<edm::TriggerResults>(theTriggerResultsCollection);
 
  // Jets
  jetLabels_ = ps.getParameter<std::vector<edm::InputTag> >("jetLabels");
  for (std::vector<edm::InputTag>::const_iterator jetlabel = jetLabels_.begin(), jetlabelEnd = jetLabels_.end();
       jetlabel != jetlabelEnd;
       ++jetlabel) {
    jetTokens_.push_back(consumes<edm::View<reco::Jet> >(*jetlabel));
  }
  sdjetLabel_ = ps.getParameter<edm::InputTag>("sdjetLabel");
  sdjetToken_ = consumes<edm::View<reco::BasicJet> >(sdjetLabel_);
 
  muonToken_ = consumes<edm::View<reco::Muon> >(ps.getParameter<edm::InputTag>("muonSrc"));
  electronToken_ = consumes<edm::View<reco::GsfElectron> >(ps.getParameter<edm::InputTag>("electronSrc"));
 
  jetPtMins_ = ps.getParameter<std::vector<double> >("jetPtMins");
  allHadPtCut_ = ps.getParameter<double>("allHadPtCut");
  allHadRapidityCut_ = ps.getParameter<double>("allHadRapidityCut");
  allHadDeltaPhiCut_ = ps.getParameter<double>("allHadDeltaPhiCut");
 
  semiMu_HadJetPtCut_ = ps.getParameter<double>("semiMu_HadJetPtCut");
  semiMu_LepJetPtCut_ = ps.getParameter<double>("semiMu_LepJetPtCut");
  semiMu_dphiHadCut_ = ps.getParameter<double>("semiMu_dphiHadCut");
  semiMu_dRMin_ = ps.getParameter<double>("semiMu_dRMin");
  semiMu_ptRel_ = ps.getParameter<double>("semiMu_ptRel");
  muonSelect_ = std::make_shared<StringCutObjectSelector<reco::Muon> >(
 
      ps.getParameter<std::string>("muonSelect"));
 
  semiE_HadJetPtCut_ = ps.getParameter<double>("semiE_HadJetPtCut");
  semiE_LepJetPtCut_ = ps.getParameter<double>("semiE_LepJetPtCut");
  semiE_dphiHadCut_ = ps.getParameter<double>("semiE_dphiHadCut");
  semiE_dRMin_ = ps.getParameter<double>("semiE_dRMin");
  semiE_ptRel_ = ps.getParameter<double>("semiE_ptRel");
  elecSelect_ = std::make_shared<StringCutObjectSelector<reco::GsfElectron> >(
 
      ps.getParameter<std::string>("elecSelect"));
 
  PFJetCorService_ = ps.getParameter<std::string>("PFJetCorService");
 
  // MET
  PFMETLabel_ = ps.getParameter<InputTag>("pfMETCollection");
  PFMETToken_ = consumes<std::vector<reco::PFMET> >(PFMETLabel_);
}

References edm::ParameterSet::getParameter(), and AlCaHLTBitMon_QueryRunRegistry::string.

~B2GDQM()

B2GDQM::~B2GDQM ( )

override

Definition at line 134 of file B2GDQM.cc.

                {
  edm::LogInfo("B2GDQM") << " Deleting B2GDQM "
                         << "\n";
}

Member Function Documentation

analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 229 of file B2GDQM.cc.

                                                                        {
  analyzeJets(iEvent, iSetup);
  analyzeSemiMu(iEvent, iSetup);
  analyzeSemiE(iEvent, iSetup);
  analyzeAllHad(iEvent, iSetup);
}

References iEvent.

analyzeAllHad()

void B2GDQM::analyzeAllHad ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

protectedvirtual

Definition at line 321 of file B2GDQM.cc.

                                                                           {
  edm::Handle<edm::View<reco::BasicJet> > jetCollection;
  bool validJets = iEvent.getByToken(sdjetToken_, jetCollection);
  if (!validJets)
    return;
 
  // Require two back-to-back jets at high pt with |delta y| < 1.0
  if (jetCollection->size() < 2)
    return;
  edm::Ptr<reco::BasicJet> jet0 = jetCollection->ptrAt(0);
  edm::Ptr<reco::BasicJet> jet1 = jetCollection->ptrAt(1);
  if (jet0.isAvailable() == false || jet1.isAvailable() == false)
    return;
  if (jet0->pt() < allHadPtCut_ || jet1->pt() < allHadPtCut_)
    return;
  if (std::abs(jet0->rapidity() - jet1->rapidity()) > allHadRapidityCut_)
    return;
  if (std::abs(reco::deltaPhi(jet0->phi(), jet1->phi())) < M_PI * 0.5)
    return;
 
  allHad_pt0->Fill(jet0->pt());
  allHad_y0->Fill(jet0->rapidity());
  allHad_phi0->Fill(jet0->phi());
  allHad_mass0->Fill(jet0->mass());
  if (jet0->numberOfDaughters() > 2) {
    double wMass =
        jet0->daughter(0)->mass() >= jet0->daughter(1)->mass() ? jet0->daughter(0)->mass() : jet0->daughter(1)->mass();
    allHad_wMass0->Fill(wMass);
  } else {
    allHad_wMass0->Fill(-1.0);
  }
 
  allHad_pt1->Fill(jet1->pt());
  allHad_y1->Fill(jet1->rapidity());
  allHad_phi1->Fill(jet1->phi());
  allHad_mass1->Fill(jet1->mass());
  if (jet1->numberOfDaughters() > 2) {
    double wMass =
        jet1->daughter(0)->mass() >= jet1->daughter(1)->mass() ? jet1->daughter(0)->mass() : jet1->daughter(1)->mass();
    allHad_wMass1->Fill(wMass);
  } else {
    allHad_wMass1->Fill(-1.0);
  }
 
  auto p4cand = (jet0->p4() + jet1->p4());
  allHad_mttbar->Fill(p4cand.mass());
}

References funct::abs(), reco::CompositePtrCandidate::daughter(), reco::deltaPhi(), iEvent, edm::Ptr< T >::isAvailable(), jetfilter_cfi::jetCollection, M_PI, reco::Candidate::mass(), reco::LeafCandidate::mass(), reco::CompositePtrCandidate::numberOfDaughters(), reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), reco::LeafCandidate::rapidity(), and pseudoTop_cfi::wMass.

analyzeJets()

void B2GDQM::analyzeJets ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

protectedvirtual

Definition at line 236 of file B2GDQM.cc.

                                                                         {
  // Loop over the different types of jets,
  //   Loop over the jets in that collection,
  //     fill PF jet information as well as substructure
  //     information for boosted jets.
  // Utilizes the CMS top-tagging algorithm and the "mass drop" W-tagger.
  for (unsigned int icoll = 0; icoll < jetLabels_.size(); ++icoll) {
    edm::Handle<edm::View<reco::Jet> > pfJetCollection;
    bool ValidPFJets = iEvent.getByToken(jetTokens_[icoll], pfJetCollection);
    if (!ValidPFJets)
      continue;
    edm::View<reco::Jet> const& pfjets = *pfJetCollection;
 
    // Jet Correction
    int countJet = 0;
    // const JetCorrector* pfcorrector =
    // JetCorrector::getJetCorrector(PFJetCorService_,iSetup);
 
    for (edm::View<reco::Jet>::const_iterator jet = pfjets.begin(), jetEnd = pfjets.end(); jet != jetEnd; ++jet) {
      if (jet->pt() < jetPtMins_[icoll])
        continue;
      pfJet_pt[icoll]->Fill(jet->pt());
      pfJet_y[icoll]->Fill(jet->rapidity());
      pfJet_phi[icoll]->Fill(jet->phi());
      pfJet_m[icoll]->Fill(jet->mass());
 
      // Dynamic cast the base class (reco::Jet) to the derived class (PFJet)
      // to access the PFJet information
      reco::PFJet const* pfjet = dynamic_cast<reco::PFJet const*>(&*jet);
 
      if (pfjet != nullptr) {
        pfJet_chef[icoll]->Fill(pfjet->chargedHadronEnergyFraction());
        pfJet_nhef[icoll]->Fill(pfjet->neutralHadronEnergyFraction());
        pfJet_cemf[icoll]->Fill(pfjet->chargedEmEnergyFraction());
        pfJet_nemf[icoll]->Fill(pfjet->neutralEmEnergyFraction());
      }
 
      // Dynamic cast the base class (reco::Jet) to the derived class (BasicJet)
      // to access the substructure information
      reco::BasicJet const* basicjet = dynamic_cast<reco::BasicJet const*>(&*jet);
 
      if (basicjet != nullptr) {
        boostedJet_subjetN[icoll]->Fill(jet->numberOfDaughters());
 
        for (unsigned int ida = 0; ida < jet->numberOfDaughters(); ++ida) {
          reco::Candidate const* subjet = jet->daughter(ida);
          boostedJet_subjetPt[icoll]->Fill(subjet->pt());
          boostedJet_subjetY[icoll]->Fill(subjet->rapidity());
          boostedJet_subjetPhi[icoll]->Fill(subjet->phi());
          boostedJet_subjetM[icoll]->Fill(subjet->mass());
        }
        // Check the various tagging algorithms
        if ((jetLabels_[icoll].label() == "ak8PFJetsPuppiSoftdrop")) {
          if (jet->numberOfDaughters() > 1) {
            reco::Candidate const* da0 = jet->daughter(0);
            reco::Candidate const* da1 = jet->daughter(1);
            if (da0->mass() > da1->mass()) {
              boostedJet_wMass[icoll]->Fill(da0->mass());
              boostedJet_massDrop[icoll]->Fill(da0->mass() / jet->mass());
            } else {
              boostedJet_wMass[icoll]->Fill(da1->mass());
              boostedJet_massDrop[icoll]->Fill(da1->mass() / jet->mass());
            }
 
          } else {
            boostedJet_massDrop[icoll]->Fill(-1.0);
          }
 
        }  // end if collection is AK8 PFJets Puppi soft-drop
 
      }  // end if basic jet != 0
      countJet++;
    }
  }
 
  // PFMETs
  edm::Handle<std::vector<reco::PFMET> > pfMETCollection;
  bool ValidPFMET = iEvent.getByToken(PFMETToken_, pfMETCollection);
  if (!ValidPFMET)
    return;
 
  pfMet_pt->Fill((*pfMETCollection)[0].pt());
  pfMet_phi->Fill((*pfMETCollection)[0].phi());
}

References edm::View< T >::begin(), reco::PFJet::chargedEmEnergyFraction(), reco::PFJet::chargedHadronEnergyFraction(), edm::View< T >::end(), icoll, iEvent, metsig::jet, label, reco::Candidate::mass(), reco::PFJet::neutralEmEnergyFraction(), reco::PFJet::neutralHadronEnergyFraction(), ExoticaDQM_cfi::pfJetCollection, B2GDQM_cfi::pfMETCollection, reco::Candidate::phi(), DiDispStaMuonMonitor_cfi::pt, reco::Candidate::pt(), and reco::Candidate::rapidity().

analyzeSemiE()

void B2GDQM::analyzeSemiE ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

protectedvirtual

Definition at line 451 of file B2GDQM.cc.

                                                                          {
  edm::Handle<edm::View<reco::GsfElectron> > electronCollection;
  bool validElectrons = iEvent.getByToken(electronToken_, electronCollection);
 
  if (!validElectrons)
    return;
  if (electronCollection->empty())
    return;
  reco::GsfElectron const& electron = electronCollection->at(0);
  if (!(*elecSelect_)(electron))
    return;
 
  edm::Handle<edm::View<reco::BasicJet> > jetCollection;
  bool validJets = iEvent.getByToken(sdjetToken_, jetCollection);
  if (!validJets)
    return;
  if (jetCollection->size() < 2)
    return;
 
  double pt0 = -1.0;
  double dRMin = 999.0;
  edm::Ptr<reco::BasicJet> hadJet;  // highest pt jet with dphi(lep,jet) > pi/2
  edm::Ptr<reco::BasicJet> lepJet;  // closest jet to lepton with pt > ptMin
 
  for (auto ijet = jetCollection->begin(), ijetBegin = ijet, ijetEnd = jetCollection->end(); ijet != ijetEnd; ++ijet) {
    // Hadronic jets
    if (std::abs(reco::deltaPhi(electron, *ijet)) > M_PI * 0.5) {
      if (ijet->pt() > pt0 && ijet->p() > semiE_HadJetPtCut_) {
        hadJet = jetCollection->ptrAt(ijet - ijetBegin);
        pt0 = hadJet->pt();
      }
    }
    // Leptonic jets
    else if (ijet->pt() > semiE_LepJetPtCut_) {
      auto idRMin = reco::deltaR(electron, *ijet);
      if (idRMin < dRMin) {
        lepJet = jetCollection->ptrAt(ijet - ijetBegin);
        dRMin = idRMin;
      }
    }
  }
  if (hadJet.isAvailable() == false || lepJet.isAvailable() == false)
    return;
 
  auto lepJetP4 = lepJet->p4();
  const auto& electronP4 = electron.p4();
 
  double tot = lepJetP4.mag2();
  double ss = electronP4.Dot(lepJet->p4());
  double per = electronP4.mag2();
  if (tot > 0.0)
    per -= ss * ss / tot;
  if (per < 0)
    per = 0;
  double ptRel = per;
  bool pass2D = dRMin > semiE_dRMin_ || ptRel > semiE_ptRel_;
 
  if (!pass2D)
    return;
 
  semiE_ePt->Fill(electron.pt());
  semiE_eEta->Fill(electron.eta());
  semiE_ePhi->Fill(electron.phi());
  semiE_eDRMin->Fill(dRMin);
  semiE_ePtRel->Fill(ptRel);
 
  semiE_hadJetDR->Fill(reco::deltaR(electron, *hadJet));
  semiE_mttbar->Fill(0.0);
 
  semiE_hadJetPt->Fill(hadJet->pt());
  semiE_hadJetY->Fill(hadJet->rapidity());
  semiE_hadJetPhi->Fill(hadJet->phi());
  semiE_hadJetMass->Fill(hadJet->mass());
  if (hadJet->numberOfDaughters() > 2) {
    double wMass = hadJet->daughter(0)->mass() >= hadJet->daughter(1)->mass() ? hadJet->daughter(0)->mass()
                                                                              : hadJet->daughter(1)->mass();
    semiE_hadJetWMass->Fill(wMass);
  } else {
    semiE_hadJetWMass->Fill(-1.0);
  }
}

References funct::abs(), reco::CompositePtrCandidate::daughter(), reco::deltaPhi(), reco::deltaR(), HLT_2018_cff::dRMin, metsig::electron, pdwgLeptonRecoSkim_cfi::electronCollection, iEvent, edm::Ptr< T >::isAvailable(), jetfilter_cfi::jetCollection, M_PI, reco::Candidate::mass(), reco::LeafCandidate::mass(), reco::CompositePtrCandidate::numberOfDaughters(), reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), electrons_cff::ptRel, reco::LeafCandidate::rapidity(), contentValuesCheck::ss, and pseudoTop_cfi::wMass.

analyzeSemiMu()

void B2GDQM::analyzeSemiMu ( edm::Event const &  e, edm::EventSetup const &  eSetup  )

protectedvirtual

Definition at line 369 of file B2GDQM.cc.

                                                                           {
  edm::Handle<edm::View<reco::Muon> > muonCollection;
  bool validMuons = iEvent.getByToken(muonToken_, muonCollection);
 
  if (!validMuons)
    return;
  if (muonCollection->empty())
    return;
  reco::Muon const& muon = muonCollection->at(0);
  if (!(*muonSelect_)(muon))
    return;
 
  edm::Handle<edm::View<reco::BasicJet> > jetCollection;
  bool validJets = iEvent.getByToken(sdjetToken_, jetCollection);
  if (!validJets)
    return;
  if (jetCollection->size() < 2)
    return;
 
  double pt0 = -1.0;
  double dRMin = 999.0;
  edm::Ptr<reco::BasicJet> hadJet;  // highest pt jet with dphi(lep,jet) > pi/2
  edm::Ptr<reco::BasicJet> lepJet;  // closest jet to lepton with pt > ptMin
 
  for (auto ijet = jetCollection->begin(), ijetBegin = ijet, ijetEnd = jetCollection->end(); ijet != ijetEnd; ++ijet) {
    // Hadronic jets
    if (std::abs(reco::deltaPhi(muon, *ijet)) > M_PI * 0.5) {
      if (ijet->pt() > pt0 && ijet->p() > semiMu_HadJetPtCut_) {
        hadJet = jetCollection->ptrAt(ijet - ijetBegin);
        pt0 = hadJet->pt();
      }
    }
    // Leptonic jets
    else if (ijet->pt() > semiMu_LepJetPtCut_) {
      auto idRMin = reco::deltaR(muon, *ijet);
      if (idRMin < dRMin) {
        lepJet = jetCollection->ptrAt(ijet - ijetBegin);
        dRMin = idRMin;
      }
    }
  }
  if (hadJet.isAvailable() == false || lepJet.isAvailable() == false)
    return;
 
  auto lepJetP4 = lepJet->p4();
  const auto& muonP4 = muon.p4();
 
  double tot = lepJetP4.mag2();
  double ss = muonP4.Dot(lepJet->p4());
  double per = muonP4.mag2();
  if (tot > 0.0)
    per -= ss * ss / tot;
  if (per < 0)
    per = 0;
  double ptRel = per;
  bool pass2D = dRMin > semiMu_dRMin_ || ptRel > semiMu_ptRel_;
 
  if (!pass2D)
    return;
 
  semiMu_muPt->Fill(muon.pt());
  semiMu_muEta->Fill(muon.eta());
  semiMu_muPhi->Fill(muon.phi());
  semiMu_muDRMin->Fill(dRMin);
  semiMu_muPtRel->Fill(ptRel);
 
  semiMu_hadJetDR->Fill(reco::deltaR(muon, *hadJet));
  semiMu_mttbar->Fill(0.0);
 
  semiMu_hadJetPt->Fill(hadJet->pt());
  semiMu_hadJetY->Fill(hadJet->rapidity());
  semiMu_hadJetPhi->Fill(hadJet->phi());
  semiMu_hadJetMass->Fill(hadJet->mass());
  if (hadJet->numberOfDaughters() > 2) {
    double wMass = hadJet->daughter(0)->mass() >= hadJet->daughter(1)->mass() ? hadJet->daughter(0)->mass()
                                                                              : hadJet->daughter(1)->mass();
    semiMu_hadJetWMass->Fill(wMass);
  } else {
    semiMu_hadJetWMass->Fill(-1.0);
  }
}

References funct::abs(), reco::CompositePtrCandidate::daughter(), reco::deltaPhi(), reco::deltaR(), HLT_2018_cff::dRMin, iEvent, edm::Ptr< T >::isAvailable(), jetfilter_cfi::jetCollection, M_PI, reco::Candidate::mass(), reco::LeafCandidate::mass(), pdwgLeptonRecoSkim_cfi::muonCollection, reco::CompositePtrCandidate::numberOfDaughters(), reco::LeafCandidate::p4(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), electrons_cff::ptRel, reco::LeafCandidate::rapidity(), contentValuesCheck::ss, and pseudoTop_cfi::wMass.

bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 142 of file B2GDQM.cc.

                                                                                     {
  bei.setCurrentFolder("Physics/B2G");
 
  //--- Jets
 
  for (unsigned int icoll = 0; icoll < jetLabels_.size(); ++icoll) {
    std::stringstream ss;
    ss << "Physics/B2G/" << jetLabels_[icoll].label();
    bei.setCurrentFolder(ss.str());
    pfJet_pt.push_back(bei.book1D("pfJet_pt", "Pt of PFJet (GeV)", 50, 0.0, 1000));
    pfJet_y.push_back(bei.book1D("pfJet_y", "Rapidity of PFJet", 60, -6.0, 6.0));
    pfJet_phi.push_back(bei.book1D("pfJet_phi", "#phi of PFJet (radians)", 60, -3.14159, 3.14159));
    pfJet_m.push_back(bei.book1D("pfJet_m", "Mass of PFJet (GeV)", 50, 0.0, 500));
    pfJet_chef.push_back(bei.book1D("pfJet_pfchef", "PFJetID CHEF", 50, 0.0, 1.0));
    pfJet_nhef.push_back(bei.book1D("pfJet_pfnhef", "PFJetID NHEF", 50, 0.0, 1.0));
    pfJet_cemf.push_back(bei.book1D("pfJet_pfcemf", "PFJetID CEMF", 50, 0.0, 1.0));
    pfJet_nemf.push_back(bei.book1D("pfJet_pfnemf", "PFJetID NEMF", 50, 0.0, 1.0));
 
    boostedJet_subjetPt.push_back(bei.book1D("boostedJet_subjetPt", "Pt of subjets (GeV)", 50, 0.0, 500));
    boostedJet_subjetY.push_back(bei.book1D("boostedJet_subjetY", "Rapidity of subjets", 60, -6.0, 6.0));
    boostedJet_subjetPhi.push_back(
        bei.book1D("boostedJet_subjetPhi", "#phi of subjets (radians)", 60, -3.14159, 3.14159));
    boostedJet_subjetM.push_back(bei.book1D("boostedJet_subjetM", "Mass of subjets (GeV)", 50, 0.0, 250.));
    boostedJet_subjetN.push_back(bei.book1D("boostedJet_subjetN", "Number of subjets", 10, 0, 10));
    boostedJet_massDrop.push_back(bei.book1D("boostedJet_massDrop", "Mass drop for W-like jets", 50, 0.0, 1.0));
    boostedJet_wMass.push_back(bei.book1D("boostedJet_wMass", "W Mass for top-like jets", 50, 0.0, 250.0));
  }
 
  bei.setCurrentFolder("Physics/B2G/MET");
  pfMet_pt = bei.book1D("pfMet_pt", "Pf Missing p_{T}; GeV", 50, 0.0, 500);
  pfMet_phi = bei.book1D("pfMet_phi", "Pf Missing p_{T} #phi;#phi (radians)", 35, -3.5, 3.5);
 
  //--- Mu+Jets
  bei.setCurrentFolder("Physics/B2G/SemiMu");
  semiMu_muPt = bei.book1D("semiMu_muPt", "Pt of Muon in #mu+Jets Channel (GeV)", 50, 0.0, 1000);
  semiMu_muEta = bei.book1D("semiMu_muEta", "#eta of Muon in #mu+Jets Channel", 60, -6.0, 6.0);
  semiMu_muPhi = bei.book1D("semiMu_muPhi", "#phi of Muon in #mu+Jets Channel (radians)", 60, -3.14159, 3.14159);
  semiMu_muDRMin = bei.book1D("semiMu_muDRMin", "#Delta R(E,nearest jet) in #mu+Jets Channel", 50, 0, 10.0);
  semiMu_muPtRel = bei.book1D("semiMu_muPtRel", "p_{T}^{REL} in #mu+Jets Channel", 60, 0, 300.);
  semiMu_hadJetDR = bei.book1D("semiMu_hadJetDR", "#Delta R(E,had jet) in #mu+Jets Channel", 50, 0, 10.0);
  semiMu_hadJetPt =
      bei.book1D("semiMu_hadJetPt", "Pt of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 1000);
  semiMu_hadJetY = bei.book1D("semiMu_hadJetY", "Rapidity of Leading Hadronic Jet in #mu+Jets Channel", 60, -6.0, 6.0);
  semiMu_hadJetPhi = bei.book1D(
      "semiMu_hadJetPhi", "#phi of Leading Hadronic Jet in #mu+Jets Channel (radians)", 60, -3.14159, 3.14159);
  semiMu_hadJetMass =
      bei.book1D("semiMu_hadJetMass", "Mass of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 500);
  semiMu_hadJetWMass =
      bei.book1D("semiMu_hadJetwMass", "W Mass for Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 250.0);
  semiMu_mttbar = bei.book1D("semiMu_mttbar", "Mass of #mu+Jets ttbar Candidate", 100, 0., 5000.);
 
  //--- E+Jets
  bei.setCurrentFolder("Physics/B2G/SemiE");
  semiE_ePt = bei.book1D("semiE_ePt", "Pt of Electron in e+Jets Channel (GeV)", 50, 0.0, 1000);
  semiE_eEta = bei.book1D("semiE_eEta", "#eta of Electron in e+Jets Channel", 60, -6.0, 6.0);
  semiE_ePhi = bei.book1D("semiE_ePhi", "#phi of Electron in e+Jets Channel (radians)", 60, -3.14159, 3.14159);
  semiE_eDRMin = bei.book1D("semiE_eDRMin", "#Delta R(E,nearest jet) in e+Jets Channel", 50, 0, 10.0);
  semiE_ePtRel = bei.book1D("semiE_ePtRel", "p_{T}^{REL} in e+Jets Channel", 60, 0, 300.);
  semiE_hadJetDR = bei.book1D("semiE_hadJetDR", "#Delta R(E,had jet) in e+Jets Channel", 50, 0, 10.0);
  semiE_hadJetPt = bei.book1D("semiE_hadJetPt", "Pt of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 1000);
  semiE_hadJetY = bei.book1D("semiE_hadJetY", "Rapidity of Leading Hadronic Jet in e+Jets Channel", 60, -6.0, 6.0);
  semiE_hadJetPhi =
      bei.book1D("semiE_hadJetPhi", "#phi of Leading Hadronic Jet in e+Jets Channel (radians)", 60, -3.14159, 3.14159);
  semiE_hadJetMass =
      bei.book1D("semiE_hadJetMass", "Mass of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 500);
  semiE_hadJetWMass =
      bei.book1D("semiE_hadJetwMass", "W Mass for Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 250.0);
  semiE_mttbar = bei.book1D("semiE_mttbar", "Mass of e+Jets ttbar Candidate", 100, 0., 5000.);
 
  //--- All-hadronic
  bei.setCurrentFolder("Physics/B2G/AllHad");
  allHad_pt0 = bei.book1D("allHad_pt0", "Pt of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 1000);
  allHad_y0 = bei.book1D("allHad_y0", "Rapidity of Leading All-Hadronic PFJet", 60, -6.0, 6.0);
  allHad_phi0 = bei.book1D("allHad_phi0", "#phi of Leading All-Hadronic PFJet (radians)", 60, -3.14159, 3.14159);
  allHad_mass0 = bei.book1D("allHad_mass0", "Mass of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 500);
  allHad_wMass0 = bei.book1D("allHad_wMass0", "W Mass for Leading All-Hadronic PFJet (GeV)", 50, 0.0, 250.0);
  allHad_pt1 = bei.book1D("allHad_pt1", "Pt of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 1000);
  allHad_y1 = bei.book1D("allHad_y1", "Rapidity of Subleading All-Hadronic PFJet", 60, -6.0, 6.0);
  allHad_phi1 = bei.book1D("allHad_phi1", "#phi of Subleading All-Hadronic PFJet (radians)", 60, -3.14159, 3.14159);
  allHad_mass1 = bei.book1D("allHad_mass1", "Mass of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 500);
  allHad_wMass1 = bei.book1D("allHad_wMass1", "W Mass for Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 250.0);
  allHad_mttbar = bei.book1D("allHad_mttbar", "Mass of All-Hadronic ttbar Candidate", 100, 0., 5000.);
}

References dqm::implementation::IBooker::book1D(), icoll, dqm::implementation::NavigatorBase::setCurrentFolder(), and contentValuesCheck::ss.

Member Data Documentation

allHad_mass0

MonitorElement* B2GDQM::allHad_mass0

private

Definition at line 198 of file B2GDQM.h.

allHad_mass1

MonitorElement* B2GDQM::allHad_mass1

private

Definition at line 203 of file B2GDQM.h.

allHad_mttbar

MonitorElement* B2GDQM::allHad_mttbar

private

Definition at line 205 of file B2GDQM.h.

allHad_phi0

MonitorElement* B2GDQM::allHad_phi0

private

Definition at line 197 of file B2GDQM.h.

allHad_phi1

MonitorElement* B2GDQM::allHad_phi1

private

Definition at line 202 of file B2GDQM.h.

allHad_pt0

MonitorElement* B2GDQM::allHad_pt0

private

Definition at line 195 of file B2GDQM.h.

allHad_pt1

MonitorElement* B2GDQM::allHad_pt1

private

Definition at line 200 of file B2GDQM.h.

allHad_wMass0

MonitorElement* B2GDQM::allHad_wMass0

private

Definition at line 199 of file B2GDQM.h.

allHad_wMass1

MonitorElement* B2GDQM::allHad_wMass1

private

Definition at line 204 of file B2GDQM.h.

allHad_y0

MonitorElement* B2GDQM::allHad_y0

private

Definition at line 196 of file B2GDQM.h.

allHad_y1

MonitorElement* B2GDQM::allHad_y1

private

Definition at line 201 of file B2GDQM.h.

allHadDeltaPhiCut_

double B2GDQM::allHadDeltaPhiCut_

private

Definition at line 129 of file B2GDQM.h.

allHadPtCut_

double B2GDQM::allHadPtCut_

private

Definition at line 127 of file B2GDQM.h.

allHadRapidityCut_

double B2GDQM::allHadRapidityCut_

private

Definition at line 128 of file B2GDQM.h.

boostedJet_massDrop

std::vector<MonitorElement*> B2GDQM::boostedJet_massDrop

private

Definition at line 164 of file B2GDQM.h.

boostedJet_subjetM

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetM

private

Definition at line 162 of file B2GDQM.h.

boostedJet_subjetN

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetN

private

Definition at line 163 of file B2GDQM.h.

boostedJet_subjetPhi

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetPhi

private

Definition at line 161 of file B2GDQM.h.

boostedJet_subjetPt

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetPt

private

Definition at line 159 of file B2GDQM.h.

boostedJet_subjetY

std::vector<MonitorElement*> B2GDQM::boostedJet_subjetY

private

Definition at line 160 of file B2GDQM.h.

boostedJet_wMass

std::vector<MonitorElement*> B2GDQM::boostedJet_wMass

private

Definition at line 165 of file B2GDQM.h.

elecSelect_

std::shared_ptr<StringCutObjectSelector<reco::GsfElectron> > B2GDQM::elecSelect_

private

Definition at line 145 of file B2GDQM.h.

electronLabel_

edm::InputTag B2GDQM::electronLabel_

private

Definition at line 118 of file B2GDQM.h.

electronToken_

edm::EDGetTokenT<edm::View<reco::GsfElectron> > B2GDQM::electronToken_

private

Definition at line 119 of file B2GDQM.h.

hltConfigProvider_

HLTConfigProvider B2GDQM::hltConfigProvider_

private

Definition at line 96 of file B2GDQM.h.

ievt

int B2GDQM::ievt

private

Definition at line 94 of file B2GDQM.h.

irun

int B2GDQM::irun

private

Definition at line 94 of file B2GDQM.h.

isValidHltConfig_

bool B2GDQM::isValidHltConfig_

private

Definition at line 97 of file B2GDQM.h.

jetLabels_

std::vector<edm::InputTag> B2GDQM::jetLabels_

private

Definition at line 105 of file B2GDQM.h.

jetPtMins_

std::vector<double> B2GDQM::jetPtMins_

private

Definition at line 125 of file B2GDQM.h.

jetTokens_

std::vector<edm::EDGetTokenT<edm::View<reco::Jet> > > B2GDQM::jetTokens_

private

Definition at line 106 of file B2GDQM.h.

muonLabel_

edm::InputTag B2GDQM::muonLabel_

private

Definition at line 115 of file B2GDQM.h.

muonSelect_

std::shared_ptr<StringCutObjectSelector<reco::Muon> > B2GDQM::muonSelect_

private

Definition at line 136 of file B2GDQM.h.

muonToken_

edm::EDGetTokenT<edm::View<reco::Muon> > B2GDQM::muonToken_

private

Definition at line 116 of file B2GDQM.h.

nEvents_

int B2GDQM::nEvents_

private

Definition at line 94 of file B2GDQM.h.

nLumiSecs_

int B2GDQM::nLumiSecs_

private

Definition at line 93 of file B2GDQM.h.

pfJet_cemf

std::vector<MonitorElement*> B2GDQM::pfJet_cemf

private

Definition at line 157 of file B2GDQM.h.

pfJet_chef

std::vector<MonitorElement*> B2GDQM::pfJet_chef

private

Definition at line 155 of file B2GDQM.h.

pfJet_m

std::vector<MonitorElement*> B2GDQM::pfJet_m

private

Definition at line 154 of file B2GDQM.h.

pfJet_nemf

std::vector<MonitorElement*> B2GDQM::pfJet_nemf

private

Definition at line 158 of file B2GDQM.h.

pfJet_nhef

std::vector<MonitorElement*> B2GDQM::pfJet_nhef

private

Definition at line 156 of file B2GDQM.h.

pfJet_phi

std::vector<MonitorElement*> B2GDQM::pfJet_phi

private

Definition at line 153 of file B2GDQM.h.

pfJet_pt

std::vector<MonitorElement*> B2GDQM::pfJet_pt

private

Definition at line 151 of file B2GDQM.h.

pfJet_y

std::vector<MonitorElement*> B2GDQM::pfJet_y

private

Definition at line 152 of file B2GDQM.h.

PFJetCorService_

std::string B2GDQM::PFJetCorService_

private

Definition at line 147 of file B2GDQM.h.

pfMet_phi

MonitorElement* B2GDQM::pfMet_phi

private

Definition at line 167 of file B2GDQM.h.

pfMet_pt

MonitorElement* B2GDQM::pfMet_pt

private

Definition at line 166 of file B2GDQM.h.

PFMETLabel_

edm::InputTag B2GDQM::PFMETLabel_

private

Definition at line 107 of file B2GDQM.h.

PFMETToken_

edm::EDGetTokenT<std::vector<reco::PFMET> > B2GDQM::PFMETToken_

private

Definition at line 108 of file B2GDQM.h.

sdjetLabel_

edm::InputTag B2GDQM::sdjetLabel_

private

Definition at line 112 of file B2GDQM.h.

sdjetToken_

edm::EDGetTokenT<edm::View<reco::BasicJet> > B2GDQM::sdjetToken_

private

Definition at line 113 of file B2GDQM.h.

semiE_dphiHadCut_

double B2GDQM::semiE_dphiHadCut_

private

Definition at line 140 of file B2GDQM.h.

semiE_dRMin_

double B2GDQM::semiE_dRMin_

private

Definition at line 142 of file B2GDQM.h.

semiE_eDRMin

MonitorElement* B2GDQM::semiE_eDRMin

private

Definition at line 185 of file B2GDQM.h.

semiE_eEta

MonitorElement* B2GDQM::semiE_eEta

private

Definition at line 183 of file B2GDQM.h.

semiE_ePhi

MonitorElement* B2GDQM::semiE_ePhi

private

Definition at line 184 of file B2GDQM.h.

semiE_ePt

MonitorElement* B2GDQM::semiE_ePt

private

Definition at line 182 of file B2GDQM.h.

semiE_ePtRel

MonitorElement* B2GDQM::semiE_ePtRel

private

Definition at line 186 of file B2GDQM.h.

semiE_hadJetDR

MonitorElement* B2GDQM::semiE_hadJetDR

private

Definition at line 187 of file B2GDQM.h.

semiE_hadJetMass

MonitorElement* B2GDQM::semiE_hadJetMass

private

Definition at line 191 of file B2GDQM.h.

semiE_hadJetPhi

MonitorElement* B2GDQM::semiE_hadJetPhi

private

Definition at line 190 of file B2GDQM.h.

semiE_hadJetPt

MonitorElement* B2GDQM::semiE_hadJetPt

private

Definition at line 188 of file B2GDQM.h.

semiE_HadJetPtCut_

double B2GDQM::semiE_HadJetPtCut_

private

Definition at line 138 of file B2GDQM.h.

semiE_hadJetWMass

MonitorElement* B2GDQM::semiE_hadJetWMass

private

Definition at line 192 of file B2GDQM.h.

semiE_hadJetY

MonitorElement* B2GDQM::semiE_hadJetY

private

Definition at line 189 of file B2GDQM.h.

semiE_LepJetPtCut_

double B2GDQM::semiE_LepJetPtCut_

private

Definition at line 139 of file B2GDQM.h.

semiE_mttbar

MonitorElement* B2GDQM::semiE_mttbar

private

Definition at line 193 of file B2GDQM.h.

semiE_ptRel_

double B2GDQM::semiE_ptRel_

private

Definition at line 144 of file B2GDQM.h.

semiMu_dphiHadCut_

double B2GDQM::semiMu_dphiHadCut_

private

Definition at line 133 of file B2GDQM.h.

semiMu_dRMin_

double B2GDQM::semiMu_dRMin_

private

Definition at line 134 of file B2GDQM.h.

semiMu_hadJetDR

MonitorElement* B2GDQM::semiMu_hadJetDR

private

Definition at line 174 of file B2GDQM.h.

semiMu_hadJetMass

MonitorElement* B2GDQM::semiMu_hadJetMass

private

Definition at line 178 of file B2GDQM.h.

semiMu_hadJetPhi

MonitorElement* B2GDQM::semiMu_hadJetPhi

private

Definition at line 177 of file B2GDQM.h.

semiMu_hadJetPt

MonitorElement* B2GDQM::semiMu_hadJetPt

private

Definition at line 175 of file B2GDQM.h.

semiMu_HadJetPtCut_

double B2GDQM::semiMu_HadJetPtCut_

private

Definition at line 131 of file B2GDQM.h.

semiMu_hadJetWMass

MonitorElement* B2GDQM::semiMu_hadJetWMass

private

Definition at line 179 of file B2GDQM.h.

semiMu_hadJetY

MonitorElement* B2GDQM::semiMu_hadJetY

private

Definition at line 176 of file B2GDQM.h.

semiMu_LepJetPtCut_

double B2GDQM::semiMu_LepJetPtCut_

private

Definition at line 132 of file B2GDQM.h.

semiMu_mttbar

MonitorElement* B2GDQM::semiMu_mttbar

private

Definition at line 180 of file B2GDQM.h.

semiMu_muDRMin

MonitorElement* B2GDQM::semiMu_muDRMin

private

Definition at line 172 of file B2GDQM.h.

semiMu_muEta

MonitorElement* B2GDQM::semiMu_muEta

private

Definition at line 170 of file B2GDQM.h.

semiMu_muPhi

MonitorElement* B2GDQM::semiMu_muPhi

private

Definition at line 171 of file B2GDQM.h.

semiMu_muPt

MonitorElement* B2GDQM::semiMu_muPt

private

Definition at line 169 of file B2GDQM.h.

semiMu_muPtRel

MonitorElement* B2GDQM::semiMu_muPtRel

private

Definition at line 173 of file B2GDQM.h.

semiMu_ptRel_

double B2GDQM::semiMu_ptRel_

private

Definition at line 135 of file B2GDQM.h.

theTriggerResultsCollection

edm::InputTag B2GDQM::theTriggerResultsCollection

private

Definition at line 100 of file B2GDQM.h.

triggerResults_

edm::Handle<edm::TriggerResults> B2GDQM::triggerResults_

private

Definition at line 103 of file B2GDQM.h.

triggerToken_

edm::EDGetTokenT<edm::TriggerResults> B2GDQM::triggerToken_

private

Definition at line 101 of file B2GDQM.h.