|
|
#include <B2GDQM.h>
Definition at line 78 of file B2GDQM.h.
◆ B2GDQM()
◆ ~B2GDQM()
◆ analyze()
◆ analyzeAllHad()
Definition at line 321 of file B2GDQM.cc.
365 auto p4cand = (jet0->
p4() + jet1->
p4());
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()
Definition at line 236 of file B2GDQM.cc.
264 reco::PFJet const* pfjet = dynamic_cast<reco::PFJet const*>(&*
jet);
266 if (pfjet !=
nullptr) {
277 if (basicjet !=
nullptr) {
280 for (
unsigned int ida = 0; ida <
jet->numberOfDaughters(); ++ida) {
289 if (
jet->numberOfDaughters() > 1) {
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()
Definition at line 451 of file B2GDQM.cc.
471 double dRMin = 999.0;
486 if (idRMin <
dRMin) {
495 auto lepJetP4 = lepJet->
p4();
496 const auto& electronP4 =
electron.p4();
498 double tot = lepJetP4.mag2();
499 double ss = electronP4.Dot(lepJet->
p4());
500 double per = electronP4.mag2();
502 per -=
ss *
ss / tot;
References funct::abs(), reco::CompositePtrCandidate::daughter(), reco::deltaPhi(), reco::deltaR(), HLT_FULL_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()
Definition at line 369 of file B2GDQM.cc.
389 double dRMin = 999.0;
404 if (idRMin <
dRMin) {
413 auto lepJetP4 = lepJet->
p4();
414 const auto& muonP4 =
muon.p4();
416 double tot = lepJetP4.mag2();
417 double ss = muonP4.Dot(lepJet->
p4());
418 double per = muonP4.mag2();
420 per -=
ss *
ss / tot;
References funct::abs(), reco::CompositePtrCandidate::daughter(), reco::deltaPhi(), reco::deltaR(), HLT_FULL_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()
Implements DQMEDAnalyzer.
Definition at line 142 of file B2GDQM.cc.
148 std::stringstream
ss;
151 pfJet_pt.push_back(bei.
book1D(
"pfJet_pt",
"Pt of PFJet (GeV)", 50, 0.0, 1000));
152 pfJet_y.push_back(bei.
book1D(
"pfJet_y",
"Rapidity of PFJet", 60, -6.0, 6.0));
153 pfJet_phi.push_back(bei.
book1D(
"pfJet_phi",
"#phi of PFJet (radians)", 60, -3.14159, 3.14159));
154 pfJet_m.push_back(bei.
book1D(
"pfJet_m",
"Mass of PFJet (GeV)", 50, 0.0, 500));
155 pfJet_chef.push_back(bei.
book1D(
"pfJet_pfchef",
"PFJetID CHEF", 50, 0.0, 1.0));
156 pfJet_nhef.push_back(bei.
book1D(
"pfJet_pfnhef",
"PFJetID NHEF", 50, 0.0, 1.0));
157 pfJet_cemf.push_back(bei.
book1D(
"pfJet_pfcemf",
"PFJetID CEMF", 50, 0.0, 1.0));
158 pfJet_nemf.push_back(bei.
book1D(
"pfJet_pfnemf",
"PFJetID NEMF", 50, 0.0, 1.0));
163 bei.
book1D(
"boostedJet_subjetPhi",
"#phi of subjets (radians)", 60, -3.14159, 3.14159));
171 pfMet_pt = bei.
book1D(
"pfMet_pt",
"Pf Missing p_{T}; GeV", 50, 0.0, 500);
172 pfMet_phi = bei.
book1D(
"pfMet_phi",
"Pf Missing p_{T} #phi;#phi (radians)", 35, -3.5, 3.5);
176 semiMu_muPt = bei.
book1D(
"semiMu_muPt",
"Pt of Muon in #mu+Jets Channel (GeV)", 50, 0.0, 1000);
177 semiMu_muEta = bei.
book1D(
"semiMu_muEta",
"#eta of Muon in #mu+Jets Channel", 60, -6.0, 6.0);
178 semiMu_muPhi = bei.
book1D(
"semiMu_muPhi",
"#phi of Muon in #mu+Jets Channel (radians)", 60, -3.14159, 3.14159);
179 semiMu_muDRMin = bei.
book1D(
"semiMu_muDRMin",
"#Delta R(E,nearest jet) in #mu+Jets Channel", 50, 0, 10.0);
181 semiMu_hadJetDR = bei.
book1D(
"semiMu_hadJetDR",
"#Delta R(E,had jet) in #mu+Jets Channel", 50, 0, 10.0);
183 bei.
book1D(
"semiMu_hadJetPt",
"Pt of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 1000);
184 semiMu_hadJetY = bei.
book1D(
"semiMu_hadJetY",
"Rapidity of Leading Hadronic Jet in #mu+Jets Channel", 60, -6.0, 6.0);
186 "semiMu_hadJetPhi",
"#phi of Leading Hadronic Jet in #mu+Jets Channel (radians)", 60, -3.14159, 3.14159);
188 bei.
book1D(
"semiMu_hadJetMass",
"Mass of Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 500);
190 bei.
book1D(
"semiMu_hadJetwMass",
"W Mass for Leading Hadronic Jet in #mu+Jets Channel (GeV)", 50, 0.0, 250.0);
191 semiMu_mttbar = bei.
book1D(
"semiMu_mttbar",
"Mass of #mu+Jets ttbar Candidate", 100, 0., 5000.);
195 semiE_ePt = bei.
book1D(
"semiE_ePt",
"Pt of Electron in e+Jets Channel (GeV)", 50, 0.0, 1000);
196 semiE_eEta = bei.
book1D(
"semiE_eEta",
"#eta of Electron in e+Jets Channel", 60, -6.0, 6.0);
197 semiE_ePhi = bei.
book1D(
"semiE_ePhi",
"#phi of Electron in e+Jets Channel (radians)", 60, -3.14159, 3.14159);
198 semiE_eDRMin = bei.
book1D(
"semiE_eDRMin",
"#Delta R(E,nearest jet) in e+Jets Channel", 50, 0, 10.0);
199 semiE_ePtRel = bei.
book1D(
"semiE_ePtRel",
"p_{T}^{REL} in e+Jets Channel", 60, 0, 300.);
200 semiE_hadJetDR = bei.
book1D(
"semiE_hadJetDR",
"#Delta R(E,had jet) in e+Jets Channel", 50, 0, 10.0);
201 semiE_hadJetPt = bei.
book1D(
"semiE_hadJetPt",
"Pt of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 1000);
202 semiE_hadJetY = bei.
book1D(
"semiE_hadJetY",
"Rapidity of Leading Hadronic Jet in e+Jets Channel", 60, -6.0, 6.0);
204 bei.
book1D(
"semiE_hadJetPhi",
"#phi of Leading Hadronic Jet in e+Jets Channel (radians)", 60, -3.14159, 3.14159);
206 bei.
book1D(
"semiE_hadJetMass",
"Mass of Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 500);
208 bei.
book1D(
"semiE_hadJetwMass",
"W Mass for Leading Hadronic Jet in e+Jets Channel (GeV)", 50, 0.0, 250.0);
209 semiE_mttbar = bei.
book1D(
"semiE_mttbar",
"Mass of e+Jets ttbar Candidate", 100, 0., 5000.);
213 allHad_pt0 = bei.
book1D(
"allHad_pt0",
"Pt of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 1000);
214 allHad_y0 = bei.
book1D(
"allHad_y0",
"Rapidity of Leading All-Hadronic PFJet", 60, -6.0, 6.0);
215 allHad_phi0 = bei.
book1D(
"allHad_phi0",
"#phi of Leading All-Hadronic PFJet (radians)", 60, -3.14159, 3.14159);
216 allHad_mass0 = bei.
book1D(
"allHad_mass0",
"Mass of Leading All-Hadronic PFJet (GeV)", 50, 0.0, 500);
217 allHad_wMass0 = bei.
book1D(
"allHad_wMass0",
"W Mass for Leading All-Hadronic PFJet (GeV)", 50, 0.0, 250.0);
218 allHad_pt1 = bei.
book1D(
"allHad_pt1",
"Pt of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 1000);
219 allHad_y1 = bei.
book1D(
"allHad_y1",
"Rapidity of Subleading All-Hadronic PFJet", 60, -6.0, 6.0);
220 allHad_phi1 = bei.
book1D(
"allHad_phi1",
"#phi of Subleading All-Hadronic PFJet (radians)", 60, -3.14159, 3.14159);
221 allHad_mass1 = bei.
book1D(
"allHad_mass1",
"Mass of Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 500);
222 allHad_wMass1 = bei.
book1D(
"allHad_wMass1",
"W Mass for Subleading All-Hadronic PFJet (GeV)", 50, 0.0, 250.0);
223 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.
◆ allHad_mass0
◆ allHad_mass1
◆ allHad_mttbar
◆ allHad_phi0
◆ allHad_phi1
◆ allHad_pt0
◆ allHad_pt1
◆ allHad_wMass0
◆ allHad_wMass1
◆ allHad_y0
◆ allHad_y1
◆ allHadDeltaPhiCut_
double B2GDQM::allHadDeltaPhiCut_ |
|
private |
◆ allHadPtCut_
double B2GDQM::allHadPtCut_ |
|
private |
◆ allHadRapidityCut_
double B2GDQM::allHadRapidityCut_ |
|
private |
◆ boostedJet_massDrop
◆ boostedJet_subjetM
◆ boostedJet_subjetN
◆ boostedJet_subjetPhi
◆ boostedJet_subjetPt
◆ boostedJet_subjetY
◆ boostedJet_wMass
◆ elecSelect_
◆ electronLabel_
◆ electronToken_
◆ hltConfigProvider_
◆ ievt
◆ irun
◆ isValidHltConfig_
bool B2GDQM::isValidHltConfig_ |
|
private |
◆ jetLabels_
◆ jetPtMins_
std::vector<double> B2GDQM::jetPtMins_ |
|
private |
◆ jetTokens_
◆ muonLabel_
◆ muonSelect_
◆ muonToken_
◆ nEvents_
◆ nLumiSecs_
◆ pfJet_cemf
◆ pfJet_chef
◆ pfJet_m
◆ pfJet_nemf
◆ pfJet_nhef
◆ pfJet_phi
◆ pfJet_pt
◆ pfJet_y
◆ PFJetCorService_
std::string B2GDQM::PFJetCorService_ |
|
private |
◆ pfMet_phi
◆ pfMet_pt
◆ PFMETLabel_
◆ PFMETToken_
◆ sdjetLabel_
◆ sdjetToken_
◆ semiE_dphiHadCut_
double B2GDQM::semiE_dphiHadCut_ |
|
private |
◆ semiE_dRMin_
double B2GDQM::semiE_dRMin_ |
|
private |
◆ semiE_eDRMin
◆ semiE_eEta
◆ semiE_ePhi
◆ semiE_ePt
◆ semiE_ePtRel
◆ semiE_hadJetDR
◆ semiE_hadJetMass
◆ semiE_hadJetPhi
◆ semiE_hadJetPt
◆ semiE_HadJetPtCut_
double B2GDQM::semiE_HadJetPtCut_ |
|
private |
◆ semiE_hadJetWMass
◆ semiE_hadJetY
◆ semiE_LepJetPtCut_
double B2GDQM::semiE_LepJetPtCut_ |
|
private |
◆ semiE_mttbar
◆ semiE_ptRel_
double B2GDQM::semiE_ptRel_ |
|
private |
◆ semiMu_dphiHadCut_
double B2GDQM::semiMu_dphiHadCut_ |
|
private |
◆ semiMu_dRMin_
double B2GDQM::semiMu_dRMin_ |
|
private |
◆ semiMu_hadJetDR
◆ semiMu_hadJetMass
◆ semiMu_hadJetPhi
◆ semiMu_hadJetPt
◆ semiMu_HadJetPtCut_
double B2GDQM::semiMu_HadJetPtCut_ |
|
private |
◆ semiMu_hadJetWMass
◆ semiMu_hadJetY
◆ semiMu_LepJetPtCut_
double B2GDQM::semiMu_LepJetPtCut_ |
|
private |
◆ semiMu_mttbar
◆ semiMu_muDRMin
◆ semiMu_muEta
◆ semiMu_muPhi
◆ semiMu_muPt
◆ semiMu_muPtRel
◆ semiMu_ptRel_
double B2GDQM::semiMu_ptRel_ |
|
private |
◆ theTriggerResultsCollection
◆ triggerResults_
◆ triggerToken_
MonitorElement * allHad_pt0
MonitorElement * allHad_wMass0
const_iterator begin() const
std::vector< MonitorElement * > pfJet_chef
std::vector< MonitorElement * > boostedJet_wMass
MonitorElement * allHad_y1
MonitorElement * pfMet_phi
virtual double mass() const =0
mass
MonitorElement * semiE_hadJetWMass
constexpr double deltaPhi(double phi1, double phi2)
double semiMu_dphiHadCut_
MonitorElement * allHad_phi1
virtual void setCurrentFolder(std::string const &fullpath)
std::vector< MonitorElement * > boostedJet_subjetY
MonitorElement * semiE_eDRMin
edm::InputTag PFMETLabel_
virtual double pt() const =0
transverse momentum
double pt() const final
transverse momentum
float chargedHadronEnergyFraction() const
chargedHadronEnergyFraction
MonitorElement * pfMet_pt
Log< level::Info, false > LogInfo
MonitorElement * semiMu_hadJetY
MonitorElement * semiE_eEta
MonitorElement * semiMu_muDRMin
MonitorElement * semiMu_hadJetDR
std::string PFJetCorService_
MonitorElement * allHad_wMass1
virtual double rapidity() const =0
rapidity
virtual void analyzeAllHad(edm::Event const &e, edm::EventSetup const &eSetup)
MonitorElement * allHad_y0
std::vector< MonitorElement * > boostedJet_subjetPt
Jets made from CaloTowers.
std::vector< MonitorElement * > pfJet_cemf
std::vector< double > jetPtMins_
std::shared_ptr< StringCutObjectSelector< reco::Muon > > muonSelect_
float chargedEmEnergyFraction() const
chargedEmEnergyFraction
std::vector< MonitorElement * > boostedJet_subjetPhi
MonitorElement * semiE_hadJetMass
MonitorElement * semiE_ePt
MonitorElement * semiE_mttbar
MonitorElement * semiMu_mttbar
std::vector< MonitorElement * > pfJet_y
edm::InputTag theTriggerResultsCollection
MonitorElement * allHad_phi0
virtual void analyzeSemiE(edm::Event const &e, edm::EventSetup const &eSetup)
std::vector< edm::InputTag > jetLabels_
double semiE_LepJetPtCut_
edm::EDGetTokenT< edm::View< reco::Muon > > muonToken_
MonitorElement * allHad_pt1
double allHadDeltaPhiCut_
MonitorElement * semiMu_muPhi
std::vector< edm::EDGetTokenT< edm::View< reco::Jet > > > jetTokens_
MonitorElement * semiE_hadJetY
edm::EDGetTokenT< edm::TriggerResults > triggerToken_
double mass() const final
mass
float neutralHadronEnergyFraction() const
neutralHadronEnergyFraction
std::vector< MonitorElement * > pfJet_m
MonitorElement * semiE_hadJetPt
const LorentzVector & p4() const final
four-momentum Lorentz vector
std::vector< MonitorElement * > pfJet_pt
std::vector< MonitorElement * > boostedJet_massDrop
edm::EDGetTokenT< edm::View< reco::BasicJet > > sdjetToken_
std::vector< MonitorElement * > boostedJet_subjetN
edm::InputTag sdjetLabel_
float neutralEmEnergyFraction() const
neutralEmEnergyFraction
MonitorElement * allHad_mass0
MonitorElement * semiMu_hadJetPhi
size_t numberOfDaughters() const override
number of daughters
edm::EDGetTokenT< std::vector< reco::PFMET > > PFMETToken_
double semiMu_LepJetPtCut_
MonitorElement * semiMu_hadJetWMass
double phi() const final
momentum azimuthal angle
MonitorElement * semiE_hadJetDR
std::vector< MonitorElement * > boostedJet_subjetM
MonitorElement * allHad_mass1
std::vector< MonitorElement * > pfJet_nhef
Jets made from PFObjects.
MonitorElement * semiMu_hadJetMass
const Candidate * daughter(size_type) const override
return daughter at a given position, i = 0, ... numberOfDaughters() - 1 (read only mode)
constexpr auto deltaR(const T1 &t1, const T2 &t2) -> decltype(t1.eta())
T getParameter(std::string const &) const
edm::EDGetTokenT< edm::View< reco::GsfElectron > > electronToken_
boost::indirect_iterator< typename seq_t::const_iterator > const_iterator
MonitorElement * semiMu_muEta
MonitorElement * semiMu_muPtRel
MonitorElement * semiE_hadJetPhi
virtual void analyzeSemiMu(edm::Event const &e, edm::EventSetup const &eSetup)
double semiMu_HadJetPtCut_
std::vector< MonitorElement * > pfJet_phi
Abs< T >::type abs(const T &t)
double rapidity() const final
rapidity
const_iterator end() const
double semiE_HadJetPtCut_
MonitorElement * allHad_mttbar
MonitorElement * semiMu_muPt
MonitorElement * semiE_ePtRel
virtual void analyzeJets(edm::Event const &e, edm::EventSetup const &eSetup)
virtual double phi() const =0
momentum azimuthal angle
MonitorElement * semiMu_hadJetPt
std::shared_ptr< StringCutObjectSelector< reco::GsfElectron > > elecSelect_
double allHadRapidityCut_
MonitorElement * semiE_ePhi
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())
std::vector< MonitorElement * > pfJet_nemf