CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EEDataCertificationTask.cc
Go to the documentation of this file.
1 /*
2  * \file EEDataCertificationTask.cc
3  *
4  * \author E. Di Marco
5  *
6 */
7 
8 #include <iostream>
9 #include <algorithm>
10 
11 #include "FWCore/ServiceRegistry/interface/Service.h"
12 #include "FWCore/MessageLogger/interface/MessageLogger.h"
13 #include "FWCore/Framework/interface/ESHandle.h"
14 #include "FWCore/Framework/interface/EventSetup.h"
15 
16 #include "DataFormats/EcalDetId/interface/EEDetId.h"
17 
18 #include "Geometry/EcalMapping/interface/EcalMappingRcd.h"
19 #include "Geometry/EcalMapping/interface/EcalElectronicsMapping.h"
20 
21 #include "DQMServices/Core/interface/MonitorElement.h"
22 #include "DQMServices/Core/interface/DQMStore.h"
23 
24 #include "DQM/EcalCommon/interface/Numbers.h"
25 #include "DQM/EcalCommon/interface/UtilsClient.h"
26 
27 #include "DQM/EcalEndcapMonitorTasks/interface/EEDataCertificationTask.h"
28 
29 EEDataCertificationTask::EEDataCertificationTask(const edm::ParameterSet& ps) {
30 
31  // cloneME switch
32  cloneME_ = ps.getUntrackedParameter<bool>("cloneME", true);
33 
34  dqmStore_ = edm::Service<DQMStore>().operator->();
35 
36  prefixME_ = ps.getUntrackedParameter<std::string>("prefixME", "");
37 
38  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup", false);
39 
40  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns", false);
41 
42  meEEDataCertificationSummary_ = 0;
43  meEEDataCertificationSummaryMap_ = 0;
44  for (int i = 0; i < 18; i++) {
45  meEEDataCertification_[i] = 0;
46  }
47 
48  hDQM_ = 0;
49  hDAQ_ = 0;
50  hIntegrityByLumi_ = 0;
51  hFrontendByLumi_ = 0;
52  hSynchronizationByLumi_ = 0;
53 
54 }
55 
56 EEDataCertificationTask::~EEDataCertificationTask() {
57 
58 }
59 
60 void EEDataCertificationTask::beginJob(void){
61 
62  if ( dqmStore_ ) {
63 
64  std::string name;
65 
66  dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
67 
68  meEEDataCertificationSummary_ = dqmStore_->bookFloat( "CertificationSummary" );
69  meEEDataCertificationSummary_->Fill(-1.0);
70 
71  name = "CertificationSummaryMap";
72  meEEDataCertificationSummaryMap_ = dqmStore_->book2D(name, name, 40, 0., 200., 20, 0., 100.);
73  meEEDataCertificationSummaryMap_->setAxisTitle("ix / ix+100", 1);
74  meEEDataCertificationSummaryMap_->setAxisTitle("iy", 2);
75 
76  dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
77 
78  for (int i = 0; i < 18; i++) {
79  meEEDataCertification_[i] = dqmStore_->bookFloat( "EcalEndcap_" + Numbers::sEE(i+1) );
80  meEEDataCertification_[i]->Fill(-1.0);
81  }
82 
83  }
84 
85 }
86 
87 void EEDataCertificationTask::endJob(void) {
88 
89  if ( enableCleanup_ ) this->cleanup();
90 
91 }
92 
93 void EEDataCertificationTask::beginLuminosityBlock(const edm::LuminosityBlock& lumiBlock, const edm::EventSetup& iSetup){
94 
95 }
96 
97 void EEDataCertificationTask::endLuminosityBlock(const edm::LuminosityBlock& lumiBlock, const edm::EventSetup& iSetup) {
98 
99  this->reset();
100 
101  MonitorElement* me;
102 
103  // evaluate the DQM quality of observables checked by lumi
104  float DQMVal[18];
105  for (int i = 0; i < 18; i++) {
106  DQMVal[i] = -1.;
107  }
108 
109  edm::ESHandle< EcalElectronicsMapping > handle;
110  iSetup.get< EcalMappingRcd >().get(handle);
111  const EcalElectronicsMapping *map = handle.product();
112  if( ! map ) edm::LogWarning("EEDaqInfoTask") << "EcalElectronicsMapping not available";
113 
114  me = dqmStore_->get( prefixME_ + "/EEIntegrityTask/EEIT weighted integrity errors by lumi" );
115  hIntegrityByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hIntegrityByLumi_ );
116 
117  me = dqmStore_->get( prefixME_ + "/EEStatusFlagsTask/FEStatus/EESFT weighted frontend errors by lumi" );
118  hFrontendByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hFrontendByLumi_ );
119 
120  me = dqmStore_->get( prefixME_ + "/EERawDataTask/EERDT FE synchronization errors by lumi" );
121  hSynchronizationByLumi_ = UtilsClient::getHisto<TH1F*>( me, cloneME_, hSynchronizationByLumi_ );
122 
123  if( hIntegrityByLumi_ && hFrontendByLumi_ && hSynchronizationByLumi_ && map) {
124 
125  float integrityErrSum = 0.;
126  float integrityQual = 1.0;
127  float frontendErrSum = 0.;
128  float frontendQual = 1.0;
129  float synchronizationErrSum = 0.;
130  float synchronizationQual = 1.0;
131 
132  for ( int i=0; i<18; i++) {
133  float ismIntegrityQual = 1.0;
134  if( hIntegrityByLumi_->GetBinContent(0) > 0 ) {
135  float errors = hIntegrityByLumi_->GetBinContent(i+1);
136  ismIntegrityQual = 1.0 - errors/hIntegrityByLumi_->GetBinContent(0);
137  integrityErrSum += errors;
138  }
139  float ismFrontendQual = 1.0;
140  if( hFrontendByLumi_->GetBinContent(0) > 0 ) {
141  float errors = hFrontendByLumi_->GetBinContent(i+1);
142  ismFrontendQual = 1.0 - errors/hFrontendByLumi_->GetBinContent(0);
143  frontendErrSum += errors;
144  }
145  float ismSynchronizationQual = 1.0;
146  if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) {
147  float errors = hSynchronizationByLumi_->GetBinContent(i+1);
148  ismSynchronizationQual = 1.0 - errors/hSynchronizationByLumi_->GetBinContent(0);
149  synchronizationErrSum += errors;
150  }
151  float minVal= std::min(ismIntegrityQual,ismFrontendQual);
152  DQMVal[i] = std::min(minVal,ismSynchronizationQual);
153  }
154 
155  if( hIntegrityByLumi_->GetBinContent(0) > 0 ) integrityQual = 1.0 - integrityErrSum/hIntegrityByLumi_->GetBinContent(0)/18.;
156  if( hFrontendByLumi_->GetBinContent(0) > 0 ) frontendQual = 1.0 - frontendErrSum/hFrontendByLumi_->GetBinContent(0)/18.;
157  if( hSynchronizationByLumi_->GetBinContent(0) > 0 ) synchronizationQual = 1.0 - synchronizationErrSum/hSynchronizationByLumi_->GetBinContent(0)/36.;
158  float minVal = std::min(integrityQual,frontendQual);
159  float totDQMVal = std::min(minVal,synchronizationQual);
160 
161  me = dqmStore_->get( prefixME_ + "/EventInfo/reportSummary" );
162  if( me ) me->Fill(totDQMVal);
163 
164  for ( int i=0; i<18; i++) {
165  me = dqmStore_->get( prefixME_ + "/EventInfo/reportSummaryContents/EcalEndcap_" + Numbers::sEE(i+1) );
166  if( me ) me->Fill(DQMVal[i]);
167 
168  me = dqmStore_->get( prefixME_ + "/EventInfo/reportSummaryMap" );
169  if( me ) {
170  for(int t=1 ; t<=nTowerMax_ ; t++){
171  if(! map->dccTowerConstituents(DccId_[i], t).size() ) continue;
172  std::vector<EcalScDetId> scs = map->getEcalScDetId(DccId_[i], t, false);
173  for(unsigned u=0 ; u<scs.size() ; u++){
174  int jx = scs[u].ix() + (scs[u].zside()<0 ? 0 : 20);
175  int jy = scs[u].iy();
176  me->setBinContent(jx,jy, DQMVal[i]);
177  }
178  }
179  }
180  }
181 
182  }
183 
184  // now combine reduced DQM with DCS and DAQ
185  me = dqmStore_->get( prefixME_ + "/EventInfo/DAQSummaryMap" );
186  hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
187 
188  me = dqmStore_->get( prefixME_ + "/EventInfo/DCSSummaryMap" );
189  hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
190 
191  float sumCert = 0.;
192  float sumCertEE[18];
193  int nValidChannels = 0;
194  int nValidChannelsEE[18];
195 
196  if ( meEEDataCertificationSummaryMap_ ){
197  for(int ix=1 ; ix<=meEEDataCertificationSummaryMap_->getNbinsX() ; ix++){
198  for(int iy=1 ; iy<=meEEDataCertificationSummaryMap_->getNbinsY() ; iy++){
199  meEEDataCertificationSummaryMap_->setBinContent( ix, iy, -1.0 );
200  }
201  }
202  }
203 
204  for (int i = 0; i < 18; i++) {
205  sumCertEE[i] = 0;
206  nValidChannelsEE[i] = 0;
207 
208  for(int t=1 ; t<=nTowerMax_ ; t++){
209 
210  std::vector<DetId> crystals = map->dccTowerConstituents(DccId_[i], t);
211  if(!crystals.size()) continue; // getEcalScDetId throws an exception when no crystal is found
212 
213  std::vector<EcalScDetId> scs = map->getEcalScDetId(DccId_[i], t, false);
214  for(unsigned u=0 ; u<scs.size() ; u++){
215 
216  int jx = scs[u].ix() + (scs[u].zside()<0 ? 0 : 20);
217  int jy = scs[u].iy();
218 
219  float xvalDQM = DQMVal[i];
220 
221  float xvalDAQ, xvalDCS;
222  xvalDAQ = xvalDCS = -1.;
223  float xcert = -1.;
224 
225  if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( jx, jy );
226  if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( jx, jy );
227 
228  if ( xvalDQM == -1 || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
229  // problems: DQM empty or DAQ and DCS not available
230  xcert = 0.0;
231  } else {
232  // do not consider the white value of DAQ and DCS (problems with DB)
233  xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
234  }
235 
236  if ( meEEDataCertificationSummaryMap_ ) meEEDataCertificationSummaryMap_->setBinContent( jx, jy, xcert );
237 
238  int ncrystals = 0;
239 
240  for(std::vector<DetId>::const_iterator it=crystals.begin() ; it!=crystals.end() ; ++it){
241  EEDetId id(*it);
242  if( id.zside() == scs[u].zside() && (id.ix()-1)/5+1 == scs[u].ix() && (id.iy()-1)/5+1 == scs[u].iy() ) ncrystals++;
243  }
244 
245  sumCertEE[i] += xcert * ncrystals;
246  nValidChannelsEE[i] += ncrystals;
247 
248  sumCert += xcert * ncrystals;
249  nValidChannels += ncrystals;
250 
251  }
252  }
253 
254  if( meEEDataCertification_[i] ) {
255  if( nValidChannelsEE[i]>0 ) meEEDataCertification_[i]->Fill( sumCertEE[i]/nValidChannelsEE[i] );
256  else meEEDataCertification_[i]->Fill( 0.0 );
257  }
258  }
259 
260  if( meEEDataCertificationSummary_ ) {
261  if( nValidChannels>0 ) meEEDataCertificationSummary_->Fill( sumCert/nValidChannels );
262  else meEEDataCertificationSummary_->Fill( 0.0 );
263  }
264 
265 }
266 
267 void EEDataCertificationTask::beginRun(const edm::Run& r, const edm::EventSetup& c) {
268 
269  if ( ! mergeRuns_ ) this->reset();
270 
271 }
272 
273 void EEDataCertificationTask::endRun(const edm::Run& r, const edm::EventSetup& c) {
274 
275  this->reset();
276 
277  edm::ESHandle< EcalElectronicsMapping > handle;
278  c.get< EcalMappingRcd >().get(handle);
279  const EcalElectronicsMapping *map = handle.product();
280  if( ! map ){
281  edm::LogWarning("EEDaqInfoTask") << "EcalElectronicsMapping not available";
282  return;
283  }
284 
285  MonitorElement* me;
286 
287  me = dqmStore_->get( prefixME_ + "/EventInfo/reportSummaryMap" );
288  hDQM_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDQM_ );
289 
290  me = dqmStore_->get( prefixME_ + "/EventInfo/DAQSummaryMap" );
291  hDAQ_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDAQ_ );
292 
293  me = dqmStore_->get( prefixME_ + "/EventInfo/DCSSummaryMap" );
294  hDCS_ = UtilsClient::getHisto<TH2F*>( me, cloneME_, hDCS_ );
295 
296  float sumCert = 0.;
297  float sumCertEE[18];
298  int nValidChannels = 0;
299  int nValidChannelsEE[18];
300 
301  if ( meEEDataCertificationSummaryMap_ ){
302  for(int ix=1 ; ix<=meEEDataCertificationSummaryMap_->getNbinsX() ; ix++){
303  for(int iy=1 ; iy<=meEEDataCertificationSummaryMap_->getNbinsY() ; iy++){
304  meEEDataCertificationSummaryMap_->setBinContent( ix, iy, -1.0 );
305  }
306  }
307  }
308 
309  for (int i = 0; i < 18; i++) {
310  sumCertEE[i] = 0;
311  nValidChannelsEE[i] = 0;
312 
313  for(int t=1 ; t<=nTowerMax_ ; t++){
314 
315  std::vector<DetId> crystals = map->dccTowerConstituents(DccId_[i], t);
316  if(!crystals.size()) continue;
317 
318  std::vector<EcalScDetId> scs = map->getEcalScDetId(DccId_[i], t, false);
319  for(unsigned u=0 ; u<scs.size() ; u++){
320 
321  int jx = scs[u].ix() + (scs[u].zside()<0 ? 0 : 20);
322  int jy = scs[u].iy();
323 
324  float xvalDQM, xvalDAQ, xvalDCS;
325  xvalDQM = xvalDAQ = xvalDCS = -1.;
326  float xcert = -1.;
327 
328  if ( hDQM_ ) xvalDQM = hDQM_->GetBinContent( jx, jy );
329  if ( hDAQ_ ) xvalDAQ = hDAQ_->GetBinContent( jx, jy );
330  if ( hDCS_ ) xvalDCS = hDCS_->GetBinContent( jx, jy );
331 
332  if ( xvalDQM == -1 || ( xvalDAQ == -1 && xvalDCS == -1 ) ) {
333  // problems: DQM empty or DAQ and DCS not available
334  xcert = 0.0;
335  } else {
336  // do not consider the white value of DAQ and DCS (problems with DB)
337  xcert = std::abs(xvalDQM) * std::abs(xvalDAQ) * std::abs(xvalDCS);
338  }
339 
340  if ( meEEDataCertificationSummaryMap_ ) meEEDataCertificationSummaryMap_->setBinContent( jx, jy, xcert );
341 
342  int ncrystals = 0;
343 
344  for(std::vector<DetId>::const_iterator it=crystals.begin() ; it!=crystals.end() ; ++it){
345  EEDetId id(*it);
346  if( id.zside() == scs[u].zside() && (id.ix()-1)/5+1 == scs[u].ix() && (id.iy()-1)/5+1 == scs[u].iy() ) ncrystals++;
347  }
348 
349  sumCertEE[i] += xcert * ncrystals;
350  nValidChannelsEE[i] += ncrystals;
351 
352  sumCert += xcert * ncrystals;
353  nValidChannels += ncrystals;
354 
355  }
356  }
357 
358  if( meEEDataCertification_[i] ) {
359  if( nValidChannelsEE[i]>0 ) meEEDataCertification_[i]->Fill( sumCertEE[i]/nValidChannelsEE[i] );
360  else meEEDataCertification_[i]->Fill( 0.0 );
361  }
362  }
363 
364  if( meEEDataCertificationSummary_ ) {
365  if( nValidChannels>0 ) meEEDataCertificationSummary_->Fill( sumCert/nValidChannels );
366  else meEEDataCertificationSummary_->Fill( 0.0 );
367  }
368 
369 }
370 
371 void EEDataCertificationTask::reset(void) {
372 
373  if ( meEEDataCertificationSummary_ ) meEEDataCertificationSummary_->Reset();
374 
375  for (int i = 0; i < 18; i++) {
376  if ( meEEDataCertification_[i] ) meEEDataCertification_[i]->Reset();
377  }
378 
379  if ( meEEDataCertificationSummaryMap_ ) meEEDataCertificationSummaryMap_->Reset();
380 
381  if ( meEEDataCertificationSummaryMap_ ) meEEDataCertificationSummaryMap_->Reset();
382 
383 }
384 
385 
386 void EEDataCertificationTask::cleanup(void){
387 
388  if ( cloneME_ ) {
389  if( hDQM_ ) delete hDQM_;
390  if( hDAQ_ ) delete hDAQ_;
391  if( hDCS_ ) delete hDCS_;
392  if( hIntegrityByLumi_ ) delete hIntegrityByLumi_;
393  if( hFrontendByLumi_ ) delete hFrontendByLumi_;
394  if( hSynchronizationByLumi_ ) delete hSynchronizationByLumi_;
395  }
396  hDQM_ = 0;
397  hDAQ_ = 0;
398  hDCS_ = 0;
399  hIntegrityByLumi_ = 0;
400  hFrontendByLumi_ = 0;
401  hSynchronizationByLumi_ = 0;
402 
403  if ( dqmStore_ ) {
404  dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo");
405  if ( meEEDataCertificationSummary_ ) dqmStore_->removeElement( meEEDataCertificationSummary_->getName() );
406  if ( meEEDataCertificationSummaryMap_ ) dqmStore_->removeElement( meEEDataCertificationSummaryMap_->getName() );
407 
408  dqmStore_->setCurrentFolder(prefixME_ + "/EventInfo/CertificationContents");
409  for (int i = 0; i < 18; i++) {
410  if ( meEEDataCertification_[i] ) dqmStore_->removeElement( meEEDataCertification_[i]->getName() );
411  }
412  }
413 
414 }
415 
416 void EEDataCertificationTask::analyze(const edm::Event& e, const edm::EventSetup& c){
417 
418 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:108
Numbers::sEE
static std::string sEE(const unsigned ism)
Definition: Numbers.cc:223
i
int i
Definition: DBlmapReader.cc:9
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:825
benchmark_cfg.errors
tuple errors
Definition: benchmark_cfg.py:136
EEDataCertificationTask::meEEDataCertification_
MonitorElement * meEEDataCertification_[18]
Definition: EEDataCertificationTask.h:79
Numbers.h
Some &quot;id&quot; conversions.
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DQMStore::bookFloat
MonitorElement * bookFloat(const char *name)
Book float.
Definition: DQMStore.cc:809
EEDataCertificationTask::~EEDataCertificationTask
virtual ~EEDataCertificationTask()
Destructor.
Definition: EEDataCertificationTask.cc:56
MonitorElement::getNbinsY
int getNbinsY(void) const
get # of bins in Y-axis
Definition: MonitorElement.cc:712
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
EEDataCertificationTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: EEDataCertificationTask.cc:416
EEDataCertificationTask::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: EEDataCertificationTask.cc:267
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:2772
EEDataCertificationTask::meEEDataCertificationSummaryMap_
MonitorElement * meEEDataCertificationSummaryMap_
Definition: EEDataCertificationTask.h:80
patZpeak.handle
tuple handle
Definition: patZpeak.py:22
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
EEDataCertificationTask::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
EndLuminosityBlock.
Definition: EEDataCertificationTask.cc:97
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1623
EcalElectronicsMapping::getEcalScDetId
std::vector< EcalScDetId > getEcalScDetId(int DCCid, int DCC_Channel, bool ignoreSingleCrystal=true) const
Definition: EcalElectronicsMapping.cc:840
UtilsClient.h
Ecal Monitor Utils for Client.
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:62
EEDataCertificationTask::EEDataCertificationTask
EEDataCertificationTask(const edm::ParameterSet &ps)
Constructor.
Definition: EEDataCertificationTask.cc:29
EEDataCertificationTask::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &lumiBlock, const edm::EventSetup &iSetup)
BeginLuminosityBlock.
Definition: EEDataCertificationTask.cc:93
MonitorElement::getNbinsX
int getNbinsX(void) const
get # of bins in X-axis
Definition: MonitorElement.cc:706
EEDataCertificationTask::endRun
void endRun(const edm::Run &r, const edm::EventSetup &c)
EndRun.
Definition: EEDataCertificationTask.cc:273
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1000
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:928
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:434
EEDataCertificationTask::meEEDataCertificationSummary_
MonitorElement * meEEDataCertificationSummary_
Definition: EEDataCertificationTask.h:78
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
edm::Run
Definition: Run.h:41
EcalElectronicsMapping::dccTowerConstituents
std::vector< DetId > dccTowerConstituents(int dccId, int tower) const
Get the constituent detids for this dccId.
Definition: EcalElectronicsMapping.cc:611
Generated for CMSSW Reference Manual by doxygen 1.8.5