CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HCALRecHitAnalyzer.cc
Go to the documentation of this file.
1 #include "DQMOffline/JetMET/interface/HCALRecHitAnalyzer.h"
2 // author: Bobby Scurlock, University of Florida
3 // first version 12/7/2006
4 // modified: Mike Schmitt
5 // date: 03.05.2007
6 // note: 1) code rewrite. 2.) changed to loop over all hcal detids;
7 // not only those within calotowers.
8 
9 #include "FWCore/Framework/interface/MakerMacros.h"
10 #include "FWCore/MessageLogger/interface/MessageLogger.h"
11 #include "FWCore/ServiceRegistry/interface/Service.h"
12 
13 //#include "PluginManager/ModuleDef.h"
14 #include "FWCore/PluginManager/interface/ModuleDef.h"
15 #include "FWCore/Framework/interface/Event.h"
16 #include "FWCore/Framework/interface/ESHandle.h"
17 #include "DataFormats/Common/interface/Handle.h"
18 //#include "FWCore/Framework/interface/Handle.h"
19 #include "FWCore/Framework/interface/EventSetup.h"
20 #include "FWCore/ParameterSet/interface/ParameterSet.h"
21 
22 //#include "Geometry/Vector/interface/GlobalPoint.h"
23 #include "DataFormats/GeometryVector/interface/GlobalPoint.h"
24 #include "Geometry/CaloGeometry/interface/CaloSubdetectorGeometry.h"
25 #include "Geometry/CaloGeometry/interface/CaloGeometry.h"
26 #include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
27 #include "Geometry/Records/interface/CaloGeometryRecord.h"
28 
29 #include "DataFormats/DetId/interface/DetId.h"
30 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
31 #include "DataFormats/HcalDetId/interface/HcalSubdetector.h"
32 
33 #include <memory>
34 #include <vector>
35 #include <utility>
36 #include <ostream>
37 #include <fstream>
38 #include <string>
39 #include <algorithm>
40 #include <cmath>
41 #include <TLorentzVector.h>
42 #include "DQMServices/Core/interface/DQMStore.h"
43 
44 #define DEBUG(X) { if (debug_) { std::cout << X << std::endl; } }
45 
46 HCALRecHitAnalyzer::HCALRecHitAnalyzer(const edm::ParameterSet& iConfig)
47 {
48 
49  // Retrieve Information from the Configuration File
50  hBHERecHitsLabel_ = consumes<HBHERecHitCollection>(iConfig.getParameter<edm::InputTag>("HBHERecHitsLabel"));
51  hORecHitsLabel_ = consumes<HORecHitCollection>(iConfig.getParameter<edm::InputTag>("HORecHitsLabel"));
52  hFRecHitsLabel_ = consumes<HFRecHitCollection>(iConfig.getParameter<edm::InputTag>("HFRecHitsLabel"));
53 
54  debug_ = iConfig.getParameter<bool>("Debug");
55  finebinning_ = iConfig.getUntrackedParameter<bool>("FineBinning");
56  FolderName_ = iConfig.getUntrackedParameter<std::string>("FolderName");
57 
58 }
59 
60 void HCALRecHitAnalyzer::dqmBeginRun(const edm::Run& iRun, const edm::EventSetup& iSetup)
61 {
62  Nevents = 0;
63  FillGeometry(iSetup);
64 }
65 
66 void HCALRecHitAnalyzer::bookHistograms(DQMStore::IBooker & ibooker,
67  edm::Run const & iRun,
68  edm::EventSetup const & ) {
69  // get ahold of back-end interface
70  // dbe_ = edm::Service<DQMStore>().operator->();
71 
72  ibooker.setCurrentFolder(FolderName_+"/geometry");
73  hHCAL_ieta_iphi_HBMap = ibooker.book2D("METTask_HCAL_ieta_iphi_HBMap","",83,-41,42,72,1,73);
74  hHCAL_ieta_iphi_HEMap = ibooker.book2D("METTask_HCAL_ieta_iphi_HEMap","",83,-41,42,72,1,73);
75  hHCAL_ieta_iphi_HFMap = ibooker.book2D("METTask_HCAL_ieta_iphi_HFMap","",83,-41,42,72,1,73);
76  hHCAL_ieta_iphi_HOMap = ibooker.book2D("METTask_HCAL_ieta_iphi_HOMap","",83,-41,42,72,1,73);
77  hHCAL_ieta_iphi_etaMap = ibooker.book2D("METTask_HCAL_ieta_iphi_etaMap","",83,-41,42,72,1,73);
78  hHCAL_ieta_iphi_phiMap = ibooker.book2D("METTask_HCAL_ieta_iphi_phiMap","",83,-41,42,72,1,73);
79  hHCAL_ieta_detaMap = ibooker.book1D("METTask_HCAL_ieta_detaMap","",83,-41,42);
80  hHCAL_ieta_dphiMap = ibooker.book1D("METTask_HCAL_ieta_dphiMap","",83,-41,42);
81 
82  // Initialize bins for geometry to -999 because z = 0 is a valid entry
83  for (int i = 1; i <= 83; i++) {
84 
85  hHCAL_ieta_detaMap->setBinContent(i,-999);
86  hHCAL_ieta_dphiMap->setBinContent(i,-999);
87 
88  for (int j = 1; j <= 72; j++) {
89 
90  hHCAL_ieta_iphi_HBMap->setBinContent(i,j,0);
91  hHCAL_ieta_iphi_HEMap->setBinContent(i,j,0);
92  hHCAL_ieta_iphi_HFMap->setBinContent(i,j,0);
93  hHCAL_ieta_iphi_HOMap->setBinContent(i,j,0);
94  hHCAL_ieta_iphi_etaMap->setBinContent(i,j,-999);
95  hHCAL_ieta_iphi_phiMap->setBinContent(i,j,-999);
96 
97  }
98 
99  }
100 
101  // ibooker.setCurrentFolder("RecoMETV/MET_HCAL/data");
102  ibooker.setCurrentFolder(FolderName_);
103  //--Store number of events used
104  hHCAL_Nevents = ibooker.book1D("METTask_HCAL_Nevents","",1,0,1);
105  //--Data integrated over all events and stored by HCAL(ieta,iphi)
106 
107  hHCAL_D1_energy_ieta_iphi = ibooker.book2D("METTask_HCAL_D1_energy_ieta_iphi","",83,-41,42,72,1,73);
108  hHCAL_D2_energy_ieta_iphi = ibooker.book2D("METTask_HCAL_D2_energy_ieta_iphi","",83,-41,42,72,1,73);
109  hHCAL_D3_energy_ieta_iphi = ibooker.book2D("METTask_HCAL_D3_energy_ieta_iphi","",83,-41,42,72,1,73);
110  hHCAL_D4_energy_ieta_iphi = ibooker.book2D("METTask_HCAL_D4_energy_ieta_iphi","",83,-41,42,72,1,73);
111 
112  hHCAL_D1_Minenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D1_Minenergy_ieta_iphi","",83,-41,42,72,1,73);
113  hHCAL_D2_Minenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D2_Minenergy_ieta_iphi","",83,-41,42,72,1,73);
114  hHCAL_D3_Minenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D3_Minenergy_ieta_iphi","",83,-41,42,72,1,73);
115  hHCAL_D4_Minenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D4_Minenergy_ieta_iphi","",83,-41,42,72,1,73);
116 
117  hHCAL_D1_Maxenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D1_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);
118  hHCAL_D2_Maxenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D2_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);
119  hHCAL_D3_Maxenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D3_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);
120  hHCAL_D4_Maxenergy_ieta_iphi = ibooker.book2D("METTask_HCAL_D4_Maxenergy_ieta_iphi","",83,-41,42,72,1,73);
121 
122  // need to initialize those
123  for (int i=1; i<=83; i++)
124  for (int j=1; j<=73; j++)
125  {
126  hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
127  hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
128  hHCAL_D3_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
129  hHCAL_D4_Maxenergy_ieta_iphi->setBinContent(i,j,-999);
130 
131  hHCAL_D1_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
132  hHCAL_D2_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
133  hHCAL_D3_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
134  hHCAL_D4_Minenergy_ieta_iphi->setBinContent(i,j, 14000);
135  }
136 
137  hHCAL_D1_Occ_ieta_iphi = ibooker.book2D("METTask_HCAL_D1_Occ_ieta_iphi","",83,-41,42,72,1,73);
138  hHCAL_D2_Occ_ieta_iphi = ibooker.book2D("METTask_HCAL_D2_Occ_ieta_iphi","",83,-41,42,72,1,73);
139  hHCAL_D3_Occ_ieta_iphi = ibooker.book2D("METTask_HCAL_D3_Occ_ieta_iphi","",83,-41,42,72,1,73);
140  hHCAL_D4_Occ_ieta_iphi = ibooker.book2D("METTask_HCAL_D4_Occ_ieta_iphi","",83,-41,42,72,1,73);
141  //--Data over eta-rings
142 
143 
144  // CaloTower values
145 
146  if(finebinning_)
147  {
148  hHCAL_D1_energyvsieta = ibooker.book2D("METTask_HCAL_D1_energyvsieta","",83,-41,42,20101,-100,2001);
149  hHCAL_D2_energyvsieta = ibooker.book2D("METTask_HCAL_D2_energyvsieta","",83,-41,42,20101,-100,2001);
150  hHCAL_D3_energyvsieta = ibooker.book2D("METTask_HCAL_D3_energyvsieta","",83,-41,42,20101,-100,2001);
151  hHCAL_D4_energyvsieta = ibooker.book2D("METTask_HCAL_D4_energyvsieta","",83,-41,42,20101,-100,2001);
152 
153  hHCAL_D1_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D1_Minenergyvsieta","",83,-41,42,20101,-100,2001);
154  hHCAL_D2_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D2_Minenergyvsieta","",83,-41,42,20101,-100,2001);
155  hHCAL_D3_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D3_Minenergyvsieta","",83,-41,42,20101,-100,2001);
156  hHCAL_D4_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D4_Minenergyvsieta","",83,-41,42,20101,-100,2001);
157 
158  hHCAL_D1_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D1_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
159  hHCAL_D2_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D2_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
160  hHCAL_D3_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D3_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
161  hHCAL_D4_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D4_Maxenergyvsieta","",83,-41,42,20101,-100,2001);
162 
163  // Integrated over phi
164  hHCAL_D1_Occvsieta = ibooker.book2D("METTask_HCAL_D1_Occvsieta","",83,-41,42,73,0,73);
165  hHCAL_D2_Occvsieta = ibooker.book2D("METTask_HCAL_D2_Occvsieta","",83,-41,42,73,0,73);
166  hHCAL_D3_Occvsieta = ibooker.book2D("METTask_HCAL_D3_Occvsieta","",83,-41,42,73,0,73);
167  hHCAL_D4_Occvsieta = ibooker.book2D("METTask_HCAL_D4_Occvsieta","",83,-41,42,73,0,73);
168 
169  hHCAL_D1_SETvsieta = ibooker.book2D("METTask_HCAL_D1_SETvsieta","",83,-41,42,20001,0,2001);
170  hHCAL_D2_SETvsieta = ibooker.book2D("METTask_HCAL_D2_SETvsieta","",83,-41,42,20001,0,2001);
171  hHCAL_D3_SETvsieta = ibooker.book2D("METTask_HCAL_D3_SETvsieta","",83,-41,42,20001,0,2001);
172  hHCAL_D4_SETvsieta = ibooker.book2D("METTask_HCAL_D4_SETvsieta","",83,-41,42,20001,0,2001);
173 
174  hHCAL_D1_METvsieta = ibooker.book2D("METTask_HCAL_D1_METvsieta","",83,-41,42,20001,0,2001);
175  hHCAL_D2_METvsieta = ibooker.book2D("METTask_HCAL_D2_METvsieta","",83,-41,42,20001,0,2001);
176  hHCAL_D3_METvsieta = ibooker.book2D("METTask_HCAL_D3_METvsieta","",83,-41,42,20001,0,2001);
177  hHCAL_D4_METvsieta = ibooker.book2D("METTask_HCAL_D4_METvsieta","",83,-41,42,20001,0,2001);
178 
179  hHCAL_D1_METPhivsieta = ibooker.book2D("METTask_HCAL_D1_METPhivsieta","",83,-41,42, 80, -4,4);
180  hHCAL_D2_METPhivsieta = ibooker.book2D("METTask_HCAL_D2_METPhivsieta","",83,-41,42, 80, -4,4);
181  hHCAL_D3_METPhivsieta = ibooker.book2D("METTask_HCAL_D3_METPhivsieta","",83,-41,42, 80, -4,4);
182  hHCAL_D4_METPhivsieta = ibooker.book2D("METTask_HCAL_D4_METPhivsieta","",83,-41,42, 80, -4,4);
183 
184  hHCAL_D1_MExvsieta = ibooker.book2D("METTask_HCAL_D1_MExvsieta","",83,-41,42, 10001,-500,501);
185  hHCAL_D2_MExvsieta = ibooker.book2D("METTask_HCAL_D2_MExvsieta","",83,-41,42, 10001,-500,501);
186  hHCAL_D3_MExvsieta = ibooker.book2D("METTask_HCAL_D3_MExvsieta","",83,-41,42, 10001,-500,501);
187  hHCAL_D4_MExvsieta = ibooker.book2D("METTask_HCAL_D4_MExvsieta","",83,-41,42, 10001,-500,501);
188 
189  hHCAL_D1_MEyvsieta = ibooker.book2D("METTask_HCAL_D1_MEyvsieta","",83,-41,42, 10001,-500,501);
190  hHCAL_D2_MEyvsieta = ibooker.book2D("METTask_HCAL_D2_MEyvsieta","",83,-41,42, 10001,-500,501);
191  hHCAL_D3_MEyvsieta = ibooker.book2D("METTask_HCAL_D3_MEyvsieta","",83,-41,42, 10001,-500,501);
192  hHCAL_D4_MEyvsieta = ibooker.book2D("METTask_HCAL_D4_MEyvsieta","",83,-41,42, 10001,-500,501);
193  }
194  else
195  {
196  hHCAL_D1_energyvsieta = ibooker.book2D("METTask_HCAL_D1_energyvsieta","",83,-41,42,1000,-10,1990);
197  hHCAL_D2_energyvsieta = ibooker.book2D("METTask_HCAL_D2_energyvsieta","",83,-41,42,1000,-10,1990);
198  hHCAL_D3_energyvsieta = ibooker.book2D("METTask_HCAL_D3_energyvsieta","",83,-41,42,1000,-10,1990);
199  hHCAL_D4_energyvsieta = ibooker.book2D("METTask_HCAL_D4_energyvsieta","",83,-41,42,1000,-10,1990);
200 
201  hHCAL_D1_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D1_Minenergyvsieta","",83,-41,42,1000,-10,1990);
202  hHCAL_D2_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D2_Minenergyvsieta","",83,-41,42,1000,-10,1990);
203  hHCAL_D3_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D3_Minenergyvsieta","",83,-41,42,1000,-10,1990);
204  hHCAL_D4_Minenergyvsieta = ibooker.book2D("METTask_HCAL_D4_Minenergyvsieta","",83,-41,42,1000,-10,1990);
205 
206  hHCAL_D1_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D1_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
207  hHCAL_D2_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D2_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
208  hHCAL_D3_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D3_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
209  hHCAL_D4_Maxenergyvsieta = ibooker.book2D("METTask_HCAL_D4_Maxenergyvsieta","",83,-41,42,1000,-10,1990);
210 
211  // Integrated over phi
212  hHCAL_D1_Occvsieta = ibooker.book2D("METTask_HCAL_D1_Occvsieta","",83,-41,42,73,0,73);
213  hHCAL_D2_Occvsieta = ibooker.book2D("METTask_HCAL_D2_Occvsieta","",83,-41,42,73,0,73);
214  hHCAL_D3_Occvsieta = ibooker.book2D("METTask_HCAL_D3_Occvsieta","",83,-41,42,73,0,73);
215  hHCAL_D4_Occvsieta = ibooker.book2D("METTask_HCAL_D4_Occvsieta","",83,-41,42,73,0,73);
216 
217  hHCAL_D1_SETvsieta = ibooker.book2D("METTask_HCAL_D1_SETvsieta","",83,-41,42,1000,0,2000);
218  hHCAL_D2_SETvsieta = ibooker.book2D("METTask_HCAL_D2_SETvsieta","",83,-41,42,1000,0,2000);
219  hHCAL_D3_SETvsieta = ibooker.book2D("METTask_HCAL_D3_SETvsieta","",83,-41,42,1000,0,2000);
220  hHCAL_D4_SETvsieta = ibooker.book2D("METTask_HCAL_D4_SETvsieta","",83,-41,42,1000,0,2000);
221 
222  hHCAL_D1_METvsieta = ibooker.book2D("METTask_HCAL_D1_METvsieta","",83,-41,42,1000,0,2000);
223  hHCAL_D2_METvsieta = ibooker.book2D("METTask_HCAL_D2_METvsieta","",83,-41,42,1000,0,2000);
224  hHCAL_D3_METvsieta = ibooker.book2D("METTask_HCAL_D3_METvsieta","",83,-41,42,1000,0,2000);
225  hHCAL_D4_METvsieta = ibooker.book2D("METTask_HCAL_D4_METvsieta","",83,-41,42,1000,0,2000);
226 
227  hHCAL_D1_METPhivsieta = ibooker.book2D("METTask_HCAL_D1_METPhivsieta","",83,-41,42, 80, -4,4);
228  hHCAL_D2_METPhivsieta = ibooker.book2D("METTask_HCAL_D2_METPhivsieta","",83,-41,42, 80, -4,4);
229  hHCAL_D3_METPhivsieta = ibooker.book2D("METTask_HCAL_D3_METPhivsieta","",83,-41,42, 80, -4,4);
230  hHCAL_D4_METPhivsieta = ibooker.book2D("METTask_HCAL_D4_METPhivsieta","",83,-41,42, 80, -4,4);
231 
232  hHCAL_D1_MExvsieta = ibooker.book2D("METTask_HCAL_D1_MExvsieta","",83,-41,42, 500,-500,500);
233  hHCAL_D2_MExvsieta = ibooker.book2D("METTask_HCAL_D2_MExvsieta","",83,-41,42, 500,-500,500);
234  hHCAL_D3_MExvsieta = ibooker.book2D("METTask_HCAL_D3_MExvsieta","",83,-41,42, 500,-500,500);
235  hHCAL_D4_MExvsieta = ibooker.book2D("METTask_HCAL_D4_MExvsieta","",83,-41,42, 500,-500,500);
236 
237  hHCAL_D1_MEyvsieta = ibooker.book2D("METTask_HCAL_D1_MEyvsieta","",83,-41,42, 500,-500,500);
238  hHCAL_D2_MEyvsieta = ibooker.book2D("METTask_HCAL_D2_MEyvsieta","",83,-41,42, 500,-500,500);
239  hHCAL_D3_MEyvsieta = ibooker.book2D("METTask_HCAL_D3_MEyvsieta","",83,-41,42, 500,-500,500);
240  hHCAL_D4_MEyvsieta = ibooker.book2D("METTask_HCAL_D4_MEyvsieta","",83,-41,42, 500,-500,500);
241 
242  }
243  // Inspect Setup for CaloTower Geometry
244  // FillGeometry(iSetup);
245 }
246 
247 void HCALRecHitAnalyzer::FillGeometry(const edm::EventSetup& iSetup)
248 {
249 
250  // ==========================================================
251  // Retrieve!
252  // ==========================================================
253 
254  const CaloSubdetectorGeometry* HBgeom;
255  const CaloSubdetectorGeometry* HEgeom;
256  const CaloSubdetectorGeometry* HOgeom;
257  const CaloSubdetectorGeometry* HFgeom;
258 
259  edm::ESHandle<CaloGeometry> pG;
260  iSetup.get<CaloGeometryRecord>().get(pG);
261 
262  if (!pG.isValid() )
263  {
264  edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Setup Handle, Aborting Task "
265  << "HCALRecHitAnalyzer::FillGeometry!\n";
266  return;
267  }
268 
269  const CaloGeometry cG = *pG;
270 
271  HBgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalBarrel);
272  HEgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalEndcap);
273  HOgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalOuter);
274  HFgeom = cG.getSubdetectorGeometry(DetId::Hcal,HcalForward);
275 
276 
277  // ==========================================================
278  // Fill Histograms!
279  // ==========================================================
280 
281  std::vector<DetId>::iterator i;
282 
283  int HBmin_ieta = 99, HBmax_ieta = -99;
284  int HBmin_iphi = 99, HBmax_iphi = -99;
285 
286  // Loop Over all Hcal Barrel DetId's
287  int nHBdetid = 0;
288  std::vector<DetId> HBids = HBgeom->getValidDetIds(DetId::Hcal,HcalBarrel);
289 
290  for (i = HBids.begin(); i != HBids.end(); i++) {
291 
292  nHBdetid++;
293 
294  const CaloCellGeometry* cell = HBgeom->getGeometry(*i);
295  HcalDetId HcalID(i->rawId());
296  //GlobalPoint p = cell->getPosition();
297 
298  int Calo_ieta = 42 + HcalID.ieta();
299  int Calo_iphi = HcalID.iphi();
300  double Calo_eta = cell->getPosition().eta();
301  double Calo_phi = cell->getPosition().phi();
302 
303  if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
304 
305  hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
306  hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
307 
308  }
309 
310  if (Calo_ieta > HBmax_ieta) HBmax_ieta = Calo_ieta;
311  if (Calo_ieta < HBmin_ieta) HBmin_ieta = Calo_ieta;
312  if (Calo_iphi > HBmax_iphi) HBmax_iphi = Calo_iphi;
313  if (Calo_iphi > HBmax_iphi) HBmin_iphi = Calo_iphi;
314 
315  }
316 
317  int HEmin_ieta = 99, HEmax_ieta = -99;
318  int HEmin_iphi = 99, HEmax_iphi = -99;
319 
320  // Loop Over all Hcal Endcap DetId's
321  int nHEdetid = 0;
322  std::vector<DetId> HEids = HEgeom->getValidDetIds(DetId::Hcal,HcalEndcap);
323 
324  for (i = HEids.begin(); i != HEids.end(); i++) {
325 
326  nHEdetid++;
327 
328  const CaloCellGeometry* cell = HEgeom->getGeometry(*i);
329  HcalDetId HcalID(i->rawId());
330  //GlobalPoint p = cell->getPosition();
331 
332  int Calo_ieta = 42 + HcalID.ieta();
333  int Calo_iphi = HcalID.iphi();
334  double Calo_eta = cell->getPosition().eta();
335  double Calo_phi = cell->getPosition().phi();
336 
337  // HCAL to HE eta, phi map comparison
338  if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
339  hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
340  hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
341  }
342 
343  if (Calo_ieta > HEmax_ieta) HEmax_ieta = Calo_ieta;
344  if (Calo_ieta < HEmin_ieta) HEmin_ieta = Calo_ieta;
345  if (Calo_iphi > HEmax_iphi) HEmax_iphi = Calo_iphi;
346  if (Calo_iphi > HEmax_iphi) HEmin_iphi = Calo_iphi;
347 
348  }
349 
350  int HFmin_ieta = 99, HFmax_ieta = -99;
351  int HFmin_iphi = 99, HFmax_iphi = -99;
352 
353  // Loop Over all Hcal Forward DetId's
354  int nHFdetid = 0;
355  std::vector<DetId> HFids = HFgeom->getValidDetIds(DetId::Hcal,HcalForward);
356 
357  for (i = HFids.begin(); i != HFids.end(); i++) {
358 
359  nHFdetid++;
360 
361  const CaloCellGeometry* cell = HFgeom->getGeometry(*i);
362  HcalDetId HcalID(i->rawId());
363  //GlobalPoint p = cell->getPosition();
364 
365  int Calo_ieta = 42 + HcalID.ieta();
366  int Calo_iphi = HcalID.iphi();
367  double Calo_eta = cell->getPosition().eta();
368  double Calo_phi = cell->getPosition().phi();
369 
370  // HCAL to HF eta, phi map comparison
371  if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
372  hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
373  hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
374  }
375 
376  if (Calo_ieta > HFmax_ieta) HFmax_ieta = Calo_ieta;
377  if (Calo_ieta < HFmin_ieta) HFmin_ieta = Calo_ieta;
378  if (Calo_iphi > HFmax_iphi) HFmax_iphi = Calo_iphi;
379  if (Calo_iphi > HFmax_iphi) HFmin_iphi = Calo_iphi;
380 
381  }
382 
383  int HOmin_ieta = 99, HOmax_ieta = -99;
384  int HOmin_iphi = 99, HOmax_iphi = -99;
385 
386  // Loop Over all Hcal Outer DetId's
387  int nHOdetid = 0;
388  std::vector<DetId> HOids = HOgeom->getValidDetIds(DetId::Hcal,HcalOuter);
389 
390  for (i = HOids.begin(); i != HOids.end(); i++) {
391 
392  nHOdetid++;
393 
394  const CaloCellGeometry* cell = HOgeom->getGeometry(*i);
395  HcalDetId HcalID(i->rawId());
396  //GlobalPoint p = cell->getPosition();
397 
398  int Calo_ieta = 42 + HcalID.ieta();
399  int Calo_iphi = HcalID.iphi();
400  double Calo_eta = cell->getPosition().eta();
401  double Calo_phi = cell->getPosition().phi();
402 
403  // HCAL to HO eta, phi map comparison
404  if (hHCAL_ieta_iphi_etaMap->getBinContent(Calo_ieta,Calo_iphi) == -999) {
405  hHCAL_ieta_iphi_etaMap->setBinContent(Calo_ieta,Calo_iphi,Calo_eta);
406  hHCAL_ieta_iphi_phiMap->setBinContent(Calo_ieta,Calo_iphi,Calo_phi*180.0/M_PI);
407  }
408 
409  if (Calo_ieta > HOmax_ieta) HOmax_ieta = Calo_ieta;
410  if (Calo_ieta < HOmin_ieta) HOmin_ieta = Calo_ieta;
411  if (Calo_iphi > HOmax_iphi) HOmax_iphi = Calo_iphi;
412  if (Calo_iphi > HOmax_iphi) HOmin_iphi = Calo_iphi;
413 
414  }
415 
416  // Set the Cell Size for each (ieta, iphi) Bin
417  double currentLowEdge_eta = 0; //double currentHighEdge_eta = 0;
418  for (int ieta = 1; ieta < 42 ; ieta++) {
419 
420  int ieta_ = 42 + ieta;
421  double eta = hHCAL_ieta_iphi_etaMap->getBinContent(ieta_,3);
422  double phi = hHCAL_ieta_iphi_phiMap->getBinContent(ieta_,3);
423  double deta = 2.0*(eta-currentLowEdge_eta);
424  deta = ((float)((int)(1.0E3*deta + 0.5)))/1.0E3;
425  double dphi = 2.0*phi;
426  if (ieta==40 || ieta==41) dphi = 20;
427  if (ieta<=39 && ieta>=21) dphi = 10;
428  if (ieta<=20) dphi = 5;
429  // BS: This is WRONG...need to correct overlap
430  if (ieta == 28) deta = 0.218;
431  if (ieta == 29) deta= 0.096;
432  currentLowEdge_eta += deta;
433 
434  // BS: This is WRONG...need to correct overlap
435  if (ieta == 29) currentLowEdge_eta = 2.964;
436 
437  hHCAL_ieta_detaMap->setBinContent(ieta_,deta); // positive rings
438  hHCAL_ieta_dphiMap->setBinContent(ieta_,dphi); // positive rings
439  hHCAL_ieta_detaMap->setBinContent(42-ieta,deta); // negative rings
440  hHCAL_ieta_dphiMap->setBinContent(42-ieta,dphi); // negative rings
441 
442  } // end loop over ieta
443 
444  edm::LogInfo("OutputInfo") << "HB ieta range: " << HBmin_ieta << " " << HBmax_ieta;
445  edm::LogInfo("OutputInfo") << "HB iphi range: " << HBmin_iphi << " " << HBmax_iphi;
446  edm::LogInfo("OutputInfo") << "HE ieta range: " << HEmin_ieta << " " << HEmax_ieta;
447  edm::LogInfo("OutputInfo") << "HE iphi range: " << HEmin_iphi << " " << HEmax_iphi;
448  edm::LogInfo("OutputInfo") << "HF ieta range: " << HFmin_ieta << " " << HFmax_ieta;
449  edm::LogInfo("OutputInfo") << "HF iphi range: " << HFmin_iphi << " " << HFmax_iphi;
450  edm::LogInfo("OutputInfo") << "HO ieta range: " << HOmin_ieta << " " << HOmax_ieta;
451  edm::LogInfo("OutputInfo") << "HO iphi range: " << HOmin_iphi << " " << HOmax_iphi;
452 
453 }
454 
455 void HCALRecHitAnalyzer::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
456 {
457  Nevents++;
458  hHCAL_Nevents->Fill(0);
459  // ==========================================================
460  // Retrieve!
461  // ==========================================================
462 
463  const HBHERecHitCollection *HBHERecHits;
464  const HORecHitCollection *HORecHits;
465  const HFRecHitCollection *HFRecHits;
466 
467  edm::Handle<HBHERecHitCollection> HBHERecHitsHandle;
468  iEvent.getByToken(hBHERecHitsLabel_,HBHERecHitsHandle);
469  if (!HBHERecHitsHandle.isValid()) {
470  edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
471  << "HCALRecHitAnalyzer::analyze!\n"; return;
472  } else {
473  HBHERecHits = HBHERecHitsHandle.product();
474  }
475  edm::Handle<HORecHitCollection> HORecHitsHandle;
476  iEvent.getByToken(hORecHitsLabel_,HORecHitsHandle);
477  if (!HORecHitsHandle.isValid()) {
478  edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
479  << "HCALRecHitAnalyzer::analyze!\n"; return;
480  } else {
481  HORecHits = HORecHitsHandle.product();
482  }
483  edm::Handle<HFRecHitCollection> HFRecHitsHandle;
484  iEvent.getByToken(hFRecHitsLabel_,HFRecHitsHandle);
485  if (!HFRecHitsHandle.isValid()) {
486  edm::LogInfo("OutputInfo") << "Failed to retrieve an Event Handle, Aborting Task "
487  << "HCALRecHitAnalyzer::analyze!\n"; return;
488  } else {
489  HFRecHits = HFRecHitsHandle.product();
490  }
491 
492  // ==========================================================
493  // Fill Histograms!
494  // ==========================================================
495 
496  TLorentzVector vHBHEMET_EtaRing[83][4];
497  int HBHEActiveRing[83][4];
498  int HBHENActiveCells[83][4];
499  double HBHESET_EtaRing[83][4];
500  double HBHEMinEnergy_EtaRing[83][4];
501  double HBHEMaxEnergy_EtaRing[83][4];
502 
503  for (int i = 0; i < 83; i++) {
504  for (int j = 0; j < 4; j++) {
505 
506  HBHEActiveRing[i][j] = 0;
507  HBHENActiveCells[i][j] = 0;
508  HBHESET_EtaRing[i][j] = 0;
509  HBHEMinEnergy_EtaRing[i][j] = 14E3;
510  HBHEMaxEnergy_EtaRing[i][j] = -999;
511  }
512  }
513 
514  // Loop over HBHERecHit's
515  HBHERecHitCollection::const_iterator hbherechit;
516  int nHBrechit = 0, nHErechit = 0;
517 
518  for (hbherechit = HBHERecHits->begin(); hbherechit != HBHERecHits->end(); hbherechit++) {
519 
520  HcalDetId det = hbherechit->id();
521  double Energy = hbherechit->energy();
522  Int_t depth = det.depth();
523  Int_t ieta = det.ieta();
524  Int_t iphi = det.iphi();
525  int EtaRing = 41 + ieta; // this counts from 0
526  double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
527  double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
528  double theta = 2*TMath::ATan(exp(-1*eta));
529  double ET = Energy*TMath::Sin(theta);
530  HcalSubdetector HcalNum = det.subdet();
531  TLorentzVector v_;
532 
533 
534 
535  if (Energy>0) // zero suppress
536  {
537  HBHEActiveRing[EtaRing][depth-1] = 1;
538  HBHENActiveCells[EtaRing][depth-1]++;
539  HBHESET_EtaRing[EtaRing][depth-1]+=ET;
540  v_.SetPtEtaPhiE(ET, 0, phi, ET);
541  vHBHEMET_EtaRing[EtaRing][depth-1]-=v_;
542 
543 
544  DEBUG( EtaRing << " " << Energy << " " << ET << " " << theta << " " << eta << " " << phi << " : " << vHBHEMET_EtaRing[EtaRing][depth-1].Phi() << " " << vHBHEMET_EtaRing[EtaRing][depth-1].Pt() );
545 
546  switch (depth) {
547  case 1:
548  hHCAL_D1_Occ_ieta_iphi->Fill(ieta,iphi);
549  break;
550  case 2:
551  hHCAL_D2_Occ_ieta_iphi->Fill(ieta,iphi);
552  break;
553  case 3:
554  hHCAL_D3_Occ_ieta_iphi->Fill(ieta,iphi);
555  break;
556  } // end switch
557  }
558 
559  if (Energy > HBHEMaxEnergy_EtaRing[EtaRing][depth-1])
560  HBHEMaxEnergy_EtaRing[EtaRing][depth-1] = Energy;
561  if (Energy < HBHEMinEnergy_EtaRing[EtaRing][depth-1])
562  HBHEMinEnergy_EtaRing[EtaRing][depth-1] = Energy;
563 
564  switch (depth) {
565  case 1:
566  hHCAL_D1_energy_ieta_iphi->Fill(ieta,iphi,Energy);
567  hHCAL_D1_energyvsieta->Fill(ieta,Energy);
568  if (Energy>hHCAL_D1_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
569  hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
570  if (Energy<hHCAL_D1_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
571  hHCAL_D1_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
572 
573  break;
574  case 2:
575  hHCAL_D2_energy_ieta_iphi->Fill(ieta,iphi,Energy);
576 
577  hHCAL_D2_energyvsieta->Fill(ieta,Energy);
578  if (Energy>hHCAL_D2_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
579  hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
580  if (Energy<hHCAL_D2_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
581  hHCAL_D2_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
582  break;
583  case 3:
584  hHCAL_D3_energy_ieta_iphi->Fill(ieta,iphi,Energy);
585 
586  hHCAL_D3_energyvsieta->Fill(ieta,Energy);
587  if (Energy>hHCAL_D3_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
588  hHCAL_D3_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
589  if (Energy<hHCAL_D3_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
590  hHCAL_D3_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
591  break;
592  } // end switch
593 
594 
595  if (HcalNum == HcalBarrel) {
596  nHBrechit++;
597  } else { // HcalEndcap
598  nHErechit++;
599  }
600  } // end loop over HBHERecHit's
601 
602  // Fill eta-ring MET quantities
603  for (int iEtaRing=0; iEtaRing < 83; iEtaRing++) {
604  for (int jDepth=0; jDepth < 3; jDepth++) {
605 
606  switch (jDepth+1) {
607  case 1:
608  hHCAL_D1_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
609  hHCAL_D1_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
610  break;
611  case 2:
612  hHCAL_D2_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
613  hHCAL_D2_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
614  break;
615  case 3:
616  hHCAL_D3_Maxenergyvsieta->Fill(iEtaRing-41, HBHEMaxEnergy_EtaRing[iEtaRing][jDepth]);
617  hHCAL_D3_Minenergyvsieta->Fill(iEtaRing-41, HBHEMinEnergy_EtaRing[iEtaRing][jDepth]);
618  break;
619  }
620 
621  if (HBHEActiveRing[iEtaRing][jDepth]) {
622 
623  switch (jDepth+1) {
624  case 1:
625  hHCAL_D1_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
626  hHCAL_D1_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
627  hHCAL_D1_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
628  hHCAL_D1_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
629  hHCAL_D1_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
630  hHCAL_D1_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
631  break;
632  case 2:
633  hHCAL_D2_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
634  hHCAL_D2_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
635  hHCAL_D2_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
636  hHCAL_D2_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
637  hHCAL_D2_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
638  hHCAL_D2_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
639  break;
640  case 3:
641  hHCAL_D3_METPhivsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Phi());
642  hHCAL_D3_MExvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Px());
643  hHCAL_D3_MEyvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Py());
644  hHCAL_D3_METvsieta->Fill(iEtaRing-41, vHBHEMET_EtaRing[iEtaRing][jDepth].Pt());
645  hHCAL_D3_SETvsieta->Fill(iEtaRing-41, HBHESET_EtaRing[iEtaRing][jDepth]);
646  hHCAL_D3_Occvsieta->Fill(iEtaRing-41, HBHENActiveCells[iEtaRing][jDepth]);
647  break;
648  } //switch
649  } //activering
650  }//depth
651  } //etaring
652 
653 
654  //-------------------------------------------------HO
655  TLorentzVector vHOMET_EtaRing[83];
656  int HOActiveRing[83];
657  int HONActiveCells[83];
658  double HOSET_EtaRing[83];
659  double HOMinEnergy_EtaRing[83];
660  double HOMaxEnergy_EtaRing[83];
661 
662  for (int i=0;i<83; i++) {
663 
664  HOActiveRing[i] = 0;
665  HONActiveCells[i] = 0;
666  HOSET_EtaRing[i] = 0;
667  HOMinEnergy_EtaRing[i] = 14E3;
668  HOMaxEnergy_EtaRing[i] = -999;
669  }
670 
671  // Loop over HORecHit's
672  HORecHitCollection::const_iterator horechit;
673  int nHOrechit = 0;
674 
675  for (horechit = HORecHits->begin(); horechit != HORecHits->end(); horechit++) {
676 
677  nHOrechit++;
678 
679  HcalDetId det = horechit->id();
680  double Energy = horechit->energy();
682  Int_t ieta = det.ieta();
683  Int_t iphi = det.iphi();
684  int EtaRing = 41+ieta; // this counts from 0
685  double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
686  double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
687  double theta = 2*TMath::ATan(exp(-1*eta));
688  double ET = Energy*TMath::Sin(theta);
689  TLorentzVector v_;
690 
691 
692  if (Energy>0) // zero suppress
693  {
694  HOActiveRing[EtaRing] = 1;
695  HONActiveCells[EtaRing]++;
696  HOSET_EtaRing[EtaRing]+=ET;
697  v_.SetPtEtaPhiE(ET, 0, phi, ET);
698  vHOMET_EtaRing[EtaRing]-=v_;
699 
700  hHCAL_D4_Occ_ieta_iphi->Fill(ieta,iphi);
701  }
702 
703  if (Energy > HOMaxEnergy_EtaRing[EtaRing])
704  HOMaxEnergy_EtaRing[EtaRing] = Energy;
705  if (Energy < HOMinEnergy_EtaRing[EtaRing])
706  HOMinEnergy_EtaRing[EtaRing] = Energy;
707 
708  hHCAL_D4_energy_ieta_iphi->Fill(ieta,iphi,Energy);
709 
710  hHCAL_D4_energyvsieta->Fill(ieta,Energy);
711  if (Energy>hHCAL_D4_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
712  hHCAL_D4_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
713  if (Energy<hHCAL_D4_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
714  hHCAL_D4_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
715 
716  } // end loop over HORecHit's
717 
718  // Fill eta-ring MET quantities
719  for (int iEtaRing=0; iEtaRing<83; iEtaRing++) {
720  hHCAL_D4_Maxenergyvsieta->Fill(iEtaRing-41, HOMaxEnergy_EtaRing[iEtaRing]);
721  hHCAL_D4_Minenergyvsieta->Fill(iEtaRing-41, HOMinEnergy_EtaRing[iEtaRing]);
722 
723  if (HOActiveRing[iEtaRing]) {
724 
725  hHCAL_D4_METPhivsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Phi());
726  hHCAL_D4_MExvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Px());
727  hHCAL_D4_MEyvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Py());
728  hHCAL_D4_METvsieta->Fill(iEtaRing-41, vHOMET_EtaRing[iEtaRing].Pt());
729  hHCAL_D4_SETvsieta->Fill(iEtaRing-41, HOSET_EtaRing[iEtaRing]);
730  hHCAL_D4_Occvsieta->Fill(iEtaRing-41, HONActiveCells[iEtaRing]);
731 
732  }
733 
734  }
735 
736  //------------------------------------------------HF
737 
738  TLorentzVector vHFMET_EtaRing[83][2];
739  int HFActiveRing[83][2];
740  int HFNActiveCells[83][2];
741  double HFSET_EtaRing[83][2];
742  double HFMinEnergy_EtaRing[83][2];
743  double HFMaxEnergy_EtaRing[83][2];
744 
745  for (int i=0;i<83; i++) {
746  for (int j=0;j<2; j++) {
747 
748  HFActiveRing[i][j] = 0;
749  HFNActiveCells[i][j] = 0;
750  HFSET_EtaRing[i][j] = 0;
751  HFMinEnergy_EtaRing[i][j] = 14E3;
752  HFMaxEnergy_EtaRing[i][j] = -999;
753 
754  }
755  }
756 
757  // Loop over HFRecHit's
758  HFRecHitCollection::const_iterator hfrechit;
759  int nHFrechit = 0;
760 
761  for (hfrechit = HFRecHits->begin(); hfrechit != HFRecHits->end(); hfrechit++) {
762 
763  nHFrechit++;
764 
765  HcalDetId det = hfrechit->id();
766  double Energy = hfrechit->energy();
767  Int_t depth = det.depth();
768  Int_t ieta = det.ieta();
769  Int_t iphi = det.iphi();
770  int EtaRing = 41+ieta; // this counts from 0
771  double eta = hHCAL_ieta_iphi_etaMap->getBinContent(EtaRing+1,iphi);
772  double phi = hHCAL_ieta_iphi_phiMap->getBinContent(EtaRing+1,iphi);
773  double theta = 2*TMath::ATan(exp(-1*eta));
774  double ET = Energy*TMath::Sin(theta);
775 
776  TLorentzVector v_;
777  if (Energy>0) // zero suppress
778  {
779  HFActiveRing[EtaRing][depth-1] = 1;
780  HFNActiveCells[EtaRing][depth-1]++;
781  HFSET_EtaRing[EtaRing][depth-1]+=ET;
782  v_.SetPtEtaPhiE(ET, 0, phi, ET);
783  vHFMET_EtaRing[EtaRing][depth-1]-=v_;
784 
785  switch (depth) {
786  case 1:
787  hHCAL_D1_Occ_ieta_iphi->Fill(ieta,iphi);
788  break;
789  case 2:
790  hHCAL_D2_Occ_ieta_iphi->Fill(ieta,iphi);
791  break;
792 
793  }
794  }
795 
796  if (Energy > HFMaxEnergy_EtaRing[EtaRing][depth-1])
797  HFMaxEnergy_EtaRing[EtaRing][depth-1] = Energy;
798  if (Energy < HFMinEnergy_EtaRing[EtaRing][depth-1])
799  HFMinEnergy_EtaRing[EtaRing][depth-1] = Energy;
800 
801 
802  switch (depth) {
803  case 1:
804  hHCAL_D1_energy_ieta_iphi->Fill(ieta,iphi,Energy);
805 
806  hHCAL_D1_energyvsieta->Fill(ieta,Energy);
807  if (Energy>hHCAL_D1_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
808  hHCAL_D1_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
809  if (Energy<hHCAL_D1_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
810  hHCAL_D1_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
811  break;
812  case 2:
813  hHCAL_D2_energy_ieta_iphi->Fill(ieta,iphi,Energy);
814 
815  hHCAL_D2_energyvsieta->Fill(ieta,Energy);
816  if (Energy>hHCAL_D2_Maxenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
817  hHCAL_D2_Maxenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
818  if (Energy<hHCAL_D2_Minenergy_ieta_iphi->getBinContent(EtaRing+1, iphi))
819  hHCAL_D2_Minenergy_ieta_iphi->setBinContent(EtaRing+1, iphi, Energy);
820  break;
821  }
822 
823  } // end loop over HFRecHit's
824 
825  // Fill eta-ring MET quantities
826  for (int iEtaRing=0; iEtaRing<83; iEtaRing++) {
827  for (int jDepth=0; jDepth<2; jDepth++) {
828  switch (jDepth+1) {
829  case 1:
830  hHCAL_D1_Maxenergyvsieta->Fill(iEtaRing-41, HFMaxEnergy_EtaRing[iEtaRing][jDepth]);
831  hHCAL_D1_Minenergyvsieta->Fill(iEtaRing-41, HFMinEnergy_EtaRing[iEtaRing][jDepth]);
832  break;
833  case 2:
834  hHCAL_D2_Maxenergyvsieta->Fill(iEtaRing-41, HFMaxEnergy_EtaRing[iEtaRing][jDepth]);
835  hHCAL_D2_Minenergyvsieta->Fill(iEtaRing-41, HFMinEnergy_EtaRing[iEtaRing][jDepth]);
836  break;
837  }
838 
839 
840  if (HFActiveRing[iEtaRing][jDepth]) {
841 
842  switch (jDepth+1) {
843  case 1:
844 
845  hHCAL_D1_METPhivsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Phi());
846  hHCAL_D1_MExvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Px());
847  hHCAL_D1_MEyvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Py());
848  hHCAL_D1_METvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Pt());
849  hHCAL_D1_SETvsieta->Fill(iEtaRing-41, HFSET_EtaRing[iEtaRing][jDepth]);
850  hHCAL_D1_Occvsieta->Fill(iEtaRing-41, HFNActiveCells[iEtaRing][jDepth]);
851  break;
852 
853  case 2:
854 
855  hHCAL_D2_METPhivsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Phi());
856  hHCAL_D2_MExvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Px());
857  hHCAL_D2_MEyvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Py());
858  hHCAL_D2_METvsieta->Fill(iEtaRing-41, vHFMET_EtaRing[iEtaRing][jDepth].Pt());
859  hHCAL_D2_SETvsieta->Fill(iEtaRing-41, HFSET_EtaRing[iEtaRing][jDepth]);
860  hHCAL_D2_Occvsieta->Fill(iEtaRing-41, HFNActiveCells[iEtaRing][jDepth]);
861  break;
862 
863  }
864  }
865  }
866  }
867 
868 
869 
870 }
871 
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
CaloGeometry::getSubdetectorGeometry
const CaloSubdetectorGeometry * getSubdetectorGeometry(const DetId &id) const
access the subdetector geometry for the given subdetector directly
Definition: CaloGeometry.cc:43
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
HCALRecHitAnalyzer::hHCAL_D4_Maxenergy_ieta_iphi
MonitorElement * hHCAL_D4_Maxenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:76
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:844
HCALRecHitAnalyzer::hHCAL_ieta_iphi_HFMap
MonitorElement * hHCAL_ieta_iphi_HFMap
Definition: HCALRecHitAnalyzer.h:54
HCALRecHitAnalyzer::hHCAL_D3_Maxenergy_ieta_iphi
MonitorElement * hHCAL_D3_Maxenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:75
create_public_lumi_plots.exp
tuple exp
Definition: create_public_lumi_plots.py:1088
HCALRecHitAnalyzer::hHCAL_D3_Maxenergyvsieta
MonitorElement * hHCAL_D3_Maxenergyvsieta
Definition: HCALRecHitAnalyzer.h:95
HcalDetId::subdet
HcalSubdetector subdet() const
get the subdetector
Definition: HcalDetId.h:49
HCALRecHitAnalyzer::hHCAL_ieta_iphi_phiMap
MonitorElement * hHCAL_ieta_iphi_phiMap
Definition: HCALRecHitAnalyzer.h:57
HCALRecHitAnalyzer::hHCAL_D1_SETvsieta
MonitorElement * hHCAL_D1_SETvsieta
Definition: HCALRecHitAnalyzer.h:103
HCALRecHitAnalyzer::hHCAL_D3_Occvsieta
MonitorElement * hHCAL_D3_Occvsieta
Definition: HCALRecHitAnalyzer.h:100
HCALRecHitAnalyzer::HCALRecHitAnalyzer
HCALRecHitAnalyzer(const edm::ParameterSet &)
Definition: HCALRecHitAnalyzer.cc:46
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
HCALRecHitAnalyzer::hHCAL_D2_Minenergy_ieta_iphi
MonitorElement * hHCAL_D2_Minenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:69
HCALRecHitAnalyzer::hHCAL_D4_METvsieta
MonitorElement * hHCAL_D4_METvsieta
Definition: HCALRecHitAnalyzer.h:111
HCALRecHitAnalyzer::hHCAL_D4_energyvsieta
MonitorElement * hHCAL_D4_energyvsieta
Definition: HCALRecHitAnalyzer.h:86
HCALRecHitAnalyzer::dqmBeginRun
void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
Definition: HCALRecHitAnalyzer.cc:60
HCALRecHitAnalyzer::hHCAL_D4_Maxenergyvsieta
MonitorElement * hHCAL_D4_Maxenergyvsieta
Definition: HCALRecHitAnalyzer.h:96
PV3DBase::phi
Geom::Phi< T > phi() const
Definition: PV3DBase.h:69
HCALRecHitAnalyzer::hHCAL_D1_Minenergy_ieta_iphi
MonitorElement * hHCAL_D1_Minenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:68
edm::SortedCollection< HBHERecHit >::const_iterator
std::vector< HBHERecHit >::const_iterator const_iterator
Definition: SortedCollection.h:82
theta
Geom::Theta< T > theta() const
Definition: Basic3DVectorLD.h:170
HCALRecHitAnalyzer::hHCAL_D4_Occvsieta
MonitorElement * hHCAL_D4_Occvsieta
Definition: HCALRecHitAnalyzer.h:101
HCALRecHitAnalyzer::FillGeometry
void FillGeometry(const edm::EventSetup &)
Definition: HCALRecHitAnalyzer.cc:247
HCALRecHitAnalyzer::hHCAL_D3_Minenergyvsieta
MonitorElement * hHCAL_D3_Minenergyvsieta
Definition: HCALRecHitAnalyzer.h:90
HCALRecHitAnalyzer::hHCAL_D3_Occ_ieta_iphi
MonitorElement * hHCAL_D3_Occ_ieta_iphi
Definition: HCALRecHitAnalyzer.h:80
HCALRecHitAnalyzer::hHCAL_D4_Occ_ieta_iphi
MonitorElement * hHCAL_D4_Occ_ieta_iphi
Definition: HCALRecHitAnalyzer.h:81
HLT_25ns10e33_v2_cff.depth
tuple depth
Definition: HLT_25ns10e33_v2_cff.py:10449
HCALRecHitAnalyzer::hHCAL_D2_Minenergyvsieta
MonitorElement * hHCAL_D2_Minenergyvsieta
Definition: HCALRecHitAnalyzer.h:89
HCALRecHitAnalyzer::hHCAL_ieta_iphi_etaMap
MonitorElement * hHCAL_ieta_iphi_etaMap
Definition: HCALRecHitAnalyzer.h:56
HCALRecHitAnalyzer::hHCAL_D4_MEyvsieta
MonitorElement * hHCAL_D4_MEyvsieta
Definition: HCALRecHitAnalyzer.h:126
HCALRecHitAnalyzer::hHCAL_D4_SETvsieta
MonitorElement * hHCAL_D4_SETvsieta
Definition: HCALRecHitAnalyzer.h:106
HCALRecHitAnalyzer::hHCAL_ieta_iphi_HBMap
MonitorElement * hHCAL_ieta_iphi_HBMap
Definition: HCALRecHitAnalyzer.h:52
HCALRecHitAnalyzer::hHCAL_D3_METvsieta
MonitorElement * hHCAL_D3_METvsieta
Definition: HCALRecHitAnalyzer.h:110
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
HCALRecHitAnalyzer::hHCAL_ieta_iphi_HOMap
MonitorElement * hHCAL_ieta_iphi_HOMap
Definition: HCALRecHitAnalyzer.h:55
CaloSubdetectorGeometry::getValidDetIds
virtual const std::vector< DetId > & getValidDetIds(DetId::Detector det=DetId::Detector(0), int subdet=0) const
Get a list of valid detector ids (for the given subdetector)
Definition: CaloSubdetectorGeometry.cc:37
HCALRecHitAnalyzer::hHCAL_D3_MExvsieta
MonitorElement * hHCAL_D3_MExvsieta
Definition: HCALRecHitAnalyzer.h:120
CaloSubdetectorGeometry::getGeometry
virtual const CaloCellGeometry * getGeometry(const DetId &id) const
Get the cell geometry of a given detector id. Should return false if not found.
Definition: CaloSubdetectorGeometry.cc:44
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
HCALRecHitAnalyzer::hHCAL_D1_Maxenergy_ieta_iphi
MonitorElement * hHCAL_D1_Maxenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:73
HCALRecHitAnalyzer::hHCAL_D1_Occ_ieta_iphi
MonitorElement * hHCAL_D1_Occ_ieta_iphi
Definition: HCALRecHitAnalyzer.h:78
HCALRecHitAnalyzer::hHCAL_D4_METPhivsieta
MonitorElement * hHCAL_D4_METPhivsieta
Definition: HCALRecHitAnalyzer.h:116
HcalDetId::depth
int depth() const
get the tower depth
Definition: HcalDetId.cc:106
HCALRecHitAnalyzer::hHCAL_D4_Minenergyvsieta
MonitorElement * hHCAL_D4_Minenergyvsieta
Definition: HCALRecHitAnalyzer.h:91
iEvent
int iEvent
Definition: GenABIO.cc:230
HCALRecHitAnalyzer::hHCAL_D2_METvsieta
MonitorElement * hHCAL_D2_METvsieta
Definition: HCALRecHitAnalyzer.h:109
HCALRecHitAnalyzer::hHCAL_D1_Occvsieta
MonitorElement * hHCAL_D1_Occvsieta
Definition: HCALRecHitAnalyzer.h:98
HCALRecHitAnalyzer::hHCAL_D1_Minenergyvsieta
MonitorElement * hHCAL_D1_Minenergyvsieta
Definition: HCALRecHitAnalyzer.h:88
HCALRecHitAnalyzer::hHCAL_D1_METPhivsieta
MonitorElement * hHCAL_D1_METPhivsieta
Definition: HCALRecHitAnalyzer.h:113
HCALRecHitAnalyzer::hHCAL_D1_METvsieta
MonitorElement * hHCAL_D1_METvsieta
Definition: HCALRecHitAnalyzer.h:108
HCALRecHitAnalyzer::hHCAL_D1_Maxenergyvsieta
MonitorElement * hHCAL_D1_Maxenergyvsieta
Definition: HCALRecHitAnalyzer.h:93
HCALRecHitAnalyzer::hHCAL_D3_energyvsieta
MonitorElement * hHCAL_D3_energyvsieta
Definition: HCALRecHitAnalyzer.h:85
HCALRecHitAnalyzer::hHCAL_D2_Maxenergyvsieta
MonitorElement * hHCAL_D2_Maxenergyvsieta
Definition: HCALRecHitAnalyzer.h:94
HCALRecHitAnalyzer::hHCAL_D2_energyvsieta
MonitorElement * hHCAL_D2_energyvsieta
Definition: HCALRecHitAnalyzer.h:84
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:56
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
HcalSubdetector
HcalSubdetector
Definition: HcalAssistant.h:31
HCALRecHitAnalyzer::hFRecHitsLabel_
edm::EDGetTokenT< HFRecHitCollection > hFRecHitsLabel_
Definition: HCALRecHitAnalyzer.h:41
j
int j
Definition: DBlmapReader.cc:9
HCALRecHitAnalyzer::hHCAL_D2_METPhivsieta
MonitorElement * hHCAL_D2_METPhivsieta
Definition: HCALRecHitAnalyzer.h:114
HCALRecHitAnalyzer::hHCAL_D2_Occvsieta
MonitorElement * hHCAL_D2_Occvsieta
Definition: HCALRecHitAnalyzer.h:99
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
HCALRecHitAnalyzer::hHCAL_D2_Maxenergy_ieta_iphi
MonitorElement * hHCAL_D2_Maxenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:74
HCALRecHitAnalyzer::hHCAL_D3_METPhivsieta
MonitorElement * hHCAL_D3_METPhivsieta
Definition: HCALRecHitAnalyzer.h:115
HCALRecHitAnalyzer::hHCAL_D1_energyvsieta
MonitorElement * hHCAL_D1_energyvsieta
Definition: HCALRecHitAnalyzer.h:83
HCALRecHitAnalyzer::hHCAL_ieta_detaMap
MonitorElement * hHCAL_ieta_detaMap
Definition: HCALRecHitAnalyzer.h:58
HCALRecHitAnalyzer::hHCAL_D2_energy_ieta_iphi
MonitorElement * hHCAL_D2_energy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:64
M_PI
#define M_PI
Definition: BXVectorInputProducer.cc:51
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:310
HCALRecHitAnalyzer::hHCAL_D1_energy_ieta_iphi
MonitorElement * hHCAL_D1_energy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:63
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.cc:101
HCALRecHitAnalyzer::hHCAL_D4_Minenergy_ieta_iphi
MonitorElement * hHCAL_D4_Minenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:71
HCALRecHitAnalyzer::hHCAL_D3_SETvsieta
MonitorElement * hHCAL_D3_SETvsieta
Definition: HCALRecHitAnalyzer.h:105
HCALRecHitAnalyzer::hORecHitsLabel_
edm::EDGetTokenT< HORecHitCollection > hORecHitsLabel_
Definition: HCALRecHitAnalyzer.h:42
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
HCALRecHitAnalyzer::hHCAL_D1_MExvsieta
MonitorElement * hHCAL_D1_MExvsieta
Definition: HCALRecHitAnalyzer.h:118
HCALRecHitAnalyzer::hHCAL_D4_MExvsieta
MonitorElement * hHCAL_D4_MExvsieta
Definition: HCALRecHitAnalyzer.h:121
HCALRecHitAnalyzer::hHCAL_D2_MEyvsieta
MonitorElement * hHCAL_D2_MEyvsieta
Definition: HCALRecHitAnalyzer.h:124
edm::Handle::product
T const * product() const
Definition: Handle.h:81
DEBUG
#define DEBUG(X)
Definition: HCALRecHitAnalyzer.cc:44
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
HCALRecHitAnalyzer::hHCAL_D1_MEyvsieta
MonitorElement * hHCAL_D1_MEyvsieta
Definition: HCALRecHitAnalyzer.h:123
HCALRecHitAnalyzer::hBHERecHitsLabel_
edm::EDGetTokenT< HBHERecHitCollection > hBHERecHitsLabel_
Definition: HCALRecHitAnalyzer.h:40
PV3DBase::eta
T eta() const
Definition: PV3DBase.h:76
HCALRecHitAnalyzer::hHCAL_D4_energy_ieta_iphi
MonitorElement * hHCAL_D4_energy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:66
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:743
HCALRecHitAnalyzer::hHCAL_D3_energy_ieta_iphi
MonitorElement * hHCAL_D3_energy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:65
HCALRecHitAnalyzer::hHCAL_D2_Occ_ieta_iphi
MonitorElement * hHCAL_D2_Occ_ieta_iphi
Definition: HCALRecHitAnalyzer.h:79
HCALRecHitAnalyzer::hHCAL_Nevents
MonitorElement * hHCAL_Nevents
Definition: HCALRecHitAnalyzer.h:61
HCALRecHitAnalyzer::hHCAL_D3_Minenergy_ieta_iphi
MonitorElement * hHCAL_D3_Minenergy_ieta_iphi
Definition: HCALRecHitAnalyzer.h:70
HCALRecHitAnalyzer::hHCAL_D2_SETvsieta
MonitorElement * hHCAL_D2_SETvsieta
Definition: HCALRecHitAnalyzer.h:104
ET
#define ET
Definition: GenericBenchmark.cc:31
HCALRecHitAnalyzer::hHCAL_D2_MExvsieta
MonitorElement * hHCAL_D2_MExvsieta
Definition: HCALRecHitAnalyzer.h:119
HCALRecHitAnalyzer::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: HCALRecHitAnalyzer.cc:455
edm::ESHandleBase::isValid
bool isValid() const
Definition: ESHandle.h:47
CaloCellGeometry::getPosition
const GlobalPoint & getPosition() const
Returns the position of reference for this cell.
Definition: CaloCellGeometry.h:84
HCALRecHitAnalyzer::hHCAL_D3_MEyvsieta
MonitorElement * hHCAL_D3_MEyvsieta
Definition: HCALRecHitAnalyzer.h:125
HCALRecHitAnalyzer::hHCAL_ieta_iphi_HEMap
MonitorElement * hHCAL_ieta_iphi_HEMap
Definition: HCALRecHitAnalyzer.h:53
HCALRecHitAnalyzer::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: HCALRecHitAnalyzer.cc:66
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:303
edm::Run
Definition: Run.h:43
HCALRecHitAnalyzer::hHCAL_ieta_dphiMap
MonitorElement * hHCAL_ieta_dphiMap
Definition: HCALRecHitAnalyzer.h:59
Generated for CMSSW Reference Manual by doxygen 1.8.5