CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripMonitorDigi.cc
Go to the documentation of this file.
1 // -*- c++ -*-
4 // Original Author: Dorian Kcira
5 // Created: Sat Feb 4 20:49:10 CET 2006
6 // $Id: SiStripMonitorDigi.cc,v 1.66 2011/10/27 19:49:20 fiori Exp $
7 #include<fstream>
8 #include "TNamed.h"
9 #include "FWCore/Framework/interface/ESHandle.h"
10 #include "FWCore/MessageLogger/interface/MessageLogger.h"
11 #include "FWCore/ServiceRegistry/interface/Service.h"
12 #include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
13 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
14 #include "DataFormats/Common/interface/DetSetNew.h"
15 #include "DataFormats/Common/interface/DetSetVectorNew.h"
16 #include "DataFormats/Common/interface/DetSetVector.h"
17 #include "DataFormats/Common/interface/Handle.h"
18 #include "DataFormats/SiStripDetId/interface/SiStripSubStructure.h"
19 #include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
20 #include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
21 #include "DQM/SiStripCommon/interface/SiStripHistoId.h"
22 #include "DQM/SiStripMonitorDigi/interface/SiStripMonitorDigi.h"
23 #include "DQMServices/Core/interface/DQMStore.h"
24 #include "DQMServices/Core/interface/MonitorElement.h"
25 #include "DPGAnalysis/SiStripTools/interface/APVCyclePhaseCollection.h"
26 #include "DPGAnalysis/SiStripTools/interface/EventWithHistory.h"
27 #include "CalibTracker/SiStripCommon/interface/SiStripDCSStatus.h"
28 
29 #include "CondFormats/SiStripObjects/interface/FedChannelConnection.h"
30 #include "DataFormats/SiStripCommon/interface/SiStripConstants.h"
31 
32 #include "TMath.h"
33 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
34 #include "DataFormats/SiStripDetId/interface/TECDetId.h"
35 #include "DataFormats/SiStripDetId/interface/TIBDetId.h"
36 #include "DataFormats/SiStripDetId/interface/TIDDetId.h"
37 #include "DataFormats/SiStripDetId/interface/TOBDetId.h"
38 
39 #include "iostream"
40 
41 
42 //--------------------------------------------------------------------------------------------
43 SiStripMonitorDigi::SiStripMonitorDigi(const edm::ParameterSet& iConfig) : dqmStore_(edm::Service<DQMStore>().operator->()), conf_(iConfig), show_mechanical_structure_view(true), show_readout_view(false), show_control_view(false), select_all_detectors(true), reset_each_run(false), folder_organizer(), m_cacheID_(0)
44 {
45  firstEvent = -1;
46  eventNb = 0;
47 
48  // Detector Partitions
49  SubDetPhasePartMap["TIB"] = "TI";
50  SubDetPhasePartMap["TID__side__1"] = "TI";
51  SubDetPhasePartMap["TID__side__2"] = "TI";
52  SubDetPhasePartMap["TOB"] = "TO";
53  SubDetPhasePartMap["TEC__side__1"] = "TM";
54  SubDetPhasePartMap["TEC__side__2"] = "TP";
55 
56  // get Digi Producer List
57  digiProducerList = conf_.getParameter<std::vector<edm::InputTag> >("DigiProducersList");
58 
59  //get on/off option for every cluster from cfi
60  edm::ParameterSet ParametersNumberOfDigis = conf_.getParameter<edm::ParameterSet>("TH1NumberOfDigis");
61  layerswitchnumdigison = ParametersNumberOfDigis.getParameter<bool>("layerswitchon");
62  moduleswitchnumdigison = ParametersNumberOfDigis.getParameter<bool>("moduleswitchon");
63 
64  edm::ParameterSet ParametersNumberOfDigisPerStrip = conf_.getParameter<edm::ParameterSet>("TH1NumberOfDigisPerStrip");
65  moduleswitchnumdigispstripon = ParametersNumberOfDigisPerStrip.getParameter<bool>("moduleswitchon");
66 
67  edm::ParameterSet ParametersADCsHottestStrip = conf_.getParameter<edm::ParameterSet>("TH1ADCsHottestStrip");
68  layerswitchadchotteston = ParametersADCsHottestStrip.getParameter<bool>("layerswitchon");
69  moduleswitchadchotteston = ParametersADCsHottestStrip.getParameter<bool>("moduleswitchon");
70 
71  edm::ParameterSet ParametersADCsCoolestStrip = conf_.getParameter<edm::ParameterSet>("TH1ADCsCoolestStrip");
72  layerswitchadccooleston = ParametersADCsCoolestStrip.getParameter<bool>("layerswitchon");
73  moduleswitchadccooleston = ParametersADCsCoolestStrip.getParameter<bool>("moduleswitchon");
74 
75  edm::ParameterSet ParametersDigiADCs = conf_.getParameter<edm::ParameterSet>("TH1DigiADCs");
76  layerswitchdigiadcson = ParametersDigiADCs.getParameter<bool>("layerswitchon");
77  moduleswitchdigiadcson = ParametersDigiADCs.getParameter<bool>("moduleswitchon");
78 
79  edm::ParameterSet ParametersStripOccupancy = conf_.getParameter<edm::ParameterSet>("TH1StripOccupancy");
80  layerswitchstripoccupancyon = ParametersStripOccupancy.getParameter<bool>("layerswitchon");
81  moduleswitchstripoccupancyon = ParametersStripOccupancy.getParameter<bool>("moduleswitchon");
82 
83  edm::ParameterSet ParametersDigiProf = conf_.getParameter<edm::ParameterSet>("TProfNumberOfDigi");
84  layerswitchnumdigisprofon = ParametersDigiProf.getParameter<bool>("layerswitchon");
85 
86  edm::ParameterSet ParametersDigiADC = conf_.getParameter<edm::ParameterSet>("TProfDigiADC");
87  layerswitchdigiadcprofon = ParametersDigiProf.getParameter<bool>("layerswitchon");
88 
89  edm::ParameterSet ParametersTotDigiProf = conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfDigis");
90  subdetswitchtotdigiprofon = ParametersTotDigiProf.getParameter<bool>("subdetswitchon");
91 
92  edm::ParameterSet ParametersTotDigisProfVsLS = conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfDigisVsLS");
93  subdetswitchtotdigiproflson = ParametersTotDigisProfVsLS.getParameter<bool>("subdetswitchon");
94 
95  edm::ParameterSet ParametersTotDigiFailure = conf_.getParameter<edm::ParameterSet>("TotalNumberOfDigisFailure");
96  subdetswitchtotdigifailureon = ParametersTotDigiFailure.getParameter<bool>("subdetswitchon");
97 
98  edm::ParameterSet ParametersDigiApvProf = conf_.getParameter<edm::ParameterSet>("TProfDigiApvCycle");
99  subdetswitchapvcycleprofon = ParametersDigiApvProf.getParameter<bool>("subdetswitchon");
100 
101  edm::ParameterSet ParametersDigiApvTH2 = conf_.getParameter<edm::ParameterSet>("TH2DigiApvCycle");
102  subdetswitchapvcycleth2on = ParametersDigiApvTH2.getParameter<bool>("subdetswitchon");
103 
104  edm::ParameterSet ParametersNApvShots = conf_.getParameter<edm::ParameterSet>("TH1NApvShots");
105  subdetswitchnapvshotson = ParametersNApvShots.getParameter<bool>("subdetswitchon");
106 
107  edm::ParameterSet ParametersNStripApvShots = conf_.getParameter<edm::ParameterSet>("TH1NStripsApvShots");
108  subdetswitchnstripsapvshotson = ParametersNStripApvShots.getParameter<bool>("subdetswitchon");
109 
110  edm::ParameterSet ParametersChargeMedianApvShots = conf_.getParameter<edm::ParameterSet>("TH1ChargeMedianApvShots");
111  subdetswitchchargemedianapvshotson = ParametersChargeMedianApvShots.getParameter<bool>("subdetswitchon");
112 
113  edm::ParameterSet ParametersApvNumberApvShots = conf_.getParameter<edm::ParameterSet>("TH1ApvNumApvShots");
114  subdetswitchapvshotsApvon = ParametersChargeMedianApvShots.getParameter<bool>("subdetswitchon");
115 
116  edm::ParameterSet ParametersNApvShotsProf = conf_.getParameter<edm::ParameterSet>("TProfNShotsVsTime");
117  subdetswitchapvshotsonprof = ParametersNApvShotsProf.getParameter<bool>("subdetswitchon");
118 
119  //Global Histograms
120 
121  edm::ParameterSet ParametersGlobalNApvShots = conf_.getParameter<edm::ParameterSet>("TH1NApvShots");
122  globalswitchnapvshotson = ParametersGlobalNApvShots.getParameter<bool>("globalswitchon");
123 
124  edm::ParameterSet ParametersGlobalNApvShotsProf = conf_.getParameter<edm::ParameterSet>("TProfGlobalNShots");
125  globalsummaryapvshotson = ParametersGlobalNApvShotsProf.getParameter<bool>("globalswitchon");
126 
127  edm::ParameterSet ParametersGlobalNStripApvShots = conf_.getParameter<edm::ParameterSet>("TH1NStripsApvShots");
128  globalswitchnstripsapvshotson = ParametersGlobalNStripApvShots.getParameter<bool>("globalswitchon");
129 
130  edm::ParameterSet ParametersGlobalApvNumApvShots = conf_.getParameter<edm::ParameterSet>("TH1ApvNumApvShots");
131  globalswitchapvshotsApvon = ParametersGlobalApvNumApvShots.getParameter<bool>("globalswitchon");
132 
133  edm::ParameterSet ParametersGlobalChargeMedianApvShots = conf_.getParameter<edm::ParameterSet>("TH1ChargeMedianApvShots");
134  globalswitchchargemedianapvshotson = ParametersGlobalChargeMedianApvShots.getParameter<bool>("globalswitchon");
135 
136  edm::ParameterSet ParametersGlobalNApvShotsTimeProf = conf_.getParameter<edm::ParameterSet>("TProfNShotsVsTime");
137  globalswitchapvshotsonprof = ParametersGlobalNApvShotsTimeProf.getParameter<bool>("globalswitchon");
138 
139  //Digi and APV Shots Maps
140 
141  digitkhistomapon = conf_.getParameter<bool>("TkHistoMap_On");
142 
143  shotshistomapon= conf_.getParameter<bool>("TkHistoMapNApvShots_On");
144  shotsstripshistomapon= conf_.getParameter<bool>("TkHistoMapNStripApvShots_On");
145  shotschargehistomapon= conf_.getParameter<bool>("TkHistoMapMedianChargeApvShots_On");
146 
147  createTrendMEs = conf_.getParameter<bool>("CreateTrendMEs");
148  Mod_On_ = conf_.getParameter<bool>("Mod_On");
149  xLumiProf = conf_.getParameter<int>("xLumiProf");
150 
151  // Event History Producer
152  historyProducer_ = conf_.getParameter<edm::InputTag>("HistoryProducer");
153 
154  // Apv Phase Producer
155  apvPhaseProducer_ = conf_.getParameter<edm::InputTag>("ApvPhaseProducer");
156 
157  // Create DCS Status
158  bool checkDCS = conf_.getParameter<bool>("UseDCSFiltering");
159  if (checkDCS) dcsStatus_ = new SiStripDCSStatus();
160  else dcsStatus_ = 0;
161 }
162 //------------------------------------------------------------------------------------------
163 
164 SiStripMonitorDigi::~SiStripMonitorDigi() {
165  if (dcsStatus_) delete dcsStatus_;
166 }
167 
168 //--------------------------------------------------------------------------------------------
169 void SiStripMonitorDigi::beginRun(const edm::Run& run, const edm::EventSetup& es){
170 
171  if (show_mechanical_structure_view) {
172  unsigned long long cacheID = es.get<SiStripDetCablingRcd>().cacheIdentifier();
173  if (m_cacheID_ != cacheID) {
174  m_cacheID_ = cacheID;
175  edm::LogInfo("SiStripMonitorDigi") <<"SiStripMonitorDigi::beginRun: "
176  << " Creating MEs for new Cabling ";
177  createMEs(es);
178  }
179  } else if (reset_each_run) {
180  edm::LogInfo("SiStripMonitorDigi") <<"SiStripMonitorDigi::beginRun: "
181  << " Resetting MEs ";
182  for (std::map<uint32_t, ModMEs >::const_iterator idet = DigiMEs.begin() ; idet!=DigiMEs.end() ; idet++) {
183  ResetModuleMEs(idet->first);
184  }
185  }
186 }
187 
188 //--------------------------------------------------------------------------------------------
189 void SiStripMonitorDigi::endRun(const edm::Run&, const edm::EventSetup&){
190 }
191 
192 //--------------------------------------------------------------------------------------------
193 void SiStripMonitorDigi::beginLuminosityBlock(const edm::LuminosityBlock& lb, const edm::EventSetup& es){
194  if (lb.id().luminosityBlock() % xLumiProf == 0) {
195  if (subdetswitchtotdigiproflson){
196  if (digiFailureMEs.SubDetTotDigiProfLS) digiFailureMEs.SubDetTotDigiProfLS->Reset();
197  }
198  }
199 }
200 
201 //--------------------------------------------------------------------------------------------
202 void SiStripMonitorDigi::endLuminosityBlock(const edm::LuminosityBlock& lb, const edm::EventSetup& es) {
203  if(lb.id().luminosityBlock() % xLumiProf == 0){
204  if (subdetswitchtotdigifailureon){
205  if (digiFailureMEs.SubDetDigiFailures) digiFailureMEs.SubDetDigiFailures->Reset();
206  MonitorElement * me = dqmStore_->get("SiStrip/MechanicalView/NumberOfDigisInLastLS");
207  if (me){
208  for (int ibin = 1; ibin<7;ibin++){
209  float value = me->getBinContent(ibin);
210  if (value < 50.){
211  digiFailureMEs.SubDetDigiFailures -> Fill(ibin,0);
212  }else{
213  digiFailureMEs.SubDetDigiFailures -> Fill(ibin,1);
214  }
215  }
216  }
217  }
218  }
219 }
220 //--------------------------------------------------------------------------------------------
221 void SiStripMonitorDigi::beginJob(){
222 }
223 
224 
225 //--------------------------------------------------------------------------------------------
226 void SiStripMonitorDigi::createMEs(const edm::EventSetup& es){
227 
228  if ( show_mechanical_structure_view ){
229  // take from eventSetup the SiStripDetCabling object - here will use SiStripDetControl later on
230  es.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
231 
232  // get list of active detectors from SiStripDetCabling
233  std::vector<uint32_t> activeDets;
234  activeDets.clear(); // just in case
235  SiStripDetCabling_->addActiveDetectorsRawIds(activeDets);
236 
237  SiStripSubStructure substructure;
238 
239  // remove any eventual zero elements - there should be none, but just in case
240  for(std::vector<uint32_t>::iterator idets = activeDets.begin(); idets != activeDets.end(); idets++){
241  if(*idets == 0) activeDets.erase(idets);
242  }
243 
244  // create SiStripFolderOrganizer
245  SiStripFolderOrganizer folder_organizer;
246 
247  // Create TkHistoMap for Digi and APV shots properies
248 
249  if (digitkhistomapon) tkmapdigi = new TkHistoMap("SiStrip/TkHisto","TkHMap_NumberOfDigi",0.0,1);
250  if (shotshistomapon) tkmapNApvshots = new TkHistoMap("SiStrip/TkHisto","TkHMap_NApvShots",0.0,1);
251  if (shotsstripshistomapon) tkmapNstripApvshot= new TkHistoMap("SiStrip/TkHisto","TkHMap_NStripApvShots",0.0,1);
252  if (shotschargehistomapon) tkmapMedianChargeApvshots= new TkHistoMap("SiStrip/TkHisto","TkHMap_MedianChargeApvShots",0.0,1);
253 
254  std::vector<uint32_t> tibDetIds;
255 
256  // loop over detectors and book MEs
257  edm::LogInfo("SiStripTkDQM|SiStripMonitorDigi")<<"nr. of activeDets: "<<activeDets.size();
258  for(std::vector<uint32_t>::const_iterator detid_iterator = activeDets.begin(); detid_iterator!=activeDets.end(); detid_iterator++){
259 
260  uint32_t detid = (*detid_iterator);
261 
262  ModMEs local_modmes;
263 
264  local_modmes.NumberOfDigis = 0;
265  local_modmes.NumberOfDigisPerStrip = 0;
266  local_modmes.ADCsHottestStrip = 0;
267  local_modmes.ADCsCoolestStrip = 0;
268  local_modmes.DigiADCs = 0;
269  local_modmes.StripOccupancy = 0;
270 
271  if (Mod_On_) {
272 
273  // set appropriate folder using SiStripFolderOrganizer
274  folder_organizer.setDetectorFolder(detid); // pass the detid to this method
275  if (reset_each_run) ResetModuleMEs(detid);
276  createModuleMEs(local_modmes, detid);
277 
278  // append to DigiMEs
279  DigiMEs.insert( std::make_pair(detid, local_modmes));
280  }
281 
282  // Create Layer Level MEs if they are not created already
283  std::pair<std::string,int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid);
284  SiStripHistoId hidmanager;
285  std::string label = hidmanager.getSubdetid(detid,false);
286 
287  // get detids for the layer
288  std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(label);
289 
290  if(iLayerME==LayerMEsMap.end()) {
291  int32_t lnumber = det_layer_pair.second;
292  std::vector<uint32_t> layerDetIds;
293  if (det_layer_pair.first == "TIB") {
294  substructure.getTIBDetectors(activeDets,layerDetIds,lnumber,0,0,0);
295  } else if (det_layer_pair.first == "TOB") {
296  substructure.getTOBDetectors(activeDets,layerDetIds,lnumber,0,0);
297  } else if (det_layer_pair.first == "TID" && lnumber > 0) {
298  substructure.getTIDDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0);
299  } else if (det_layer_pair.first == "TID" && lnumber < 0) {
300  substructure.getTIDDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0);
301  } else if (det_layer_pair.first == "TEC" && lnumber > 0) {
302  substructure.getTECDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0,0,0);
303  } else if (det_layer_pair.first == "TEC" && lnumber < 0) {
304  substructure.getTECDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0,0,0);
305  }
306 
307  LayerDetMap[label] = layerDetIds;
308 
309  // book Layer plots
310  folder_organizer.setLayerFolder(detid,det_layer_pair.second);
311  createLayerMEs(label, layerDetIds.size());
312  }
313 
314  // book sub-detector plots
315  std::pair<std::string,std::string> sdet_pair = folder_organizer.getSubDetFolderAndTag(detid);
316  if (SubDetMEsMap.find(sdet_pair.second) == SubDetMEsMap.end()){
317  dqmStore_->setCurrentFolder(sdet_pair.first);
318  createSubDetMEs(sdet_pair.second);
319  }
320 
321  }//end of loop over detectors
322 
323  //book some Summary histograms on APV shots in the MechanicalView
324 
325  dqmStore_->setCurrentFolder("SiStrip/MechanicalView/");
326 
327  if (globalsummaryapvshotson) {
328  const char* HistoName = "Summary Mean Apv shots for SubDets";
329  NApvShotsGlobalProf= dqmStore_->bookProfile(HistoName, HistoName,
330  6,0.5,6.5,
331  100, 0., 0., "" );
332  NApvShotsGlobalProf->setBinLabel(1, std::string("TECB"));
333  NApvShotsGlobalProf->setBinLabel(2, std::string("TECF"));
334  NApvShotsGlobalProf->setBinLabel(3, std::string("TIB"));
335  NApvShotsGlobalProf->setBinLabel(4, std::string("TIDB"));
336  NApvShotsGlobalProf->setBinLabel(5, std::string("TIDF"));
337  NApvShotsGlobalProf->setBinLabel(6, std::string("TOB"));
338  NApvShotsGlobalProf->setAxisTitle(" mean APV shots # / evt",2);
339  }
340 
341  //cumulative number of APV shots Vs Time
342  if (globalswitchapvshotsonprof){
343  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfNShotsVsTime");
344  const char* HistoName = "NApv_Shots_vs_Time";
345  ShotsVsTimeApvShotsGlobal=dqmStore_->bookProfile(HistoName,HistoName,
346  Parameters.getParameter<int32_t>("Nbins"),
347  Parameters.getParameter<double>("xmin"),
348  Parameters.getParameter<double>("xmax"),
349  200, //that parameter should not be there !?
350  Parameters.getParameter<double>("ymin"),
351  Parameters.getParameter<double>("ymax"),
352  "" );
353  ShotsVsTimeApvShotsGlobal->setAxisTitle("Time (s)",1);
354  ShotsVsTimeApvShotsGlobal->setAxisTitle("# Apv Shots",2);
355  if (ShotsVsTimeApvShotsGlobal->kind() == MonitorElement::DQM_KIND_TPROFILE) ShotsVsTimeApvShotsGlobal->getTH1()->SetBit(TH1::kCanRebin);
356  }
357 
358  //cumulative number of Strips in APV shots
359  if (globalswitchnstripsapvshotson){
360  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1NStripsApvShots");
361  const char* HistoName = "Number_of_Strips_in_Apv_Shots";
362  StripMultiplicityApvShotsGlobal=dqmStore_->book1D(HistoName,HistoName,
363  Parameters.getParameter<int32_t>("Nbins"),
364  Parameters.getParameter<double>("xmin"),
365  Parameters.getParameter<double>("xmax"));
366  StripMultiplicityApvShotsGlobal->setAxisTitle("# strips in Apv Shots",1);
367  }
368 
369  //cumulative number of APV shots
370  if (globalswitchnapvshotson){
371  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1NApvShots");
372  const char* HistoName = "Number_of_Apv_Shots";
373  NApvShotsGlobal=dqmStore_->book1D(HistoName,HistoName,
374  Parameters.getParameter<int32_t>("Nbins"),
375  Parameters.getParameter<double>("xmin"),
376  Parameters.getParameter<double>("xmax"));
377  NApvShotsGlobal->setAxisTitle("# Apv Shots",1);
378  }
379 
380  //cumulative Median Charge in APV shots
381  if (globalswitchchargemedianapvshotson){
382  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1ChargeMedianApvShots");
383  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
384  const char* HistoName = "Apv_Shots_Charge_Median";
385  MedianChargeApvShotsGlobal=dqmStore_->book1D(HistoName,HistoName,
386  Parameters.getParameter<int32_t>("Nbins"),
387  Parameters.getParameter<double>("xmin"),
388  Parameters.getParameter<double>("xmax"));
389  MedianChargeApvShotsGlobal->setAxisTitle("Apv Shots Charge Median (ADC)",1);
390  }
391 
392  //cmulative APV number with shots
393  if (globalswitchapvshotsApvon){
394  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1ApvNumApvShots");
395  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
396  const char* HistoName = "Apv_Shots_Apv_Number";
397  NApvApvShotsGlobal=dqmStore_->book1D(HistoName,HistoName,
398  Parameters.getParameter<int32_t>("Nbins"),
399  Parameters.getParameter<double>("xmin"),
400  Parameters.getParameter<double>("xmax"));
401  NApvApvShotsGlobal->setAxisTitle("Apv Number",1);
402  }
403 
404  //
405  // Book new histogram to monitor digi in last LS
406  //
407 
408  digiFailureMEs.SubDetTotDigiProfLS = 0;
409  digiFailureMEs.SubDetDigiFailures = 0;
410 
411  std::stringstream ss;
412 
413  folder_organizer.getLayerFolderName(ss, 0);
414  dqmStore_->setCurrentFolder(ss.str().c_str());
415 
416  if (subdetswitchtotdigiproflson) {
417  const char* HistoName = "NumberOfDigisInLastLS";
418  digiFailureMEs.SubDetTotDigiProfLS= dqmStore_->bookProfile(HistoName, HistoName,
419  6,0.5,6.5,
420  100, 0., 10000., "" );
421  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(1, std::string("TECB"));
422  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(2, std::string("TECF"));
423  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(3, std::string("TIB"));
424  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(4, std::string("TIDB"));
425  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(5, std::string("TIDF"));
426  digiFailureMEs.SubDetTotDigiProfLS->setBinLabel(6, std::string("TOB"));
427  }
428 
429  if (subdetswitchtotdigifailureon) {
430  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TotalNumberOfDigisFailure");
431  std::string HistoName = "DataPresentInLastLS";
432  digiFailureMEs.SubDetDigiFailures = dqmStore_->book2D(HistoName,HistoName,
433  Parameters.getParameter<int32_t>("Nbins"),
434  Parameters.getParameter<double>("xmin"),
435  Parameters.getParameter<double>("xmax"),
436  Parameters.getParameter<int32_t>("Nbinsy"),
437  Parameters.getParameter<double>("ymin"),
438  Parameters.getParameter<double>("ymax"));
439  digiFailureMEs.SubDetDigiFailures->setBinLabel(1, std::string("TECB"));
440  digiFailureMEs.SubDetDigiFailures->setBinLabel(2, std::string("TECF"));
441  digiFailureMEs.SubDetDigiFailures->setBinLabel(3, std::string("TIB"));
442  digiFailureMEs.SubDetDigiFailures->setBinLabel(4, std::string("TIDB"));
443  digiFailureMEs.SubDetDigiFailures->setBinLabel(5, std::string("TIDF"));
444  digiFailureMEs.SubDetDigiFailures->setBinLabel(6, std::string("TOB"));
445  }
446  }//end of if
447 
448 }//end of method
449 
450 //--------------------------------------------------------------------------------------------
451 void SiStripMonitorDigi::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){
452 
453  // Filter out events if DCS Event if requested
454  if (dcsStatus_ && !dcsStatus_->getStatus(iEvent, iSetup)) return;
455 
456  TotalNShots=0;
457 
458  runNb = iEvent.id().run();
459  eventNb++;
460 
461  float iOrbitSec = iEvent.orbitNumber()/11223.0;
462 
463  digi_detset_handles.clear();
464  for(std::vector<edm::InputTag>::iterator itDigiProducerList = digiProducerList.begin();
465  itDigiProducerList != digiProducerList.end(); ++itDigiProducerList ) {
466 
467  edm::Handle< edm::DetSetVector<SiStripDigi> > digi_handle;
468 
469  iEvent.getByLabel((*itDigiProducerList),digi_handle);
470 
471  if (digi_handle.isValid()) digi_detset_handles.push_back(digi_handle.product());
472  }
473 
474  // initialise # of clusters to zero
475  for (std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.begin();
476  iSubdet != SubDetMEsMap.end(); iSubdet++) {
477  iSubdet->second.totNDigis = 0;
478  iSubdet->second.SubDetApvShots.clear();
479  }
480 
481  for (std::map<std::string, std::vector< uint32_t > >::const_iterator iterLayer = LayerDetMap.begin();
482  iterLayer != LayerDetMap.end(); iterLayer++) {
483 
484  std::string layer_label = iterLayer->first;
485 
486  std::vector< uint32_t > layer_dets = iterLayer->second;
487  std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(layer_label);
488 
489  //get Layer MEs
490  LayerMEs local_layermes;
491 
492  if(iLayerME == LayerMEsMap.end()) continue;
493  else local_layermes = iLayerME->second;
494 
495  int largest_adc_layer= 0;
496  int smallest_adc_layer= 99999;
497 
498  int ndigi_layer = 0;
499 
500  uint16_t iDet = 0;
501 
502  std::string subdet_label = "";
503 
504  // loop over all modules in the layer
505  for (std::vector< uint32_t >::const_iterator iterDets = layer_dets.begin() ;
506  iterDets != layer_dets.end() ; iterDets++) {
507  iDet++;
508 
509  // detid and type of ME
510  uint32_t detid = (*iterDets);
511 
512  // Get SubDet label once
513  if (subdet_label.size() == 0) subdet_label = folder_organizer.getSubDetFolderAndTag(detid).second;
514 
515  // DetId and corresponding set of MEs
516 
517  std::map<uint32_t, ModMEs >::iterator pos = DigiMEs.find(detid);
518  ModMEs local_modmes = pos->second;
519 
520  // search digis of detid
521  int loc = getDigiSourceIndex(detid);
522 
523  int ndigi_det = 0;
524 
525  if (loc > -1) {
526  ndigi_det = (*(digi_detset_handles[loc]))[detid].size();
527  APVShotFinder theShotFinder = APVShotFinder((*(digi_detset_handles[loc]))[detid]);
528  const std::vector<APVShot>& shots = theShotFinder.getShots();
529  AddApvShotsToSubDet(shots,SubDetMEsMap[subdet_label].SubDetApvShots);
530  if (shotshistomapon) tkmapNApvshots->fill(detid,shots.size());
531  if (shotsstripshistomapon) FillApvShotsMap(tkmapNstripApvshot,shots,detid,1);
532  if (shotschargehistomapon) FillApvShotsMap(tkmapMedianChargeApvshots,shots,detid,2);
533  }
534 
535  if(Mod_On_ && moduleswitchnumdigison && (local_modmes.NumberOfDigis != NULL))
536  (local_modmes.NumberOfDigis)->Fill(ndigi_det);
537 
538  if (layerswitchnumdigisprofon)
539  local_layermes.LayerNumberOfDigisProfile->Fill(iDet*1.0,ndigi_det);
540 
541  if (digitkhistomapon) tkmapdigi->fill(detid,ndigi_det);
542 
543  if (ndigi_det == 0) continue; // no digis for this detid => jump to next step of loop
544 
545  const edm::DetSet<SiStripDigi> & digi_detset = (*(digi_detset_handles[loc]))[detid];
546 
547  ndigi_layer += ndigi_det;
548 
549  // ADCs
550  int largest_adc=(digi_detset.data.begin())->adc();
551  int smallest_adc=(digi_detset.data.begin())->adc();
552 
553 
554  // Check if these parameters are really needed
555  float det_occupancy = 0.0;
556 
557  for(edm::DetSet<SiStripDigi>::const_iterator digiIter = digi_detset.data.begin();
558  digiIter!= digi_detset.data.end(); digiIter++ ){
559 
560  int this_adc = digiIter->adc();
561 
562  if (this_adc > 0.0) det_occupancy++;
563 
564  if(this_adc>largest_adc) largest_adc = this_adc;
565  if(this_adc<smallest_adc) smallest_adc = this_adc;
566 
567  if(Mod_On_ && moduleswitchnumdigispstripon && (local_modmes.NumberOfDigisPerStrip != NULL) && (this_adc > 0.0) )
568  (local_modmes.NumberOfDigisPerStrip)->Fill(digiIter->strip());
569 
570  if(Mod_On_ && moduleswitchdigiadcson && (local_modmes.DigiADCs != NULL) )
571  (local_modmes.DigiADCs)->Fill(static_cast<float>(this_adc));
572 
573  //Fill #ADCs for this digi at layer level
574  if(layerswitchdigiadcson) {
575  fillME(local_layermes.LayerDigiADCs , this_adc);
576  if (createTrendMEs) fillTrend(local_layermes.LayerDigiADCsTrend, this_adc, iOrbitSec);
577  }
578 
579  if (layerswitchdigiadcprofon)
580  local_layermes.LayerDigiADCProfile->Fill(iDet*1.0,this_adc);
581 
582  }//end of loop over digis in this det
583 
584  // Occupancy
585  short nstrips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128;
586  if (nstrips > 0 && det_occupancy > 0 ) {
587  det_occupancy = det_occupancy/nstrips;
588  if (Mod_On_ && moduleswitchstripoccupancyon && (local_modmes.StripOccupancy != NULL))
589  (local_modmes.StripOccupancy)->Fill(det_occupancy);
590  if (layerswitchstripoccupancyon) {
591  fillME(local_layermes.LayerStripOccupancy, det_occupancy);
592  if (createTrendMEs) fillTrend(local_layermes.LayerStripOccupancyTrend, det_occupancy, iOrbitSec);
593  }
594  }
595 
596  if (largest_adc > largest_adc_layer) largest_adc_layer = largest_adc;
597  if (smallest_adc < smallest_adc_layer) smallest_adc_layer = smallest_adc;
598 
599  // nr. of adcs for hottest strip
600  if( Mod_On_ && moduleswitchadchotteston && (local_modmes.ADCsHottestStrip != NULL))
601  (local_modmes.ADCsHottestStrip)->Fill(static_cast<float>(largest_adc));
602 
603  // nr. of adcs for coolest strip
604  if(Mod_On_ && moduleswitchadccooleston && (local_modmes.ADCsCoolestStrip != NULL))
605  (local_modmes.ADCsCoolestStrip)->Fill(static_cast<float>(smallest_adc));
606 
607  }//end of loop over DetIds
608 
609  if(layerswitchnumdigison) {
610  fillME(local_layermes.LayerNumberOfDigis,ndigi_layer);
611  if (createTrendMEs) fillTrend(local_layermes.LayerNumberOfDigisTrend, ndigi_layer, iOrbitSec);
612  }
613  if(layerswitchadchotteston) {
614  fillME(local_layermes.LayerADCsHottestStrip,largest_adc_layer);
615  if (createTrendMEs) fillTrend(local_layermes.LayerADCsHottestStripTrend, largest_adc_layer, iOrbitSec);
616  }
617  if(layerswitchadccooleston) {
618  fillME(local_layermes.LayerADCsCoolestStrip ,smallest_adc_layer);
619  if (createTrendMEs) fillTrend(local_layermes.LayerADCsCoolestStripTrend, smallest_adc_layer, iOrbitSec);
620  }
621 
622  std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(subdet_label);
623  if(iSubdet != SubDetMEsMap.end()) {
624  iSubdet->second.totNDigis += ndigi_layer;
625  //std::cout << " totDigis" << iSubdet->second.totNDigis << " in " << subdet_label << std::endl;
626  }
627  }
628 
629  for (std::map<std::string, SubDetMEs>::iterator it = SubDetMEsMap.begin();
630  it != SubDetMEsMap.end(); it++) {
631 
632  if (subdetswitchtotdigiproflson) {
633  if (strcmp(it->first.c_str(),"TEC__side__1")==0){
634  digiFailureMEs.SubDetTotDigiProfLS->Fill(1, it->second.totNDigis);
635  }else if (strcmp(it->first.c_str(),"TEC__side__2")==0){
636  digiFailureMEs.SubDetTotDigiProfLS->Fill(2, it->second.totNDigis);
637  }else if (strcmp(it->first.c_str(),"TIB")==0){
638  digiFailureMEs.SubDetTotDigiProfLS->Fill(3, it->second.totNDigis);
639  }else if (strcmp(it->first.c_str(),"TID__side__1")==0){
640  digiFailureMEs.SubDetTotDigiProfLS->Fill(4, it->second.totNDigis);
641  }else if (strcmp(it->first.c_str(),"TID__side__2")==0){
642  digiFailureMEs.SubDetTotDigiProfLS->Fill(5, it->second.totNDigis);
643  }else if (strcmp(it->first.c_str(),"TOB")==0){
644  digiFailureMEs.SubDetTotDigiProfLS->Fill(6, it->second.totNDigis);
645  }
646  }
647 
648  if (globalsummaryapvshotson) {
649  if (strcmp(it->first.c_str(),"TEC__side__1")==0){
650  NApvShotsGlobalProf->Fill(1,it->second.SubDetApvShots.size());
651  }else if (strcmp(it->first.c_str(),"TEC__side__2")==0){
652  NApvShotsGlobalProf->Fill(2,it->second.SubDetApvShots.size());
653  }else if (strcmp(it->first.c_str(),"TIB")==0){
654  NApvShotsGlobalProf->Fill(3,it->second.SubDetApvShots.size());
655  }else if (strcmp(it->first.c_str(),"TID__side__1")==0){
656  NApvShotsGlobalProf->Fill(4,it->second.SubDetApvShots.size());
657  }else if (strcmp(it->first.c_str(),"TID__side__2")==0){
658  NApvShotsGlobalProf->Fill(5,it->second.SubDetApvShots.size());
659  }else if (strcmp(it->first.c_str(),"TOB")==0){
660  NApvShotsGlobalProf->Fill(6,it->second.SubDetApvShots.size());
661  }
662  }
663 
664  SubDetMEs subdetmes= it->second;
665  std::string subdet = it->first;
666 
667  // Fill APV shots histograms for SubDet
668 
669  uint ShotsSize=subdetmes.SubDetApvShots.size();
670  TotalNShots+=ShotsSize; //Counter for total Shots in the SiStrip Tracker
671 
672  if (subdetswitchnapvshotson ) subdetmes.SubDetNApvShotsTH1->Fill(ShotsSize);// N shots
673  if (subdetswitchapvshotsonprof) subdetmes.SubDetNApvShotsProf ->Fill(iOrbitSec,ShotsSize); //N shots vs time
674 
675  for (uint i=0; i< ShotsSize; ++i){ // Strip multiplicity, charge median and APV number distributions for APV shots
676 
677  if (subdetswitchapvshotsApvon) subdetmes.SubDetNApvShotsNApvTH1->Fill((subdetmes.SubDetApvShots[i].apvNumber()+1));//APV are defined by 0 to 5 I want 1 to 6
678  if (globalswitchapvshotsApvon) NApvApvShotsGlobal->Fill((subdetmes.SubDetApvShots[i].apvNumber()+1));
679 
680  if (subdetswitchnstripsapvshotson) subdetmes.SubDetNStripsApvShotsTH1->Fill(subdetmes.SubDetApvShots[i].nStrips());
681  if (globalswitchnstripsapvshotson) StripMultiplicityApvShotsGlobal->Fill(subdetmes.SubDetApvShots[i].nStrips());
682 
683  if (subdetswitchchargemedianapvshotson) subdetmes.SubDetChargeMedianApvShotsTH1->Fill(subdetmes.SubDetApvShots[i].median());
684  if (globalswitchchargemedianapvshotson) MedianChargeApvShotsGlobal->Fill(subdetmes.SubDetApvShots[i].median());
685 
686  }
687 
688  if (subdetswitchtotdigiprofon)subdetmes.SubDetTotDigiProf->Fill(iOrbitSec,subdetmes.totNDigis);
689  }
690 
691  if (globalswitchnapvshotson) NApvShotsGlobal->Fill(TotalNShots);
692  if (globalswitchapvshotsonprof) ShotsVsTimeApvShotsGlobal->Fill(iOrbitSec,TotalNShots);
693 
694  // get EventHistory
695 
696  edm::Handle<EventWithHistory> event_history;
697  iEvent.getByLabel(historyProducer_,event_history);
698 
699  // get Phase of APV
700  edm::Handle<APVCyclePhaseCollection> apv_phase_collection;
701  iEvent.getByLabel(apvPhaseProducer_,apv_phase_collection);
702 
703  if (event_history.isValid()
704  && !event_history.failedToGet()
705  && apv_phase_collection.isValid()
706  && !apv_phase_collection.failedToGet()) {
707 
708 
709  long long tbx = event_history->absoluteBX();
710 
711 
712  for (std::map<std::string, SubDetMEs>::iterator it = SubDetMEsMap.begin();
713  it != SubDetMEsMap.end(); it++) {
714 
715  SubDetMEs subdetmes;
716  std::string subdet = it->first;
717  subdetmes = it->second;
718 
719  int the_phase = APVCyclePhaseCollection::invalid;
720  long long tbx_corr = tbx;
721 
722  if (SubDetPhasePartMap.find(subdet) != SubDetPhasePartMap.end()) the_phase = apv_phase_collection->getPhase(SubDetPhasePartMap[subdet]);
723  if(the_phase==APVCyclePhaseCollection::nopartition ||
724  the_phase==APVCyclePhaseCollection::multiphase ||
725  the_phase==APVCyclePhaseCollection::invalid) the_phase=30;
726  tbx_corr -= the_phase;
727 
728  if (subdetswitchapvcycleprofon)subdetmes.SubDetDigiApvProf->Fill(tbx_corr%70,subdetmes.totNDigis);
729  if (subdetswitchapvcycleth2on) subdetmes.SubDetDigiApvTH2->Fill(tbx_corr%70,subdetmes.totNDigis);
730  }
731  }
732 }//end of method analyze
733 //--------------------------------------------------------------------------------------------
734 
735 void SiStripMonitorDigi::endJob(void){
736  bool outputMEsInRootFile = conf_.getParameter<bool>("OutputMEsInRootFile");
737  std::string outputFileName = conf_.getParameter<std::string>("OutputFileName");
738 
739  // save histograms in a file
740  if(outputMEsInRootFile) dqmStore_->save(outputFileName);
741 
742 }//end of method
743 //--------------------------------------------------------------------------------------------
744 void SiStripMonitorDigi::ResetModuleMEs(uint32_t idet){
745  std::map<uint32_t, ModMEs >::iterator pos = DigiMEs.find(idet);
746  ModMEs mod_me = pos->second;
747 
748  if(Mod_On_ && moduleswitchnumdigison) mod_me.NumberOfDigis->Reset();
749  if(Mod_On_ && moduleswitchnumdigispstripon) mod_me.NumberOfDigisPerStrip ->Reset();
750  if(Mod_On_ && moduleswitchadchotteston) mod_me.ADCsHottestStrip->Reset();
751  if(Mod_On_ && moduleswitchadccooleston) mod_me.ADCsCoolestStrip->Reset();
752  if(Mod_On_ && moduleswitchdigiadcson) mod_me.DigiADCs->Reset();
753  if(Mod_On_ && moduleswitchstripoccupancyon) mod_me.StripOccupancy->Reset();
754 
755 }
756 //------------------------------------------------------------------------------------------
757 MonitorElement* SiStripMonitorDigi::bookMETrend(const char* ParameterSetLabel, const char* HistoName)
758 {
759  edm::ParameterSet ParametersTrend = conf_.getParameter<edm::ParameterSet>("Trending");
760  MonitorElement* me = dqmStore_->bookProfile(HistoName,HistoName,
761  ParametersTrend.getParameter<int32_t>("Nbins"),
762  // 0,
763  ParametersTrend.getParameter<double>("xmin"),
764  ParametersTrend.getParameter<double>("xmax"),
765  // ParametersTrend.getParameter<int32_t>("Nbins"),
766  100, //that parameter should not be there !?
767  ParametersTrend.getParameter<double>("ymin"),
768  ParametersTrend.getParameter<double>("ymax"),
769  "" );
770  if(!me) return me;
771 
772  me->setAxisTitle("Event Time in Seconds",1);
773  if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) me->getTH1()->SetBit(TH1::kCanRebin);
774  return me;
775 }
776 
777 //------------------------------------------------------------------------------------------
778 MonitorElement* SiStripMonitorDigi::bookME1D(const char* ParameterSetLabel, const char* HistoName)
779 {
780  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
781  return dqmStore_->book1D(HistoName,HistoName,
782  Parameters.getParameter<int32_t>("Nbinx"),
783  Parameters.getParameter<double>("xmin"),
784  Parameters.getParameter<double>("xmax")
785  );
786 }
787 
788 //--------------------------------------------------------------------------------
789 void SiStripMonitorDigi::fillTrend(MonitorElement* me ,float value, float timeinorbit)
790 {
791  if(!me) return;
792  me->Fill(timeinorbit,value);
793 }
794 
795 //
796 // -- Create Module Level MEs
797 //
798 void SiStripMonitorDigi::createModuleMEs(ModMEs& mod_single, uint32_t detid) {
799 
800  // use SistripHistoId for producing histogram id (and title)
801  SiStripHistoId hidmanager;
802  std::string hid;
803 
804  //nr. of digis per module
805  if(moduleswitchnumdigison) {
806  hid = hidmanager.createHistoId("NumberOfDigis","det",detid);
807  mod_single.NumberOfDigis = dqmStore_->book1D(hid, hid, 21, -0.5, 20.5);
808  dqmStore_->tag(mod_single.NumberOfDigis, detid);
809  mod_single.NumberOfDigis->setAxisTitle("number of digis in one detector module");
810  mod_single.NumberOfDigis->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
811  }
812 
813  //nr. of digis per strip in module
814  if(moduleswitchnumdigispstripon){
815  hid = hidmanager.createHistoId("NumberOfDigisPerStrip","det",detid);
816  short nstrips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128;
817  mod_single.NumberOfDigisPerStrip = dqmStore_->book1D(hid, hid, nstrips, -0.5, nstrips+0.5);
818  dqmStore_->tag(mod_single.NumberOfDigisPerStrip, detid);
819  mod_single.NumberOfDigisPerStrip->setAxisTitle("number of (digis > 0) per strip");
820  mod_single.NumberOfDigisPerStrip->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
821  }
822  //#ADCs for hottest strip
823  if(moduleswitchadchotteston) {
824  hid = hidmanager.createHistoId("ADCsHottestStrip","det",detid);
825  mod_single.ADCsHottestStrip = bookME1D("TH1ADCsHottestStrip", hid.c_str());
826  dqmStore_->tag(mod_single.ADCsHottestStrip, detid); // 6 APVs -> 768 strips
827  mod_single.ADCsHottestStrip->setAxisTitle("number of ADCs for hottest strip");
828  }
829 
830  //#ADCs for coolest strip
831  if(moduleswitchadccooleston) {
832  hid = hidmanager.createHistoId("ADCsCoolestStrip","det",detid);
833  mod_single.ADCsCoolestStrip = bookME1D("TH1ADCsCoolestStrip", hid.c_str());
834  dqmStore_->tag(mod_single.ADCsCoolestStrip, detid);
835  mod_single.ADCsCoolestStrip->setAxisTitle("number of ADCs for coolest strip");
836  }
837 
838  //#ADCs for each digi
839  if(moduleswitchdigiadcson) {
840  hid = hidmanager.createHistoId("DigiADCs","det",detid);
841  mod_single.DigiADCs = bookME1D("TH1DigiADCs", hid.c_str());
842  dqmStore_->tag(mod_single.DigiADCs, detid);
843  mod_single.DigiADCs->setAxisTitle("number of ADCs for each digi");
844  }
845 
846  //Strip occupancy
847  if(moduleswitchstripoccupancyon) {
848  hid = hidmanager.createHistoId("StripOccupancy","det",detid);
849  mod_single.StripOccupancy = bookME1D("TH1StripOccupancy", hid.c_str());
850  dqmStore_->tag(mod_single.StripOccupancy, detid);
851  mod_single.StripOccupancy->setAxisTitle("strip occupancy");
852  }
853 
854 }
855 
856 //
857 // -- Create Module Level MEs
858 //
859 
860 void SiStripMonitorDigi::createLayerMEs(std::string label, int ndets) {
861 
862  std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(label);
863  if(iLayerME==LayerMEsMap.end()){
864  SiStripHistoId hidmanager;
865  LayerMEs layerMEs;
866  layerMEs.LayerNumberOfDigis = 0;
867  layerMEs.LayerNumberOfDigisTrend = 0;
868  layerMEs.LayerADCsHottestStrip = 0;
869  layerMEs.LayerADCsHottestStripTrend = 0;
870  layerMEs.LayerADCsCoolestStrip = 0;
871  layerMEs.LayerADCsCoolestStripTrend = 0;
872  layerMEs.LayerDigiADCs = 0;
873  layerMEs.LayerDigiADCsTrend = 0;
874  layerMEs.LayerStripOccupancy = 0;
875  layerMEs.LayerStripOccupancyTrend = 0;
876  layerMEs.LayerNumberOfDigisProfile = 0;
877  layerMEs.LayerDigiADCProfile = 0;
878 
879 
880  //#Digis
881  if(layerswitchnumdigison) {
882  layerMEs.LayerNumberOfDigis=bookME1D("TH1NumberOfDigis", hidmanager.createHistoLayer("Summary_TotalNumberOfDigis","layer",label,"").c_str());
883  if (createTrendMEs) layerMEs.LayerNumberOfDigisTrend=bookMETrend("TH1NumberOfDigis", hidmanager.createHistoLayer("Trend_NumberOfDigis","layer",label,"").c_str());
884  }
885 
886  //#ADCs for hottest strip
887  if(layerswitchadchotteston) {
888  layerMEs.LayerADCsHottestStrip=bookME1D("TH1ADCsHottestStrip", hidmanager.createHistoLayer("Summary_ADCsHottestStrip","layer",label,"").c_str());
889  if (createTrendMEs) layerMEs.LayerADCsHottestStripTrend=bookMETrend("TH1ADCsHottestStrip", hidmanager.createHistoLayer("Trend_ADCsHottestStrip","layer",label,"").c_str());
890  }
891 
892  //#ADCs for coolest strip
893  if(layerswitchadccooleston) {
894  layerMEs.LayerADCsCoolestStrip=bookME1D("TH1ADCsCoolestStrip", hidmanager.createHistoLayer("Summary_ADCsCoolestStrip","layer",label,"").c_str());
895  if (createTrendMEs) layerMEs.LayerADCsCoolestStripTrend=bookMETrend("TH1ADCsCoolestStrip", hidmanager.createHistoLayer("Trend_ADCsCoolestStrip","layer",label,"").c_str());
896  }
897 
898  //#ADCs for each digi
899  if(layerswitchdigiadcson) {
900  layerMEs.LayerDigiADCs=bookME1D("TH1DigiADCs", hidmanager.createHistoLayer("Summary_DigiADCs","layer",label,"").c_str());
901  if (createTrendMEs) layerMEs.LayerDigiADCsTrend=bookMETrend("TH1DigiADCs", hidmanager.createHistoLayer("Trend_DigiADCs","layer",label,"").c_str());
902  }
903 
904  //Strip Occupancy
905  if(layerswitchstripoccupancyon) {
906  layerMEs.LayerStripOccupancy=bookME1D("TH1StripOccupancy", hidmanager.createHistoLayer("Summary_StripOccupancy","layer",label,"").c_str());
907  if (createTrendMEs) layerMEs.LayerStripOccupancyTrend=bookMETrend("TH1StripOccupancy", hidmanager.createHistoLayer("Trend_StripOccupancy","layer",label,"").c_str());
908 
909  }
910  // # of Digis
911  if(layerswitchnumdigisprofon) {
912  std::string hid = hidmanager.createHistoLayer("NumberOfDigiProfile","layer",label,"");
913  layerMEs.LayerNumberOfDigisProfile = dqmStore_->bookProfile(hid, hid, ndets, 0.5, ndets+0.5,21, -0.5, 200.5);
914  }
915 
916  // # of Digis
917  if(layerswitchdigiadcprofon) {
918  std::string hid = hidmanager.createHistoLayer("DigiADCProfile","layer",label,"");
919  layerMEs.LayerDigiADCProfile = dqmStore_->bookProfile(hid, hid, ndets, 0.5, ndets+0.5, 64, -0.5, 255.5);
920  }
921 
922  LayerMEsMap[label]=layerMEs;
923  }
924 }
925 //
926 // -- Create SubDetector MEs
927 //
928 void SiStripMonitorDigi::createSubDetMEs(std::string label) {
929 
930  SubDetMEs subdetMEs;
931  subdetMEs.totNDigis = 0;
932  subdetMEs.SubDetTotDigiProf = 0;
933  subdetMEs.SubDetDigiApvProf = 0;
934  subdetMEs.SubDetDigiApvTH2 = 0;
935 
936  subdetMEs.SubDetApvShots.clear();
937  subdetMEs.SubDetNApvShotsTH1 = 0;
938  subdetMEs.SubDetChargeMedianApvShotsTH1 = 0;
939  subdetMEs.SubDetNStripsApvShotsTH1 = 0;
940  subdetMEs.SubDetNApvShotsProf = 0;
941 
942  std::string HistoName;
943 
944  // Total Number of Digi - Profile
945  if(subdetswitchtotdigiprofon){
946  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfDigis");
947  HistoName = "TotalNumberOfDigiProfile__" + label;
948  subdetMEs.SubDetTotDigiProf=dqmStore_->bookProfile(HistoName,HistoName,
949  Parameters.getParameter<int32_t>("Nbins"),
950  Parameters.getParameter<double>("xmin"),
951  Parameters.getParameter<double>("xmax"),
952  100, //that parameter should not be there !?
953  Parameters.getParameter<double>("ymin"),
954  Parameters.getParameter<double>("ymax"),
955  "" );
956  subdetMEs.SubDetTotDigiProf->setAxisTitle("Event Time in Seconds",1);
957  if (subdetMEs.SubDetTotDigiProf->kind() == MonitorElement::DQM_KIND_TPROFILE) subdetMEs.SubDetTotDigiProf->getTH1()->SetBit(TH1::kCanRebin);
958  }
959 
960  // Number of Digi vs Bx - Profile
961  if(subdetswitchapvcycleprofon){
962  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfDigiApvCycle");
963  HistoName = "Digi_vs_ApvCycle__" + label;
964  subdetMEs.SubDetDigiApvProf=dqmStore_->bookProfile(HistoName,HistoName,
965  Parameters.getParameter<int32_t>("Nbins"),
966  Parameters.getParameter<double>("xmin"),
967  Parameters.getParameter<double>("xmax"),
968  200, //that parameter should not be there !?
969  Parameters.getParameter<double>("ymin"),
970  Parameters.getParameter<double>("ymax"),
971  "" );
972  subdetMEs.SubDetDigiApvProf->setAxisTitle("ApvCycle (Corrected Absolute Bx % 70)",1);
973  }
974 
975  // Number of Digi vs Bx - TH2
976  if(subdetswitchapvcycleth2on){
977  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH2DigiApvCycle");
978  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
979  HistoName = "Digi_vs_ApvCycle_2D__" + label;
980  // Adjusting the scale for 2D histogram
981  double h2ymax = 9999.0;
982  double yfact = Parameters.getParameter<double>("yfactor");
983  if(label.find("TIB") != std::string::npos) h2ymax = (6984.*256.)*yfact;
984  else if (label.find("TID") != std::string::npos) h2ymax = (2208.*256.)*yfact;
985  else if (label.find("TOB") != std::string::npos) h2ymax = (12906.*256.)*yfact;
986  else if (label.find("TEC") != std::string::npos) h2ymax = (7552.*2.*256.)*yfact;
987  subdetMEs.SubDetDigiApvTH2=dqmStore_->book2D(HistoName,HistoName,
988  Parameters.getParameter<int32_t>("Nbins"),
989  Parameters.getParameter<double>("xmin"),
990  Parameters.getParameter<double>("xmax"),
991  Parameters.getParameter<int32_t>("Nbinsy"), //it was 100 that parameter should not be there !?
992  Parameters.getParameter<double>("ymin"),
993  h2ymax);
994  subdetMEs.SubDetDigiApvTH2->setAxisTitle("absolute Bx mod(70)",1);
995  }
996 
997  //Number of APV Shots
998  if (subdetswitchnapvshotson){
999  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1NApvShots");
1000  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
1001  HistoName = "Number_of_Apv_Shots_" + label;
1002  subdetMEs.SubDetNApvShotsTH1=dqmStore_->book1D(HistoName,HistoName,
1003  Parameters.getParameter<int32_t>("Nbins"),
1004  Parameters.getParameter<double>("xmin"),
1005  Parameters.getParameter<double>("xmax"));
1006  subdetMEs.SubDetNApvShotsTH1->setAxisTitle("# Apv Shots",1);
1007  }
1008 
1009  //Strip multiplicity of APV Shots
1010  if (subdetswitchnstripsapvshotson){
1011  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1NStripsApvShots");
1012  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
1013  HistoName = "Number_of_Strips_in_Apv_Shots_" + label;
1014  subdetMEs.SubDetNStripsApvShotsTH1=dqmStore_->book1D(HistoName,HistoName,
1015  Parameters.getParameter<int32_t>("Nbins"),
1016  Parameters.getParameter<double>("xmin"),
1017  Parameters.getParameter<double>("xmax"));
1018  subdetMEs.SubDetNStripsApvShotsTH1->setAxisTitle("# strips in Apv Shots",1);
1019  }
1020 
1021  //Charge median of APV Shots
1022  if (subdetswitchchargemedianapvshotson){
1023  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1ChargeMedianApvShots");
1024  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
1025  HistoName = "Apv_Shots_Charge_Median_" + label;
1026  subdetMEs.SubDetChargeMedianApvShotsTH1=dqmStore_->book1D(HistoName,HistoName,
1027  Parameters.getParameter<int32_t>("Nbins"),
1028  Parameters.getParameter<double>("xmin"),
1029  Parameters.getParameter<double>("xmax"));
1030  subdetMEs.SubDetChargeMedianApvShotsTH1->setAxisTitle("Apv Shots Charge Median (ADC)",1);
1031  }
1032 
1033 
1034  if (subdetswitchchargemedianapvshotson){
1035  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH1ApvNumApvShots");
1036  //dqmStore_->setCurrentFolder("SiStrip/MechanicalView/"+label);
1037  HistoName = "Apv_Shots_Apv_Number_" + label;
1038  subdetMEs.SubDetNApvShotsNApvTH1=dqmStore_->book1D(HistoName,HistoName,
1039  Parameters.getParameter<int32_t>("Nbins"),
1040  Parameters.getParameter<double>("xmin"),
1041  Parameters.getParameter<double>("xmax"));
1042  subdetMEs.SubDetNApvShotsNApvTH1->setAxisTitle("Apv Number",1);
1043  }
1044 
1045 
1046 
1047 
1048 
1049  //APV Shots number Vs time
1050  if(subdetswitchapvshotsonprof){
1051  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfNShotsVsTime");
1052  HistoName = "NApv_Shots_vs_Time_" + label;
1053  subdetMEs.SubDetNApvShotsProf=dqmStore_->bookProfile(HistoName,HistoName,
1054  Parameters.getParameter<int32_t>("Nbins"),
1055  Parameters.getParameter<double>("xmin"),
1056  Parameters.getParameter<double>("xmax"),
1057  200, //that parameter should not be there !?
1058  Parameters.getParameter<double>("ymin"),
1059  Parameters.getParameter<double>("ymax"),
1060  "" );
1061  subdetMEs.SubDetNApvShotsProf->setAxisTitle("Time (s)",1);
1062  subdetMEs.SubDetNApvShotsProf->setAxisTitle("# Apv Shots",2);
1063  if (subdetMEs.SubDetNApvShotsProf->kind() == MonitorElement::DQM_KIND_TPROFILE) subdetMEs.SubDetNApvShotsProf->getTH1()->SetBit(TH1::kCanRebin);
1064  }
1065 
1066 
1067 
1068 
1069 
1070 
1071 
1072  SubDetMEsMap[label]=subdetMEs;
1073 }
1074 //
1075 // -- Get DetSet vector for a given Detector
1076 //
1077 int SiStripMonitorDigi::getDigiSourceIndex(uint32_t id) {
1078  int location = -1;
1079  for (unsigned int ival = 0; ival < digi_detset_handles.size(); ++ival){
1080  edm::DetSetVector<SiStripDigi>::const_iterator isearch = digi_detset_handles[ival]->find(id);
1081  if(isearch != digi_detset_handles[ival]->end()) {
1082  location = ival;
1083  break;
1084  }
1085  }
1086  return location;
1087 }
1088 
1089 void SiStripMonitorDigi::AddApvShotsToSubDet(const std::vector<APVShot> & moduleShots, std::vector<APVShot> & subdetShots){
1090 
1091  for (uint i=0; i<moduleShots.size(); i++){
1092  subdetShots.push_back(moduleShots[i]);
1093  }
1094 }
1095 
1096 void SiStripMonitorDigi::FillApvShotsMap(TkHistoMap* the_map, const std::vector<APVShot> & shots, uint32_t id ,int mode){
1097 
1098  for (uint i=0; i<shots.size(); i++){
1099  if (mode==1) the_map->fill(id,shots[i].nStrips()); //mode == 1 fill with strip multiplicity
1100  if (mode==2) the_map->fill(id,shots[i].median()); // mode == 2 fill with charge median
1101  }
1102 }
1103 
1104 
1105 //define this as a plug-in
1106 DEFINE_FWK_MODULE(SiStripMonitorDigi);
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
edm::LuminosityBlockBase::id
LuminosityBlockID id() const
Definition: LuminosityBlockBase.h:53
SiStripMonitorDigi::ResetModuleMEs
void ResetModuleMEs(uint32_t idet)
Definition: SiStripMonitorDigi.cc:744
SiStripMonitorDigi::createSubDetMEs
void createSubDetMEs(std::string label)
Definition: SiStripMonitorDigi.cc:928
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
i
int i
Definition: DBlmapReader.cc:9
SiStripMonitorDigi::LayerMEsMap
std::map< std::string, LayerMEs > LayerMEsMap
Definition: SiStripMonitorDigi.h:122
SiStripMonitorDigi::LayerMEs::LayerNumberOfDigisTrend
MonitorElement * LayerNumberOfDigisTrend
Definition: SiStripMonitorDigi.h:54
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:717
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
SiStripMonitorDigi::LayerMEs::LayerStripOccupancy
MonitorElement * LayerStripOccupancy
Definition: SiStripMonitorDigi.h:61
SiStripMonitorDigi::LayerMEs::LayerADCsHottestStrip
MonitorElement * LayerADCsHottestStrip
Definition: SiStripMonitorDigi.h:55
DQMStore::save
void save(const std::string &filename, const std::string &path="", const std::string &pattern="", const std::string &rewrite="", SaveReferenceTag ref=SaveWithReference, int minStatus=dqm::qstatus::STATUS_OK, const std::string &fileupdate="RECREATE")
Definition: DQMStore.cc:2113
SiStripMonitorDigi::historyProducer_
edm::InputTag historyProducer_
Definition: SiStripMonitorDigi.h:195
SiStripMonitorDigi::createModuleMEs
void createModuleMEs(ModMEs &mod_single, uint32_t detid)
Definition: SiStripMonitorDigi.cc:798
SiStripMonitorDigi::LayerMEs::LayerADCsHottestStripTrend
MonitorElement * LayerADCsHottestStripTrend
Definition: SiStripMonitorDigi.h:56
MonitorElement::setBinLabel
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:889
abs
#define abs(x)
Definition: mlp_lapack.h:159
SiStripMonitorDigi::NApvShotsGlobalProf
MonitorElement * NApvShotsGlobalProf
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::DigiMEs
std::map< uint32_t, ModMEs > DigiMEs
Definition: SiStripMonitorDigi.h:118
NULL
#define NULL
Definition: scimark2.h:8
SiStripMonitorDigi::SubDetPhasePartMap
std::map< std::string, std::string > SubDetPhasePartMap
Definition: SiStripMonitorDigi.h:124
SiStripMonitorDigi::MedianChargeApvShotsGlobal
MonitorElement * MedianChargeApvShotsGlobal
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::StripMultiplicityApvShotsGlobal
MonitorElement * StripMultiplicityApvShotsGlobal
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::tkmapNstripApvshot
TkHistoMap * tkmapNstripApvshot
Definition: SiStripMonitorDigi.h:139
alignBH_cfg.mode
list mode
Definition: alignBH_cfg.py:7
SiStripMonitorDigi::LayerMEs::LayerNumberOfDigisProfile
MonitorElement * LayerNumberOfDigisProfile
Definition: SiStripMonitorDigi.h:63
Association.map
dictionary map
Definition: Association.py:205
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:152
SiStripMonitorDigi::SubDetMEs::SubDetApvShots
std::vector< APVShot > SubDetApvShots
Definition: SiStripMonitorDigi.h:76
SiStripMonitorDigi::beginLuminosityBlock
virtual void beginLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
Definition: SiStripMonitorDigi.cc:193
DQMStore::tag
void tag(MonitorElement *me, unsigned int myTag)
Definition: DQMStore.cc:1354
SiStripMonitorDigi::m_cacheID_
unsigned long long m_cacheID_
Definition: SiStripMonitorDigi.h:132
SiStripMonitorDigi::SubDetMEs::SubDetNApvShotsNApvTH1
MonitorElement * SubDetNApvShotsNApvTH1
Definition: SiStripMonitorDigi.h:81
SiStripDCSStatus::getStatus
bool getStatus(edm::Event const &e, edm::EventSetup const &eSetup)
Definition: SiStripDCSStatus.cc:46
dumpDBToFile_GT_ttrig_cfg.outputFileName
tuple outputFileName
Definition: dumpDBToFile_GT_ttrig_cfg.py:31
SiStripMonitorDigi::ModMEs::NumberOfDigisPerStrip
MonitorElement * NumberOfDigisPerStrip
Definition: SiStripMonitorDigi.h:44
APVShotFinder::getShots
const std::vector< APVShot > & getShots() const
Definition: APVShotFinder.cc:84
iEvent
int iEvent
Definition: GenABIO.cc:243
SiStripHistoId::getSubdetid
std::string getSubdetid(uint32_t id, bool flag_ring)
Definition: SiStripHistoId.cc:87
SiStripMonitorDigi::ModMEs::ADCsCoolestStrip
MonitorElement * ADCsCoolestStrip
Definition: SiStripMonitorDigi.h:46
Parameters
vector< ParameterSet > Parameters
Definition: HLTMuonPlotter.cc:34
SiStripFolderOrganizer::getSubDetFolderAndTag
std::pair< std::string, std::string > getSubDetFolderAndTag(const uint32_t &detid)
Definition: SiStripFolderOrganizer.cc:355
SiStripMonitorDigi::dcsStatus_
SiStripDCSStatus * dcsStatus_
Definition: SiStripMonitorDigi.h:198
SiStripMonitorDigi::SiStripDetCabling_
edm::ESHandle< SiStripDetCabling > SiStripDetCabling_
Definition: SiStripMonitorDigi.h:133
SiStripSubStructure::getTOBDetectors
void getTOBDetectors(const std::vector< uint32_t > &inputDetRawIds, std::vector< uint32_t > &tobDetRawIds, uint32_t layer=0, uint32_t bkw_frw=0, uint32_t rod=0) const
Definition: SiStripSubStructure.cc:78
SiStripMonitorDigi::bookMETrend
MonitorElement * bookMETrend(const char *ParameterSetLabel, const char *HistoName)
Definition: SiStripMonitorDigi.cc:757
SiStripMonitorDigi::LayerMEs::LayerADCsCoolestStripTrend
MonitorElement * LayerADCsCoolestStripTrend
Definition: SiStripMonitorDigi.h:58
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
SiStripMonitorDigi::SubDetMEs::SubDetNStripsApvShotsTH1
MonitorElement * SubDetNStripsApvShotsTH1
Definition: SiStripMonitorDigi.h:79
SiStripMonitorDigi::digiFailureMEs
DigiFailureMEs digiFailureMEs
Definition: SiStripMonitorDigi.h:125
TkHistoMap::fill
void fill(uint32_t &detid, float value)
Definition: TkHistoMap.cc:130
SiStripMonitorDigi::fillME
void fillME(MonitorElement *ME, float value1)
Definition: SiStripMonitorDigi.h:98
SiStripMonitorDigi::NApvApvShotsGlobal
MonitorElement * NApvApvShotsGlobal
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::digi_detset_handles
std::vector< const edm::DetSetVector< SiStripDigi > * > digi_detset_handles
Definition: SiStripMonitorDigi.h:130
SiStripMonitorDigi::createLayerMEs
void createLayerMEs(std::string label, int ndet)
Definition: SiStripMonitorDigi.cc:860
SiStripFolderOrganizer::getLayerFolderName
void getLayerFolderName(std::stringstream &ss, uint32_t rawdetid, bool ring_flag=0)
Definition: SiStripFolderOrganizer.cc:315
SiStripMonitorDigi::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: SiStripMonitorDigi.cc:451
MonitorElement::getTH1
TH1 * getTH1(void) const
Definition: MonitorElement.cc:1479
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:1031
end
#define end
Definition: vmac.h:38
edm::EventBase::orbitNumber
int orbitNumber() const
Definition: EventBase.h:63
MonitorElement::kind
Kind kind(void) const
Get the type of the monitor element.
Definition: MonitorElement.h:91
SiStripFolderOrganizer::setDetectorFolder
void setDetectorFolder(uint32_t rawdetid=0)
Definition: SiStripFolderOrganizer.cc:154
SiStripMonitorDigi::SubDetMEsMap
std::map< std::string, SubDetMEs > SubDetMEsMap
Definition: SiStripMonitorDigi.h:123
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:1468
SiStripMonitorDigi::ModMEs::ADCsHottestStrip
MonitorElement * ADCsHottestStrip
Definition: SiStripMonitorDigi.h:45
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
edm::Event::getByLabel
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:356
SiStripMonitorDigi::fillTrend
void fillTrend(MonitorElement *me, float value, float timeinorbit)
Definition: SiStripMonitorDigi.cc:789
SiStripMonitorDigi::NApvShotsGlobal
MonitorElement * NApvShotsGlobal
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::LayerMEs::LayerDigiADCProfile
MonitorElement * LayerDigiADCProfile
Definition: SiStripMonitorDigi.h:64
edm::HandleBase::failedToGet
bool failedToGet() const
Definition: HandleBase.h:80
SiStripMonitorDigi::getDigiSourceIndex
int getDigiSourceIndex(uint32_t id)
Definition: SiStripMonitorDigi.cc:1077
SiStripMonitorDigi::SiStripMonitorDigi
SiStripMonitorDigi(const edm::ParameterSet &)
Definition: SiStripMonitorDigi.cc:43
SiStripMonitorDigi::apvPhaseProducer_
edm::InputTag apvPhaseProducer_
Definition: SiStripMonitorDigi.h:196
SiStripMonitorDigi::LayerDetMap
std::map< std::string, std::vector< uint32_t > > LayerDetMap
Definition: SiStripMonitorDigi.h:121
SiStripHistoId::createHistoId
std::string createHistoId(std::string description, std::string id_type, uint32_t component_id)
Definition: SiStripHistoId.cc:36
SiStripMonitorDigi::SubDetMEs::SubDetChargeMedianApvShotsTH1
MonitorElement * SubDetChargeMedianApvShotsTH1
Definition: SiStripMonitorDigi.h:78
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
edm::LuminosityBlockID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockID.h:43
edm::Handle::product
T const * product() const
Definition: Handle.h:74
SiStripMonitorDigi::conf_
edm::ParameterSet conf_
Definition: SiStripMonitorDigi.h:116
cscdqm::HistoName
std::string HistoName
Definition: CSCDQM_HistoDef.h:32
SiStripSubStructure::getTECDetectors
void getTECDetectors(const std::vector< uint32_t > &inputDetRawIds, std::vector< uint32_t > &tecDetRawIds, uint32_t side=0, uint32_t wheel=0, uint32_t petal_bkw_frw=0, uint32_t petal=0, uint32_t ring=0, uint32_t ster=0) const
Definition: SiStripSubStructure.cc:100
SiStripMonitorDigi::digiProducerList
std::vector< edm::InputTag > digiProducerList
Definition: SiStripMonitorDigi.h:117
SiStripMonitorDigi::tkmapMedianChargeApvshots
TkHistoMap * tkmapMedianChargeApvshots
Definition: SiStripMonitorDigi.h:139
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:712
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
edm::DetSet::data
collection_type data
Definition: DetSet.h:79
SiStripMonitorDigi::tkmapNApvshots
TkHistoMap * tkmapNApvshots
Definition: SiStripMonitorDigi.h:139
SiStripMonitorDigi::LayerMEs::LayerStripOccupancyTrend
MonitorElement * LayerStripOccupancyTrend
Definition: SiStripMonitorDigi.h:62
SiStripSubStructure::getTIBDetectors
void getTIBDetectors(const std::vector< uint32_t > &inputDetRawIds, std::vector< uint32_t > &tibDetRawIds, uint32_t layer=0, uint32_t bkw_frw=0, uint32_t int_ext=0, uint32_t string=0) const
Definition: SiStripSubStructure.cc:30
SiStripMonitorDigi::endRun
virtual void endRun(const edm::Run &, const edm::EventSetup &)
Definition: SiStripMonitorDigi.cc:189
SiStripMonitorDigi::LayerMEs::LayerADCsCoolestStrip
MonitorElement * LayerADCsCoolestStrip
Definition: SiStripMonitorDigi.h:57
SiStripFolderOrganizer::GetSubDetAndLayer
std::pair< std::string, int32_t > GetSubDetAndLayer(const uint32_t &detid, bool ring_flag=0)
Definition: SiStripFolderOrganizer.cc:120
SiStripMonitorDigi::AddApvShotsToSubDet
void AddApvShotsToSubDet(const std::vector< APVShot > &, std::vector< APVShot > &)
Definition: SiStripMonitorDigi.cc:1089
SiStripMonitorDigi::endLuminosityBlock
virtual void endLuminosityBlock(const edm::LuminosityBlock &, const edm::EventSetup &)
Definition: SiStripMonitorDigi.cc:202
SiStripMonitorDigi::folder_organizer
SiStripFolderOrganizer folder_organizer
Definition: SiStripMonitorDigi.h:128
SiStripSubStructure::getTIDDetectors
void getTIDDetectors(const std::vector< uint32_t > &inputDetRawIds, std::vector< uint32_t > &tidDetRawIds, uint32_t side=0, uint32_t wheel=0, uint32_t ring=0, uint32_t ster=0) const
Definition: SiStripSubStructure.cc:54
SiStripMonitorDigi::ShotsVsTimeApvShotsGlobal
MonitorElement * ShotsVsTimeApvShotsGlobal
Definition: SiStripMonitorDigi.h:137
SiStripMonitorDigi::createMEs
void createMEs(const edm::EventSetup &es)
Definition: SiStripMonitorDigi.cc:226
SiStripHistoId::createHistoLayer
std::string createHistoLayer(std::string description, std::string id_type, std::string path, std::string flag)
Definition: SiStripHistoId.cc:61
edm::DetSet::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSet.h:34
SiStripMonitorDigi::beginRun
virtual void beginRun(const edm::Run &, const edm::EventSetup &)
Definition: SiStripMonitorDigi.cc:169
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:845
edm::DetSetVector::const_iterator
collection_type::const_iterator const_iterator
Definition: DetSetVector.h:106
SiStripFolderOrganizer::setLayerFolder
void setLayerFolder(uint32_t rawdetid=0, int32_t layer=0, bool ring_flag=0)
Definition: SiStripFolderOrganizer.cc:234
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:916
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:422
SiStripMonitorDigi::FillApvShotsMap
void FillApvShotsMap(TkHistoMap *, const std::vector< APVShot > &, uint32_t id, int)
Definition: SiStripMonitorDigi.cc:1096
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:429
edm::Run
Definition: Run.h:33
SiStripMonitorDigi::bookME1D
MonitorElement * bookME1D(const char *ParameterSetLabel, const char *HistoName)
Definition: SiStripMonitorDigi.cc:778
Generated for CMSSW Reference Manual by doxygen 1.8.5