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