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