|
|
Go to the documentation of this file.
41 #include <TLorentzVector.h>
111 verbose_ =
pset.getUntrackedParameter<
bool>(
"verbose",
false);
128 edm::LogError(
"HLTPhotonJetVal") <<
"Initialization of HLTConfigProvider failed!!";
140 nvertices_ = iBooker.
book1D(
"nvertices",
"Number of vertices", 100, 0, 100);
142 nphotons_ = iBooker.
book1D(
"nphotons",
"Number of photons", 100, 0, 10);
148 pfmet_ = iBooker.
book1D(
"pfmet",
"PF MET", 100, 0, 250);
150 pfmetphi_ = iBooker.
book1D(
"pfmetphi",
"PF MET phi;#phi_{PFMET}", 100, -4, 4);
152 iBooker.
book1D(
"delphiphomet_reco",
"Reco: #Delta#phi(photon, MET);#Delta#phi(#gamma,MET)", 100, 0, 4);
153 delphiphomet_ = iBooker.
book1D(
"delphiphomet",
"#Delta#phi(photon, MET);#Delta#phi(#gamma,MET)", 100, 0, 4);
155 npfjets_ = iBooker.
book1D(
"npfjets",
"Number of PF Jets", 100, 0, 20);
157 iBooker.
book1D(
"delphijetmet_reco",
"Reco: #Delta#phi(PFJet, MET);#Delta#phi(Jet,MET)", 100, 0, 4);
158 delphijetmet_ = iBooker.
book1D(
"delphijetmet",
"#Delta#phi(PFJet, MET);#Delta#phi(Jet,MET)", 100, 0, 4);
159 invmassjj_reco_ = iBooker.
book1D(
"invmassjj_reco",
"Reco: Inv mass two leading jets;M_{jj}[GeV]", 100, 0, 2000);
160 invmassjj_ = iBooker.
book1D(
"invmassjj",
"Inv mass two leading jets;M_{jj}[GeV]", 100, 0, 2000);
161 deletajj_reco_ = iBooker.
book1D(
"deletajj_reco",
"Reco: #Delta#eta(jj);|#Delta#eta_{jj}|", 100, 0, 6);
162 deletajj_ = iBooker.
book1D(
"deletajj",
"#Delta#eta(jj);|#Delta#eta_{jj}|", 100, 0, 6);
167 "Reco: Trigger vs. # vertices;N_{vertices};Trigger",
175 "Trigger vs. # vertices;N_{vertices};Trigger",
191 edm::LogError(
"HigPhotonJetHLT") <<
"Missing triggerResults collection" << std::endl;
218 for (
unsigned int itrig = 0; itrig <
triggerResults->size(); itrig++) {
246 std::cout <<
"xshi:: number of pfmets: " << pfmets->size() << std::endl;
258 for (
auto const& phoIter : *
photons) {
268 double tmp_delphiphomet = fabs(
deltaPhi(phoIter.phi(), pfmet.
phi()));
283 std::cout <<
"xshi:: N pfjets : " << pfjets->size() << std::endl;
285 double min_delphijetmet = 6.0;
286 TLorentzVector p4jet1, p4jet2, p4jj;
288 double etajet1(0), etajet2(0);
290 for (
auto const& jetIter : *pfjets) {
295 double tmp_delphijetmet = fabs(
deltaPhi(jetIter.phi(), pfmet.
phi()));
296 if (tmp_delphijetmet < min_delphijetmet)
297 min_delphijetmet = tmp_delphijetmet;
300 p4jet1.SetXYZM(jetIter.px(), jetIter.py(), jetIter.pz(), jetIter.mass());
301 etajet1 = jetIter.eta();
304 p4jet2.SetXYZM(jetIter.px(), jetIter.py(), jetIter.pz(), jetIter.mass());
305 etajet2 = jetIter.eta();
315 p4jj = p4jet1 + p4jet2;
316 double deletajj = etajet1 - etajet2;
318 std::cout <<
"xshi:: invmass jj " << p4jj.M() << std::endl;
335 std::cout <<
"xshi:: endRun total number of events: " <<
evtsrun_ <<
", integral = " <<
h->Integral()
336 <<
", norm = " << norm << std::endl;
342 for (
unsigned int itrig = 0; itrig <
triggerResults->size(); itrig++) {
348 if (triggername.find(
i) != std::string::npos) {
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
MonitorElement * trigvsnvtx_
MonitorElement * nphotons_reco_
edm::EDGetTokenT< reco::PFJetCollection > pfJetsToken_
MonitorElement * deletajj_reco_
void dqmEndRun(const edm::Run &, const edm::EventSetup &) override
edm::EDGetTokenT< edm::TriggerResults > triggerResultsToken_
MonitorElement * photonpt_
MonitorElement * pfmetphi_reco_
void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override
static const std::string triggerResults
virtual void setCurrentFolder(std::string const &fullpath)
edm::EDGetTokenT< reco::VertexCollection > pvToken_
MonitorElement * nvertices_reco_
MonitorElement * pfmetphi_
bool accept(const edm::Event &event, const edm::TriggerResults &triggerTable, const std::string &triggerPath)
MonitorElement * triggers_reco_
std::string hltProcessName_
MonitorElement * nvertices_
#define DEFINE_FWK_MODULE(type)
Integral< F, X >::type integral(const F &f)
bool isMonitoredTriggerAccepted(const edm::TriggerNames &, const edm::Handle< edm::TriggerResults > &)
edm::EDGetTokenT< reco::PhotonCollection > photonsToken_
MonitorElement * delphiphomet_
void analyze(const edm::Event &, const edm::EventSetup &) override
MonitorElement * pfmet_reco_
MonitorElement * invmassjj_reco_
MonitorElement * triggers_
MonitorElement * npfjets_reco_
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
HigPhotonJetHLTOfflineSource(const edm::ParameterSet &)
MonitorElement * delphijetmet_reco_
MonitorElement * nphotons_
Log< level::Error, false > LogError
MonitorElement * delphiphomet_reco_
MonitorElement * delphijetmet_
MonitorElement * trigvsnvtx_reco_
double et() const final
transverse energy
double phi() const final
momentum azimuthal angle
MonitorElement * photonrapidity_reco_
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
MonitorElement * photonpt_reco_
MonitorElement * photonrapidity_
MonitorElement * invmassjj_
MonitorElement * deletajj_
std::vector< std::string > hltPathsToCheck_
MonitorElement * npfjets_
edm::EDGetTokenT< reco::PFMETCollection > pfMetToken_
MonitorElement * book1D(TString const &name, TString const &title, int const nchX, double const lowX, double const highX, FUNC onbooking=NOOP())