CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Groups Pages
L1TCaloParamsViewer.cc
Go to the documentation of this file.
5 
8 
13 #include <iomanip>
14 
16 private:
28  bool printHfSF;
39 
41 
42  std::string hash(void* buf, size_t len) const;
43 
44 public:
45  void analyze(const edm::Event&, const edm::EventSetup&) override;
46 
48  printPUSParams = pset.getUntrackedParameter<bool>("printPUSParams", false);
49  printTauCalibLUT = pset.getUntrackedParameter<bool>("printTauCalibLUT", false);
50  printTauCompressLUT = pset.getUntrackedParameter<bool>("printTauCompressLUT", false);
51  printJetCalibLUT = pset.getUntrackedParameter<bool>("printJetCalibLUT", false);
52  printJetCalibPar = pset.getUntrackedParameter<bool>("printJetCalibParams", false);
53  printJetPUSPar = pset.getUntrackedParameter<bool>("printJetPUSPar", false);
54  printJetCompressPtLUT = pset.getUntrackedParameter<bool>("printJetCompressPtLUT", false);
55  printJetCompressEtaLUT = pset.getUntrackedParameter<bool>("printJetCompressEtaLUT", false);
56  printEgCalibLUT = pset.getUntrackedParameter<bool>("printEgCalibLUT", false);
57  printEgIsoLUT = pset.getUntrackedParameter<bool>("printEgIsoLUT", false);
58  printEtSumMetPUSLUT = pset.getUntrackedParameter<bool>("printEtSumMetPUSLUT", false);
59  printHfSF = pset.getUntrackedParameter<bool>("printHfSF", false);
60  printHcalSF = pset.getUntrackedParameter<bool>("printHcalSF", false);
61  printEcalSF = pset.getUntrackedParameter<bool>("printEcalSF", false);
62  printEtSumEttPUSLUT = pset.getUntrackedParameter<bool>("printEtSumEttPUSLUT", false);
63  printEtSumEcalSumPUSLUT = pset.getUntrackedParameter<bool>("printEtSumEcalSumPUSLUT", false);
64  printMetCalibrationLUT = pset.getUntrackedParameter<bool>("printMetCalibrationLUT", false);
65  printMetHFCalibrationLUT = pset.getUntrackedParameter<bool>("printMetHFCalibrationLUT", false);
66  printEtSumEttCalibrationLUT = pset.getUntrackedParameter<bool>("printEtSumEttCalibrationLUT", false);
67  printEtSumEcalSumCalibrationLUT = pset.getUntrackedParameter<bool>("printEtSumEcalSumCalibrationLUT", false);
68 
69  useStage2Rcd = pset.getUntrackedParameter<bool>("useStage2Rcd", false);
70  }
71 
72  ~L1TCaloParamsViewer(void) override {}
73 };
74 
75 #include <openssl/sha.h>
76 #include <cmath>
77 #include <iostream>
78 using namespace std;
79 
80 std::string L1TCaloParamsViewer::hash(void* buf, size_t len) const {
81  char tmp[SHA_DIGEST_LENGTH * 2 + 1];
82  bzero(tmp, sizeof(tmp));
83  SHA_CTX ctx;
84  if (!SHA1_Init(&ctx))
85  throw cms::Exception("L1TCaloParamsViewer::hash") << "SHA1 initialization error";
86 
87  if (!SHA1_Update(&ctx, buf, len))
88  throw cms::Exception("L1TCaloParamsViewer::hash") << "SHA1 processing error";
89 
90  unsigned char hash[SHA_DIGEST_LENGTH];
91  if (!SHA1_Final(hash, &ctx))
92  throw cms::Exception("L1TCaloParamsViewer::hash") << "SHA1 finalization error";
93 
94  // re-write bytes in hex
95  for (unsigned int i = 0; i < 20; i++)
96  ::sprintf(&tmp[i * 2], "%02x", hash[i]);
97 
98  tmp[20 * 2] = 0;
99  return std::string(tmp);
100 }
101 
104  if (useStage2Rcd)
105  evSetup.get<L1TCaloStage2ParamsRcd>().get(handle1);
106  else
107  evSetup.get<L1TCaloParamsRcd>().get(handle1);
108 
109  std::shared_ptr<l1t::CaloParams> ptr(new l1t::CaloParams(*(handle1.product())));
110 
111  l1t::CaloParamsHelper* ptr1 = nullptr;
112  ptr1 = (l1t::CaloParamsHelper*)(&(*ptr));
113 
114  edm::LogInfo("") << "L1TCaloParamsViewer:";
115 
116  cout << endl << " Towers: " << endl;
117  cout << " towerLsbH= " << ptr1->towerLsbH() << endl;
118  cout << " towerLsbE= " << ptr1->towerLsbE() << endl;
119  cout << " towerLsbSum= " << ptr1->towerLsbSum() << endl;
120  cout << " towerNBitsH= " << ptr1->towerNBitsH() << endl;
121  cout << " towerNBitsE= " << ptr1->towerNBitsE() << endl;
122  cout << " towerNBitsSum= " << ptr1->towerNBitsSum() << endl;
123  cout << " towerNBitsRatio= " << ptr1->towerNBitsRatio() << endl;
124  cout << " towerMaskE= " << ptr1->towerMaskE() << endl;
125  cout << " towerMaskH= " << ptr1->towerMaskH() << endl;
126  cout << " towerMaskSum= " << ptr1->towerMaskSum() << endl;
127  cout << " towerEncoding= " << ptr1->doTowerEncoding() << endl;
128 
129  cout << endl << " Regions: " << endl;
130  cout << " regionLsb= " << ptr1->regionLsb() << endl;
131  cout << " regionPUSType= " << ptr1->regionPUSType() << endl;
132  cout << " regionPUSParams= [" << ptr1->regionPUSParams().size() << "] ";
133  float pusParams[ptr1->regionPUSParams().size()];
134  for (unsigned int i = 0; i < ptr1->regionPUSParams().size(); i++) {
135  pusParams[i] = ceil(2 * ptr1->regionPUSParams()[i]);
136  if (printPUSParams)
137  cout << " " << ceil(2 * pusParams[i]) << endl;
138  }
139 
140  if (!ptr1->regionPUSParams().empty())
141  cout << hash(pusParams, sizeof(float) * ptr1->regionPUSParams().size()) << endl;
142  else
143  cout << endl;
144 
145  if (!ptr1->regionPUSLUT()->empty()) {
146  cout << " regionPUSLUT= [" << ptr1->regionPUSLUT()->maxSize() << "] ";
147  int regionPUSLUT[ptr1->regionPUSLUT()->maxSize()];
148  for (unsigned int i = 0; i < ptr1->regionPUSLUT()->maxSize(); i++)
149  regionPUSLUT[i] = ptr1->regionPUSLUT()->data(i);
150  cout << hash(regionPUSLUT, sizeof(int) * ptr1->regionPUSLUT()->maxSize()) << endl;
151  } else {
152  cout << " regionPUSLUT= [0]" << endl;
153  }
154 
155  cout << " pileUpTowerThreshold= " << ptr1->pileUpTowerThreshold() << endl;
156 
157  cout << endl << " EG: " << endl;
158  cout << " egLsb= " << ptr1->egLsb() << endl;
159  cout << " egSeedThreshold= " << ptr1->egSeedThreshold() << endl;
160  cout << " egNeighbourThreshold= " << ptr1->egNeighbourThreshold() << endl;
161  cout << " egHcalThreshold= " << ptr1->egHcalThreshold() << endl;
162 
163  if (!ptr1->egTrimmingLUT()->empty()) {
164  cout << " egTrimmingLUT= [" << ptr1->egTrimmingLUT()->maxSize() << "] " << flush;
165  int egTrimming[ptr1->egTrimmingLUT()->maxSize()];
166  for (unsigned int i = 0; i < ptr1->egTrimmingLUT()->maxSize(); i++)
167  egTrimming[i] = ptr1->egTrimmingLUT()->data(i);
168  cout << hash(egTrimming, sizeof(int) * ptr1->egTrimmingLUT()->maxSize()) << endl;
169  } else {
170  cout << " egTrimmingLUT= [0] " << endl;
171  }
172 
173  cout << " egMaxHcalEt= " << ptr1->egMaxHcalEt() << endl;
174  cout << " egMaxPtHOverE= " << ptr1->egMaxPtHOverE() << endl;
175  cout << " egMinPtJetIsolation= " << ptr1->egMinPtJetIsolation() << endl;
176  cout << " egMaxPtJetIsolation= " << ptr1->egMaxPtJetIsolation() << endl;
177  cout << " egMinPtHOverEIsolation= " << ptr1->egMinPtHOverEIsolation() << endl;
178  cout << " egMaxPtHOverEIsolation= " << ptr1->egMaxPtHOverEIsolation() << endl;
179 
180  if (!ptr1->egMaxHOverELUT()->empty()) {
181  cout << " egMaxHOverELUT= [" << ptr1->egMaxHOverELUT()->maxSize() << "] ";
182  int egMaxHOverE[ptr1->egMaxHOverELUT()->maxSize()];
183  for (unsigned int i = 0; i < ptr1->egMaxHOverELUT()->maxSize(); i++)
184  egMaxHOverE[i] = ptr1->egMaxHOverELUT()->data(i);
185  cout << hash(egMaxHOverE, sizeof(int) * ptr1->egMaxHOverELUT()->maxSize()) << endl;
186  } else {
187  cout << " egMaxHOverELUT= [0]" << endl;
188  }
189 
190  if (!ptr1->egCompressShapesLUT()->empty()) {
191  cout << " egCompressShapesLUT= [" << ptr1->egCompressShapesLUT()->maxSize() << "] ";
192  int egCompressShapes[ptr1->egCompressShapesLUT()->maxSize()];
193  for (unsigned int i = 0; i < ptr1->egCompressShapesLUT()->maxSize(); i++)
194  egCompressShapes[i] = ptr1->egCompressShapesLUT()->data(i);
195  cout << hash(egCompressShapes, sizeof(int) * ptr1->egCompressShapesLUT()->maxSize()) << endl;
196  } else {
197  cout << " egCompressShapesLUT= [0]" << endl;
198  }
199 
200  cout << " egShapeIdType= " << ptr1->egShapeIdType() << endl;
201  cout << " egShapeIdVersion= " << ptr1->egShapeIdVersion() << endl;
202  if (!ptr1->egShapeIdLUT()->empty()) {
203  cout << " egShapeIdLUT= [" << ptr1->egShapeIdLUT()->maxSize() << "] " << flush;
204  int egShapeId[ptr1->egShapeIdLUT()->maxSize()];
205  for (unsigned int i = 0; i < ptr1->egShapeIdLUT()->maxSize(); i++)
206  egShapeId[i] = ptr1->egShapeIdLUT()->data(i);
207  cout << hash(egShapeId, sizeof(int) * ptr1->egShapeIdLUT()->maxSize()) << endl;
208  } else {
209  cout << " egShapeIdLUT= [0]" << endl;
210  }
211 
212  cout << " egBypassEGVetos= " << ptr1->egBypassEGVetos() << endl;
213  cout << " egBypassShape= " << ptr1->egBypassShape() << endl;
214  cout << " egBypassExtHoverE= " << ptr1->egBypassExtHOverE() << endl;
215  cout << " egBypassECALFG= " << ptr1->egBypassECALFG() << endl;
216  cout << " egHOverEcutBarrel= " << ptr1->egHOverEcutBarrel() << endl;
217  cout << " egHOverEcutEndcap= " << ptr1->egHOverEcutEndcap() << endl;
218 
219  cout << " egPUSType= " << ptr1->egPUSType() << endl;
220 
221  cout << " egIsolationType= " << ptr1->egIsolationType() << endl;
222  if (!ptr1->egIsolationLUT()->empty()) {
223  cout << " egIsoLUT= [" << ptr1->egIsolationLUT()->maxSize() << "] " << flush;
224  int egIsolation[ptr1->egIsolationLUT()->maxSize()];
225  for (unsigned int i = 0; i < ptr1->egIsolationLUT()->maxSize(); i++)
226  egIsolation[i] = ptr1->egIsolationLUT()->data(i);
227  cout << hash(egIsolation, sizeof(int) * ptr1->egIsolationLUT()->maxSize()) << endl;
228  if (printEgIsoLUT)
229  for (unsigned int i = 0; i < ptr1->egIsolationLUT()->maxSize(); i++)
230  cout << i << " " << egIsolation[i] << endl;
231  } else {
232  cout << " egIsoLUT= [0]" << endl;
233  }
234  if (!ptr1->egIsolationLUT2()->empty()) {
235  cout << " egIsoLUT2= [" << ptr1->egIsolationLUT2()->maxSize() << "] " << flush;
236  int egIsolation2[ptr1->egIsolationLUT2()->maxSize()];
237  for (unsigned int i = 0; i < ptr1->egIsolationLUT2()->maxSize(); i++)
238  egIsolation2[i] = ptr1->egIsolationLUT2()->data(i);
239  cout << hash(egIsolation2, sizeof(int) * ptr1->egIsolationLUT2()->maxSize()) << endl;
240  if (printEgIsoLUT)
241  for (unsigned int i = 0; i < ptr1->egIsolationLUT2()->maxSize(); i++)
242  cout << i << " " << egIsolation2[i] << endl;
243  } else {
244  cout << " egIsoLUT2= [0]" << endl;
245  }
246 
247  cout << " egIsoAreaNrTowersEta= " << ptr1->egIsoAreaNrTowersEta() << endl;
248  cout << " egIsoAreaNrTowersPhi= " << ptr1->egIsoAreaNrTowersPhi() << endl;
249  cout << " egIsoVetoNrTowersPhi= " << ptr1->egIsoVetoNrTowersPhi() << endl;
250  cout << " egPUSParams= [" << ptr1->egPUSParams().size() << "] " << flush;
251  float egPUSParams[ptr1->egPUSParams().size()];
252  for (unsigned int i = 0; i < ptr1->egPUSParams().size(); i++)
253  egPUSParams[i] = ptr1->egPUSParams()[i];
254 
255  if (!ptr1->egPUSParams().empty())
256  cout << hash(egPUSParams, sizeof(float) * ptr1->egPUSParams().size()) << endl;
257  else
258  cout << endl;
259 
260  cout << " egCalibrationParams= [" << ptr1->egCalibrationParams().size() << "] " << flush;
261  double egCalibrationParams[ptr1->egCalibrationParams().size()];
262  for (unsigned int i = 0; i < ptr1->egCalibrationParams().size(); i++)
263  egCalibrationParams[i] = ptr1->egCalibrationParams()[i];
264 
265  if (!ptr1->egCalibrationParams().empty())
266  cout << hash(egCalibrationParams, sizeof(double) * ptr1->egCalibrationParams().size()) << endl;
267  else
268  cout << endl;
269 
270  cout << " egCalibrationType= " << ptr1->egCalibrationType() << endl;
271  cout << " egCalibrationVersion= " << ptr1->egCalibrationVersion() << endl;
272  if (!ptr1->egCalibrationLUT()->empty()) {
273  cout << " egCalibrationLUT= [" << ptr1->egCalibrationLUT()->maxSize() << "] " << flush;
274  int egCalibration[ptr1->egCalibrationLUT()->maxSize()];
275  for (unsigned int i = 0; i < ptr1->egCalibrationLUT()->maxSize(); i++)
276  egCalibration[i] = ptr1->egCalibrationLUT()->data(i);
277  cout << hash(egCalibration, sizeof(int) * ptr1->egCalibrationLUT()->maxSize()) << endl;
278  if (printEgCalibLUT)
279  for (unsigned int i = 0; i < ptr1->egCalibrationLUT()->maxSize(); i++)
280  cout << i << " " << egCalibration[i] << endl;
281  } else {
282  cout << " egCalibrationLUT= [0]" << endl;
283  }
284 
285  cout << endl << " Tau: " << endl;
286  cout << " tauLsb= " << ptr1->tauLsb() << endl;
287  //cout<<" tauSeedThreshold= "<<ptr1->tauSeedThreshold()<<endl;
288  //cout<<" tauNeighbourThreshold= "<<ptr1->tauNeighbourThreshold()<<endl;
289  cout << " tauMaxPtTauVeto= " << ptr1->tauMaxPtTauVeto() << endl;
290  cout << " tauMinPtJetIsolationB= " << ptr1->tauMinPtJetIsolationB() << endl;
291  cout << " tauPUSType= " << ptr1->tauPUSType() << endl;
292  cout << " tauMaxJetIsolationB= " << ptr1->tauMaxJetIsolationB() << endl;
293  cout << " tauMaxJetIsolationA= " << ptr1->tauMaxJetIsolationA() << endl;
294  cout << " tauIsoAreaNrTowersEta= " << ptr1->tauIsoAreaNrTowersEta() << endl;
295  cout << " tauIsoAreaNrTowersPhi= " << ptr1->tauIsoAreaNrTowersPhi() << endl;
296  cout << " tauIsoVetoNrTowersPhi= " << ptr1->tauIsoVetoNrTowersPhi() << endl;
297  if (!ptr1->tauIsolationLUT()->empty()) {
298  cout << " tauIsoLUT= [" << ptr1->tauIsolationLUT()->maxSize() << "] " << flush;
299  int tauIsolation[ptr1->tauIsolationLUT()->maxSize()];
300  for (unsigned int i = 0; i < ptr1->tauIsolationLUT()->maxSize(); i++)
301  tauIsolation[i] = ptr1->tauIsolationLUT()->data(i);
302  cout << hash(tauIsolation, sizeof(int) * ptr1->tauIsolationLUT()->maxSize()) << endl;
303  } else {
304  cout << " tauIsoLUT= [0]" << endl;
305  }
306  if (!ptr1->tauIsolationLUT2()->empty()) {
307  cout << " tauIsoLUT2= [" << ptr1->tauIsolationLUT2()->maxSize() << "] " << flush;
308  int tauIsolation2[ptr1->tauIsolationLUT2()->maxSize()];
309  for (unsigned int i = 0; i < ptr1->tauIsolationLUT2()->maxSize(); i++)
310  tauIsolation2[i] = ptr1->tauIsolationLUT2()->data(i);
311  cout << hash(tauIsolation2, sizeof(int) * ptr1->tauIsolationLUT2()->maxSize()) << endl;
312  } else {
313  cout << " tauIsoLUT2= [0]" << endl;
314  }
315  if (!ptr1->tauTrimmingShapeVetoLUT()->empty()) {
316  cout << " tauTrimmingShapeVetoLUT=[" << ptr1->tauTrimmingShapeVetoLUT()->maxSize() << "] " << flush;
317  int tauTrimmingShapeVetoLUT[ptr1->tauTrimmingShapeVetoLUT()->maxSize()];
318  for (unsigned int i = 0; i < ptr1->tauTrimmingShapeVetoLUT()->maxSize(); i++)
319  tauTrimmingShapeVetoLUT[i] = ptr1->tauTrimmingShapeVetoLUT()->data(i);
320  cout << hash(tauTrimmingShapeVetoLUT, sizeof(int) * ptr1->tauTrimmingShapeVetoLUT()->maxSize()) << endl;
321  } else {
322  cout << " tauTrimmingShapeVetoLUT=[0]" << endl;
323  }
324 
325  if (!ptr1->tauCalibrationLUT()->empty()) {
326  cout << " tauCalibrationLUT= [" << ptr1->tauCalibrationLUT()->maxSize() << "] " << flush;
327  int tauCalibration[512]; //ptr1->tauCalibrationLUT()->maxSize()];
328  for (unsigned int i = 0; i < 512; i++) //ptr1->tauCalibrationLUT()->maxSize(); i++)
329  tauCalibration[i] = ptr1->tauCalibrationLUT()->data(i);
330  cout << hash(tauCalibration, sizeof(int) * 512 /*ptr1->tauCalibrationLUT()->maxSize() */) << endl;
331 
332  if (printTauCalibLUT)
333  for (unsigned int i = 0; i < 512 /*ptr1->tauCalibrationLUT()->maxSize()*/; i++)
334  cout << i << " " << tauCalibration[i] << endl;
335 
336  } else {
337  cout << " tauCalibrationLUT= [0]" << endl;
338  }
339 
340  cout << " tauCalibrationType= " << ptr1->tauCalibrationType() << endl;
341 
342  cout << " tauCalibrationParams= [" << ptr1->tauCalibrationParams().size() << "] " << flush;
343  double tauCalibrationParams[ptr1->tauCalibrationParams().size()];
344  for (unsigned int i = 0; i < ptr1->tauCalibrationParams().size(); i++)
345  tauCalibrationParams[i] = ptr1->tauCalibrationParams()[i];
346 
347  if (!ptr1->tauCalibrationParams().empty())
348  cout << hash(tauCalibrationParams, sizeof(double) * ptr1->tauCalibrationParams().size()) << endl;
349  else
350  cout << endl;
351 
352  if (!ptr1->tauCompressLUT()->empty()) {
353  cout << " tauCompressLUT= [" << ptr1->tauCompressLUT()->maxSize() << "] " << flush;
354  int tauCompress[ptr1->tauCompressLUT()->maxSize()];
355  for (unsigned int i = 0; i < ptr1->tauCompressLUT()->maxSize(); i++)
356  tauCompress[i] = ptr1->tauCompressLUT()->data(i);
357  cout << hash(tauCompress, sizeof(int) * ptr1->tauCompressLUT()->maxSize()) << endl;
358 
359  if (printTauCompressLUT)
360  for (unsigned int i = 0; i < ptr1->tauCompressLUT()->maxSize(); i++)
361  cout << i << " " << tauCompress[i] << endl;
362 
363  } else {
364  cout << " tauCompressLUT= [0]" << endl;
365  }
366 
367  if (!ptr1->tauEtToHFRingEtLUT()->empty()) {
368  cout << " tauEtToHFRingEtLUT= [" << ptr1->tauEtToHFRingEtLUT()->maxSize() << "] " << flush;
369  int tauEtToHFRingEt[ptr1->tauEtToHFRingEtLUT()->maxSize()];
370  for (unsigned int i = 0; i < ptr1->tauEtToHFRingEtLUT()->maxSize(); i++)
371  tauEtToHFRingEt[i] = ptr1->tauEtToHFRingEtLUT()->data(i);
372 
373  cout << hash(tauEtToHFRingEt, sizeof(int) * ptr1->tauEtToHFRingEtLUT()->maxSize()) << endl;
374  } else {
375  cout << " tauEtToHFRingEtLUT= [0]" << endl;
376  }
377 
378  cout << " isoTauEtaMin= " << ptr1->isoTauEtaMin() << endl;
379  cout << " isoTauEtaMax= " << ptr1->isoTauEtaMax() << endl;
380  cout << " tauPUSParams= [" << ptr1->tauPUSParams().size() << "] " << flush;
381  float tauPUSParams[ptr1->tauPUSParams().size()];
382  for (unsigned int i = 0; i < ptr1->tauPUSParams().size(); i++)
383  tauPUSParams[i] = ptr1->tauPUSParams()[i];
384 
385  if (!ptr1->tauPUSParams().empty())
386  cout << hash(tauPUSParams, sizeof(float) * ptr1->tauPUSParams().size()) << endl;
387  else
388  cout << endl;
389 
390  cout << endl << " Jets: " << endl;
391  cout << " jetLsb= " << ptr1->jetLsb() << endl;
392  cout << " jetSeedThreshold= " << ptr1->jetSeedThreshold() << endl;
393  cout << " jetNeighbourThreshold= " << ptr1->jetNeighbourThreshold() << endl;
394  cout << " jetRegionMask= " << ptr1->jetRegionMask() << endl;
395  cout << " jetBypassPUS= " << ptr1->jetBypassPUS() << endl;
396  //cout<<" jetPUSType= "<<ptr1->jetPUSType()<<endl;
397  cout << " jetPUSUsePhiRing= " << ptr1->jetPUSUsePhiRing() << endl;
398  cout << " jetCalibrationType= " << ptr1->jetCalibrationType() << endl;
399  //cout<<" jetCalibrationParams= ["<<ptr1->jetCalibrationParams().size()<<"] "<<flush;
400  //float jetCalibrationParams[ptr1->jetCalibrationParams().size()]; // deliberately drop double precision
401  //for(unsigned int i=0; i<ptr1->jetCalibrationParams().size(); i++) jetCalibrationParams[i] = ptr1->jetCalibrationParams()[i];
402 
403  /*if( !ptr1->jetCalibrationParams().empty() ){
404  cout << hash( jetCalibrationParams, sizeof(float)*ptr1->jetCalibrationParams().size() ) << endl;
405  if( printJetCalibPar )
406  for(unsigned int i=0; i<ptr1->jetCalibrationParams().size(); i++)
407  cout<<i<<" " << std::setprecision(14) << jetCalibrationParams[i]<<endl;
408 
409  } else cout<<endl;
410 
411  cout<<" jetPUSParams= ["<<ptr1->jetPUSParams().size()<<"] "<<flush;
412  float jetPUSParams[ptr1->jetPUSParams().size()]; // deliberately drop double precision
413  for(unsigned int i=0; i<ptr1->jetPUSParams().size(); i++) jetPUSParams[i] = ptr1->jetPUSParams()[i];
414  if( !ptr1->jetPUSParams().empty() ){
415  cout << hash( jetPUSParams, sizeof(float)*ptr1->jetPUSParams().size() ) << endl;
416  if( printJetPUSPar )
417  for(unsigned int i=0; i<ptr1->jetPUSParams().size(); i++)
418  cout<<i<<" " << std::setprecision(14) << jetPUSParams[i]<<endl;
419 
420  } else cout<<endl;
421  */
422 
423  if (!ptr1->jetCalibrationLUT()->empty()) {
424  cout << " jetCalibrationLUT= [" << ptr1->jetCalibrationLUT()->maxSize() << "] " << flush;
425  int jetCalibration[ptr1->jetCalibrationLUT()->maxSize()];
426  for (unsigned int i = 0; i < ptr1->jetCalibrationLUT()->maxSize(); i++)
427  jetCalibration[i] = ptr1->jetCalibrationLUT()->data(i);
428 
429  cout << hash(jetCalibration, sizeof(int) * ptr1->jetCalibrationLUT()->maxSize()) << endl;
430 
431  if (printJetCalibLUT)
432  for (unsigned int i = 0; i < ptr1->jetCalibrationLUT()->maxSize(); i++)
433  cout << i << " " << jetCalibration[i] << endl;
434 
435  } else {
436  cout << " jetCalibrationLUT= [0]" << endl;
437  }
438 
439  if (!ptr1->jetCompressPtLUT()->empty()) {
440  cout << " jetCompressPtLUT= [" << ptr1->jetCompressPtLUT()->maxSize() << "] " << flush;
441  int jetCompressPt[ptr1->jetCompressPtLUT()->maxSize()];
442  for (unsigned int i = 0; i < ptr1->jetCompressPtLUT()->maxSize(); i++)
443  jetCompressPt[i] = ptr1->jetCompressPtLUT()->data(i);
444 
445  cout << hash(jetCompressPt, sizeof(int) * ptr1->jetCompressPtLUT()->maxSize()) << endl;
446 
447  if (printJetCompressPtLUT)
448  for (unsigned int i = 0; i < ptr1->jetCompressPtLUT()->maxSize(); i++)
449  cout << i << " " << jetCompressPt[i] << endl;
450 
451  } else {
452  cout << " jetCompressPtLUT= [0]" << endl;
453  }
454 
455  if (!ptr1->jetCompressEtaLUT()->empty()) {
456  cout << " jetCompressEtaLUT= [" << ptr1->jetCompressEtaLUT()->maxSize() << "] " << flush;
457  int jetCompressEta[ptr1->jetCompressEtaLUT()->maxSize()];
458  for (unsigned int i = 0; i < ptr1->jetCompressEtaLUT()->maxSize(); i++)
459  jetCompressEta[i] = ptr1->jetCompressEtaLUT()->data(i);
460 
461  cout << hash(jetCompressEta, sizeof(int) * ptr1->jetCompressEtaLUT()->maxSize()) << endl;
462 
463  if (printJetCompressEtaLUT)
464  for (unsigned int i = 0; i < ptr1->jetCompressEtaLUT()->maxSize(); i++)
465  cout << i << " " << jetCompressEta[i] << endl;
466 
467  } else {
468  cout << " jetCompressEtaLUT= [0]" << endl;
469  }
470 
471  cout << endl << " Sums: " << endl;
472  unsigned int nEntities = 0;
473  cout << " etSumLsb= " << ptr1->etSumLsb() << endl;
474  cout << " etSumEtaMin= [";
475  for (unsigned int i = 0; ptr1->etSumEtaMin(i) > 0.001; i++)
476  cout << (i == 0 ? "" : ",") << ptr1->etSumEtaMin(i);
477  cout << "]" << endl;
478  cout << " etSumEtaMax= [";
479  for (unsigned int i = 0; ptr1->etSumEtaMax(i) > 0.001; i++, nEntities++)
480  cout << (i == 0 ? "" : ",") << ptr1->etSumEtaMax(i);
481  cout << "]" << endl;
482  cout << " etSumEtThreshold= [";
483  for (unsigned int i = 0; i < nEntities; i++)
484  cout << (i == 0 ? "" : ",") << ptr1->etSumEtThreshold(i);
485  cout << "]" << endl;
486 
487  cout << " etSumBypassMetPUS= " << ptr1->etSumBypassMetPUS() << endl;
488  cout << " etSumBypassEttPUS= " << ptr1->etSumBypassEttPUS() << endl;
489  cout << " etSumBypassEcalSumPUS " << ptr1->etSumBypassEcalSumPUS() << endl;
490 
491  cout << " etSumMetPUSType= " << ptr1->etSumMetPUSType() << endl;
492  cout << " etSumEttPUSType= " << ptr1->etSumEttPUSType() << endl;
493  cout << " etSumEcalSumPUSType= " << ptr1->etSumEcalSumPUSType() << endl;
494 
495  cout << " etSumCentralityUpper= [";
496  for (unsigned int i = 0; ptr1->etSumCentUpper(i) > 0.001; i++)
497  cout << (i == 0 ? "" : ",") << ptr1->etSumCentUpper(i);
498  cout << "]" << endl;
499  cout << " etSumCentralityLower= [";
500  for (unsigned int i = 0; ptr1->etSumCentLower(i) > 0.001; i++)
501  cout << (i == 0 ? "" : ",") << ptr1->etSumCentLower(i);
502  cout << "]" << endl;
503 
504  cout << " metCalibrationType= " << ptr1->metCalibrationType() << endl;
505  cout << " metHFCalibrationType= " << ptr1->metHFCalibrationType() << endl;
506  cout << " etSumEttCalibrationType=" << ptr1->etSumEttCalibrationType() << endl;
507  cout << " etSumEcalSumCalibrationType=" << ptr1->etSumEcalSumCalibrationType() << endl;
508 
509  if (!ptr1->etSumMetPUSLUT()->empty()) {
510  cout << " etSumMetPUSLUT= [" << ptr1->etSumMetPUSLUT()->maxSize() << "] " << flush;
511  int etSumMetPUSLUT[ptr1->etSumMetPUSLUT()->maxSize()];
512  for (unsigned int i = 0; i < ptr1->etSumMetPUSLUT()->maxSize(); i++)
513  etSumMetPUSLUT[i] = ptr1->etSumMetPUSLUT()->data(i);
514 
515  cout << hash(etSumMetPUSLUT, sizeof(int) * ptr1->etSumMetPUSLUT()->maxSize()) << endl;
516 
517  if (printEtSumMetPUSLUT)
518  for (unsigned int i = 0; i < ptr1->etSumMetPUSLUT()->maxSize(); i++)
519  cout << i << " " << etSumMetPUSLUT[i] << endl;
520 
521  } else {
522  cout << " etSumMetPUSLUT= [0]" << endl;
523  }
524 
525  if (!ptr1->etSumEttPUSLUT()->empty()) {
526  cout << " etSumEttPUSLUT= [" << ptr1->etSumEttPUSLUT()->maxSize() << "] " << flush;
527  int etSumEttPUSLUT[ptr1->etSumEttPUSLUT()->maxSize()];
528  for (unsigned int i = 0; i < ptr1->etSumEttPUSLUT()->maxSize(); i++)
529  etSumEttPUSLUT[i] = ptr1->etSumEttPUSLUT()->data(i);
530 
531  cout << hash(etSumEttPUSLUT, sizeof(int) * ptr1->etSumEttPUSLUT()->maxSize()) << endl;
532 
533  if (printEtSumEttPUSLUT)
534  for (unsigned int i = 0; i < ptr1->etSumEttPUSLUT()->maxSize(); i++)
535  cout << i << " " << etSumEttPUSLUT[i] << endl;
536 
537  } else {
538  cout << " etSumEttPUSLUT= [0]" << endl;
539  }
540 
541  if (!ptr1->etSumEcalSumPUSLUT()->empty()) {
542  cout << " etSumEcalSumPUSLUT= [" << ptr1->etSumEcalSumPUSLUT()->maxSize() << "] " << flush;
543  int etSumEcalSumPUSLUT[ptr1->etSumEcalSumPUSLUT()->maxSize()];
544  for (unsigned int i = 0; i < ptr1->etSumEcalSumPUSLUT()->maxSize(); i++)
545  etSumEcalSumPUSLUT[i] = ptr1->etSumEcalSumPUSLUT()->data(i);
546 
547  cout << hash(etSumEcalSumPUSLUT, sizeof(int) * ptr1->etSumEcalSumPUSLUT()->maxSize()) << endl;
548 
549  if (printEtSumEcalSumPUSLUT)
550  for (unsigned int i = 0; i < ptr1->etSumEcalSumPUSLUT()->maxSize(); i++)
551  cout << i << " " << etSumEcalSumPUSLUT[i] << endl;
552 
553  } else {
554  cout << " etSumEcalSumPUSLUT= [0]" << endl;
555  }
556 
557  if (!ptr1->metCalibrationLUT()->empty()) {
558  cout << " metCalibrationLUT= [" << ptr1->metCalibrationLUT()->maxSize() << "] " << flush;
559  int metCalibrationLUT[ptr1->metCalibrationLUT()->maxSize()];
560  for (unsigned int i = 0; i < ptr1->metCalibrationLUT()->maxSize(); i++)
561  metCalibrationLUT[i] = ptr1->metCalibrationLUT()->data(i);
562 
563  cout << hash(metCalibrationLUT, sizeof(int) * ptr1->metCalibrationLUT()->maxSize()) << endl;
564 
565  if (printMetCalibrationLUT)
566  for (unsigned int i = 0; i < ptr1->metCalibrationLUT()->maxSize(); i++)
567  cout << i << " " << metCalibrationLUT[i] << endl;
568 
569  } else {
570  cout << " metCalibrationLUT= [0]" << endl;
571  }
572 
573  if (!ptr1->metHFCalibrationLUT()->empty()) {
574  cout << " metHFCalibrationLUT= [" << ptr1->metHFCalibrationLUT()->maxSize() << "] " << flush;
575  int metHFCalibrationLUT[ptr1->metHFCalibrationLUT()->maxSize()];
576  for (unsigned int i = 0; i < ptr1->metHFCalibrationLUT()->maxSize(); i++)
577  metHFCalibrationLUT[i] = ptr1->metHFCalibrationLUT()->data(i);
578 
579  cout << hash(metHFCalibrationLUT, sizeof(int) * ptr1->metHFCalibrationLUT()->maxSize()) << endl;
580 
581  if (printMetHFCalibrationLUT)
582  for (unsigned int i = 0; i < ptr1->metHFCalibrationLUT()->maxSize(); i++)
583  cout << i << " " << metHFCalibrationLUT[i] << endl;
584 
585  } else {
586  cout << " metHFCalibrationLUT= [0]" << endl;
587  }
588 
589  if (!ptr1->metPhiCalibrationLUT()->empty()) {
590  cout << " metPhiCalibrationLUT= [" << ptr1->metPhiCalibrationLUT()->maxSize() << "] " << flush;
591  int metPhiCalibrationLUT[ptr1->metPhiCalibrationLUT()->maxSize()];
592  for (unsigned int i = 0; i < ptr1->metPhiCalibrationLUT()->maxSize(); i++)
593  metPhiCalibrationLUT[i] = ptr1->metPhiCalibrationLUT()->data(i);
594 
595  cout << hash(metPhiCalibrationLUT, sizeof(int) * ptr1->metPhiCalibrationLUT()->maxSize()) << endl;
596 
597  if (printMetPhiCalibrationLUT)
598  for (unsigned int i = 0; i < ptr1->metPhiCalibrationLUT()->maxSize(); i++)
599  cout << i << " " << metPhiCalibrationLUT[i] << endl;
600 
601  } else {
602  cout << " metPhiCalibrationLUT= [0]" << endl;
603  }
604 
605  if (!ptr1->metHFPhiCalibrationLUT()->empty()) {
606  cout << " metHFPhiCalibrationLUT= [" << ptr1->metHFPhiCalibrationLUT()->maxSize() << "] " << flush;
607  int metHFPhiCalibrationLUT[ptr1->metHFPhiCalibrationLUT()->maxSize()];
608  for (unsigned int i = 0; i < ptr1->metHFPhiCalibrationLUT()->maxSize(); i++)
609  metHFPhiCalibrationLUT[i] = ptr1->metHFPhiCalibrationLUT()->data(i);
610 
611  cout << hash(metHFPhiCalibrationLUT, sizeof(int) * ptr1->metHFPhiCalibrationLUT()->maxSize()) << endl;
612 
613  if (printMetHFCalibrationLUT)
614  for (unsigned int i = 0; i < ptr1->metHFPhiCalibrationLUT()->maxSize(); i++)
615  cout << i << " " << metHFPhiCalibrationLUT[i] << endl;
616 
617  } else {
618  cout << " metHFPhiCalibrationLUT= [0]" << endl;
619  }
620 
621  if (!ptr1->etSumEttCalibrationLUT()->empty()) {
622  cout << " etSumEttCalibrationLUT= [" << ptr1->etSumEttCalibrationLUT()->maxSize() << "] " << flush;
623  int etSumEttCalibrationLUT[ptr1->etSumEttCalibrationLUT()->maxSize()];
624  for (unsigned int i = 0; i < ptr1->etSumEttCalibrationLUT()->maxSize(); i++)
625  etSumEttCalibrationLUT[i] = ptr1->etSumEttCalibrationLUT()->data(i);
626 
627  cout << hash(etSumEttCalibrationLUT, sizeof(int) * ptr1->etSumEttCalibrationLUT()->maxSize()) << endl;
628 
629  if (printEtSumEttCalibrationLUT)
630  for (unsigned int i = 0; i < ptr1->etSumEttCalibrationLUT()->maxSize(); i++)
631  cout << i << " " << etSumEttCalibrationLUT[i] << endl;
632 
633  } else {
634  cout << " etSumEttCalibrationLUT= [0]" << endl;
635  }
636 
637  if (!ptr1->etSumEcalSumCalibrationLUT()->empty()) {
638  cout << " etSumEcalSumCalibrationLUT=[" << ptr1->etSumEttCalibrationLUT()->maxSize() << "] " << flush;
639  int etSumEcalSumCalibrationLUT[ptr1->etSumEcalSumCalibrationLUT()->maxSize()];
640  for (unsigned int i = 0; i < ptr1->etSumEcalSumCalibrationLUT()->maxSize(); i++)
641  etSumEcalSumCalibrationLUT[i] = ptr1->etSumEcalSumCalibrationLUT()->data(i);
642 
643  cout << hash(etSumEcalSumCalibrationLUT, sizeof(int) * ptr1->etSumEcalSumCalibrationLUT()->maxSize()) << endl;
644 
645  if (printEtSumEcalSumCalibrationLUT)
646  for (unsigned int i = 0; i < ptr1->etSumEcalSumCalibrationLUT()->maxSize(); i++)
647  cout << i << " " << etSumEcalSumCalibrationLUT[i] << endl;
648 
649  } else {
650  cout << " etSumEcalSumCalibrationLUT=[0]" << endl;
651  }
652 
653  cout << endl << " HI centrality trigger: " << endl;
654  cout << " centralityLUT= [";
655  for (unsigned int i = 0; i < ptr1->centralityLUT()->maxSize(); i++)
656  cout << (i == 0 ? "" : ",") << ptr1->centralityLUT()->data(i);
657  cout << "]" << endl;
658 
659  std::vector<int> mbt = ptr1->minimumBiasThresholds();
660  cout << " minimumBiasThresholds= [";
661  for (unsigned int i = 0; i < mbt.size(); i++)
662  cout << mbt[i];
663  cout << "]" << endl;
664 
665  cout << endl << "centralityRegionMask() = " << ptr1->centralityRegionMask() << endl;
666  cout << endl << "jetRegionMask() = " << ptr1->jetRegionMask() << endl;
667  cout << endl << "tauRegionMask() = " << ptr1->tauRegionMask() << endl;
668 
669  cout << endl << " HI Q2 trigger: " << endl;
670  cout << " q2LUT= [";
671  for (unsigned int i = 0; i < ptr1->q2LUT()->maxSize(); i++)
672  cout << (i == 0 ? "" : ",") << ptr1->q2LUT()->data(i);
673  cout << "]" << endl;
674 
675  cout << endl << " Layer1: " << endl;
676  std::vector<double> ecalSF = ptr1->layer1ECalScaleFactors();
677  cout << " layer1ECalScaleFactors= [" << ecalSF.size() << "] " << flush;
678  int _ecalSF[ecalSF.size()];
679  for (unsigned int i = 0; i < ecalSF.size(); i++)
680  _ecalSF[i] = int(ecalSF[i] * 100000.);
681  cout << hash(_ecalSF, sizeof(int) * ecalSF.size()) << endl;
682  if (printEcalSF) {
683  cout << endl << " [" << endl;
684  for (unsigned int i = 0; i < ecalSF.size(); i++)
685  cout << (i == 0 ? "" : ",") << int(ecalSF[i] * 1000.) / 1000.;
686  cout << "]" << endl;
687  }
688  std::vector<double> hcalSF = ptr1->layer1HCalScaleFactors();
689  cout << " layer1HCalScaleFactors= [" << hcalSF.size() << "] " << flush;
690  int _hcalSF[hcalSF.size()];
691  for (unsigned int i = 0; i < hcalSF.size(); i++) {
692  // round false precision
693  // double significand;
694  // int exponent;
695  // significand = frexp( hcalSF[i], &exponent );
696  // _hcalSF[i] = ldexp( int(significand*10000)/10000., exponent );
697  _hcalSF[i] = int(hcalSF[i] * 100000.);
698  }
699  cout << hash(_hcalSF, sizeof(int) * hcalSF.size()) << endl;
700  if (printHcalSF) {
701  cout << endl << " [" << endl;
702  for (unsigned int i = 0; i < hcalSF.size(); i++)
703  cout << (i == 0 ? "" : ",") << int(hcalSF[i] * 1000.) / 1000.;
704  cout << "]" << endl;
705  }
706  std::vector<double> hfSF = ptr1->layer1HFScaleFactors();
707  cout << " layer1HFScaleFactors= [" << hfSF.size() << "] " << flush;
708  int _hfSF[hfSF.size()];
709  for (unsigned int i = 0; i < hfSF.size(); i++)
710  _hfSF[i] = int(hfSF[i] * 100000.);
711  cout << hash(_hfSF, sizeof(int) * hfSF.size()) << endl;
712  if (printHfSF) {
713  cout << endl << " [" << endl;
714  for (unsigned int i = 0; i < hfSF.size(); i++)
715  cout << (i == 0 ? "" : ",") << int(hfSF[i] * 1000.) / 1000.;
716  cout << "]" << endl;
717  }
718 
719  std::vector<int> ecalScaleET = ptr1->layer1ECalScaleETBins();
720  cout << " layer1ECalScaleETBins= [";
721  for (unsigned int i = 0; i < ecalScaleET.size(); i++)
722  cout << (i == 0 ? "" : ",") << ecalScaleET[i];
723  cout << "]" << endl;
724  std::vector<int> hcalScaleET = ptr1->layer1HCalScaleETBins();
725  cout << " layer1HCalScaleETBins= [";
726  for (unsigned int i = 0; i < hcalScaleET.size(); i++)
727  cout << (i == 0 ? "" : ",") << hcalScaleET[i];
728  cout << "]" << endl;
729  std::vector<int> hfScaleET = ptr1->layer1HFScaleETBins();
730  cout << " layer1HFScaleETBins= [";
731  for (unsigned int i = 0; i < hfScaleET.size(); i++)
732  cout << (i == 0 ? "" : ",") << hfScaleET[i];
733  cout << "]" << endl;
734 
735  std::vector<unsigned> layer1ECalScalePhi = ptr1->layer1ECalScalePhiBins();
736  cout << " layer1ECalScalePhi= [";
737  for (unsigned int i = 0; i < layer1ECalScalePhi.size(); i++)
738  cout << (i == 0 ? "" : ",") << layer1ECalScalePhi[i];
739  cout << "]" << endl;
740  std::vector<unsigned> layer1HCalScalePhi = ptr1->layer1HCalScalePhiBins();
741  cout << " layer1HCalScalePhi= [";
742  for (unsigned int i = 0; i < layer1HCalScalePhi.size(); i++)
743  cout << (i == 0 ? "" : ",") << layer1HCalScalePhi[i];
744  cout << "]" << endl;
745  std::vector<unsigned> layer1HFScalePhiBins = ptr1->layer1HFScalePhiBins();
746  cout << " layer1HFScalePhiBins= [";
747  for (unsigned int i = 0; i < layer1HFScalePhiBins.size(); i++)
748  cout << (i == 0 ? "" : ",") << layer1HFScalePhiBins[i];
749  cout << "]" << endl;
750 
751  // std::vector<unsigned> layer1SecondStageLUT = ptr1->layer1SecondStageLUT();
752  // cout<<" layer1HFScalePhiBins= ["; for(unsigned int i=0; i<layer1SecondStageLUT.size(); i++) cout<<(i==0?"":",")<<layer1SecondStageLUT[i]; cout<<"]"<<endl;
753 }
754 
758 
constexpr int32_t ceil(float num)
T getUntrackedParameter(std::string const &, T const &) const
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
int iEvent
Definition: GenABIO.cc:224
~L1TCaloParamsViewer(void) override
L1TCaloParamsViewer(const edm::ParameterSet &pset)
Log< level::Info, false > LogInfo
T const * product() const
Definition: ESHandle.h:86
T get() const
Definition: EventSetup.h:82
tuple cout
Definition: gather_cfg.py:144
std::string hash(void *buf, size_t len) const
tmp
align.sh
Definition: createJobs.py:716
void analyze(const edm::Event &, const edm::EventSetup &) override