CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
CommissioningHistograms.cc
Go to the documentation of this file.
1 #include "DQM/SiStripCommissioningClients/interface/CommissioningHistograms.h"
2 #include "CondFormats/SiStripObjects/interface/CommissioningAnalysis.h"
3 #include "DataFormats/SiStripCommon/interface/SiStripFecKey.h"
4 #include "DataFormats/SiStripCommon/interface/SiStripFedKey.h"
5 #include "DataFormats/SiStripCommon/interface/SiStripDetKey.h"
6 #include "DQM/SiStripCommissioningSummary/interface/SummaryGenerator.h"
7 #include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
8 #include "FWCore/MessageLogger/interface/MessageLogger.h"
9 #include <iomanip>
10 
11 using namespace std;
12 using namespace sistrip;
13 
14 
15 // -----------------------------------------------------------------------------
17 CommissioningHistograms::CommissioningHistograms( const edm::ParameterSet& pset,
18  DQMStore* bei,
19  const sistrip::RunType& task )
20  : factory_(0),
21  task_(task),
22  bei_(bei),
23  data_(),
24  histos_(),
25  pset_(pset)
26 {
27  LogTrace(mlDqmClient_)
28  << "[" << __PRETTY_FUNCTION__ << "]"
29  << " Constructing object...";
30 
31  // DQMStore
32  if ( !bei_ ) {
33  edm::LogError(mlDqmClient_)
34  << "[CommissioningHistograms::" << __func__ << "]"
35  << " NULL pointer to DQMStore!";
36  }
37 
38  clearHistosMap();
39 }
40 
41 // -----------------------------------------------------------------------------
43 CommissioningHistograms::CommissioningHistograms()
44  : factory_(0),
45  task_(sistrip::UNDEFINED_RUN_TYPE),
46  bei_(0),
47  data_(),
48  histos_()
49 {
50  LogTrace(mlDqmClient_)
51  << "[" << __PRETTY_FUNCTION__ << "]"
52  << " Constructing object...";
53 }
54 
55 // -----------------------------------------------------------------------------
57 CommissioningHistograms::~CommissioningHistograms() {
58  LogTrace(mlDqmClient_)
59  << "[" << __PRETTY_FUNCTION__ << "]"
60  << " Destructing object...";
61  clearHistosMap();
62  //@@ do not delete BEI ptrs!
63 }
64 
65 // -----------------------------------------------------------------------------
67 void CommissioningHistograms::Histo::print( std::stringstream& ss ) const {
68  ss << " [Histo::" << __func__ << "]" << std::endl
69  << " Histogram title : " << title_ << std::endl
70  << " MonitorElement* : 0x"
71  << std::hex
72  << std::setw(8) << std::setfill('0') << me_ << std::endl
73  << std::dec
74  << " CollateME* : 0x"
75  << std::hex
76  << std::setw(8) << std::setfill('0') << cme_ << std::endl
77  << std::dec;
78 }
79 
80 // -----------------------------------------------------------------------------
81 //
82 uint32_t CommissioningHistograms::runNumber( DQMStore* const bei,
83  const std::vector<std::string>& contents ) {
84 
85  // Check if histograms present
86  if ( contents.empty() ) {
87  edm::LogError(mlDqmClient_)
88  << "[CommissioningHistograms::" << __func__ << "]"
89  << " Found no histograms!";
90  return 0;
91  }
92 
93  // Iterate through added contents
94  std::vector<std::string>::const_iterator istr = contents.begin();
95  while ( istr != contents.end() ) {
96 
97  // Extract source directory path
98  std::string source_dir = istr->substr( 0, istr->find(":") );
99 
100  // Generate corresponding client path (removing trailing "/")
101  SiStripFecKey path( source_dir );
102  std::string client_dir = path.path();
103  std::string slash = client_dir.substr( client_dir.size()-1, 1 );
104  if ( slash == sistrip::dir_ ) { client_dir = client_dir.substr( 0, client_dir.size()-1 ); }
105  client_dir = std::string(sistrip::collate_) + sistrip::dir_ + client_dir;
106 
107  // Iterate though MonitorElements from source directory
108  std::vector<MonitorElement*> me_list = bei->getContents( source_dir );
109  std::vector<MonitorElement*>::iterator ime = me_list.begin();
110  for ( ; ime != me_list.end(); ime++ ) {
111 
112  if ( !(*ime) ) {
113  edm::LogError(mlDqmClient_)
114  << "[CommissioningHistograms::" << __func__ << "]"
115  << " NULL pointer to MonitorElement!";
116  continue;
117  }
118 
119  // Search for run type in string
120  std::string title = (*ime)->getName();
121  std::string::size_type pos = title.find( sistrip::runNumber_ );
122 
123  // Extract run number from string
124  if ( pos != std::string::npos ) {
125  std::string value = title.substr( pos+sizeof(sistrip::runNumber_) , std::string::npos );
126  if ( !value.empty() ) {
127  LogTrace(mlDqmClient_)
128  << "[CommissioningHistograms::" << __func__ << "]"
129  << " Found string \"" << title.substr(pos,std::string::npos)
130  << "\" with value \"" << value << "\"";
131  if ( !(bei->get(client_dir+"/"+title.substr(pos,std::string::npos))) ) {
132  bei->setCurrentFolder(client_dir);
133  bei->bookString( title.substr(pos,std::string::npos), value );
134  LogTrace(mlDqmClient_)
135  << "[CommissioningHistograms::" << __func__ << "]"
136  << " Booked string \"" << title.substr(pos,std::string::npos)
137  << "\" in directory \"" << client_dir << "\"";
138  }
139  uint32_t run;
140  std::stringstream ss;
141  ss << value;
142  ss >> std::dec >> run;
143  return run;
144  }
145  }
146 
147  }
148 
149  istr++;
150 
151  }
152  return 0;
153 }
154 
155 // -----------------------------------------------------------------------------
157 sistrip::RunType CommissioningHistograms::runType( DQMStore* const bei,
158  const std::vector<std::string>& contents ) {
159 
160  // Check if histograms present
161  if ( contents.empty() ) {
162  edm::LogError(mlDqmClient_)
163  << "[CommissioningHistograms::" << __func__ << "]"
164  << " Found no histograms!";
165  return sistrip::UNKNOWN_RUN_TYPE;
166  }
167 
168  // Iterate through added contents
169  std::vector<std::string>::const_iterator istr = contents.begin();
170  while ( istr != contents.end() ) {
171 
172  // Extract source directory path
173  std::string source_dir = istr->substr( 0, istr->find(":") );
174 
175  // Generate corresponding client path (removing trailing "/")
176  SiStripFecKey path( source_dir );
177  std::string client_dir = path.path();
178  std::string slash = client_dir.substr( client_dir.size()-1, 1 );
179  if ( slash == sistrip::dir_ ) { client_dir = client_dir.substr( 0, client_dir.size()-1 ); }
180  client_dir = std::string(sistrip::collate_) + sistrip::dir_ + client_dir;
181 
182  // Iterate though MonitorElements from source directory
183  std::vector<MonitorElement*> me_list = bei->getContents( source_dir );
184 
185  if ( me_list.empty() ) {
186  edm::LogError(mlDqmClient_)
187  << "[CommissioningHistograms::" << __func__ << "]"
188  << " No MonitorElements found in dir " << source_dir;
189  return sistrip::UNKNOWN_RUN_TYPE;
190  }
191 
192  std::vector<MonitorElement*>::iterator ime = me_list.begin();
193  for ( ; ime != me_list.end(); ime++ ) {
194 
195  if ( !(*ime) ) {
196  edm::LogError(mlDqmClient_)
197  << "[CommissioningHistograms::" << __func__ << "]"
198  << " NULL pointer to MonitorElement!";
199  continue;
200  }
201 
202  // Search for run type in string
203  std::string title = (*ime)->getName();
204  std::string::size_type pos = title.find( sistrip::taskId_ );
205 
206  // Extract commissioning task from string
207  if ( pos != std::string::npos ) {
208  std::string value = title.substr( pos+sizeof(sistrip::taskId_) , std::string::npos );
209  if ( !value.empty() ) {
210  LogTrace(mlDqmClient_)
211  << "[CommissioningHistograms::" << __func__ << "]"
212  << " Found string \"" << title.substr(pos,std::string::npos)
213  << "\" with value \"" << value << "\"";
214  if ( !(bei->get(client_dir+sistrip::dir_+title.substr(pos,std::string::npos))) ) {
215  bei->setCurrentFolder(client_dir);
216  bei->bookString( title.substr(pos,std::string::npos), value );
217  LogTrace(mlDqmClient_)
218  << "[CommissioningHistograms::" << __func__ << "]"
219  << " Booked string \"" << title.substr(pos,std::string::npos)
220  << "\" in directory \"" << client_dir << "\"";
221  }
222  return SiStripEnumsAndStrings::runType( value );
223  }
224  }
225 
226  }
227 
228  istr++;
229 
230  }
231 
232  edm::LogError(mlDqmClient_)
233  << "[CommissioningHistograms::" << __func__ << "]"
234  << " Unable to extract RunType!";
235  return sistrip::UNKNOWN_RUN_TYPE;
236 
237 }
238 
239 // -----------------------------------------------------------------------------
240 // Temporary fix: builds a list of histogram directories
241 void CommissioningHistograms::getContents( DQMStore* const bei,
242  std::vector<std::string>& contents ) {
243 }
244 
245 // -----------------------------------------------------------------------------
246 //
247 void CommissioningHistograms::copyCustomInformation( DQMStore* const bei,
248  const std::vector<std::string>& contents ) {
249 
250  // Check if histograms present
251  if ( contents.empty() ) {
252  edm::LogWarning(mlDqmClient_)
253  << "[CommissioningHistograms::" << __func__ << "]"
254  << " Found no histograms!";
255  return;
256  }
257 
258  // Iterate through added contents
259  std::vector<std::string>::const_iterator istr = contents.begin();
260  while ( istr != contents.end() ) {
261 
262  // Extract source directory path
263  std::string source_dir = istr->substr( 0, istr->find(":") );
264 
265  // Generate corresponding client path (removing trailing "/")
266  SiStripFecKey path( source_dir );
267  std::string client_dir = path.path();
268  std::string slash = client_dir.substr( client_dir.size()-1, 1 );
269  if ( slash == sistrip::dir_ ) { client_dir = client_dir.substr( 0, client_dir.size()-1 ); }
270 
271  // Iterate though MonitorElements from source directory
272  std::vector<MonitorElement*> me_list = bei->getContents( source_dir );
273  std::vector<MonitorElement*>::iterator ime = me_list.begin();
274  for ( ; ime != me_list.end(); ime++ ) {
275 
276  if ( !(*ime) ) {
277  edm::LogWarning(mlDqmClient_)
278  << "[CommissioningHistograms::" << __func__ << "]"
279  << " NULL pointer to MonitorElement!";
280  continue;
281  }
282  // Search for calchan, isha or vfs
283  if((*ime)->kind()==MonitorElement::DQM_KIND_INT) {
284  std::string title = (*ime)->getName();
285  std::string::size_type pos = title.find("calchan");
286  if( pos == std::string::npos ) pos = title.find("isha");
287  if( pos == std::string::npos ) pos = title.find("vfs");
288  if( pos != std::string::npos ) {
289  int value = (*ime)->getIntValue();
290  if ( value>=0 ) {
291  edm::LogVerbatim(mlDqmClient_)
292  << "[CommissioningHistograms::" << __func__ << "]"
293  << " Found \"" << title.substr(pos,std::string::npos)
294  << "\" with value \"" << value << "\"";
295  if ( !(bei->get(client_dir+"/"+title.substr(pos,std::string::npos))) ) {
296  bei->setCurrentFolder(client_dir);
297  bei->bookInt( title.substr(pos,std::string::npos))->Fill(value);
298  edm::LogVerbatim(mlDqmClient_)
299  << "[CommissioningHistograms::" << __func__ << "]"
300  << " Booked \"" << title.substr(pos,std::string::npos)
301  << "\" in directory \"" << client_dir << "\"";
302  }
303  }
304  }
305  }
306  }
307  istr++;
308  }
309 }
310 
311 // -----------------------------------------------------------------------------
312 
314 void CommissioningHistograms::extractHistograms( const std::vector<std::string>& contents ) {
315  LogTrace(mlDqmClient_)
316  << "[CommissioningHistograms::" << __func__ << "]"
317  << " Extracting available histograms...";
318 
319  // Check pointer
320  if ( !bei_ ) {
321  edm::LogError(mlDqmClient_)
322  << "[CommissioningHistograms::" << __func__ << "]"
323  << " NULL pointer to DQMStore!";
324  return;
325  }
326 
327  // Check list of histograms
328  if ( contents.empty() ) {
329  edm::LogError(mlDqmClient_)
330  << "[CommissioningHistograms::" << __func__ << "]"
331  << " Empty contents vector!";
332  return;
333  }
334 
335  // Iterate through list of histograms
336  std::vector<std::string>::const_iterator idir;
337  for ( idir = contents.begin(); idir != contents.end(); idir++ ) {
338 
339  // Ignore "DQM source" directories if looking in client file
340  if ( idir->find(sistrip::collate_) == std::string::npos ) { continue; }
341 
342  // Extract source directory path
343  std::string source_dir = idir->substr( 0, idir->find(":") );
344 
345  // Extract view and create key
346  sistrip::View view = SiStripEnumsAndStrings::view( source_dir );
347  SiStripKey path;
348  if ( view == sistrip::CONTROL_VIEW ) { path = SiStripFecKey( source_dir ); }
349  else if ( view == sistrip::READOUT_VIEW ) { path = SiStripFedKey( source_dir ); }
350  else if ( view == sistrip::DETECTOR_VIEW ) { path = SiStripDetKey( source_dir ); }
351  else { path = SiStripKey(); }
352 
353  // Check path is valid
354  if ( path.granularity() == sistrip::UNKNOWN_GRAN ||
355  path.granularity() == sistrip::UNDEFINED_GRAN ) {
356  continue;
357  }
358 
359  // Generate corresponding client path (removing trailing "/")
360  std::string client_dir(sistrip::undefinedView_);
361  if ( view == sistrip::CONTROL_VIEW ) { client_dir = SiStripFecKey( path.key() ).path(); }
362  else if ( view == sistrip::READOUT_VIEW ) { client_dir = SiStripFedKey( path.key() ).path(); }
363  else if ( view == sistrip::DETECTOR_VIEW ) { client_dir = SiStripDetKey( path.key() ).path(); }
364  else { client_dir = SiStripKey( path.key() ).path(); }
365  std::string slash = client_dir.substr( client_dir.size()-1, 1 );
366  if ( slash == sistrip::dir_ ) { client_dir = client_dir.substr( 0, client_dir.size()-1 ); }
367  client_dir = std::string(sistrip::collate_) + sistrip::dir_ + client_dir;
368 
369  // Retrieve MonitorElements from source directory
370  std::vector<MonitorElement*> me_list = bei_->getContents( source_dir );
371 
372  // Iterate though MonitorElements and create CMEs
373  std::vector<MonitorElement*>::iterator ime = me_list.begin();
374  for ( ; ime != me_list.end(); ime++ ) {
375 
376  // Retrieve histogram title
377  SiStripHistoTitle title( (*ime)->getName() );
378 
379  // Check histogram type
380  //if ( title.histoType() != sistrip::EXPERT_HISTO ) { continue; }
381 
382  // Check granularity
383  uint16_t channel = sistrip::invalid_;
384  if ( title.granularity() == sistrip::APV ) {
385  channel = SiStripFecKey::lldChan( title.channel() );
386  } else if ( title.granularity() == sistrip::UNKNOWN_GRAN ||
387  title.granularity() == sistrip::UNDEFINED_GRAN ) {
388  std::stringstream ss;
389  ss << "[CommissioningHistograms::" << __func__ << "]"
390  << " Unexpected granularity for histogram title: "
391  << std::endl << title
392  << " found in path "
393  << std::endl << path;
394  edm::LogError(mlDqmClient_) << ss.str();
395  } else {
396  channel = title.channel();
397  }
398 
399  // Build key
400  uint32_t key = sistrip::invalid32_;
401 
402  if ( view == sistrip::CONTROL_VIEW ) {
403 
404  // for all runs except cabling
405  SiStripFecKey temp( path.key() );
406  key = SiStripFecKey( temp.fecCrate(),
407  temp.fecSlot(),
408  temp.fecRing(),
409  temp.ccuAddr(),
410  temp.ccuChan(),
411  channel ).key();
412  mapping_[title.keyValue()] = key;
413 
414  } else if ( view == sistrip::READOUT_VIEW ) {
415 
416  // for cabling run
417  key = SiStripFedKey( path.key() ).key();
418  uint32_t temp = SiStripFecKey( sistrip::invalid_,
419  sistrip::invalid_,
420  sistrip::invalid_,
421  sistrip::invalid_,
422  sistrip::invalid_,
423  channel ).key(); // just record lld channel
424  mapping_[title.keyValue()] = temp;
425 
426  } else if ( view == sistrip::DETECTOR_VIEW ) {
427 
428  SiStripDetKey temp( path.key() );
429  key = SiStripDetKey( temp.partition() ).key();
430  mapping_[title.keyValue()] = key;
431 
432  } else { key = SiStripKey( path.key() ).key(); }
433 
434  // Find CME in histos map
435  Histo* histo = 0;
436  HistosMap::iterator ihistos = histos_.find( key );
437  if ( ihistos != histos_.end() ) {
438  Histos::iterator ihis = ihistos->second.begin();
439  while ( !histo && ihis < ihistos->second.end() ) {
440  if ( (*ime)->getName() == (*ihis)->title_ ) { histo = *ihis; }
441  ihis++;
442  }
443  }
444 
445  // Create CollateME if it doesn't exist
446  if ( !histo ) {
447 
448  histos_[key].push_back( new Histo() );
449  histo = histos_[key].back();
450  histo->title_ = (*ime)->getName();
451 
452  // If histogram present in client directory, add to map
453  if ( source_dir.find(sistrip::collate_) != std::string::npos ) {
454  histo->me_ = bei_->get( client_dir +"/"+(*ime)->getName() );
455  if ( !histo->me_ ) {
456  edm::LogError(mlDqmClient_)
457  << "[CommissioningHistograms::" << __func__ << "]"
458  << " NULL pointer to MonitorElement!";
459  }
460  }
461 
462  }
463 
464  }
465 
466  }
467 
468  //printHistosMap();
469 
470  edm::LogVerbatim(mlDqmClient_)
471  << "[CommissioningHistograms::" << __func__ << "]"
472  << " Found histograms for " << histos_.size()
473  << " structures in cached histogram map!";
474 
475 }
476 
477 
478 // -----------------------------------------------------------------------------
480 void CommissioningHistograms::histoAnalysis( bool debug ) {
481  LogTrace(mlDqmClient_)
482  << "[CommissioningHistograms::" << __func__ << "]"
483  << " (Derived) implementation to come...";
484 }
485 
486 // -----------------------------------------------------------------------------
488 void CommissioningHistograms::printAnalyses() {
489  Analyses::iterator ianal = data().begin();
490  Analyses::iterator janal = data().end();
491  for ( ; ianal != janal; ++ianal ) {
492  if ( ianal->second ) {
493  std::stringstream ss;
494  ianal->second->print( ss );
495  if ( ianal->second->isValid() ) { LogTrace(mlDqmClient_) << ss.str();
496  } else { edm::LogWarning(mlDqmClient_) << ss.str(); }
497  }
498  }
499 }
500 
501 // -----------------------------------------------------------------------------
503 void CommissioningHistograms::printSummary() {
504 
505  std::stringstream good;
506  std::stringstream bad;
507 
508  Analyses::iterator ianal = data().begin();
509  Analyses::iterator janal = data().end();
510  for ( ; ianal != janal; ++ianal ) {
511  if ( ianal->second ) {
512  if ( ianal->second->isValid() ) { ianal->second->summary( good ); }
513  else { ianal->second->summary( bad ); }
514  }
515  }
516 
517  if ( good.str().empty() ) { good << "None found!"; }
518  LogTrace(mlDqmClient_)
519  << "[CommissioningHistograms::" << __func__ << "]"
520  << " Printing summary of good analyses:" << "\n"
521  << good.str();
522 
523  if ( bad.str().empty() ) { return; } //@@ bad << "None found!"; }
524  LogTrace(mlDqmClient_)
525  << "[CommissioningHistograms::" << __func__ << "]"
526  << " Printing summary of bad analyses:" << "\n"
527  << bad.str();
528 
529 }
530 
531 // -----------------------------------------------------------------------------
533 void CommissioningHistograms::printHistosMap() {
534  LogTrace(mlDqmClient_)
535  << "[CommissioningHistograms::" << __func__ << "]"
536  << " Printing histogram map, which has "
537  << histos_.size() << " entries...";
538  HistosMap::const_iterator ihistos = histos_.begin();
539  for ( ; ihistos != histos_.end(); ihistos++ ) {
540  std::stringstream ss;
541  ss << " Found " << ihistos->second.size()
542  << " histogram(s) for key: " << std::endl
543  << SiStripFedKey(ihistos->first) << std::endl;
544  Histos::const_iterator ihisto = ihistos->second.begin();
545  for ( ; ihisto != ihistos->second.end(); ihisto++ ) {
546  if ( *ihisto ) { (*ihisto)->print(ss); }
547  else { ss << " NULL pointer to Histo object!"; }
548  }
549  LogTrace(mlDqmClient_) << ss.str();
550  }
551 }
552 
553 // -----------------------------------------------------------------------------
555 void CommissioningHistograms::clearHistosMap() {
556  LogTrace(mlDqmClient_)
557  << "[CommissioningHistograms::" << __func__ << "]"
558  << " Clearing histogram map...";
559  HistosMap::iterator ihistos = histos_.begin();
560  for ( ; ihistos != histos_.end(); ihistos++ ) {
561  Histos::iterator ihisto = ihistos->second.begin();
562  for ( ; ihisto != ihistos->second.end(); ihisto++ ) {
563  if ( *ihisto ) { delete *ihisto; }
564  }
565  ihistos->second.clear();
566  }
567  histos_.clear();
568 }
569 
570 // -----------------------------------------------------------------------------
572 void CommissioningHistograms::createSummaryHisto( const sistrip::Monitorable& mon,
573  const sistrip::Presentation& pres,
574  const std::string& dir,
575  const sistrip::Granularity& gran ) {
576  LogTrace(mlDqmClient_)
577  << "[CommissioningHistograms::" << __func__ << "]";
578 
579  // Check view
580  sistrip::View view = SiStripEnumsAndStrings::view(dir);
581  if ( view == sistrip::UNKNOWN_VIEW ) { return; }
582 
583  // Analyze histograms
584  if ( data().empty() ) { histoAnalysis( false ); }
585 
586  // Check
587  if ( data().empty() ) {
588  edm::LogError(mlDqmClient_)
589  << "[CommissioningHistograms::" << __func__ << "]"
590  << " No analyses generated!";
591  return;
592  }
593 
594  // Extract data to be histogrammed
595  uint32_t xbins = factory()->init( mon, pres, view, dir, gran, data() );
596 
597  // Only create histograms if entries are found!
598  if ( !xbins ) { return; }
599 
600  // Create summary histogram (if it doesn't already exist)
601  TH1* summary = 0;
602  if ( pres != sistrip::HISTO_1D ) { summary = histogram( mon, pres, view, dir, xbins ); }
603  else { summary = histogram( mon, pres, view, dir, sistrip::FED_ADC_RANGE, 0., sistrip::FED_ADC_RANGE*1. ); }
604 
605  // Fill histogram with data
606  factory()->fill( *summary );
607 
608 }
609 
610 // -----------------------------------------------------------------------------
612 void CommissioningHistograms::remove( std::string pattern ) {
613 
614  if ( !bei_ ) {
615  edm::LogError(mlDqmClient_)
616  << "[CommissioningHistograms::" << __func__ << "]"
617  << " NULL pointer to DQMStore!";
618  return;
619  }
620 
621  bei_->setVerbose(0);
622 
623  LogTrace(mlDqmClient_)
624  << "[CommissioningHistograms::" << __func__ << "]"
625  << " Removing histograms...";
626 
627  if ( !pattern.empty() ) {
628 
629  if ( bei_->dirExists(pattern) ) {
630  bei_->rmdir(pattern);
631  }
632 
633  LogTrace(mlDqmClient_)
634  << "[CommissioningHistograms::" << __func__ << "]"
635  << " Removing directories (and MonitorElements"
636  << " therein) that match the pattern \""
637  << pattern << "\"";
638 
639  } else {
640 
641  bei_->cd();
642  bei_->removeContents();
643 
644  if( bei_->dirExists("Collector") ) {
645  bei_->rmdir("Collector");
646  }
647  if( bei_->dirExists("EvF") ) {
648  bei_->rmdir("EvF");
649  }
650  if( bei_->dirExists("SiStrip") ) {
651  bei_->rmdir("SiStrip");
652  }
653 
654  LogTrace(mlDqmClient_)
655  << "[CommissioningHistograms::" << __func__ << "]"
656  << " Removing \"DQM source\" directories (and MonitorElements therein)";
657 
658  }
659 
660  LogTrace(mlDqmClient_)
661  << "[CommissioningHistograms::" << __func__ << "]"
662  << " Removed histograms!";
663 
664  bei_->setVerbose(1);
665 
666 }
667 
668 // -----------------------------------------------------------------------------
670 void CommissioningHistograms::save( std::string& path,
671  uint32_t run_number ) {
672 
673  // Construct path and filename
674  std::stringstream ss;
675 
676  if ( !path.empty() ) {
677 
678  ss << path;
679  if ( ss.str().find(".root") == std::string::npos ) { ss << ".root"; }
680 
681  } else {
682 
683  // Retrieve SCRATCH directory
684  std::string scratch = "SCRATCH";
685  std::string dir = "";
686  if ( getenv(scratch.c_str()) != NULL ) {
687  dir = getenv(scratch.c_str());
688  }
689 
690  // Add directory path
691  if ( !dir.empty() ) { ss << dir << "/"; }
692  else { ss << "/tmp/"; }
693 
694  // Add filename with run number and ".root" extension
695  ss << sistrip::dqmClientFileName_ << "_"
696  << std::setfill('0') << std::setw(8) << run_number
697  << ".root";
698 
699  }
700 
701  // Save file with appropriate filename
702  LogTrace(mlDqmClient_)
703  << "[CommissioningHistograms::" << __func__ << "]"
704  << " Saving histograms to root file"
705  << " (This may take some time!)";
706  path = ss.str();
707  bei_->save( path, sistrip::collate_ );
708  edm::LogVerbatim(mlDqmClient_)
709  << "[CommissioningHistograms::" << __func__ << "]"
710  << " Saved histograms to root file \""
711  << ss.str() << "\"!";
712 
713 }
714 
715 // -----------------------------------------------------------------------------
716 //
717 TH1* CommissioningHistograms::histogram( const sistrip::Monitorable& mon,
718  const sistrip::Presentation& pres,
719  const sistrip::View& view,
720  const std::string& directory,
721  const uint32_t& xbins,
722  const float& xlow,
723  const float& xhigh ) {
724 
725  // Remember pwd
726  std::string pwd = bei_->pwd();
727  bei_->setCurrentFolder( std::string(sistrip::collate_) + sistrip::dir_ + directory );
728 
729  // Construct histogram name
730  std::string name = SummaryGenerator::name( task_, mon, pres, view, directory );
731 
732  // Check if summary plot already exists and remove
733  MonitorElement* me = bei_->get( bei_->pwd() + "/" + name );
734  if ( me ) {
735  bei_->removeElement( name );
736  me = 0;
737  }
738 
739  // Create summary plot
740  float high = static_cast<float>( xbins );
741  if ( pres == sistrip::HISTO_1D ) {
742  if ( xlow < 1. * sistrip::valid_ &&
743  xhigh < 1. * sistrip::valid_ ) {
744  me = bei_->book1D( name, name, xbins, xlow, xhigh );
745  } else {
746  me = bei_->book1D( name, name, xbins, 0., high );
747  }
748  } else if ( pres == sistrip::HISTO_2D_SUM ) {
749  me = bei_->book1D( name, name,
750  xbins, 0., high );
751  } else if ( pres == sistrip::HISTO_2D_SCATTER ) {
752  me = bei_->book2D( name, name, xbins, 0., high,
753  sistrip::FED_ADC_RANGE+1,
754  0.,
755  sistrip::FED_ADC_RANGE*1. );
756  } else if ( pres == sistrip::PROFILE_1D ) {
757  me = bei_->bookProfile( name, name, xbins, 0., high,
758  sistrip::FED_ADC_RANGE+1,
759  0.,
760  sistrip::FED_ADC_RANGE*1. );
761  } else {
762  me = 0;
763  edm::LogWarning(mlDqmClient_)
764  << "[CommissioningHistograms::" << __func__ << "]"
765  << " Unexpected presentation \""
766  << SiStripEnumsAndStrings::presentation( pres )
767  << "\" Unable to create summary plot!";
768  }
769 
770  // Check pointer
771  if ( me ) {
772  LogTrace(mlDqmClient_)
773  << "[CommissioningHistograms::" << __func__ << "]"
774  << " Created summary plot with name \"" << me->getName()
775  << "\" in directory \""
776  << bei_->pwd() << "\"!";
777  } else {
778  edm::LogWarning(mlDqmClient_)
779  << "[CommissioningHistograms::" << __func__ << "]"
780  << " NULL pointer to MonitorElement!"
781  << " Unable to create summary plot!";
782  }
783 
784  // Extract root object
785  TH1* summary = ExtractTObject<TH1>().extract( me );
786  if ( !summary ) {
787  edm::LogWarning(mlDqmClient_)
788  << "[CommissioningHistograms::" << __func__ << "]"
789  << " Unable to extract root object!"
790  << " Returning NULL pointer!";
791  }
792 
793  // Return to pwd
794  bei_->setCurrentFolder( pwd );
795 
796  return summary;
797 
798 }
799 
SummaryGenerator::name
static std::string name(const sistrip::RunType &, const sistrip::Monitorable &, const sistrip::Presentation &, const sistrip::View &, const std::string &directory)
Definition: SummaryGenerator.cc:57
sistrip::runNumber_
static const char runNumber_[]
Definition: ConstantsForDqm.h:33
MonitorElement::getName
const std::string & getName(void) const
get name of ME
Definition: MonitorElement.h:111
CommissioningHistograms::histogram
TH1 * histogram(const sistrip::Monitorable &, const sistrip::Presentation &, const sistrip::View &, const std::string &directory, const uint32_t &xbins, const float &xlow=1.*sistrip::invalid_, const float &xhigh=1.*sistrip::invalid_)
Definition: CommissioningHistograms.cc:717
SiStripFedKey::print
virtual void print(std::stringstream &ss) const
Definition: SiStripFedKey.cc:510
fw3dlego::xbins
const double xbins[]
Definition: fw3dlego_xbins.cc:17
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:942
sistrip::dir_
static const char dir_[]
Definition: ConstantsForDqm.h:17
DQMStore::rmdir
void rmdir(const std::string &fullpath)
Definition: DQMStore.cc:3051
CommissioningHistograms::Histo::print
void print(std::stringstream &) const
Definition: CommissioningHistograms.cc:67
SiStripHistoTitle
Utility class that holds histogram title.
Definition: SiStripHistoTitle.h:20
SiStripFedKey
A container class for generic run and event-related info, information required by the commissioning a...
Definition: SiStripFedKey.h:56
sistrip::invalid32_
static const uint32_t invalid32_
Definition: Constants.h:15
DQMStore::cd
void cd(void)
go to top directory (ie. root)
Definition: DQMStore.cc:632
sistrip::dqmClientFileName_
static const char dqmClientFileName_[]
Definition: ConstantsForDqm.h:26
CommissioningHistograms::getContents
static void getContents(DQMStore *const, std::vector< std::string > &)
Definition: CommissioningHistograms.cc:241
sistrip::mlDqmClient_
static const char mlDqmClient_[]
Definition: ConstantsForLogger.h:19
NULL
#define NULL
Definition: scimark2.h:8
SiStripFecKey::lldChan
const uint16_t & lldChan() const
Definition: SiStripFecKey.h:209
sistrip::FED_ADC_RANGE
static const uint16_t FED_ADC_RANGE
Definition: ConstantsForHardwareSystems.h:17
CommissioningHistograms::save
void save(std::string &filename, uint32_t run_number=0)
Definition: CommissioningHistograms.cc:670
sistrip::valid_
static const uint16_t valid_
Definition: Constants.h:17
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:19
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiStripKey::key
const uint32_t & key() const
Definition: SiStripKey.h:125
cond::ecalcond::bad
int bad(Items const &cont)
Definition: EcalChannelStatusPyWrapper.cc:37
sistrip::undefinedView_
static const char undefinedView_[]
Definition: ConstantsForView.h:18
SiStripFecKey
Utility class that identifies a position within the strip tracker control structure, down to the level of an APV25.
Definition: SiStripFecKey.h:45
SiStripKey::granularity
const sistrip::Granularity & granularity() const
Definition: SiStripKey.h:127
edm::second
U second(std::pair< T, U > const &p)
Definition: ParameterSet.cc:250
SiStripEnumsAndStrings::view
static std::string view(const sistrip::View &)
Definition: SiStripEnumsAndStrings.cc:6
SiStripEnumsAndStrings::runType
static std::string runType(const sistrip::RunType &)
Definition: SiStripEnumsAndStrings.cc:27
DQMStore::removeContents
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
Definition: DQMStore.cc:3085
cmsHarvester.path
tuple path
else: Piece not in the list, fine.
Definition: cmsHarvester.py:1591
CommissioningHistograms::extractHistograms
void extractHistograms(const std::vector< std::string > &)
Definition: CommissioningHistograms.cc:314
DQMStore::bookString
MonitorElement * bookString(const char *name, const char *value)
Book string.
Definition: DQMStore.cc:908
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
CommissioningHistograms::runType
static sistrip::RunType runType(DQMStore *const, const std::vector< std::string > &)
Definition: CommissioningHistograms.cc:157
DQMStore::removeElement
void removeElement(const std::string &name)
Definition: DQMStore.cc:3093
sistrip::taskId_
static const char taskId_[]
Definition: ConstantsForDqm.h:31
SummaryPlotFactory< CommissioningAnalysis * >::init
uint32_t init(const sistrip::Monitorable &, const sistrip::Presentation &, const sistrip::View &, const std::string &top_level_dir, const sistrip::Granularity &, const std::map< uint32_t, CommissioningAnalysis * > &data)
Definition: CommissioningSummaryFactory.cc:6
SiStripKey
Base utility class that identifies a position within a logical structure of the strip tracker...
Definition: SiStripKey.h:23
SiStripKey::path
const std::string & path() const
Definition: SiStripKey.h:126
CommissioningHistograms::runNumber
static uint32_t runNumber(DQMStore *const, const std::vector< std::string > &)
Definition: CommissioningHistograms.cc:82
DQMStore::bookProfile
MonitorElement * bookProfile(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, const char *option="s")
Definition: DQMStore.cc:1256
DQMStore::setVerbose
void setVerbose(unsigned level)
Definition: DQMStore.cc:619
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:1696
DQMStore::getContents
std::vector< MonitorElement * > getContents(const std::string &path) const
Definition: DQMStore.cc:1725
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
SiStripDetKey
Utility class that identifies a position within the strip tracker geometrical structure, down to the level of an APV25 chip.
Definition: SiStripDetKey.h:28
DQMStore::dirExists
bool dirExists(const std::string &path) const
true if directory exists
Definition: DQMStore.cc:720
extract
int extract(std::vector< int > *output, const std::string &dati)
Definition: InvMatrixUtils.cc:313
relmon_rootfiles_spy.contents
tuple contents
Definition: relmon_rootfiles_spy.py:128
debug
#define debug
Definition: HDRShower.cc:19
SiStripEnumsAndStrings::presentation
static std::string presentation(const sistrip::Presentation &)
Definition: SiStripEnumsAndStrings.cc:526
sistrip::invalid_
static const uint16_t invalid_
Definition: Constants.h:16
CommissioningHistograms::createSummaryHisto
virtual void createSummaryHisto(const sistrip::Monitorable &, const sistrip::Presentation &, const std::string &top_level_dir, const sistrip::Granularity &)
Definition: CommissioningHistograms.cc:572
combine.key
list key
Definition: combine.py:13
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", const uint32_t run=0, const uint32_t lumi=0, SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE", const bool resetMEsAfterWriting=false)
Definition: DQMStore.cc:2490
LaserDQM_cfi.mon
tuple mon
Definition: LaserDQM_cfi.py:3
CommissioningHistograms::copyCustomInformation
static void copyCustomInformation(DQMStore *const, const std::vector< std::string > &)
Definition: CommissioningHistograms.cc:247
CommissioningHistograms::bei
DQMStore *const bei() const
Definition: CommissioningHistograms.h:163
dir
dbl *** dir
Definition: mlp_gen.cc:35
DQMStore::bookInt
MonitorElement * bookInt(const char *name)
Book int.
Definition: DQMStore.cc:849
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:1070
sistrip::collate_
static const char collate_[]
Definition: ConstantsForDqm.h:29
CommissioningHistograms::remove
void remove(std::string pattern="")
Definition: CommissioningHistograms.cc:612
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:655
CommissioningHistograms::histoAnalysis
virtual void histoAnalysis(bool debug)
Definition: CommissioningHistograms.cc:480
DQMStore::pwd
const std::string & pwd(void) const
Definition: DQMStore.cc:627
Generated for CMSSW Reference Manual by doxygen 1.8.5