40 #include <openssl/sha.h> 46 char tmp[SHA_DIGEST_LENGTH*2+1];
47 bzero(tmp,
sizeof(tmp));
49 if( !SHA1_Init( &ctx ) )
50 throw cms::Exception(
"L1TCaloParamsViewer::hash")<<
"SHA1 initialization error";
52 if( !SHA1_Update( &ctx, buf, len ) )
53 throw cms::Exception(
"L1TCaloParamsViewer::hash")<<
"SHA1 processing error";
55 unsigned char hash[SHA_DIGEST_LENGTH];
56 if( !SHA1_Final(hash, &ctx) )
57 throw cms::Exception(
"L1TCaloParamsViewer::hash")<<
"SHA1 finalization error";
60 for(
unsigned int i=0;
i<20;
i++)
61 ::sprintf(&tmp[
i*2],
"%02x", hash[i]);
78 cout<<endl<<
" Towers: "<<endl;
79 cout<<
" towerLsbH= "<<ptr1->towerLsbH()<<endl;
80 cout<<
" towerLsbE= "<<ptr1->towerLsbE()<<endl;
81 cout<<
" towerLsbSum= "<<ptr1->towerLsbSum()<<endl;
82 cout<<
" towerNBitsH= "<<ptr1->towerNBitsH()<<endl;
83 cout<<
" towerNBitsE= "<<ptr1->towerNBitsE()<<endl;
84 cout<<
" towerNBitsSum= "<<ptr1->towerNBitsSum()<<endl;
85 cout<<
" towerNBitsRatio= "<<ptr1->towerNBitsRatio()<<endl;
86 cout<<
" towerMaskE= "<<ptr1->towerMaskE()<<endl;
87 cout<<
" towerMaskH= "<<ptr1->towerMaskH()<<endl;
88 cout<<
" towerMaskSum= "<<ptr1->towerMaskSum()<<endl;
89 cout<<
" towerEncoding= "<<ptr1->doTowerEncoding()<<endl;
91 cout<<endl<<
" Regions: "<<endl;
92 cout<<
" regionLsb= "<<ptr1->regionLsb()<<endl;
93 cout<<
" regionPUSType= "<<ptr1->regionPUSType()<<endl;
94 cout<<
" regionPUSParams= ["<<ptr1->regionPUSParams().size()<<
"] ";
95 float pusParams[ptr1->regionPUSParams().size()];
97 for(
unsigned int i=0;
i<ptr1->regionPUSParams().size();
i++){
98 pusParams[
i] = ceil(2*ptr1->regionPUSParams()[
i]);
102 if( !ptr1->regionPUSParams().empty() )
103 cout <<
hash(pusParams,
sizeof(
float)*ptr1->regionPUSParams().size()) << endl;
107 cout<<endl<<
" EG: "<<endl;
108 cout<<
" egLsb= "<<ptr1->egLsb()<<endl;
109 cout<<
" egSeedThreshold= "<<ptr1->egSeedThreshold()<<endl;
110 cout<<
" egNeighbourThreshold= "<<ptr1->egNeighbourThreshold()<<endl;
111 cout<<
" egHcalThreshold= "<<ptr1->egHcalThreshold()<<endl;
113 if( !ptr1->egTrimmingLUT()->empty() ){
114 cout<<
" egTrimmingLUT= ["<<ptr1->egTrimmingLUT()->maxSize()<<
"] "<<flush;
115 int egTrimming[ptr1->egTrimmingLUT()->maxSize()];
116 for(
unsigned int i=0;
i<ptr1->egTrimmingLUT()->maxSize();
i++) egTrimming[
i] = ptr1->egTrimmingLUT()->data(
i);
117 cout <<
hash( egTrimming,
sizeof(
int)*ptr1->egTrimmingLUT()->maxSize() ) << endl;
119 cout<<
" egTrimmingLUT= [0] "<<endl;
122 cout<<
" egMaxHcalEt= "<<ptr1->egMaxHcalEt()<<endl;
123 cout<<
" egMaxPtHOverE= "<<ptr1->egMaxPtHOverE()<<endl;
124 cout<<
" egMinPtJetIsolation= "<<ptr1->egMinPtJetIsolation()<<endl;
125 cout<<
" egMaxPtJetIsolation= "<<ptr1->egMaxPtJetIsolation()<<endl;
126 cout<<
" egMinPtHOverEIsolation= "<<ptr1->egMinPtHOverEIsolation()<<endl;
127 cout<<
" egMaxPtHOverEIsolation= "<<ptr1->egMaxPtHOverEIsolation()<<endl;
129 if( !ptr1->egMaxHOverELUT()->empty() ){
130 cout<<
" egMaxHOverELUT= ["<<ptr1->egMaxHOverELUT()->maxSize()<<
"] ";
131 int egMaxHOverE[ptr1->egMaxHOverELUT()->maxSize()];
132 for(
unsigned int i=0;
i<ptr1->egMaxHOverELUT()->maxSize();
i++) egMaxHOverE[
i] = ptr1->egMaxHOverELUT()->data(
i);
133 cout <<
hash( egMaxHOverE,
sizeof(
int)*ptr1->egMaxHOverELUT()->maxSize() ) << endl;
135 cout<<
" egMaxHOverELUT= [0]"<<endl;
138 if( !ptr1->egCompressShapesLUT()->empty() ){
139 cout<<
" egCompressShapesLUT= ["<<ptr1->egCompressShapesLUT()->maxSize()<<
"] ";
140 int egCompressShapes[ptr1->egCompressShapesLUT()->maxSize()];
141 for(
unsigned int i=0;
i<ptr1->egCompressShapesLUT()->maxSize();
i++) egCompressShapes[
i] = ptr1->egCompressShapesLUT()->data(
i);
142 cout <<
hash( egCompressShapes,
sizeof(
int)*ptr1->egCompressShapesLUT()->maxSize() ) << endl;
144 cout<<
" egCompressShapesLUT= [0]"<<endl;
149 if( !ptr1->egShapeIdLUT()->empty() ){
150 cout<<
" egShapeIdLUT= ["<<ptr1->egShapeIdLUT()->maxSize()<<
"] "<<flush;
151 int egShapeId[ptr1->egShapeIdLUT()->maxSize()];
152 for(
unsigned int i=0;
i<ptr1->egShapeIdLUT()->maxSize();
i++) egShapeId[
i] = ptr1->egShapeIdLUT()->data(
i);
153 cout <<
hash( egShapeId,
sizeof(
int)*ptr1->egShapeIdLUT()->maxSize() )<<endl;
155 cout<<
" egShapeIdLUT= [0]"<<endl;
158 cout<<
" egBypassEGVetos= "<<ptr1->egBypassEGVetos()<<endl;
159 cout<<
" egHOverEcutBarrel= "<<ptr1->egHOverEcutBarrel()<<endl;
160 cout<<
" egHOverEcutEndcap= "<<ptr1->egHOverEcutEndcap()<<endl;
162 cout<<
" egPUSType= "<<ptr1->egPUSType()<<endl;
164 cout<<
" egIsolationType= "<<ptr1->egIsolationType()<<endl;
165 if( !ptr1->egIsolationLUT()->empty() ){
166 cout<<
" egIsoLUT= ["<<ptr1->egIsolationLUT()->maxSize()<<
"] "<<flush;
167 int egIsolation[ptr1->egIsolationLUT()->maxSize()];
168 for(
unsigned int i=0;
i<ptr1->egIsolationLUT()->maxSize();
i++) egIsolation[
i] = ptr1->egIsolationLUT()->data(
i);
169 cout <<
hash( egIsolation,
sizeof(
int)*ptr1->egIsolationLUT()->maxSize() ) << endl;
171 for(
unsigned int i=0;
i<ptr1->egIsolationLUT()->maxSize();
i++)
172 cout<<
i<<
" " << egIsolation[
i]<<endl;
174 cout<<
" egIsoLUT= [0]"<<endl;
177 cout<<
" egIsoAreaNrTowersEta= "<<ptr1->egIsoAreaNrTowersEta()<<endl;
178 cout<<
" egIsoAreaNrTowersPhi= "<<ptr1->egIsoAreaNrTowersPhi()<<endl;
179 cout<<
" egIsoVetoNrTowersPhi= "<<ptr1->egIsoVetoNrTowersPhi()<<endl;
180 cout<<
" egPUSParams= ["<<ptr1->egPUSParams().size()<<
"] "<<flush;
181 float egPUSParams[ptr1->egPUSParams().size()];
182 for(
unsigned int i=0;
i<ptr1->egPUSParams().size();
i++) egPUSParams[
i] = ptr1->egPUSParams()[
i];
184 if( !ptr1->egPUSParams().empty() )
185 cout <<
hash( egPUSParams,
sizeof(
float)*ptr1->egPUSParams().size() ) << endl;
188 cout<<
" egCalibrationParams= ["<<ptr1->egCalibrationParams().size()<<
"] "<<flush;
189 double egCalibrationParams[ptr1->egCalibrationParams().size()];
190 for(
unsigned int i=0;
i<ptr1->egCalibrationParams().size();
i++) egCalibrationParams[
i] = ptr1->egCalibrationParams()[
i];
192 if( !ptr1->egCalibrationParams().empty() )
193 cout <<
hash( egCalibrationParams,
sizeof(
double)*ptr1->egCalibrationParams().size() ) << endl;
196 cout<<
" egCalibrationType= "<<ptr1->egCalibrationType()<<endl;
198 if( !ptr1->egCalibrationLUT()->empty() ){
199 cout<<
" egCalibrationLUT= ["<<ptr1->egCalibrationLUT()->maxSize()<<
"] "<<flush;
200 int egCalibration[ptr1->egCalibrationLUT()->maxSize()];
201 for(
unsigned int i=0;
i<ptr1->egCalibrationLUT()->maxSize();
i++) egCalibration[
i] = ptr1->egCalibrationLUT()->data(
i);
202 cout <<
hash( egCalibration,
sizeof(
int)*ptr1->egCalibrationLUT()->maxSize() ) << endl;
204 for(
unsigned int i=0;
i<ptr1->egCalibrationLUT()->maxSize();
i++)
205 cout<<
i<<
" " << egCalibration[
i]<<endl;
207 cout<<
" egCalibrationLUT= [0]"<<endl;
210 cout<<endl<<
" Tau: "<<endl;
211 cout<<
" tauLsb= "<<ptr1->tauLsb()<<endl;
212 cout<<
" tauSeedThreshold= "<<ptr1->tauSeedThreshold()<<endl;
213 cout<<
" tauNeighbourThreshold= "<<ptr1->tauNeighbourThreshold()<<endl;
214 cout<<
" tauMaxPtTauVeto= "<<ptr1->tauMaxPtTauVeto()<<endl;
215 cout<<
" tauMinPtJetIsolationB= "<<ptr1->tauMinPtJetIsolationB()<<endl;
216 cout<<
" tauPUSType= "<<ptr1->tauPUSType()<<endl;
217 cout<<
" tauMaxJetIsolationB= "<<ptr1->tauMaxJetIsolationB()<<endl;
218 cout<<
" tauMaxJetIsolationA= "<<ptr1->tauMaxJetIsolationA()<<endl;
219 cout<<
" tauIsoAreaNrTowersEta= "<<ptr1->tauIsoAreaNrTowersEta()<<endl;
220 cout<<
" tauIsoAreaNrTowersPhi= "<<ptr1->tauIsoAreaNrTowersPhi()<<endl;
221 cout<<
" tauIsoVetoNrTowersPhi= "<<ptr1->tauIsoVetoNrTowersPhi()<<endl;
222 if( !ptr1->tauIsolationLUT()->empty() ){
223 cout<<
" tauIsoLUT= ["<<ptr1->tauIsolationLUT()->maxSize()<<
"] "<<flush;
224 int tauIsolation[ptr1->tauIsolationLUT()->maxSize()];
225 for(
unsigned int i=0;
i<ptr1->tauIsolationLUT()->maxSize();
i++) tauIsolation[
i] = ptr1->tauIsolationLUT()->data(
i);
226 cout <<
hash( tauIsolation,
sizeof(
int)*ptr1->tauIsolationLUT()->maxSize() ) << endl;
228 cout<<
" tauIsoLUT= [0]"<<endl;
231 if( !ptr1->tauCalibrationLUT()->empty() ){
232 cout<<
" tauCalibrationLUT= ["<<ptr1->tauCalibrationLUT()->maxSize()<<
"] "<<flush;
233 int tauCalibration[512];
234 for(
unsigned int i=0;
i<512;
i++)
235 tauCalibration[
i] = ptr1->tauCalibrationLUT()->data(
i);
236 cout <<
hash( tauCalibration,
sizeof(
int)*512) << endl;
239 for(
unsigned int i=0;
i<512;
i++)
240 cout<<
i<<
" "<<tauCalibration[
i]<<endl;
243 cout<<
" tauCalibrationLUT= [0]"<<endl;
246 if( !ptr1->tauEtToHFRingEtLUT()->empty() ){
247 cout<<
" tauEtToHFRingEtLUT= ["<<ptr1->tauEtToHFRingEtLUT()->maxSize()<<
"] "<<flush;
248 int tauEtToHFRingEt[ptr1->tauEtToHFRingEtLUT()->maxSize()];
249 for(
unsigned int i=0;
i<ptr1->tauEtToHFRingEtLUT()->maxSize();
i++) tauEtToHFRingEt[
i] = ptr1->tauEtToHFRingEtLUT()->data(
i);
251 cout <<
hash( tauEtToHFRingEt,
sizeof(
int)*ptr1->tauEtToHFRingEtLUT()->maxSize() ) << endl;
253 cout<<
" tauEtToHFRingEtLUT= [0]"<<endl;
256 cout<<
" isoTauEtaMin= "<<ptr1->isoTauEtaMin()<<endl;
257 cout<<
" isoTauEtaMax= "<<ptr1->isoTauEtaMax()<<endl;
258 cout<<
" tauPUSParams= ["<<ptr1->tauPUSParams().size()<<
"] "<<flush;
259 float tauPUSParams[ptr1->tauPUSParams().size()];
260 for(
unsigned int i=0;
i<ptr1->tauPUSParams().size();
i++) tauPUSParams[
i] = ptr1->tauPUSParams()[
i];
262 if( !ptr1->tauPUSParams().empty() )
263 cout <<
hash( tauPUSParams,
sizeof(
float)*ptr1->tauPUSParams().size() ) << endl;
266 cout<<endl<<
" Jets: "<<endl;
267 cout<<
" jetLsb= "<<ptr1->jetLsb()<<endl;
268 cout<<
" jetSeedThreshold= "<<ptr1->jetSeedThreshold()<<endl;
269 cout<<
" jetNeighbourThreshold= "<<ptr1->jetNeighbourThreshold()<<endl;
270 cout<<
" jetPUSType= "<<ptr1->jetPUSType()<<endl;
271 cout<<
" jetCalibrationType= "<<ptr1->jetCalibrationType()<<endl;
272 cout<<
" jetCalibrationParams= ["<<ptr1->jetCalibrationParams().size()<<
"] "<<flush;
273 float jetCalibrationParams[ptr1->jetCalibrationParams().size()];
274 for(
unsigned int i=0;
i<ptr1->jetCalibrationParams().size();
i++) jetCalibrationParams[
i] = ptr1->jetCalibrationParams()[
i];
276 if( !ptr1->jetCalibrationParams().empty() ){
277 cout <<
hash( jetCalibrationParams,
sizeof(
float)*ptr1->jetCalibrationParams().size() ) << endl;
279 for(
unsigned int i=0;
i<ptr1->jetCalibrationParams().size();
i++)
280 cout<<
i<<
" " << std::setprecision(14) << jetCalibrationParams[
i]<<endl;
284 if( !ptr1->jetCalibrationLUT()->empty() ){
285 cout<<
" jetCalibrationLUT= ["<<ptr1->jetCalibrationLUT()->maxSize()<<
"] "<<flush;
286 int jetCalibration[ptr1->jetCalibrationLUT()->maxSize()];
287 for(
unsigned int i=0;
i<ptr1->jetCalibrationLUT()->maxSize();
i++)
288 jetCalibration[
i] = ptr1->jetCalibrationLUT()->data(
i);
290 cout <<
hash( jetCalibration,
sizeof(
int)*ptr1->jetCalibrationLUT()->maxSize() ) << endl;
293 for(
unsigned int i=0;
i<ptr1->jetCalibrationLUT()->maxSize();
i++)
294 cout<<
i<<
" "<<jetCalibration[
i]<<endl;
297 cout<<
" jetCalibrationLUT= [0]"<<endl;
300 cout<<endl<<
" Sums: "<<endl;
301 unsigned int nEntities = 0;
302 cout<<
" etSumLsb= "<<ptr1->etSumLsb()<<endl;
303 cout<<
" etSumEtaMin= [";
for(
unsigned int i=0; ptr1->etSumEtaMin(
i)>0.001;
i++)
cout<<(
i==0?
"":
",")<<ptr1->etSumEtaMin(
i);
cout<<
"]"<<endl;
304 cout<<
" etSumEtaMax= [";
for(
unsigned int i=0; ptr1->etSumEtaMax(
i)>0.001;
i++,nEntities++)
cout<<(
i==0?
"":
",")<<ptr1->etSumEtaMax(
i);
cout<<
"]"<<endl;
305 cout<<
" etSumEtThreshold= [";
for(
unsigned int i=0;
i<nEntities;
i++)
cout<<(
i==0?
"":
",")<<ptr1->etSumEtThreshold(
i);
cout<<
"]"<<endl;
307 cout<<endl<<
" HI centrality trigger: "<<endl;
308 cout<<
" centralityLUT= [";
for(
unsigned int i=0;
i<ptr1->centralityLUT()->maxSize();
i++)
cout<<(
i==0?
"":
",")<<ptr1->centralityLUT()->data(
i);
cout<<
"]"<<endl;
310 std::vector<int> mbt = ptr1->minimumBiasThresholds();
311 cout<<
" minimumBiasThresholds= [";
for(
unsigned int i=0;
i<mbt.size();
i++)
cout<<mbt[
i];
cout<<
"]"<<endl;
313 cout<<endl<<
"centralityRegionMask() = "<<ptr1->centralityRegionMask()<<endl;
314 cout<<endl<<
"jetRegionMask() = "<<ptr1->jetRegionMask()<<endl;
315 cout<<endl<<
"tauRegionMask() = "<<ptr1->tauRegionMask()<<endl;
317 cout<<endl<<
" HI Q2 trigger: "<<endl;
318 cout<<
" q2LUT= [";
for(
unsigned int i=0;
i<ptr1->q2LUT()->maxSize();
i++)
cout<<(
i==0?
"":
",")<<ptr1->q2LUT()->data(
i);
cout<<
"]"<<endl;
320 cout<<endl<<
" Layer1: "<<endl;
321 std::vector<double> ecalSF = ptr1->layer1ECalScaleFactors();
322 cout<<
" layer1ECalScaleFactors= ["<< ecalSF.size()<<
"] "<<flush;
323 float _ecalSF[ ecalSF.size() ];
324 for(
unsigned int i=0;
i<ecalSF.size();
i++) _ecalSF[
i] = ecalSF[
i];
325 cout <<
hash( _ecalSF ,
sizeof(
float)*ecalSF.size() ) << endl;
327 std::vector<double> hcalSF = ptr1->layer1HCalScaleFactors();
328 cout<<
" layer1HCalScaleFactors= ["<< hcalSF.size()<<
"] "<<flush;
329 float _hcalSF[ hcalSF.size() ];
330 for(
unsigned int i=0;
i<hcalSF.size();
i++){
336 _hcalSF[
i] =
int(hcalSF[
i]*1000.)/1000.;
338 cout <<
hash( _hcalSF,
sizeof(
float)*hcalSF.size() ) << endl;
340 std::vector<double> hfSF = ptr1->layer1HFScaleFactors();
341 cout<<
" layer1HFScaleFactors= ["<< hfSF.size()<<
"] "<<flush;
342 float _hfSF[ hfSF.size() ];
343 for(
unsigned int i=0;
i<hfSF.size();
i++) _hfSF[
i] = hfSF[
i];
344 cout <<
hash( _hfSF,
sizeof(
float)*hfSF.size() ) << endl;
347 std::vector<int> ecalScaleET = ptr1->layer1ECalScaleETBins();
348 cout<<
" layer1ECalScaleETBins= [";
for(
unsigned int i=0;
i<ecalScaleET.size();
i++)
cout<<(
i==0?
"":
",")<<ecalScaleET[
i];
cout<<
"]"<<endl;
349 std::vector<int> hcalScaleET = ptr1->layer1HCalScaleETBins();
350 cout<<
" layer1HCalScaleETBins= [";
for(
unsigned int i=0;
i<hcalScaleET.size();
i++)
cout<<(
i==0?
"":
",")<<hcalScaleET[
i];
cout<<
"]"<<endl;
351 std::vector<int> hfScaleET = ptr1->layer1HFScaleETBins();
352 cout<<
" layer1HFScaleETBins= [";
for(
unsigned int i=0;
i<hfScaleET.size();
i++)
cout<<(
i==0?
"":
",")<<hfScaleET[
i];
cout<<
"]"<<endl;
T getUntrackedParameter(std::string const &, T const &) const
#define DEFINE_FWK_MODULE(type)
~L1TCaloParamsViewer(void) override
L1TCaloParamsViewer(const edm::ParameterSet &pset)
std::vector< std::vector< double > > tmp
std::string hash(void *buf, size_t len) const
T const * product() const
void analyze(const edm::Event &, const edm::EventSetup &) override