CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
QcdPhotonsDQM.cc
Go to the documentation of this file.
1 /*
2  * See header file for a description of this class.
3  *
4  * \author Michael B. Anderson, University of Wisconsin Madison
5  */
6 
7 #include "DQM/Physics/src/QcdPhotonsDQM.h"
8 
9 #include "FWCore/Framework/interface/Event.h"
10 #include "FWCore/Framework/interface/EventSetup.h"
11 #include "FWCore/ParameterSet/interface/ParameterSet.h"
12 
13 #include "DQMServices/Core/interface/DQMStore.h"
14 #include "DQMServices/Core/interface/MonitorElement.h"
15 #include "FWCore/ServiceRegistry/interface/Service.h"
16 #include "FWCore/Common/interface/TriggerNames.h"
17 
18 #include "DataFormats/Common/interface/Handle.h"
19 
20 #include "FWCore/MessageLogger/interface/MessageLogger.h"
21 
22 // Physics Objects
23 #include "DataFormats/EgammaCandidates/interface/Photon.h"
24 #include "DataFormats/JetReco/interface/Jet.h"
25 
26 // Vertex
27 #include "DataFormats/VertexReco/interface/Vertex.h"
28 
29 // For removing ECAL Spikes
30 #include "RecoEcal/EgammaCoreTools/interface/EcalClusterLazyTools.h"
31 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgo.h"
32 #include "RecoLocalCalo/EcalRecAlgos/interface/EcalSeverityLevelAlgoRcd.h"
33 #include "FWCore/ServiceRegistry/interface/Service.h"
34 
35 #include "FWCore/Framework/interface/ESHandle.h"
36 #include "CondFormats/EcalObjects/interface/EcalCondObjectContainer.h"
37 
38 //geometry
39 //#include "Geometry/CaloGeometry/interface/CaloGeometry.h"
40 //#include "Geometry/Records/interface/IdealGeometryRecord.h"
41 //#include "Geometry/Records/interface/CaloGeometryRecord.h"
42 //#include "Geometry/CaloEventSetup/interface/CaloTopologyRecord.h"
43 //#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
44 //#include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
45 //#include "Geometry/EcalAlgo/interface/EcalPreshowerGeometry.h"
46 
47 // Math stuff
48 #include "DataFormats/Math/interface/deltaR.h"
49 #include "DataFormats/Math/interface/deltaPhi.h"
50 
51 #include <vector>
52 
53 #include <string>
54 #include <cmath>
55 using namespace std;
56 using namespace edm;
57 using namespace reco;
58 
59 
60 QcdPhotonsDQM::QcdPhotonsDQM(const ParameterSet& parameters) {
61  // Get parameters from configuration file
62  theTriggerPathToPass_ = parameters.getParameter<string>("triggerPathToPass");
63  thePlotTheseTriggersToo_ = parameters.getParameter<vector<string> >("plotTheseTriggersToo");
64  theJetCollectionLabel_ = parameters.getParameter<InputTag>("jetCollection");
65  trigTagToken_ = consumes<edm::TriggerResults>(
66  parameters.getUntrackedParameter<edm::InputTag>("trigTag"));
67  thePhotonCollectionToken_ = consumes<reco::PhotonCollection>(
68  parameters.getParameter<InputTag>("photonCollection"));
69  theJetCollectionToken_ = consumes<edm::View<reco::Jet> >(
70  parameters.getParameter<InputTag>("jetCollection"));
71  theVertexCollectionToken_ = consumes<reco::VertexCollection>(
72  parameters.getParameter<InputTag>("vertexCollection"));
73  theMinJetPt_ = parameters.getParameter<double>("minJetPt");
74  theMinPhotonEt_ = parameters.getParameter<double>("minPhotonEt");
75  theRequirePhotonFound_ = parameters.getParameter<bool>("requirePhotonFound");
76  thePlotPhotonMaxEt_ = parameters.getParameter<double>("plotPhotonMaxEt");
77  thePlotPhotonMaxEta_ = parameters.getParameter<double>("plotPhotonMaxEta");
78  thePlotJetMaxEta_ = parameters.getParameter<double>("plotJetMaxEta");
79  theBarrelRecHitTag_ = parameters.getParameter<InputTag>("barrelRecHitTag");
80  theEndcapRecHitTag_ = parameters.getParameter<InputTag>("endcapRecHitTag");
81  theBarrelRecHitToken_ = consumes<EcalRecHitCollection>(
82  parameters.getParameter<InputTag>("barrelRecHitTag"));
83  theEndcapRecHitToken_ = consumes<EcalRecHitCollection>(
84  parameters.getParameter<InputTag>("endcapRecHitTag"));
85  // just to initialize
86  isValidHltConfig_ = false;
87 
88  // coverity says...
89  h_deltaEt_photon_jet = 0;
90  h_deltaPhi_jet_jet2 = 0;
91  h_deltaPhi_photon_jet = 0;
92  h_deltaPhi_photon_jet2 = 0;
93  h_deltaR_jet_jet2 = 0;
94  h_deltaR_photon_jet2 = 0;
95  h_jet2_eta = 0;
96  h_jet2_pt = 0;
97  h_jet2_ptOverPhotonEt = 0;
98  h_jet_count = 0;
99  h_jet_eta = 0;
100  h_jet_pt = 0;
101  h_photon_count_bar = 0;
102  h_photon_count_end = 0;
103  h_photon_et = 0;
104  h_photon_et_beforeCuts = 0;
105  h_photon_et_jetco = 0;
106  h_photon_et_jetcs = 0;
107  h_photon_et_jetfo = 0;
108  h_photon_et_jetfs = 0;
109  h_photon_et_ratio_co_cs = 0;
110  h_photon_et_ratio_co_fo = 0;
111  h_photon_et_ratio_co_fs = 0;
112  h_photon_et_ratio_cs_fo = 0;
113  h_photon_et_ratio_cs_fs = 0;
114  h_photon_et_ratio_fo_fs = 0;
115  h_photon_eta = 0;
116  h_triggers_passed = 0;
117 
118  theDbe = Service<DQMStore>().operator->();
119 
120 }
121 
122 QcdPhotonsDQM::~QcdPhotonsDQM() {
123 }
124 
125 
126 void QcdPhotonsDQM::beginJob() {
127 
128  logTraceName = "QcdPhotonAnalyzer";
129 
130  LogTrace(logTraceName)<<"Parameters initialization";
131 
132  theDbe->setCurrentFolder("Physics/QcdPhotons"); // Use folder with name of PAG
133 
134  std::stringstream aStringStream;
135  std::string aString;
136  aStringStream << theMinJetPt_;
137  aString = aStringStream.str();
138 
139  // Monitor of triggers passed
140  int numOfTriggersToMonitor = thePlotTheseTriggersToo_.size();
141  h_triggers_passed = theDbe->book1D("triggers_passed", "Events passing these trigger paths", numOfTriggersToMonitor, 0, numOfTriggersToMonitor);
142  for (int i=0; i<numOfTriggersToMonitor; i++) {
143  h_triggers_passed->setBinLabel(i+1,thePlotTheseTriggersToo_[i]);
144  }
145 
146  // Keep the number of plots and number of bins to a minimum!
147  h_photon_et_beforeCuts = theDbe->book1D("photon_et_beforeCuts", "#gamma with highest E_{T};E_{T}(#gamma) (GeV)", 20, 0., thePlotPhotonMaxEt_);
148  h_photon_et = theDbe->book1D("photon_et", "#gamma with highest E_{T};E_{T}(#gamma) (GeV)", 20, 0., thePlotPhotonMaxEt_);
149  h_photon_eta = theDbe->book1D("photon_eta", "#gamma with highest E_{T};#eta(#gamma)", 40, -thePlotPhotonMaxEta_, thePlotPhotonMaxEta_);
150  h_photon_count_bar = theDbe->book1D("photon_count_bar","Number of #gamma's passing selection (Barrel);Number of #gamma's", 8, -0.5, 7.5);
151  h_photon_count_end = theDbe->book1D("photon_count_end","Number of #gamma's passing selection (Endcap);Number of #gamma's", 8, -0.5, 7.5);
152 
153  h_jet_pt = theDbe->book1D("jet_pt", "Jet with highest p_{T} (from "+theJetCollectionLabel_.label()+");p_{T}(1^{st} jet) (GeV)", 20, 0., thePlotPhotonMaxEt_);
154  h_jet_eta = theDbe->book1D("jet_eta", "Jet with highest p_{T} (from "+theJetCollectionLabel_.label()+");#eta(1^{st} jet)", 20, -thePlotJetMaxEta_, thePlotJetMaxEta_);
155  h_deltaPhi_photon_jet = theDbe->book1D("deltaPhi_photon_jet", "#Delta#phi between Highest E_{T} #gamma and jet;#Delta#phi(#gamma,1^{st} jet)", 20, 0, 3.1415926);
156  h_deltaPhi_jet_jet2 = theDbe->book1D("deltaPhi_jet_jet2", "#Delta#phi between Highest E_{T} jet and 2^{nd} jet;#Delta#phi(1^{st} jet,2^{nd} jet)", 20, 0, 3.1415926);
157  h_deltaEt_photon_jet = theDbe->book1D("deltaEt_photon_jet", "(E_{T}(#gamma)-p_{T}(jet))/E_{T}(#gamma) when #Delta#phi(#gamma,1^{st} jet) > 2.8;#DeltaE_{T}(#gamma,1^{st} jet)/E_{T}(#gamma)", 20, -1.0, 1.0);
158  h_jet_count = theDbe->book1D("jet_count", "Number of "+theJetCollectionLabel_.label()+" (p_{T} > "+aString+" GeV);Number of Jets", 8, -0.5, 7.5);
159  h_jet2_pt = theDbe->book1D("jet2_pt", "Jet with 2^{nd} highest p_{T} (from "+theJetCollectionLabel_.label()+");p_{T}(2^{nd} jet) (GeV)", 20, 0., thePlotPhotonMaxEt_);
160  h_jet2_eta = theDbe->book1D("jet2_eta", "Jet with 2^{nd} highest p_{T} (from "+theJetCollectionLabel_.label()+");#eta(2^{nd} jet)", 20, -thePlotJetMaxEta_, thePlotJetMaxEta_);
161  h_jet2_ptOverPhotonEt = theDbe->book1D("jet2_ptOverPhotonEt", "p_{T}(2^{nd} highest jet) / E_{T}(#gamma);p_{T}(2^{nd} Jet)/E_{T}(#gamma)", 20, 0.0, 4.0);
162  h_deltaPhi_photon_jet2 = theDbe->book1D("deltaPhi_photon_jet2","#Delta#phi between Highest E_{T} #gamma and 2^{nd} highest jet;#Delta#phi(#gamma,2^{nd} jet)", 20, 0, 3.1415926);
163  h_deltaR_jet_jet2 = theDbe->book1D("deltaR_jet_jet2", "#DeltaR between Highest Jet and 2^{nd} Highest;#DeltaR(1^{st} jet,2^{nd} jet)", 30, 0, 6.0);
164  h_deltaR_photon_jet2 = theDbe->book1D("deltaR_photon_jet2", "#DeltaR between Highest E_{T} #gamma and 2^{nd} jet;#DeltaR(#gamma, 2^{nd} jet)", 30, 0, 6.0);
165 
166  // Photon Et for different jet configurations
167  Float_t bins_et[] = {15,20,30,50,80};
168  int num_bins_et = 4;
169  h_photon_et_jetcs = theDbe->book1D("photon_et_jetcs", "#gamma with highest E_{T} (#eta(jet)<1.45, #eta(#gamma)#eta(jet)>0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
170  h_photon_et_jetco = theDbe->book1D("photon_et_jetco", "#gamma with highest E_{T} (#eta(jet)<1.45, #eta(#gamma)#eta(jet)<0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
171  h_photon_et_jetfs = theDbe->book1D("photon_et_jetfs", "#gamma with highest E_{T} (1.55<#eta(jet)<2.5, #eta(#gamma)#eta(jet)>0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
172  h_photon_et_jetfo = theDbe->book1D("photon_et_jetfo", "#gamma with highest E_{T} (1.55<#eta(jet)<2.5, #eta(#gamma)#eta(jet)<0);E_{T}(#gamma) (GeV)", num_bins_et, bins_et);
173  h_photon_et_jetcs->getTH1F()->Sumw2();
174  h_photon_et_jetco->getTH1F()->Sumw2();
175  h_photon_et_jetfs->getTH1F()->Sumw2();
176  h_photon_et_jetfo->getTH1F()->Sumw2();
177  // Ratio of the above Photon Et distributions
178  h_photon_et_ratio_co_cs = theDbe->book1D("photon_et_ratio_00_co_cs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
179  h_photon_et_ratio_fo_fs = theDbe->book1D("photon_et_ratio_01_fo_fs", "D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
180  h_photon_et_ratio_cs_fs = theDbe->book1D("photon_et_ratio_02_cs_fs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
181  h_photon_et_ratio_co_fs = theDbe->book1D("photon_et_ratio_03_co_fs", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)>0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
182  h_photon_et_ratio_cs_fo = theDbe->book1D("photon_et_ratio_04_cs_fo", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)>0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
183  h_photon_et_ratio_co_fo = theDbe->book1D("photon_et_ratio_05_co_fo", "D(|#eta(jet)|<1.45, #eta(jet)*#eta(#gamma)<0) / D(1.55<|#eta(jet)|<2.6, #eta(jet)*#eta(#gamma)<0);E_{T}(#gamma) (GeV); ratio", num_bins_et, bins_et);
184  h_photon_et_ratio_co_cs->getTH1F()->Sumw2();
185  h_photon_et_ratio_fo_fs->getTH1F()->Sumw2();
186  h_photon_et_ratio_cs_fs->getTH1F()->Sumw2();
187  h_photon_et_ratio_co_fs->getTH1F()->Sumw2();
188  h_photon_et_ratio_cs_fo->getTH1F()->Sumw2();
189  h_photon_et_ratio_co_fo->getTH1F()->Sumw2();
190 }
191 
192 
196 void QcdPhotonsDQM::beginRun( const edm::Run &iRun, const edm::EventSetup &iSet ) {
197 
198  // passed as parameter to HLTConfigProvider::init(), not yet used
199  bool isConfigChanged = false;
200 
201  // isValidHltConfig_ could be used to short-circuit analyze() in case of problems
202  isValidHltConfig_ = hltConfigProvider_.init( iRun, iSet, "HLT", isConfigChanged );
203 
204  num_events_in_run = 0;
205 }
206 
207 
208 void QcdPhotonsDQM::analyze(const Event& iEvent, const EventSetup& iSetup) {
209  num_events_in_run++;
210 
212  //if( ! isValidHltConfig_ ) return;
213 
214  LogTrace(logTraceName)<<"Analysis of event # ";
215 
217  // Did event pass HLT paths?
218  Handle<TriggerResults> HLTresults;
219  iEvent.getByToken(trigTagToken_, HLTresults);
220  if (!HLTresults.isValid()) {
221  //LogWarning("") << ">>> TRIGGER collection does not exist !!!";
222  return;
223  }
224  const edm::TriggerNames & trigNames = iEvent.triggerNames(*HLTresults);
225 
226  bool passed_HLT=false;
227 
228  // See if event passed trigger paths
229  // increment that bin in the trigger plot
230  for (unsigned int i=0; i<thePlotTheseTriggersToo_.size(); i++) {
231  passed_HLT = false;
232  for (unsigned int ti=0; (ti<trigNames.size()) && !passed_HLT; ++ti) {
233  size_t pos = trigNames.triggerName(ti).find(thePlotTheseTriggersToo_[i]);
234  if (pos==0) passed_HLT = HLTresults->accept(ti);
235  }
236  if (passed_HLT) h_triggers_passed->Fill(i);
237  }
238 
239 
240  // grab photons
241  Handle<PhotonCollection> photonCollection;
242  iEvent.getByToken(thePhotonCollectionToken_, photonCollection);
243 
244  // If photon collection is empty, exit
245  if (!photonCollection.isValid()) return;
246 
247 
248  // Quit if the event did not pass the HLT path we care about
249  passed_HLT = false;
250  {
251  //bool found=false;
252  for (unsigned int ti=0; ti<trigNames.size(); ++ti) {
253  size_t pos = trigNames.triggerName(ti).find(theTriggerPathToPass_);
254  if (pos==0) {
255  passed_HLT = HLTresults->accept(ti);
256  //found=true;
257  break;
258  }
259  }
260 
261  // Assumption: reco photons are ordered by Et
262  for (PhotonCollection::const_iterator recoPhoton = photonCollection->begin(); recoPhoton!=photonCollection->end(); recoPhoton++){
263  // stop looping over photons once we get to too low Et
264  if ( recoPhoton->et() < theMinPhotonEt_ ) break;
265 
266  h_photon_et_beforeCuts->Fill(recoPhoton->et());
267  break; // leading photon only
268  }
269 
270  if (!passed_HLT) {
271  return;
272  }
273  }
274 
276 
277  //std::cout << "\tpassed main trigger (" << theTriggerPathToPass_ << ")" << std::endl;
278 
280  // Does event have valid vertex?
281  // Get the primary event vertex
282  Handle<VertexCollection> vertexHandle;
283  iEvent.getByToken(theVertexCollectionToken_, vertexHandle);
284  VertexCollection vertexCollection = *(vertexHandle.product());
285  //double vtx_ndof = -1.0;
286  //double vtx_z = 0.0;
287  //bool vtx_isFake = true;
288  //if (vertexCollection.size()>0) {
289  // vtx_ndof = vertexCollection.begin()->ndof();
290  // vtx_z = vertexCollection.begin()->z();
291  // vtx_isFake = false;
292  //}
293  //if (vtx_isFake || fabs(vtx_z)>15 || vtx_ndof<4) return;
294 
295  int nvvertex = 0;
296  for (unsigned int i=0; i<vertexCollection.size(); ++i) {
297  if (vertexCollection[i].isValid()) nvvertex++;
298  }
299  if (nvvertex==0) return;
300 
302 
303  //std::cout << "\tpassed vertex selection" << std::endl;
304 
305 
307  // Did the event pass certain L1 Technical Trigger bits?
308  // It's probably beam halo
309  // TODO: ADD code
311 
312 
313  // For finding spikes
314  Handle<EcalRecHitCollection> EBReducedRecHits;
315  iEvent.getByToken(theBarrelRecHitToken_, EBReducedRecHits);
316  Handle<EcalRecHitCollection> EEReducedRecHits;
317  iEvent.getByToken(theEndcapRecHitToken_, EEReducedRecHits);
318  EcalClusterLazyTools lazyTool(iEvent, iSetup, theBarrelRecHitTag_, theEndcapRecHitTag_);
319 
320 
321  // Find the highest et "decent" photon
322  float photon_et = -9.0;
323  float photon_eta = -9.0;
324  float photon_phi = -9.0;
325  bool photon_passPhotonID = false;
326  bool found_lead_pho = false;
327  int photon_count_bar = 0;
328  int photon_count_end = 0;
329  // Assumption: reco photons are ordered by Et
330  for (PhotonCollection::const_iterator recoPhoton = photonCollection->begin(); recoPhoton!=photonCollection->end(); recoPhoton++){
331 
332  // stop looping over photons once we get to too low Et
333  if ( recoPhoton->et() < theMinPhotonEt_ ) break;
334 
335  /*
336  // Ignore ECAL Spikes
337  const reco::CaloClusterPtr seed = recoPhoton->superCluster()->seed();
338  DetId id = lazyTool.getMaximum(*seed).first; // Cluster shape variables
339  // float time = -999., outOfTimeChi2 = -999., chi2 = -999.; // UNUSED
340  int flags=-1, severity = -1;
341  const EcalRecHitCollection & rechits = ( recoPhoton->isEB() ? *EBReducedRecHits : *EEReducedRecHits);
342  EcalRecHitCollection::const_iterator it = rechits.find( id );
343  if( it != rechits.end() ) {
344  // time = it->time(); // UNUSED
345  // outOfTimeChi2 = it->outOfTimeChi2(); // UNUSED
346  // chi2 = it->chi2(); // UNUSED
347  flags = it->recoFlag();
348 
349  edm::ESHandle<EcalSeverityLevelAlgo> sevlv;
350  iSetup.get<EcalSeverityLevelAlgoRcd>().get(sevlv);
351  severity = sevlv->severityLevel( id, rechits);
352  }
353  bool isNotSpike = ((recoPhoton->isEB() && (severity!=3 && severity!=4 ) && (flags != 2) ) || recoPhoton->isEE());
354  if (!isNotSpike) continue; // move on to next photon
355  // END of determining ECAL Spikes
356  */
357 
358  bool pho_current_passPhotonID = false;
359  bool pho_current_isEB = recoPhoton->isEB();
360  bool pho_current_isEE = recoPhoton->isEE();
361 
362  if ( pho_current_isEB && (recoPhoton->sigmaIetaIeta() < 0.01 || recoPhoton->hadronicOverEm() < 0.05) ) {
363  // Photon object in barrel passes photon ID
364  pho_current_passPhotonID = true;
365  photon_count_bar++;
366  } else if ( pho_current_isEE && (recoPhoton->hadronicOverEm() < 0.05) ) {
367  // Photon object in endcap passes photon ID
368  pho_current_passPhotonID = true;
369  photon_count_end++;
370  }
371 
372  if (!found_lead_pho) {
373  found_lead_pho = true;
374  photon_passPhotonID = pho_current_passPhotonID;
375  photon_et = recoPhoton->et();
376  photon_eta = recoPhoton->eta();
377  photon_phi = recoPhoton->phi();
378  }
379  }
380 
381  // If user requires a photon to be found, but none is, return.
382  // theRequirePhotonFound should pretty much always be set to 'True'
383  // except when running on qcd monte carlo just to see the jets.
384  if ( theRequirePhotonFound_ && (!photon_passPhotonID || photon_et<theMinPhotonEt_) ) return;
385 
386 
388  // Find the highest et jet
389  Handle<View<Jet> > jetCollection;
390  iEvent.getByToken (theJetCollectionToken_,jetCollection);
391  if (!jetCollection.isValid()) return;
392 
393  float jet_pt = -8.0;
394  float jet_eta = -8.0;
395  float jet_phi = -8.0;
396  int jet_count = 0;
397  float jet2_pt = -9.0;
398  float jet2_eta = -9.0;
399  float jet2_phi = -9.0;
400  // Assumption: jets are ordered by Et
401  for (unsigned int i_jet = 0; i_jet < jetCollection->size(); i_jet++) {
402  const Jet* jet = & jetCollection->at(i_jet);
403 
404  float jet_current_pt = jet->pt();
405 
406  // don't care about jets that overlap with the lead photon
407  if ( deltaR(jet->eta(), jet->phi(), photon_eta, photon_phi) < 0.5 ) continue;
408  // stop looping over jets once we get to too low Et
409  if (jet_current_pt < theMinJetPt_) break;
410 
411  jet_count++;
412  if (jet_current_pt > jet_pt) {
413  jet2_pt = jet_pt; // 2nd highest jet get's et from current highest
414  jet2_eta = jet_eta;
415  jet2_phi = jet_phi;
416  jet_pt = jet_current_pt; // current highest jet gets et from the new highest
417  jet_eta = jet->eta();
418  jet_phi = jet->phi();
419  } else if (jet_current_pt > jet2_pt) {
420  jet2_pt = jet_current_pt;
421  jet2_eta = jet->eta();
422  jet2_phi = jet->phi();
423  }
424  }
426 
427 
429  // Fill histograms if a jet found
430  // NOTE: if a photon was required to be found, but wasn't
431  // we wouldn't have made it to this point in the code
432  if ( jet_pt > 0.0 ) {
433 
434  // Photon Plots
435  h_photon_et ->Fill( photon_et );
436  h_photon_eta ->Fill( photon_eta );
437  h_photon_count_bar->Fill( photon_count_bar );
438  h_photon_count_end->Fill( photon_count_end );
439 
440  // Photon Et hists for different orientations to the jet
441  if ( fabs(photon_eta)<1.45 && photon_passPhotonID ) { // Lead photon is in barrel
442  if (fabs(jet_eta)<1.45){ // jet is in barrel
443  if (photon_eta*jet_eta>0) {
444  h_photon_et_jetcs->Fill(photon_et);
445  } else {
446  h_photon_et_jetco->Fill(photon_et);
447  }
448  } else if (jet_eta>1.55 && jet_eta<2.5) { // jet is in endcap
449  if (photon_eta*jet_eta>0) {
450  h_photon_et_jetfs->Fill(photon_et);
451  } else {
452  h_photon_et_jetfo->Fill(photon_et);
453  }
454  }
455  } // END of Lead Photon is in Barrel
456 
457  // Jet Plots
458  h_jet_pt ->Fill( jet_pt );
459  h_jet_eta ->Fill( jet_eta );
460  h_jet_count ->Fill( jet_count );
461  h_deltaPhi_photon_jet ->Fill( abs(deltaPhi(photon_phi, jet_phi)) );
462  if ( abs(deltaPhi(photon_phi,jet_phi))>2.8 ) h_deltaEt_photon_jet->Fill( (photon_et-jet_pt)/photon_et );
463 
464  // 2nd Highest Jet Plots
465  if ( jet2_pt > 0.0 ) {
466  h_jet2_pt ->Fill( jet2_pt );
467  h_jet2_eta ->Fill( jet2_eta );
468  h_jet2_ptOverPhotonEt ->Fill( jet2_pt/photon_et );
469  h_deltaPhi_photon_jet2->Fill( abs(deltaPhi(photon_phi, jet2_phi)) );
470  h_deltaPhi_jet_jet2 ->Fill( abs(deltaPhi( jet_phi, jet2_phi)) );
471  h_deltaR_jet_jet2 ->Fill( deltaR( jet_eta, jet_phi, jet2_eta, jet2_phi) );
472  h_deltaR_photon_jet2 ->Fill( deltaR(photon_eta, photon_phi, jet2_eta, jet2_phi) );
473  }
474  }
475  // End of Filling histograms
477 }
478 
479 
480 void QcdPhotonsDQM::endJob(void) {}
481 
482 void QcdPhotonsDQM::endRun(const edm::Run& run, const edm::EventSetup& es) {
483  if (num_events_in_run>0) {
484  h_triggers_passed->getTH1F()->Scale(1.0/num_events_in_run);
485  }
486  h_photon_et_ratio_co_cs->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetcs->getTH1F() );
487  h_photon_et_ratio_fo_fs->getTH1F()->Divide( h_photon_et_jetfo->getTH1F(), h_photon_et_jetfs->getTH1F() );
488  h_photon_et_ratio_cs_fs->getTH1F()->Divide( h_photon_et_jetcs->getTH1F(), h_photon_et_jetfs->getTH1F() );
489  h_photon_et_ratio_co_fs->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetfs->getTH1F() );
490  h_photon_et_ratio_cs_fo->getTH1F()->Divide( h_photon_et_jetcs->getTH1F(), h_photon_et_jetfo->getTH1F() );
491  h_photon_et_ratio_co_fo->getTH1F()->Divide( h_photon_et_jetco->getTH1F(), h_photon_et_jetfo->getTH1F() );
492 }
493 
494 
495 // Local Variables:
496 // show-trailing-whitespace: t
497 // truncate-lines: t
498 // End:
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
Parameters.parameters
dictionary parameters
Definition: Parameters.py:2
edm::Event::triggerNames
virtual edm::TriggerNames const & triggerNames(edm::TriggerResults const &triggerResults) const
Definition: Event.cc:204
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
reco::Jet
Base class for all types of Jets.
Definition: Jet.h:20
QcdPhotonsDQM::beginJob
void beginJob()
Inizialize parameters for histo binning.
Definition: QcdPhotonsDQM.cc:126
reco::VertexCollection
std::vector< Vertex > VertexCollection
collection of Vertex objects
Definition: VertexFwd.h:9
edm::TriggerNames::size
Strings::size_type size() const
Definition: TriggerNames.cc:39
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:105
GoodVertex_cfg.vertexCollection
tuple vertexCollection
Definition: GoodVertex_cfg.py:61
QcdPhotonsDQM::endRun
void endRun(const edm::Run &, const edm::EventSetup &)
Definition: QcdPhotonsDQM.cc:482
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
QcdPhotonsDQM::QcdPhotonsDQM
QcdPhotonsDQM(const edm::ParameterSet &)
Constructor.
Definition: QcdPhotonsDQM.cc:60
reco::LeafCandidate::phi
virtual float phi() const GCC11_FINAL
momentum azimuthal angle
Definition: LeafCandidate.h:155
iEvent
int iEvent
Definition: GenABIO.cc:243
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
reco::LeafCandidate::eta
virtual float eta() const GCC11_FINAL
momentum pseudorapidity
Definition: LeafCandidate.h:159
QcdPhotonsDQM::endJob
void endJob(void)
Save the histos.
Definition: QcdPhotonsDQM.cc:480
deltaR
double deltaR(double eta1, double eta2, double phi1, double phi2)
Definition: TreeUtility.cc:17
trigNames
static const char *const trigNames[]
Definition: EcalDumpRaw.cc:74
SiPixelRawToDigiRegional_cfi.deltaPhi
tuple deltaPhi
Definition: SiPixelRawToDigiRegional_cfi.py:9
edm::TriggerNames::triggerName
std::string const & triggerName(unsigned int index) const
Definition: TriggerNames.cc:27
QcdPhotonsDQM::beginRun
void beginRun(const edm::Run &, const edm::EventSetup &)
Definition: QcdPhotonsDQM.cc:196
edm::Handle::product
T const * product() const
Definition: Handle.h:81
QcdPhotonsDQM::~QcdPhotonsDQM
virtual ~QcdPhotonsDQM()
Destructor.
Definition: QcdPhotonsDQM.cc:122
reco::LeafCandidate::pt
virtual float pt() const GCC11_FINAL
transverse momentum
Definition: LeafCandidate.h:153
QcdPhotonsDQM::analyze
void analyze(const edm::Event &, const edm::EventSetup &)
Get the analysis.
Definition: QcdPhotonsDQM.cc:208
edm::Run
Definition: Run.h:41
Generated for CMSSW Reference Manual by doxygen 1.8.5