CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
SiStripMonitorCluster.cc
Go to the documentation of this file.
1 
2 // -*- C++ -*-
3 // Package: SiStripMonitorCluster
4 // Class: SiStripMonitorCluster
7 // Original Author: Dorian Kcira
8 // Created: Wed Feb 1 16:42:34 CET 2006
9 #include <vector>
10 #include <numeric>
11 #include <fstream>
12 #include <math.h>
13 #include "TNamed.h"
14 #include "FWCore/MessageLogger/interface/MessageLogger.h"
15 #include "FWCore/ServiceRegistry/interface/Service.h"
16 
17 #include "CalibTracker/Records/interface/SiStripDetCablingRcd.h"
18 #include "CalibFormats/SiStripObjects/interface/SiStripDetCabling.h"
19 #include "CondFormats/DataRecord/interface/SiStripCondDataRecords.h"
20 #include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
21 #include "CalibFormats/SiStripObjects/interface/SiStripGain.h"
22 #include "CalibFormats/SiStripObjects/interface/SiStripQuality.h"
23 #include "DataFormats/SiStripCluster/interface/SiStripClusterCollection.h"
24 #include "DataFormats/SiStripCluster/interface/SiStripCluster.h"
25 #include "DQM/SiStripCommon/interface/SiStripFolderOrganizer.h"
26 #include "DQM/SiStripCommon/interface/SiStripHistoId.h"
27 #include "DQM/SiStripMonitorCluster/interface/SiStripMonitorCluster.h"
28 #include "DQMServices/Core/interface/DQMStore.h"
29 #include "DQMServices/Core/interface/MonitorElement.h"
30 #include "DataFormats/SiStripDetId/interface/StripSubdetector.h"
31 #include "DataFormats/SiStripDetId/interface/SiStripSubStructure.h"
32 #include "CalibTracker/SiStripCommon/interface/SiStripDCSStatus.h"
33 #include "DataFormats/SiPixelCluster/interface/SiPixelCluster.h"
34 
35 #include "DPGAnalysis/SiStripTools/interface/APVCyclePhaseCollection.h"
36 #include "DPGAnalysis/SiStripTools/interface/EventWithHistory.h"
37 
38 #include "CommonTools/TriggerUtils/interface/GenericTriggerEventFlag.h"
39 #include "Geometry/Records/interface/TrackerTopologyRcd.h"
40 
41 #include "TMath.h"
42 #include <iostream>
43 
44 //--------------------------------------------------------------------------------------------
45 SiStripMonitorCluster::SiStripMonitorCluster(const edm::ParameterSet& iConfig)
46  : conf_(iConfig), show_mechanical_structure_view(true), show_readout_view(false), show_control_view(false), select_all_detectors(false), reset_each_run(false), m_cacheID_(0)
47  // , genTriggerEventFlag_(new GenericTriggerEventFlag(iConfig, consumesCollector()))
48 {
49 
50  // initialize
51  passBPTXfilter_ = true;
52 
53  // initialize GenericTriggerEventFlag by specific configuration
54  // in this way, one can set specific selections for different MEs
55  genTriggerEventFlagBPTXfilter_ = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("BPTXfilter") , consumesCollector(), *this );
56  genTriggerEventFlagPixelDCSfilter_ = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("PixelDCSfilter"), consumesCollector(), *this );
57  genTriggerEventFlagStripDCSfilter_ = new GenericTriggerEventFlag(iConfig.getParameter<edm::ParameterSet>("StripDCSfilter"), consumesCollector(), *this );
58 
59  firstEvent = -1;
60  eventNb = 0;
61 
62  // Detector Partitions
63  SubDetPhasePartMap["TIB"] = "TI";
64  SubDetPhasePartMap["TID__MINUS"] = "TI";
65  SubDetPhasePartMap["TID__PLUS"] = "TI";
66  SubDetPhasePartMap["TOB"] = "TO";
67  SubDetPhasePartMap["TEC__MINUS"] = "TM";
68  SubDetPhasePartMap["TEC__PLUS"] = "TP";
69 
70  //get on/off option for every cluster from cfi
71  edm::ParameterSet ParametersnClusters = conf_.getParameter<edm::ParameterSet>("TH1nClusters");
72  layerswitchncluson = ParametersnClusters.getParameter<bool>("layerswitchon");
73  moduleswitchncluson = ParametersnClusters.getParameter<bool>("moduleswitchon");
74 
75  edm::ParameterSet ParametersClusterCharge = conf_.getParameter<edm::ParameterSet>("TH1ClusterCharge");
76  layerswitchcluschargeon = ParametersClusterCharge.getParameter<bool>("layerswitchon");
77  moduleswitchcluschargeon = ParametersClusterCharge.getParameter<bool>("moduleswitchon");
78  subdetswitchcluschargeon = ParametersClusterCharge.getParameter<bool>("subdetswitchon");
79 
80  edm::ParameterSet ParametersClusterStoN = conf_.getParameter<edm::ParameterSet>("TH1ClusterStoN");
81  layerswitchclusstonon = ParametersClusterStoN.getParameter<bool>("layerswitchon");
82  moduleswitchclusstonon = ParametersClusterStoN.getParameter<bool>("moduleswitchon");
83 
84  edm::ParameterSet ParametersClusterStoNVsPos = conf_.getParameter<edm::ParameterSet>("TH1ClusterStoNVsPos");
85  layerswitchclusstonVsposon = ParametersClusterStoNVsPos.getParameter<bool>("layerswitchon");
86  moduleswitchclusstonVsposon = ParametersClusterStoNVsPos.getParameter<bool>("moduleswitchon");
87 
88  edm::ParameterSet ParametersClusterPos = conf_.getParameter<edm::ParameterSet>("TH1ClusterPos");
89  layerswitchclusposon = ParametersClusterPos.getParameter<bool>("layerswitchon");
90  moduleswitchclusposon = ParametersClusterPos.getParameter<bool>("moduleswitchon");
91 
92  edm::ParameterSet ParametersClusterDigiPos = conf_.getParameter<edm::ParameterSet>("TH1ClusterDigiPos");
93  layerswitchclusdigiposon = ParametersClusterDigiPos.getParameter<bool>("layerswitchon");
94  moduleswitchclusdigiposon = ParametersClusterDigiPos.getParameter<bool>("moduleswitchon");
95 
96  edm::ParameterSet ParametersClusterNoise = conf_.getParameter<edm::ParameterSet>("TH1ClusterNoise");
97  layerswitchclusnoiseon = ParametersClusterNoise.getParameter<bool>("layerswitchon");
98  moduleswitchclusnoiseon = ParametersClusterNoise.getParameter<bool>("moduleswitchon");
99 
100  edm::ParameterSet ParametersClusterWidth = conf_.getParameter<edm::ParameterSet>("TH1ClusterWidth");
101  layerswitchcluswidthon = ParametersClusterWidth.getParameter<bool>("layerswitchon");
102  moduleswitchcluswidthon = ParametersClusterWidth.getParameter<bool>("moduleswitchon");
103  subdetswitchcluswidthon = ParametersClusterWidth.getParameter<bool>("subdetswitchon");
104 
105  edm::ParameterSet ParametersModuleLocalOccupancy = conf_.getParameter<edm::ParameterSet>("TH1ModuleLocalOccupancy");
106  layerswitchlocaloccupancy = ParametersModuleLocalOccupancy.getParameter<bool>("layerswitchon");
107  moduleswitchlocaloccupancy = ParametersModuleLocalOccupancy.getParameter<bool>("moduleswitchon");
108 
109  edm::ParameterSet ParametersNrOfClusterizedStrips = conf_.getParameter<edm::ParameterSet>("TH1NrOfClusterizedStrips");
110  layerswitchnrclusterizedstrip = ParametersNrOfClusterizedStrips.getParameter<bool>("layerswitchon");
111  moduleswitchnrclusterizedstrip = ParametersNrOfClusterizedStrips.getParameter<bool>("moduleswitchon");
112 
113  edm::ParameterSet ParametersClusterProf = conf_.getParameter<edm::ParameterSet>("TProfNumberOfCluster");
114  layerswitchnumclusterprofon = ParametersClusterProf.getParameter<bool>("layerswitchon");
115 
116  edm::ParameterSet ParametersClusterWidthProf = conf_.getParameter<edm::ParameterSet>("TProfClusterWidth");
117  layerswitchclusterwidthprofon = ParametersClusterWidthProf.getParameter<bool>("layerswitchon");
118 
119  edm::ParameterSet ParametersTotClusterProf = conf_.getParameter<edm::ParameterSet>("TProfTotalNumberOfClusters");
120  subdetswitchtotclusprofon = ParametersTotClusterProf.getParameter<bool>("subdetswitchon");
121 
122  edm::ParameterSet ParametersTotClusterTH1 = conf_.getParameter<edm::ParameterSet>("TH1TotalNumberOfClusters");
123  subdetswitchtotclusth1on = ParametersTotClusterTH1.getParameter<bool>("subdetswitchon");
124 
125  edm::ParameterSet ParametersClusterApvProf = conf_.getParameter<edm::ParameterSet>("TProfClustersApvCycle");
126  subdetswitchapvcycleprofon = ParametersClusterApvProf.getParameter<bool>("subdetswitchon");
127 
128  edm::ParameterSet ParametersClustersApvTH2 = conf_.getParameter<edm::ParameterSet>("TH2ClustersApvCycle");
129  subdetswitchapvcycleth2on = ParametersClustersApvTH2.getParameter<bool>("subdetswitchon");
130 
131  edm::ParameterSet ParametersApvCycleDBxProf2 = conf_.getParameter<edm::ParameterSet>("TProf2ApvCycleVsDBx");
132  subdetswitchapvcycledbxprof2on = ParametersApvCycleDBxProf2.getParameter<bool>("subdetswitchon");
133 
134  edm::ParameterSet ParametersDBxCycleProf = conf_.getParameter<edm::ParameterSet>("TProfClustersVsDBxCycle");
135  subdetswitchdbxcycleprofon = ParametersDBxCycleProf.getParameter<bool>("subdetswitchon");
136 
137  edm::ParameterSet ParametersCStripVsCPix = conf_.getParameter<edm::ParameterSet>("TH2CStripVsCpixel");
138  globalswitchcstripvscpix = ParametersCStripVsCPix.getParameter<bool>("globalswitchon");
139 
140  edm::ParameterSet ParametersMultiplicityRegionsTH1 = conf_.getParameter<edm::ParameterSet>("TH1MultiplicityRegions");
141  globalswitchMultiRegions = ParametersMultiplicityRegionsTH1.getParameter<bool>("globalswitchon");
142 
143  edm::ParameterSet ParametersApvCycleVsDBxGlobalTH2 = conf_.getParameter<edm::ParameterSet>("TH2ApvCycleVsDBxGlobal");
144  globalswitchapvcycledbxth2on = ParametersApvCycleVsDBxGlobalTH2.getParameter<bool>("globalswitchon");
145 
146  edm::ParameterSet ParametersNoiseStrip2ApvCycle = conf_.getParameter<edm::ParameterSet>("TH1StripNoise2ApvCycle");
147  globalswitchstripnoise2apvcycle = ParametersNoiseStrip2ApvCycle.getParameter<bool>("globalswitchon");
148 
149  edm::ParameterSet ParametersNoiseStrip3ApvCycle = conf_.getParameter<edm::ParameterSet>("TH1StripNoise3ApvCycle");
150  globalswitchstripnoise3apvcycle = ParametersNoiseStrip3ApvCycle.getParameter<bool>("globalswitchon");
151 
152  edm::ParameterSet ParametersMainDiagonalPosition = conf_.getParameter<edm::ParameterSet>("TH1MainDiagonalPosition");
153  globalswitchmaindiagonalposition= ParametersMainDiagonalPosition.getParameter<bool>("globalswitchon");
154 
155  edm::ParameterSet ClusterMultiplicityRegions = conf_.getParameter<edm::ParameterSet>("MultiplicityRegions");
156  k0 = ClusterMultiplicityRegions.getParameter<double>("k0");
157  q0 = ClusterMultiplicityRegions.getParameter<double>("q0");
158  dk0 = ClusterMultiplicityRegions.getParameter<double>("dk0");
159  maxClus = ClusterMultiplicityRegions.getParameter<double>("MaxClus");
160  minPix = ClusterMultiplicityRegions.getParameter<double>("MinPix");
161 
162  edm::ParameterSet ParametersNclusVsCycleTimeProf2D = conf_.getParameter<edm::ParameterSet>("NclusVsCycleTimeProf2D");
163  globalswitchnclusvscycletimeprof2don = ParametersNclusVsCycleTimeProf2D.getParameter<bool>("globalswitchon");
164 
165  edm::ParameterSet ParametersClusWidthVsAmpTH2 = conf_.getParameter<edm::ParameterSet>("ClusWidthVsAmpTH2");
166  clusterWidth_vs_amplitude_on = ParametersClusWidthVsAmpTH2.getParameter<bool>("globalswitchon");
167  layer_clusterWidth_vs_amplitude_on = ParametersClusWidthVsAmpTH2.getParameter<bool>("layerswitchon");
168  subdet_clusterWidth_vs_amplitude_on = ParametersClusWidthVsAmpTH2.getParameter<bool>("subdetswitchon");
169  module_clusterWidth_vs_amplitude_on = ParametersClusWidthVsAmpTH2.getParameter<bool>("moduleswitchon");
170 
171  clustertkhistomapon = conf_.getParameter<bool>("TkHistoMap_On");
172  clusterchtkhistomapon = conf_.getParameter<bool>("ClusterChTkHistoMap_On");
173  createTrendMEs = conf_.getParameter<bool>("CreateTrendMEs");
174  trendVsLs_ = conf_.getParameter<bool>("TrendVsLS");
175  Mod_On_ = conf_.getParameter<bool>("Mod_On");
176  ClusterHisto_ = conf_.getParameter<bool>("ClusterHisto");
177 
178  topFolderName_ = conf_.getParameter<std::string>("TopFolderName");
179 
180 
181  // Poducer name of input StripClusterCollection
182  clusterProducerStripToken_ = consumes<edmNew::DetSetVector<SiStripCluster> >(conf_.getParameter<edm::InputTag>("ClusterProducerStrip") );
183  clusterProducerPixToken_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(conf_.getParameter<edm::InputTag>("ClusterProducerPix") );
184  /*
185  clusterProducerStrip_ = conf_.getParameter<edm::InputTag>("ClusterProducerStrip");
186  clusterProducerPix_ = conf_.getParameter<edm::InputTag>("ClusterProducerPix");
187  */
188  // SiStrip Quality Label
189  qualityLabel_ = conf_.getParameter<std::string>("StripQualityLabel");
190  // cluster quality conditions
191  edm::ParameterSet cluster_condition = conf_.getParameter<edm::ParameterSet>("ClusterConditions");
192  applyClusterQuality_ = cluster_condition.getParameter<bool>("On");
193  sToNLowerLimit_ = cluster_condition.getParameter<double>("minStoN");
194  sToNUpperLimit_ = cluster_condition.getParameter<double>("maxStoN");
195  widthLowerLimit_ = cluster_condition.getParameter<double>("minWidth");
196  widthUpperLimit_ = cluster_condition.getParameter<double>("maxWidth");
197 
198  // Event History Producer
199  // historyProducer_ = conf_.getParameter<edm::InputTag>("HistoryProducer");
200  historyProducerToken_ = consumes<EventWithHistory>(conf_.getParameter<edm::InputTag>("HistoryProducer") );
201  // Apv Phase Producer
202  // apvPhaseProducer_ = conf_.getParameter<edm::InputTag>("ApvPhaseProducer");
203  apvPhaseProducerToken_ = consumes<APVCyclePhaseCollection>(conf_.getParameter<edm::InputTag>("ApvPhaseProducer") );
204  // Create DCS Status
205  bool checkDCS = conf_.getParameter<bool>("UseDCSFiltering");
206  if (checkDCS) dcsStatus_ = new SiStripDCSStatus(consumesCollector());
207  else dcsStatus_ = 0;
208 
209 }
210 
211 SiStripMonitorCluster::~SiStripMonitorCluster() {
212  if (dcsStatus_) delete dcsStatus_;
213  if (genTriggerEventFlagBPTXfilter_ ) delete genTriggerEventFlagBPTXfilter_;
214  if (genTriggerEventFlagPixelDCSfilter_) delete genTriggerEventFlagPixelDCSfilter_;
215  if (genTriggerEventFlagStripDCSfilter_) delete genTriggerEventFlagStripDCSfilter_;
216 }
217 
218 
219 //--------------------------------------------------------------------------------------------
220 void SiStripMonitorCluster::dqmBeginRun(const edm::Run& run, const edm::EventSetup &es ) {
221 
222  // Initialize the GenericTriggerEventFlag
223  if ( genTriggerEventFlagBPTXfilter_->on() )
224  genTriggerEventFlagBPTXfilter_->initRun( run, es );
225  if ( genTriggerEventFlagPixelDCSfilter_->on() )
226  genTriggerEventFlagPixelDCSfilter_->initRun( run, es );
227  if ( genTriggerEventFlagStripDCSfilter_->on() )
228  genTriggerEventFlagStripDCSfilter_->initRun( run, es );
229 }
230 
231 //--------------------------------------------------------------------------------------------
232 void SiStripMonitorCluster::createMEs(const edm::EventSetup& es , DQMStore::IBooker & ibooker){
233 
234  if ( show_mechanical_structure_view ){
235 
236  //Retrieve tracker topology from geometry
237  edm::ESHandle<TrackerTopology> tTopoHandle;
238  es.get<TrackerTopologyRcd>().get(tTopoHandle);
239  const TrackerTopology* const tTopo = tTopoHandle.product();
240 
241  // take from eventSetup the SiStripDetCabling object - here will use SiStripDetControl later on
242  es.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
243 
244  // get list of active detectors from SiStripDetCabling
245  std::vector<uint32_t> activeDets;
246  SiStripDetCabling_->addActiveDetectorsRawIds(activeDets);
247 
248  SiStripSubStructure substructure;
249 
250  SiStripFolderOrganizer folder_organizer;
251  folder_organizer.setSiStripFolderName(topFolderName_);
252  folder_organizer.setSiStripFolder();
253 
254 
255  // Create TkHistoMap for Cluster
256  if (clustertkhistomapon) {
257  // std::cout << "[SiStripMonitorCluster::createMEs] topFolderName_: " << topFolderName_ << " ";
258  if ( (topFolderName_ == "SiStrip") or (std::string::npos != topFolderName_.find("HLT")) )
259  tkmapcluster = new TkHistoMap(ibooker , topFolderName_,"TkHMap_NumberOfCluster",0.,true);
260  else tkmapcluster = new TkHistoMap(ibooker , topFolderName_+"/TkHistoMap","TkHMap_NumberOfCluster",0.,false);
261  }
262  if (clusterchtkhistomapon) {
263  if ( (topFolderName_ == "SiStrip") or (std::string::npos != topFolderName_.find("HLT")) )
264  tkmapclusterch = new TkHistoMap(ibooker , topFolderName_,"TkHMap_ClusterCharge",0.,true);
265  else tkmapclusterch = new TkHistoMap(ibooker , topFolderName_+"/TkHistoMap","TkHMap_ClusterCharge",0.,false);
266  }
267 
268  // loop over detectors and book MEs
269  edm::LogInfo("SiStripTkDQM|SiStripMonitorCluster")<<"nr. of activeDets: "<<activeDets.size();
270  for(std::vector<uint32_t>::iterator detid_iterator = activeDets.begin(); detid_iterator!=activeDets.end(); detid_iterator++){
271  uint32_t detid = (*detid_iterator);
272  // remove any eventual zero elements - there should be none, but just in case
273  if(detid == 0) {
274  activeDets.erase(detid_iterator);
275  continue;
276  }
277 
278  if (Mod_On_) {
279  ModMEs mod_single;
280  // set appropriate folder using SiStripFolderOrganizer
281  folder_organizer.setDetectorFolder(detid, tTopo); // pass the detid to this method
282  if (reset_each_run) ResetModuleMEs(detid);
283  createModuleMEs(mod_single, detid , ibooker);
284  // append to ModuleMEsMap
285  ModuleMEsMap.insert( std::make_pair(detid, mod_single));
286  }
287 
288  // Create Layer Level MEs if they are not created already
289  std::pair<std::string,int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo);
290  SiStripHistoId hidmanager;
291  std::string label = hidmanager.getSubdetid(detid,tTopo,false);
292 
293  std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(label);
294  if(iLayerME==LayerMEsMap.end()) {
295 
296  // get detids for the layer
297  int32_t lnumber = det_layer_pair.second;
298  std::vector<uint32_t> layerDetIds;
299  if (det_layer_pair.first == "TIB") {
300  substructure.getTIBDetectors(activeDets,layerDetIds,lnumber,0,0,0);
301  } else if (det_layer_pair.first == "TOB") {
302  substructure.getTOBDetectors(activeDets,layerDetIds,lnumber,0,0);
303  } else if (det_layer_pair.first == "TID" && lnumber > 0) {
304  substructure.getTIDDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0);
305  } else if (det_layer_pair.first == "TID" && lnumber < 0) {
306  substructure.getTIDDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0);
307  } else if (det_layer_pair.first == "TEC" && lnumber > 0) {
308  substructure.getTECDetectors(activeDets,layerDetIds,2,abs(lnumber),0,0,0,0);
309  } else if (det_layer_pair.first == "TEC" && lnumber < 0) {
310  substructure.getTECDetectors(activeDets,layerDetIds,1,abs(lnumber),0,0,0,0);
311  }
312  LayerDetMap[label] = layerDetIds;
313 
314  // book Layer MEs
315  folder_organizer.setLayerFolder(detid,tTopo,det_layer_pair.second);
316  createLayerMEs(label, layerDetIds.size() , ibooker );
317  }
318  // book sub-detector plots
319  auto sdet_pair = folder_organizer.getSubDetFolderAndTag(detid, tTopo);
320  if (SubDetMEsMap.find(sdet_pair.second) == SubDetMEsMap.end()){
321  ibooker.setCurrentFolder(sdet_pair.first);
322 
323  createSubDetMEs(sdet_pair.second , ibooker );
324  }
325  }//end of loop over detectors
326 
327  // Create Global Histogram
328  if (globalswitchapvcycledbxth2on) {
329  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
330  edm::ParameterSet GlobalTH2Parameters = conf_.getParameter<edm::ParameterSet>("TH2ApvCycleVsDBxGlobal");
331  std::string HistoName = "DeltaBx_vs_ApvCycle";
332  GlobalApvCycleDBxTH2 = ibooker.book2D(HistoName,HistoName,
333  GlobalTH2Parameters.getParameter<int32_t>("Nbinsx"),
334  GlobalTH2Parameters.getParameter<double>("xmin"),
335  GlobalTH2Parameters.getParameter<double>("xmax"),
336  GlobalTH2Parameters.getParameter<int32_t>("Nbinsy"),
337  GlobalTH2Parameters.getParameter<double>("ymin"),
338  GlobalTH2Parameters.getParameter<double>("ymax"));
339  GlobalApvCycleDBxTH2->setAxisTitle("APV Cycle (Corrected Absolute Bx % 70)",1);
340  GlobalApvCycleDBxTH2->setAxisTitle("Delta Bunch Crossing Cycle",2);
341 
342  // plot DeltaBX ***************************
343  edm::ParameterSet GlobalTH1Parameters = conf_.getParameter<edm::ParameterSet>("TH1DBxGlobal");
344  HistoName = "DeltaBx";
345  GlobalDBxTH1 = ibooker.book1D(HistoName,HistoName,
346  GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
347  GlobalTH1Parameters.getParameter<double>("xmin"),
348  GlobalTH1Parameters.getParameter<double>("xmax"));
349  GlobalDBxTH1->setAxisTitle("Delta Bunch Crossing",1);
350 
351 
352 
353  // plot DeltaBXCycle ***************************
354  edm::ParameterSet DBxCycle = conf_.getParameter<edm::ParameterSet>("TH1DBxCycleGlobal");
355  HistoName = "DeltaBxCycle";
356  GlobalDBxCycleTH1 = ibooker.book1D(HistoName,HistoName,
357  DBxCycle.getParameter<int32_t>("Nbinsx"),
358  DBxCycle.getParameter<double>("xmin"),
359  DBxCycle.getParameter<double>("xmax"));
360  GlobalDBxCycleTH1->setAxisTitle("Delta Bunch Crossing Cycle",1);
361 
362 
363  }
364 
365  if (globalswitchcstripvscpix) {
366  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
367  edm::ParameterSet GlobalTH2Parameters = conf_.getParameter<edm::ParameterSet>("TH2CStripVsCpixel");
368  std::string HistoName = "StripClusVsPixClus";
369  GlobalCStripVsCpix = ibooker.book2D(HistoName,HistoName,
370  GlobalTH2Parameters.getParameter<int32_t>("Nbinsx"),
371  GlobalTH2Parameters.getParameter<double>("xmin"),
372  GlobalTH2Parameters.getParameter<double>("xmax"),
373  GlobalTH2Parameters.getParameter<int32_t>("Nbinsy"),
374  GlobalTH2Parameters.getParameter<double>("ymin"),
375  GlobalTH2Parameters.getParameter<double>("ymax"));
376  GlobalCStripVsCpix->setAxisTitle("Strip Clusters",1);
377  GlobalCStripVsCpix->setAxisTitle("Pix Clusters",2);
378 
379  // Absolute Bunch Crossing ***********************
380  edm::ParameterSet GlobalTH1Parameters = conf_.getParameter<edm::ParameterSet>("TH1ABx_CSCP");
381  HistoName = "AbsoluteBx_CStripVsCpixel";
382  GlobalABXTH1_CSCP = ibooker.book1D(HistoName,HistoName,
383  GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
384  GlobalTH1Parameters.getParameter<double>("xmin"),
385  GlobalTH1Parameters.getParameter<double>("xmax"));
386  GlobalABXTH1_CSCP->setAxisTitle("Absolute Bunch Crossing",1);
387  }
388 
389  if (globalswitchMultiRegions){
390  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
391  edm::ParameterSet GlobalTH2Parameters = conf_.getParameter<edm::ParameterSet>("TH1MultiplicityRegions");
392  std::string HistoName = "ClusterMultiplicityRegions";
393  PixVsStripMultiplicityRegions = ibooker.book1D(HistoName,HistoName,
394  GlobalTH2Parameters.getParameter<int32_t>("Nbinx"),
395  GlobalTH2Parameters.getParameter<double>("xmin"),
396  GlobalTH2Parameters.getParameter<double>("xmax"));
397  PixVsStripMultiplicityRegions->setAxisTitle("");
398  PixVsStripMultiplicityRegions->setBinLabel(1,"Main Diagonal");
399  PixVsStripMultiplicityRegions->setBinLabel(2,"Strip Noise");
400  PixVsStripMultiplicityRegions->setBinLabel(3,"High Strip Noise");
401  PixVsStripMultiplicityRegions->setBinLabel(4,"Beam Background");
402  PixVsStripMultiplicityRegions->setBinLabel(5,"No Strip Clusters");
403  }
404 
405  if (globalswitchmaindiagonalposition){
406  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
407  edm::ParameterSet GlobalTH1Parameters = conf_.getParameter<edm::ParameterSet>("TH1MainDiagonalPosition");
408  std::string HistoName = "MainDiagonal Position";
409  GlobalMainDiagonalPosition = ibooker.book1D(HistoName,HistoName,
410  GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
411  GlobalTH1Parameters.getParameter<double>("xmin"),
412  GlobalTH1Parameters.getParameter<double>("xmax"));
413  GlobalMainDiagonalPosition->setAxisTitle("atan(NPix/(k*NStrip))");
414 
415  //PLOT MainDiagonalPosition_vs_BX ***************************
416  edm::ParameterSet GlobalTProfParameters = conf_.getParameter<edm::ParameterSet>("TProfMainDiagonalPosition");
417  HistoName = "MainDiagonalPosition_vs_BX";
418  GlobalMainDiagonalPosition_vs_BX = ibooker.bookProfile(HistoName, HistoName,
419  GlobalTProfParameters.getParameter<int32_t>("Nbinsx"),
420  GlobalTProfParameters.getParameter<double>("xmin"),
421  GlobalTProfParameters.getParameter<double>("xmax"),
422  GlobalTProfParameters.getParameter<int32_t>("Nbinsy"),
423  GlobalTProfParameters.getParameter<double>("ymin"),
424  GlobalTProfParameters.getParameter<double>("ymax"));
425 
426  GlobalMainDiagonalPosition_vs_BX->setAxisTitle("Absolute BX", 1);
427  GlobalMainDiagonalPosition_vs_BX->setAxisTitle("tan^{-1}(NPix/k*NStrip))", 2);
428 
429 
430  edm::ParameterSet GlobalTH2Parameters = conf_.getParameter<edm::ParameterSet>("TH2MainDiagonalPosition");
431  HistoName = "TH2MainDiagonalPosition_vs_BX";
432  GlobalTH2MainDiagonalPosition_vs_BX = ibooker.book2D(HistoName,HistoName,
433  GlobalTH2Parameters.getParameter<int32_t>("Nbinsx"),
434  GlobalTH2Parameters.getParameter<double>("xmin"),
435  GlobalTH2Parameters.getParameter<double>("xmax"),
436  GlobalTH2Parameters.getParameter<int32_t>("Nbinsy"),
437  GlobalTH2Parameters.getParameter<double>("ymin"),
438  GlobalTH2Parameters.getParameter<double>("ymax"));
439  GlobalTH2MainDiagonalPosition_vs_BX->setAxisTitle("Absolute BX",1);
440  GlobalTH2MainDiagonalPosition_vs_BX->setAxisTitle("tan^{-1}(NPix/k*NStrip))",2);
441 
442 
443 
444  }
445 
446  // TO BE ADDED !!!
447  /*
448  if ( globalswitchapvcycledbxth2on or globalswitchcstripvscpix or globalswitchMultiRegions or ClusterHisto_ ) {
449  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
450  std::string HistoName = "BPTX rate";
451  BPTXrateTrend = ibooker.bookProfile(HistoName,HistoName, LSBin, LSMin, LSMax, 0, 10000.,"");
452  BPTXrateTrend->getTH1()->SetCanExtend(TH1::kAllAxes);
453  BPTXrateTrend->setAxisTitle("#Lumi section",1);
454  BPTXrateTrend->setAxisTitle("Number of BPTX events per LS",2);
455  }
456  */
457 
458  if (globalswitchstripnoise2apvcycle){
459  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
460  edm::ParameterSet GlobalTH1Parameters = conf_.getParameter<edm::ParameterSet>("TH1StripNoise2ApvCycle");
461  std::string HistoName = "StripNoise_ApvCycle";
462  StripNoise2Cycle = ibooker.book1D(HistoName,HistoName,
463  GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
464  GlobalTH1Parameters.getParameter<double>("xmin"),
465  GlobalTH1Parameters.getParameter<double>("xmax"));
466  StripNoise2Cycle->setAxisTitle("APV Cycle");
467  }
468 
469  if (globalswitchstripnoise3apvcycle){
470  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
471  edm::ParameterSet GlobalTH1Parameters = conf_.getParameter<edm::ParameterSet>("TH1StripNoise3ApvCycle");
472  std::string HistoName = "HighStripNoise_ApvCycle";
473  StripNoise3Cycle = ibooker.book1D(HistoName,HistoName,
474  GlobalTH1Parameters.getParameter<int32_t>("Nbinsx"),
475  GlobalTH1Parameters.getParameter<double>("xmin"),
476  GlobalTH1Parameters.getParameter<double>("xmax"));
477  StripNoise3Cycle->setAxisTitle("APV Cycle");
478  }
479 
480  if ( globalswitchnclusvscycletimeprof2don ) {
481  const char* HistoName = "StripClusVsBXandOrbit";
482  const char* HistoTitle = "Strip cluster multiplicity vs BX mod(70) and Orbit;Event 1 BX mod(70);time [Orb#]";
483  edm::ParameterSet ParametersNclusVsCycleTimeProf2D = conf_.getParameter<edm::ParameterSet>("NclusVsCycleTimeProf2D");
484  NclusVsCycleTimeProf2D = ibooker.bookProfile2D ( HistoName , HistoTitle ,
485  ParametersNclusVsCycleTimeProf2D.getParameter<int32_t>("Nbins"),
486  ParametersNclusVsCycleTimeProf2D.getParameter<double>("xmin"),
487  ParametersNclusVsCycleTimeProf2D.getParameter<double>("xmax"),
488  ParametersNclusVsCycleTimeProf2D.getParameter<int32_t>("Nbinsy"),
489  ParametersNclusVsCycleTimeProf2D.getParameter<double>("ymin"),
490  ParametersNclusVsCycleTimeProf2D.getParameter<double>("ymax"),
491  0 , 0 );
492  if (NclusVsCycleTimeProf2D->kind() == MonitorElement::DQM_KIND_TPROFILE2D) NclusVsCycleTimeProf2D->getTH1()->SetCanExtend(TH1::kAllAxes);
493  }
494  if ( clusterWidth_vs_amplitude_on ) {
495  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
496  edm::ParameterSet ParametersClusWidthVsAmpTH2 = conf_.getParameter<edm::ParameterSet>("ClusWidthVsAmpTH2");
497  const char* HistoName = "ClusterWidths_vs_Amplitudes";
498  const char* HistoTitle = "Cluster widths vs amplitudes;Amplitudes (integrated ADC counts);Cluster widths";
499  ClusWidthVsAmpTH2 = ibooker.book2D(HistoName,HistoTitle,
500  ParametersClusWidthVsAmpTH2.getParameter<int32_t>("Nbinsx"),
501  ParametersClusWidthVsAmpTH2.getParameter<double>("xmin"),
502  ParametersClusWidthVsAmpTH2.getParameter<double>("xmax"),
503  ParametersClusWidthVsAmpTH2.getParameter<int32_t>("Nbinsy"),
504  ParametersClusWidthVsAmpTH2.getParameter<double>("ymin"),
505  ParametersClusWidthVsAmpTH2.getParameter<double>("ymax"));
506  }
507 
508  if (ClusterHisto_){
509  ibooker.setCurrentFolder(topFolderName_+"/MechanicalView/");
510  edm::ParameterSet PixelCluster = conf_.getParameter<edm::ParameterSet>("TH1NClusPx");
511  std::string HistoName = "NumberOfClustersInPixel";
512  NumberOfPixelClus = ibooker.book1D(HistoName, HistoName,
513  PixelCluster.getParameter<int32_t>("Nbinsx"),
514  PixelCluster.getParameter<double>("xmin"),
515  PixelCluster.getParameter<double>("xmax"));
516  NumberOfPixelClus->setAxisTitle("# of Clusters in Pixel", 1);
517  NumberOfPixelClus->setAxisTitle("Number of Events", 2);
518  //
519  edm::ParameterSet StripCluster = conf_.getParameter<edm::ParameterSet>("TH1NClusStrip");
520  HistoName = "NumberOfClustersInStrip";
521  NumberOfStripClus = ibooker.book1D(HistoName, HistoName,
522  StripCluster.getParameter<int32_t>("Nbinsx"),
523  StripCluster.getParameter<double>("xmin"),
524  StripCluster.getParameter<double>("xmax"));
525  NumberOfStripClus->setAxisTitle("# of Clusters in Strip", 1);
526  NumberOfStripClus->setAxisTitle("Number of Events", 2);
527 
528  // NumberOfClustersinStrip vs BX PLOT ****************************
529  edm::ParameterSet StripClusterBX = conf_.getParameter<edm::ParameterSet>("TProfNClusStrip");
530  HistoName = "NumberOfClustersInStrip_vs_BX";
531  NumberOfStripClus_vs_BX = ibooker.bookProfile(HistoName, HistoName,
532  StripClusterBX.getParameter<int32_t>("Nbinsx"),
533  StripClusterBX.getParameter<double>("xmin"),
534  StripClusterBX.getParameter<double>("xmax"),
535  StripClusterBX.getParameter<int32_t>("Nbinsy"),
536  StripClusterBX.getParameter<double>("ymin"),
537  StripClusterBX.getParameter<double>("ymax"));
538 
539  NumberOfStripClus_vs_BX->setAxisTitle("Absolute BX", 1);
540  NumberOfStripClus_vs_BX->setAxisTitle("# of Clusters in Strip", 2);
541 
542  // NumberOfClustersinStrip vs BX PLOT ****************************
543  edm::ParameterSet PixelClusterBX = conf_.getParameter<edm::ParameterSet>("TProfNClusPixel");
544  HistoName = "NumberOfClustersInPixel_vs_BX";
545  NumberOfPixelClus_vs_BX = ibooker.bookProfile(HistoName, HistoName,
546  PixelClusterBX.getParameter<int32_t>("Nbinsx"),
547  PixelClusterBX.getParameter<double>("xmin"),
548  PixelClusterBX.getParameter<double>("xmax"),
549  PixelClusterBX.getParameter<int32_t>("Nbinsy"),
550  PixelClusterBX.getParameter<double>("ymin"),
551  PixelClusterBX.getParameter<double>("ymax"));
552 
553  NumberOfPixelClus_vs_BX->setAxisTitle("Absolute BX", 1);
554  NumberOfPixelClus_vs_BX->setAxisTitle("# of Clusters in Pixel", 2);
555 
556  }
557 
558  }//end of if
559 }//end of method
560 
561 //--------------------------------------------------------------------------------------------
562 
563 void SiStripMonitorCluster::bookHistograms(DQMStore::IBooker & ibooker, const edm::Run& run, const edm::EventSetup& es)
564 {
565  if (show_mechanical_structure_view) {
566  unsigned long long cacheID = es.get<SiStripDetCablingRcd>().cacheIdentifier();
567  if (m_cacheID_ != cacheID) {
568  m_cacheID_ = cacheID;
569  edm::LogInfo("SiStripMonitorCluster") <<"SiStripMonitorCluster::bookHistograms: "
570  << " Creating MEs for new Cabling ";
571 
572  createMEs(es , ibooker);
573  }
574  } else if (reset_each_run) {
575  edm::LogInfo("SiStripMonitorCluster") <<"SiStripMonitorCluster::bookHistograms: "
576  << " Resetting MEs ";
577  for (std::map<uint32_t, ModMEs >::const_iterator idet = ModuleMEsMap.begin() ; idet!=ModuleMEsMap.end() ; idet++) {
578  ResetModuleMEs(idet->first);
579  }
580  }
581 
582 }
583 
584 //--------------------------------------------------------------------------------------------
585 void SiStripMonitorCluster::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup)
586 {
587  //Retrieve tracker topology from geometry
588  edm::ESHandle<TrackerTopology> tTopoHandle;
589  iSetup.get<TrackerTopologyRcd>().get(tTopoHandle);
590  const TrackerTopology* const tTopo = tTopoHandle.product();
591 
592  // Filter out events if Trigger Filtering is requested
593  passBPTXfilter_ = ( iEvent.isRealData() and genTriggerEventFlagBPTXfilter_->on() ) ? genTriggerEventFlagBPTXfilter_->accept( iEvent, iSetup) : true;
594  passPixelDCSfilter_ = ( iEvent.isRealData() and genTriggerEventFlagPixelDCSfilter_->on() ) ? genTriggerEventFlagPixelDCSfilter_->accept( iEvent, iSetup) : true;
595  passStripDCSfilter_ = ( iEvent.isRealData() and genTriggerEventFlagStripDCSfilter_->on() ) ? genTriggerEventFlagStripDCSfilter_->accept( iEvent, iSetup) : true;
596  // std::cout << "passBPTXfilter_ ? " << passBPTXfilter_ << std::endl;
597 
598  // Filter out events if DCS Event if requested
599  if (dcsStatus_ && !dcsStatus_->getStatus(iEvent,iSetup)) return;
600 
601  runNb = iEvent.id().run();
602  eventNb++;
603  trendVar = trendVsLs_ ? iEvent.orbitNumber()/262144.0 : iEvent.orbitNumber()/11223.0; // lumisection : seconds
604 
605  int NPixClusters=0, NStripClusters=0, MultiplicityRegion=0;
606  bool isPixValid=false;
607 
608  edm::ESHandle<SiStripNoises> noiseHandle;
609  iSetup.get<SiStripNoisesRcd>().get(noiseHandle);
610 
611  edm::ESHandle<SiStripGain> gainHandle;
612  iSetup.get<SiStripGainRcd>().get(gainHandle);
613 
614  edm::ESHandle<SiStripQuality> qualityHandle;
615  iSetup.get<SiStripQualityRcd>().get(qualityLabel_,qualityHandle);
616 
617  iSetup.get<SiStripDetCablingRcd>().get(SiStripDetCabling_);
618 
619  // get collection of DetSetVector of clusters from Event
620  edm::Handle< edmNew::DetSetVector<SiStripCluster> > cluster_detsetvektor;
621  iEvent.getByToken(clusterProducerStripToken_, cluster_detsetvektor);
622 
623  //get pixel clusters
624  edm::Handle< edmNew::DetSetVector<SiPixelCluster> > cluster_detsetvektor_pix;
625  iEvent.getByToken(clusterProducerPixToken_, cluster_detsetvektor_pix);
626 
627  if (!cluster_detsetvektor.isValid()) return;
628 
629  const edmNew::DetSetVector<SiStripCluster> * StrC= cluster_detsetvektor.product();
630  NStripClusters= StrC->data().size();
631 
632  if (cluster_detsetvektor_pix.isValid()){
633  const edmNew::DetSetVector<SiPixelCluster> * PixC= cluster_detsetvektor_pix.product();
634  NPixClusters= PixC->data().size();
635  isPixValid=true;
636  MultiplicityRegion=FindRegion(NStripClusters,NPixClusters);
637 
638  if ( passBPTXfilter_ and passPixelDCSfilter_ and passStripDCSfilter_ ) {
639  if (globalswitchcstripvscpix) GlobalCStripVsCpix->Fill(NStripClusters,NPixClusters);
640  if (globalswitchmaindiagonalposition && NStripClusters > 0) GlobalMainDiagonalPosition->Fill(atan(NPixClusters/(k0*NStripClusters)));
641 
642  if (globalswitchMultiRegions) PixVsStripMultiplicityRegions->Fill(MultiplicityRegion);
643  }
644 
645  if (ClusterHisto_){
646  if ( passBPTXfilter_ and passPixelDCSfilter_ )
647  NumberOfPixelClus->Fill(NPixClusters);
648  if ( passBPTXfilter_ and passStripDCSfilter_ )
649  NumberOfStripClus->Fill(NStripClusters);
650  }
651  }
652  // initialise # of clusters to zero
653  for (std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.begin();
654  iSubdet != SubDetMEsMap.end(); iSubdet++) {
655  iSubdet->second.totNClusters = 0;
656  }
657 
658  SiStripFolderOrganizer folder_organizer;
659  bool found_layer_me = false;
660  for (std::map<std::string, std::vector< uint32_t > >::const_iterator iterLayer = LayerDetMap.begin();
661  iterLayer != LayerDetMap.end(); iterLayer++) {
662 
663  std::string layer_label = iterLayer->first;
664 
665  int ncluster_layer = 0;
666  std::map<std::string, LayerMEs>::iterator iLayerME = LayerMEsMap.find(layer_label);
667 
668  //get Layer MEs
669  LayerMEs layer_single;
670  if(iLayerME != LayerMEsMap.end()) {
671  layer_single = iLayerME->second;
672  found_layer_me = true;
673  }
674 
675  bool found_module_me = false;
676  uint16_t iDet = 0;
677  std::string subdet_label = "";
678  // loop over all modules in the layer
679  for (std::vector< uint32_t >::const_iterator iterDets = iterLayer->second.begin() ;
680  iterDets != iterLayer->second.end() ; iterDets++) {
681  iDet++;
682  // detid and type of ME
683  uint32_t detid = (*iterDets);
684 
685  // Get SubDet label once
686  if (subdet_label.size() == 0) subdet_label = folder_organizer.getSubDetFolderAndTag(detid, tTopo).second;
687 
688  // DetId and corresponding set of MEs
689  ModMEs mod_single;
690  if (Mod_On_) {
691  std::map<uint32_t, ModMEs >::iterator imodME = ModuleMEsMap.find(detid);
692  if (imodME != ModuleMEsMap.end()) {
693  mod_single = imodME->second;
694  found_module_me = true;
695  }
696  } else found_module_me = false;
697 
698  edmNew::DetSetVector<SiStripCluster>::const_iterator isearch = cluster_detsetvektor->find(detid); // search clusters of detid
699 
700  if(isearch==cluster_detsetvektor->end()){
701  if(found_module_me && moduleswitchncluson && (mod_single.NumberOfClusters)){
702  (mod_single.NumberOfClusters)->Fill(0.); // no clusters for this detector module,fill histogram with 0
703  }
704  if(clustertkhistomapon) tkmapcluster->fill(detid,0.);
705  if(clusterchtkhistomapon) tkmapclusterch->fill(detid,0.);
706  if (found_layer_me && layerswitchnumclusterprofon) layer_single.LayerNumberOfClusterProfile->Fill(iDet, 0.0);
707  continue; // no clusters for this detid => jump to next step of loop
708  }
709 
710  //cluster_detset is a structure, cluster_detset.data is a std::vector<SiStripCluster>, cluster_detset.id is uint32_t
711  // edmNew::DetSet<SiStripCluster> cluster_detset = (*cluster_detsetvektor)[detid]; // the statement above makes sure there exists an element with 'detid'
712  edmNew::DetSet<SiStripCluster> cluster_detset = (*isearch);
713 
714  // Filling TkHistoMap with number of clusters for each module
715  if(clustertkhistomapon) {
716  tkmapcluster->fill(detid,static_cast<float>(cluster_detset.size()));
717  }
718 
719  if(moduleswitchncluson && found_module_me && (mod_single.NumberOfClusters != NULL)){ // nr. of clusters per module
720  (mod_single.NumberOfClusters)->Fill(static_cast<float>(cluster_detset.size()));
721  }
722 
723  if (found_layer_me && layerswitchnumclusterprofon)
724  layer_single.LayerNumberOfClusterProfile->Fill(iDet, static_cast<float>(cluster_detset.size()));
725  ncluster_layer += cluster_detset.size();
726 
727  short total_clusterized_strips = 0;
728 
729  SiStripNoises::Range detNoiseRange = noiseHandle->getRange(detid);
730  SiStripApvGain::Range detGainRange = gainHandle->getRange(detid);
731  SiStripQuality::Range qualityRange = qualityHandle->getRange(detid);
732 
733  for(edmNew::DetSet<SiStripCluster>::const_iterator clusterIter = cluster_detset.begin(); clusterIter!= cluster_detset.end(); clusterIter++){
734 
735  const auto & ampls = clusterIter->amplitudes();
736  // cluster position
737  float cluster_position = clusterIter->barycenter();
738  // start defined as nr. of first strip beloning to the cluster
739  short cluster_start = clusterIter->firstStrip();
740  // width defined as nr. of strips that belong to cluster
741  short cluster_width = ampls.size();
742  // add nr of strips of this cluster to total nr. of clusterized strips
743  total_clusterized_strips = total_clusterized_strips + cluster_width;
744 
745  if(clusterchtkhistomapon) tkmapclusterch->fill(detid,static_cast<float>(clusterIter->charge()));
746 
747  // cluster signal and noise from the amplitudes
748  float cluster_signal = 0.0;
749  float cluster_noise = 0.0;
750  int nrnonzeroamplitudes = 0;
751  float noise2 = 0.0;
752  float noise = 0.0;
753  for(uint iamp=0; iamp<ampls.size(); iamp++){
754  if(ampls[iamp]>0){ // nonzero amplitude
755  cluster_signal += ampls[iamp];
756  if(!qualityHandle->IsStripBad(qualityRange, clusterIter->firstStrip()+iamp)){
757  noise = noiseHandle->getNoise(clusterIter->firstStrip()+iamp,detNoiseRange)/gainHandle->getStripGain(clusterIter->firstStrip()+iamp, detGainRange);
758  }
759  noise2 += noise*noise;
760  nrnonzeroamplitudes++;
761  }
762  } // End loop over cluster amplitude
763 
764  if (nrnonzeroamplitudes > 0) cluster_noise = sqrt(noise2/nrnonzeroamplitudes);
765 
766  if( applyClusterQuality_ &&
767  (cluster_signal/cluster_noise < sToNLowerLimit_ ||
768  cluster_signal/cluster_noise > sToNUpperLimit_ ||
769  cluster_width < widthLowerLimit_ ||
770  cluster_width > widthUpperLimit_) ) continue;
771 
772  ClusterProperties cluster_properties;
773  cluster_properties.charge = cluster_signal;
774  cluster_properties.position = cluster_position;
775  cluster_properties.start = cluster_start;
776  cluster_properties.width = cluster_width;
777  cluster_properties.noise = cluster_noise;
778 
779  // Fill Module Level MEs
780  if (found_module_me) fillModuleMEs(mod_single, cluster_properties);
781 
782  // Fill Layer Level MEs
783  if (found_layer_me) {
784  fillLayerMEs(layer_single, cluster_properties);
785  if (layerswitchclusterwidthprofon)
786  layer_single.LayerClusterWidthProfile->Fill(iDet, cluster_width);
787  }
788 
789  if (subdetswitchcluschargeon || subdetswitchcluswidthon)
790  {
791  std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(subdet_label);
792  if(iSubdet != SubDetMEsMap.end())
793  {
794  if (subdetswitchcluschargeon) iSubdet->second.SubDetClusterChargeTH1->Fill(cluster_signal);
795  if (subdetswitchcluswidthon) iSubdet->second.SubDetClusterWidthTH1->Fill(cluster_width);
796  }
797  }
798 
799 
800  if (subdet_clusterWidth_vs_amplitude_on){
801  std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(subdet_label);
802  if(iSubdet != SubDetMEsMap.end()) iSubdet->second.SubDetClusWidthVsAmpTH2->Fill(cluster_signal, cluster_width);
803  }
804 
805  if ( clusterWidth_vs_amplitude_on ) {
806  ClusWidthVsAmpTH2->Fill(cluster_signal, cluster_width);
807  }
808 
809 
810  if (subdetswitchtotclusprofon) {
811  std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(subdet_label);
812  std::pair<std::string,int32_t> det_layer_pair = folder_organizer.GetSubDetAndLayer(detid, tTopo);
813  if(iSubdet != SubDetMEsMap.end()) iSubdet->second.SubDetNumberOfClusterPerLayerTrend->Fill(trendVar,abs(det_layer_pair.second));
814  }
815 
816  if (subdetswitchtotclusprofon && (subdet_label.find("TID")!=std::string::npos || subdet_label.find("TEC")!=std::string::npos)){
817  std::pair<std::string,int32_t> det_ring_pair = folder_organizer.GetSubDetAndLayer(detid,tTopo,true);
818  fillME(layer_single.LayerNumberOfClusterPerRingTrend, trendVar, abs(det_ring_pair.second));
819  }
820 
821  } // end loop over clusters
822 
823  short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs
824  float local_occupancy = static_cast<float>(total_clusterized_strips)/static_cast<float>(total_nr_strips);
825  if (found_module_me) {
826  if(moduleswitchnrclusterizedstrip && mod_single.NrOfClusterizedStrips ){ // nr of clusterized strips
827  mod_single.NrOfClusterizedStrips->Fill(static_cast<float>(total_clusterized_strips));
828  }
829 
830  if(moduleswitchlocaloccupancy && mod_single.ModuleLocalOccupancy ){ // Occupancy
831  mod_single.ModuleLocalOccupancy->Fill(local_occupancy);
832  }
833  }
834  if (layerswitchlocaloccupancy && found_layer_me && layer_single.LayerLocalOccupancy) {
835  fillME(layer_single.LayerLocalOccupancy,local_occupancy);
836  if (createTrendMEs)
837  fillME(layer_single.LayerLocalOccupancyTrend,trendVar,local_occupancy);
838  }
839  }
840 
841  if (subdetswitchtotclusprofon)
842  fillME(layer_single.LayerNumberOfClusterTrend,trendVar,ncluster_layer);
843 
844  std::map<std::string, SubDetMEs>::iterator iSubdet = SubDetMEsMap.find(subdet_label);
845  if(iSubdet != SubDetMEsMap.end()) iSubdet->second.totNClusters += ncluster_layer;
846  }
847 
848  // EventHistory
849  edm::Handle<EventWithHistory> event_history;
850  iEvent.getByToken(historyProducerToken_,event_history);
851 
852  // Phase of APV
853  edm::Handle<APVCyclePhaseCollection> apv_phase_collection;
854  iEvent.getByToken(apvPhaseProducerToken_,apv_phase_collection);
855 
856  if (event_history.isValid()
857  && !event_history.failedToGet()
858  && apv_phase_collection.isValid()
859  && !apv_phase_collection.failedToGet()) {
860 
861 
862  long long dbx = event_history->deltaBX();
863  long long tbx = event_history->absoluteBX();
864 
865  bool global_histo_filled = false;
866  bool MultiplicityRegion_Vs_APVcycle_filled=false;
867 
868  // plot n 2
869  if ( passBPTXfilter_ and passPixelDCSfilter_ and passStripDCSfilter_ ) {
870  if (globalswitchcstripvscpix) GlobalABXTH1_CSCP->Fill(tbx%3564);
871  }
872  // plot n 2
873 
874  for (std::map<std::string, SubDetMEs>::iterator it = SubDetMEsMap.begin();
875  it != SubDetMEsMap.end(); it++) {
876  std::string sdet = it->first;
877  //std::string sdet = sdet_tag.substr(0,sdet_tag.find_first_of("_"));
878  SubDetMEs sdetmes = it->second;
879 
880  int the_phase = APVCyclePhaseCollection::invalid;
881  long long tbx_corr = tbx;
882 
883 
884  if (SubDetPhasePartMap.find(sdet) != SubDetPhasePartMap.end()) the_phase = apv_phase_collection->getPhase(SubDetPhasePartMap[sdet]);
885  if(the_phase==APVCyclePhaseCollection::nopartition ||
886  the_phase==APVCyclePhaseCollection::multiphase ||
887  the_phase==APVCyclePhaseCollection::invalid) {
888  the_phase=30;
889  //std::cout << " subdet " << it->first << " not valid" << " MR " << MultiplicityRegion <<std::endl;
890  }
891  tbx_corr -= the_phase;
892 
893  long long dbxincycle = event_history->deltaBXinCycle(the_phase);
894  //std::cout<<"dbx in cycle: " << dbxincycle << std::endl;
895  if (globalswitchapvcycledbxth2on && !global_histo_filled) {
896  GlobalApvCycleDBxTH2->Fill(tbx_corr%70,dbx);
897  GlobalDBxTH1->Fill(dbx);
898  GlobalDBxCycleTH1->Fill(dbxincycle);
899  global_histo_filled = true;
900  }
901 
902 
903  // Fill MainDiagonalPosition plots ***************
904  if (cluster_detsetvektor_pix.isValid()){
905  if (ClusterHisto_){
906  if ( passBPTXfilter_ and passStripDCSfilter_ ){
907  NumberOfStripClus_vs_BX->Fill(tbx%3564,NStripClusters);
908  if(passPixelDCSfilter_ ){
909  NumberOfPixelClus_vs_BX->Fill(tbx%3564,NPixClusters);
910  if (globalswitchmaindiagonalposition && NStripClusters > 0){
911  GlobalMainDiagonalPosition_vs_BX->Fill(tbx%3564,atan(NPixClusters/(k0*NStripClusters)));
912  GlobalTH2MainDiagonalPosition_vs_BX->Fill(tbx%3564,atan(NPixClusters/(k0*NStripClusters)));
913  }
914  }
915 
916  }
917  }
918  }
919 
920  if (isPixValid && !MultiplicityRegion_Vs_APVcycle_filled){
921  if (globalswitchstripnoise2apvcycle && MultiplicityRegion==2) {StripNoise2Cycle->Fill(tbx_corr%70);}
922  if (globalswitchstripnoise3apvcycle && MultiplicityRegion==3) {StripNoise3Cycle->Fill(tbx_corr%70);}
923  MultiplicityRegion_Vs_APVcycle_filled=true;
924  }
925 
926  if (subdetswitchtotclusth1on)
927  sdetmes.SubDetTotClusterTH1->Fill(sdetmes.totNClusters);
928  if (subdetswitchtotclusprofon)
929  sdetmes.SubDetTotClusterProf->Fill(trendVar,sdetmes.totNClusters);
930  if (subdetswitchapvcycleprofon)
931  sdetmes.SubDetClusterApvProf->Fill(tbx_corr%70,sdetmes.totNClusters);
932  if (subdetswitchapvcycleth2on)
933  sdetmes.SubDetClusterApvTH2->Fill(tbx_corr%70,sdetmes.totNClusters);
934  if (subdetswitchdbxcycleprofon) {
935  sdetmes.SubDetClusterDBxCycleProf->Fill(dbxincycle,sdetmes.totNClusters);
936  }
937  if (subdetswitchapvcycledbxprof2on)
938  sdetmes.SubDetApvDBxProf2->Fill(tbx_corr%70,dbx,sdetmes.totNClusters);
939  }
940 
941  if ( globalswitchnclusvscycletimeprof2don )
942  {
943  long long tbx_corr = tbx;
944  int the_phase = apv_phase_collection->getPhase("All");
945 
946  if( the_phase == APVCyclePhaseCollection::nopartition ||
947  the_phase == APVCyclePhaseCollection::multiphase ||
948  the_phase == APVCyclePhaseCollection::invalid )
949  the_phase=30;
950 
951  tbx_corr -= the_phase;
952 
953  NclusVsCycleTimeProf2D->Fill( tbx_corr%70 , (int)event_history->_orbit , NStripClusters );
954  }
955  }
956 }
957 //
958 // -- Reset MEs
959 //------------------------------------------------------------------------------
960 void SiStripMonitorCluster::ResetModuleMEs(uint32_t idet){
961  std::map<uint32_t, ModMEs >::iterator pos = ModuleMEsMap.find(idet);
962  ModMEs mod_me = pos->second;
963 
964  if (moduleswitchncluson) mod_me.NumberOfClusters->Reset();
965  if (moduleswitchclusposon) mod_me.ClusterPosition->Reset();
966  if (moduleswitchclusdigiposon) mod_me.ClusterDigiPosition->Reset();
967  if (moduleswitchclusstonVsposon) mod_me.ClusterSignalOverNoiseVsPos->Reset();
968  if (moduleswitchcluswidthon) mod_me.ClusterWidth->Reset();
969  if (moduleswitchcluschargeon) mod_me.ClusterCharge->Reset();
970  if (moduleswitchclusnoiseon) mod_me.ClusterNoise->Reset();
971  if (moduleswitchclusstonon) mod_me.ClusterSignalOverNoise->Reset();
972  if (moduleswitchlocaloccupancy) mod_me.ModuleLocalOccupancy->Reset();
973  if (moduleswitchnrclusterizedstrip) mod_me.NrOfClusterizedStrips->Reset();
974  if (module_clusterWidth_vs_amplitude_on) mod_me.Module_ClusWidthVsAmpTH2->Reset();
975 }
976 //
977 // -- Create Module Level MEs
978 //
979 void SiStripMonitorCluster::createModuleMEs(ModMEs& mod_single, uint32_t detid , DQMStore::IBooker & ibooker) {
980 
981  // use SistripHistoId for producing histogram id (and title)
982  SiStripHistoId hidmanager;
983  std::string hid;
984 
985  //nr. of clusters per module
986  if(moduleswitchncluson) {
987  hid = hidmanager.createHistoId("NumberOfClusters","det",detid);
988  mod_single.NumberOfClusters = bookME1D("TH1nClusters", hid.c_str() , ibooker);
989  ibooker.tag(mod_single.NumberOfClusters, detid);
990  mod_single.NumberOfClusters->setAxisTitle("number of clusters in one detector module");
991  mod_single.NumberOfClusters->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
992  }
993 
994  //ClusterPosition
995  if(moduleswitchclusposon) {
996  short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs
997  hid = hidmanager.createHistoId("ClusterPosition","det",detid);
998  mod_single.ClusterPosition = ibooker.book1D(hid, hid, total_nr_strips, 0.5, total_nr_strips+0.5);
999  ibooker.tag(mod_single.ClusterPosition, detid);
1000  mod_single.ClusterPosition->setAxisTitle("cluster position [strip number +0.5]");
1001  }
1002 
1003  //ClusterDigiPosition
1004  if(moduleswitchclusdigiposon) {
1005  short total_nr_strips = SiStripDetCabling_->nApvPairs(detid) * 2 * 128; // get correct # of avp pairs
1006  hid = hidmanager.createHistoId("ClusterDigiPosition","det",detid);
1007  mod_single.ClusterDigiPosition = ibooker.book1D(hid, hid, total_nr_strips, 0.5, total_nr_strips+0.5);
1008  ibooker.tag(mod_single.ClusterDigiPosition, detid);
1009  mod_single.ClusterDigiPosition->setAxisTitle("digi in cluster position [strip number +0.5]");
1010  }
1011 
1012  //ClusterWidth
1013  if(moduleswitchcluswidthon) {
1014  hid = hidmanager.createHistoId("ClusterWidth","det",detid);
1015  mod_single.ClusterWidth = bookME1D("TH1ClusterWidth", hid.c_str() , ibooker);
1016  ibooker.tag(mod_single.ClusterWidth, detid);
1017  mod_single.ClusterWidth->setAxisTitle("cluster width [nr strips]");
1018  }
1019 
1020  //ClusterCharge
1021  if(moduleswitchcluschargeon) {
1022  hid = hidmanager.createHistoId("ClusterCharge","det",detid);
1023  mod_single.ClusterCharge = bookME1D("TH1ClusterCharge", hid.c_str() , ibooker );
1024  ibooker.tag(mod_single.ClusterCharge, detid);
1025  mod_single.ClusterCharge->setAxisTitle("cluster charge [ADC]");
1026  }
1027 
1028  //ClusterNoise
1029  if(moduleswitchclusnoiseon) {
1030  hid = hidmanager.createHistoId("ClusterNoise","det",detid);
1031  mod_single.ClusterNoise = bookME1D("TH1ClusterNoise", hid.c_str() , ibooker );
1032  ibooker.tag(mod_single.ClusterNoise, detid);
1033  mod_single.ClusterNoise->setAxisTitle("cluster noise");
1034  }
1035 
1036  //ClusterSignalOverNoise
1037  if(moduleswitchclusstonon) {
1038  hid = hidmanager.createHistoId("ClusterSignalOverNoise","det",detid);
1039  mod_single.ClusterSignalOverNoise = bookME1D("TH1ClusterStoN", hid.c_str() , ibooker);
1040  ibooker.tag(mod_single.ClusterSignalOverNoise, detid);
1041  mod_single.ClusterSignalOverNoise->setAxisTitle("ratio of signal to noise for each cluster");
1042  }
1043 
1044  //ClusterSignalOverNoiseVsPos
1045  if(moduleswitchclusstonVsposon) {
1046  hid = hidmanager.createHistoId("ClusterSignalOverNoiseVsPos","det",detid);
1047  Parameters = conf_.getParameter<edm::ParameterSet>("TH1ClusterStoNVsPos");
1048  mod_single.ClusterSignalOverNoiseVsPos= ibooker.bookProfile(hid.c_str(),hid.c_str(),
1049  Parameters.getParameter<int32_t>("Nbinx"),
1050  Parameters.getParameter<double>("xmin"),
1051  Parameters.getParameter<double>("xmax"),
1052  Parameters.getParameter<int32_t>("Nbiny"),
1053  Parameters.getParameter<double>("ymin"),
1054  Parameters.getParameter<double>("ymax")
1055  );
1056  ibooker.tag(mod_single.ClusterSignalOverNoiseVsPos, detid);
1057  mod_single.ClusterSignalOverNoiseVsPos->setAxisTitle("pos");
1058  }
1059 
1060  //ModuleLocalOccupancy
1061  if (moduleswitchlocaloccupancy) {
1062  hid = hidmanager.createHistoId("ClusterLocalOccupancy","det",detid);
1063  mod_single.ModuleLocalOccupancy = bookME1D("TH1ModuleLocalOccupancy", hid.c_str() , ibooker);
1064  ibooker.tag(mod_single.ModuleLocalOccupancy, detid);
1065  mod_single.ModuleLocalOccupancy->setAxisTitle("module local occupancy [% of clusterized strips]");
1066  }
1067 
1068  //NrOfClusterizedStrips
1069  if (moduleswitchnrclusterizedstrip) {
1070  hid = hidmanager.createHistoId("NrOfClusterizedStrips","det",detid);
1071  mod_single.NrOfClusterizedStrips = bookME1D("TH1NrOfClusterizedStrips", hid.c_str() , ibooker );
1072  ibooker.tag(mod_single.NrOfClusterizedStrips, detid);
1073  mod_single.NrOfClusterizedStrips->setAxisTitle("number of clusterized strips");
1074  }
1075 
1076  if (module_clusterWidth_vs_amplitude_on){
1077  hid = hidmanager.createHistoId("ClusterWidths_vs_Amplitudes","det",detid);
1078  Parameters = conf_.getParameter<edm::ParameterSet>("ClusWidthVsAmpTH2");
1079  int32_t Nbinsx = Parameters.getParameter<int32_t>("Nbinsx");
1080  Nbinsx = Nbinsx/2; //Without this "rebinning" the job is killed on lxplus. We think it's due to the high memory needed to create all those histograms.
1081  mod_single.Module_ClusWidthVsAmpTH2 = ibooker.book2D(hid.c_str() , hid.c_str(),
1082  Nbinsx,
1083  Parameters.getParameter<double>("xmin"),
1084  Parameters.getParameter<double>("xmax"),
1085  Parameters.getParameter<int32_t>("Nbinsy"),
1086  Parameters.getParameter<double>("ymin"),
1087  Parameters.getParameter<double>("ymax")
1088  );
1089  ibooker.tag(mod_single.Module_ClusWidthVsAmpTH2, detid);
1090  mod_single.Module_ClusWidthVsAmpTH2->setAxisTitle("Amplitudes (integrated ADC counts)",1);
1091  mod_single.Module_ClusWidthVsAmpTH2->setAxisTitle("Cluster widths",2);
1092  }
1093 
1094 }
1095 //
1096 // -- Create Module Level MEs
1097 //
1098 void SiStripMonitorCluster::createLayerMEs(std::string label, int ndets , DQMStore::IBooker & ibooker) {
1099 
1100  SiStripHistoId hidmanager;
1101 
1102  LayerMEs layerMEs;
1103  layerMEs.LayerClusterStoN = 0;
1104  layerMEs.LayerClusterStoNTrend = 0;
1105  layerMEs.LayerClusterCharge = 0;
1106  layerMEs.LayerClusterChargeTrend = 0;
1107  layerMEs.LayerClusterNoise = 0;
1108  layerMEs.LayerClusterNoiseTrend = 0;
1109  layerMEs.LayerClusterWidth = 0;
1110  layerMEs.LayerClusterWidthTrend = 0;
1111  layerMEs.LayerLocalOccupancy = 0;
1112  layerMEs.LayerLocalOccupancyTrend = 0;
1113  layerMEs.LayerNumberOfClusterProfile = 0;
1114  layerMEs.LayerNumberOfClusterTrend = 0;
1115  layerMEs.LayerNumberOfClusterPerRingTrend = 0;
1116  layerMEs.LayerClusterWidthProfile = 0;
1117  layerMEs.LayerClusWidthVsAmpTH2 = 0;
1118  layerMEs.LayerClusterPosition = 0;
1119 
1120  //Cluster Width
1121  if(layerswitchcluswidthon) {
1122  layerMEs.LayerClusterWidth=bookME1D("TH1ClusterWidth", hidmanager.createHistoLayer("Summary_ClusterWidth","layer",label,"").c_str() , ibooker );
1123  if (createTrendMEs)
1124  layerMEs.LayerClusterWidthTrend=bookMETrend(hidmanager.createHistoLayer("Trend_ClusterWidth","layer",label,"").c_str() , ibooker );
1125  }
1126 
1127  //Cluster Noise
1128  if(layerswitchclusnoiseon) {
1129  layerMEs.LayerClusterNoise=bookME1D("TH1ClusterNoise", hidmanager.createHistoLayer("Summary_ClusterNoise","layer",label,"").c_str() , ibooker );
1130  if (createTrendMEs)
1131  layerMEs.LayerClusterNoiseTrend=bookMETrend(hidmanager.createHistoLayer("Trend_ClusterNoise","layer",label,"").c_str() , ibooker);
1132  }
1133 
1134  //Cluster Charge
1135  if(layerswitchcluschargeon) {
1136  layerMEs.LayerClusterCharge=bookME1D("TH1ClusterCharge", hidmanager.createHistoLayer("Summary_ClusterCharge","layer",label,"").c_str() , ibooker );
1137  if (createTrendMEs)
1138  layerMEs.LayerClusterChargeTrend=bookMETrend(hidmanager.createHistoLayer("Trend_ClusterCharge","layer",label,"").c_str(),ibooker);
1139  }
1140 
1141  //Cluster StoN
1142  if(layerswitchclusstonon) {
1143  layerMEs.LayerClusterStoN=bookME1D("TH1ClusterStoN", hidmanager.createHistoLayer("Summary_ClusterSignalOverNoise","layer",label,"").c_str() , ibooker );
1144  if (createTrendMEs)
1145  layerMEs.LayerClusterStoNTrend=bookMETrend(hidmanager.createHistoLayer("Trend_ClusterSignalOverNoise","layer",label,"").c_str(),ibooker);
1146  }
1147 
1148  //Cluster Occupancy
1149  if(layerswitchlocaloccupancy) {
1150  layerMEs.LayerLocalOccupancy=bookME1D("TH1ModuleLocalOccupancy", hidmanager.createHistoLayer("Summary_ClusterLocalOccupancy","layer",label,"").c_str() , ibooker );
1151  if (createTrendMEs)
1152  layerMEs.LayerLocalOccupancyTrend=bookMETrend(hidmanager.createHistoLayer("Trend_ClusterLocalOccupancy","layer",label,"").c_str(),ibooker);
1153  }
1154 
1155  // # of Cluster Trend
1156  if (subdetswitchtotclusprofon){
1157  layerMEs.LayerNumberOfClusterTrend = bookMETrend(hidmanager.createHistoLayer("NumberOfClusterTrend","layer",label,"").c_str(),ibooker);
1158 
1159  if (label.find("TID")!=std::string::npos || label.find("TEC")!=std::string::npos){
1160  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("NumberOfClusterPerRingVsTrendVarTH2");
1161  layerMEs.LayerNumberOfClusterPerRingTrend = bookME2D("NumberOfClusterPerRingVsTrendVarTH2", hidmanager.createHistoLayer("NumberOfClusterPerRing_vs_TrendVar","layer",label,"").c_str() , ibooker );
1162  }
1163  }
1164 
1165  // # of Cluster Profile
1166  if(layerswitchnumclusterprofon) {
1167  std::string hid = hidmanager.createHistoLayer("NumberOfClusterProfile","layer",label,"");
1168  layerMEs.LayerNumberOfClusterProfile = ibooker.bookProfile(hid, hid, ndets, 0.5, ndets+0.5,21, -0.5, 20.5);
1169  }
1170 
1171  // Cluster Width Profile
1172  if(layerswitchclusterwidthprofon) {
1173  std::string hid = hidmanager.createHistoLayer("ClusterWidthProfile","layer",label,"");
1174  layerMEs.LayerClusterWidthProfile = ibooker.bookProfile(hid, hid, ndets, 0.5, ndets+0.5, 20, -0.5, 19.5);
1175  }
1176 
1177  if(layer_clusterWidth_vs_amplitude_on) {
1178  layerMEs.LayerClusWidthVsAmpTH2=bookME2D("ClusWidthVsAmpTH2", hidmanager.createHistoLayer("ClusterWidths_vs_Amplitudes","layer",label,"").c_str() , ibooker );
1179  }
1180 
1181  // Cluster Position
1182  if (layerswitchclusposon) {
1183  std::string hid = hidmanager.createHistoLayer("ClusterPosition","layer",label,"");
1184  layerMEs.LayerClusterPosition=ibooker.book1D(hid, hid, 128*6, 0.5, 128*6+0.5);
1185  }
1186 
1187  LayerMEsMap[label]=layerMEs;
1188 }
1189 //
1190 // -- Create SubDetector MEs
1191 //
1192 void SiStripMonitorCluster::createSubDetMEs(std::string label , DQMStore::IBooker & ibooker ) {
1193 
1194  SubDetMEs subdetMEs;
1195  subdetMEs.totNClusters = 0;
1196  subdetMEs.SubDetTotClusterTH1 = 0;
1197  subdetMEs.SubDetTotClusterProf = 0;
1198  subdetMEs.SubDetClusterApvProf = 0;
1199  subdetMEs.SubDetClusterApvTH2 = 0;
1200  subdetMEs.SubDetClusterDBxCycleProf = 0;
1201  subdetMEs.SubDetApvDBxProf2 = 0;
1202  subdetMEs.SubDetClusterChargeTH1 = 0;
1203  subdetMEs.SubDetClusterWidthTH1 = 0;
1204  subdetMEs.SubDetClusWidthVsAmpTH2 = 0;
1205  subdetMEs.SubDetNumberOfClusterPerLayerTrend = 0;
1206 
1207 
1208  std::string HistoName;
1209  // cluster charge
1210  if (subdetswitchcluschargeon){
1211  HistoName = "ClusterCharge__" + label;
1212  subdetMEs.SubDetClusterChargeTH1 = bookME1D("TH1ClusterCharge",HistoName.c_str() , ibooker);
1213  subdetMEs.SubDetClusterChargeTH1->setAxisTitle("Cluster charge [ADC counts]");
1214  subdetMEs.SubDetClusterChargeTH1->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
1215  }
1216  // cluster width
1217  if (subdetswitchcluswidthon){
1218  HistoName = "ClusterWidth__" + label;
1219  subdetMEs.SubDetClusterWidthTH1 = bookME1D("TH1ClusterWidth",HistoName.c_str() , ibooker);
1220  subdetMEs.SubDetClusterWidthTH1->setAxisTitle("Cluster width [strips]");
1221  subdetMEs.SubDetClusterWidthTH1->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
1222  }
1223  // Total Number of Cluster - 1D
1224  if (subdetswitchtotclusth1on){
1225  HistoName = "TotalNumberOfCluster__" + label;
1226  subdetMEs.SubDetTotClusterTH1 = bookME1D("TH1TotalNumberOfClusters",HistoName.c_str() , ibooker);
1227  subdetMEs.SubDetTotClusterTH1->setAxisTitle("Total number of clusters in subdetector");
1228  subdetMEs.SubDetTotClusterTH1->getTH1()->StatOverflows(kTRUE); // over/underflows in Mean calculation
1229  }
1230  // Total Number of Cluster vs Time - Profile
1231  if (subdetswitchtotclusprofon){
1232  edm::ParameterSet Parameters = trendVsLs_ ? conf_.getParameter<edm::ParameterSet>("TrendingLS") : conf_.getParameter<edm::ParameterSet>("Trending");
1233  HistoName = "TotalNumberOfClusterProfile__" + label;
1234  subdetMEs.SubDetTotClusterProf = ibooker.bookProfile(HistoName,HistoName,
1235  Parameters.getParameter<int32_t>("Nbins"),
1236  Parameters.getParameter<double>("xmin"),
1237  Parameters.getParameter<double>("xmax"),
1238  0 , 0 , "" );
1239  subdetMEs.SubDetTotClusterProf->setAxisTitle(Parameters.getParameter<std::string>("xaxis"),1);
1240  if (subdetMEs.SubDetTotClusterProf->kind() == MonitorElement::DQM_KIND_TPROFILE) subdetMEs.SubDetTotClusterProf->getTH1()->SetCanExtend(TH1::kAllAxes);
1241 
1242  HistoName = "TotalNumberOfClusterPerLayer__" + label;
1243  subdetMEs.SubDetNumberOfClusterPerLayerTrend = bookME2D("NumberOfClusterPerLayerTrendVarTH2", HistoName.c_str() , ibooker );
1244  subdetMEs.SubDetNumberOfClusterPerLayerTrend->setAxisTitle("Lumisection",1);
1245  subdetMEs.SubDetNumberOfClusterPerLayerTrend->setAxisTitle("Layer Number",2);
1246 
1247  }
1248 
1249  // Total Number of Cluster vs APV cycle - Profile
1250  if(subdetswitchapvcycleprofon){
1251  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfClustersApvCycle");
1252  HistoName = "Cluster_vs_ApvCycle__" + label;
1253  subdetMEs.SubDetClusterApvProf=ibooker.bookProfile(HistoName,HistoName,
1254  Parameters.getParameter<int32_t>("Nbins"),
1255  Parameters.getParameter<double>("xmin"),
1256  Parameters.getParameter<double>("xmax"),
1257  200, //that parameter should not be there !?
1258  Parameters.getParameter<double>("ymin"),
1259  Parameters.getParameter<double>("ymax"),
1260  "" );
1261  subdetMEs.SubDetClusterApvProf->setAxisTitle("Apv Cycle (Corrected Absolute Bx % 70)",1);
1262  }
1263 
1264  // Total Number of Clusters vs ApvCycle - 2D
1265  if(subdetswitchapvcycleth2on){
1266  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TH2ClustersApvCycle");
1267  HistoName = "Cluster_vs_ApvCycle_2D__" + label;
1268  // Adjusting the scale for 2D histogram
1269  double h2ymax = 9999.0;
1270  double yfact = Parameters.getParameter<double>("yfactor");
1271  if(label.find("TIB") != std::string::npos) h2ymax = (6984.*256.)*yfact;
1272  else if (label.find("TID") != std::string::npos) h2ymax = (2208.*256.)*yfact;
1273  else if (label.find("TOB") != std::string::npos) h2ymax = (12906.*256.)*yfact;
1274  else if (label.find("TEC") != std::string::npos) h2ymax = (7552.*2.*256.)*yfact;
1275 
1276  subdetMEs.SubDetClusterApvTH2=ibooker.book2D(HistoName,HistoName,
1277  Parameters.getParameter<int32_t>("Nbinsx"),
1278  Parameters.getParameter<double>("xmin"),
1279  Parameters.getParameter<double>("xmax"),
1280  Parameters.getParameter<int32_t>("Nbinsy"),
1281  Parameters.getParameter<double>("ymin"),
1282  h2ymax);
1283  subdetMEs.SubDetClusterApvTH2->setAxisTitle("Apv Cycle (Corrected Absolute Bx % 70))",1);
1284  subdetMEs.SubDetClusterApvTH2->setAxisTitle("Total # of Clusters",2);
1285 
1286  }
1287 
1288  // Cluster widths vs amplitudes - 2D
1289  if(subdet_clusterWidth_vs_amplitude_on){
1290  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("ClusWidthVsAmpTH2");
1291  HistoName = "ClusterWidths_vs_Amplitudes__" + label;
1292  subdetMEs.SubDetClusWidthVsAmpTH2=ibooker.book2D(HistoName,HistoName,
1293  Parameters.getParameter<int32_t>("Nbinsx"),
1294  Parameters.getParameter<double>("xmin"),
1295  Parameters.getParameter<double>("xmax"),
1296  Parameters.getParameter<int32_t>("Nbinsy"),
1297  Parameters.getParameter<double>("ymin"),
1298  Parameters.getParameter<double>("ymax"));
1299  subdetMEs.SubDetClusWidthVsAmpTH2->setAxisTitle("Amplitudes (integrated ADC counts)",1);
1300  subdetMEs.SubDetClusWidthVsAmpTH2->setAxisTitle("Cluster widths",2);
1301 
1302  }
1303 
1304 
1305 
1306  // Total Number of Cluster vs DeltaBxCycle - Profile
1307  if(subdetswitchdbxcycleprofon){
1308  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProfClustersVsDBxCycle");
1309  HistoName = "Cluster_vs_DeltaBxCycle__" + label;
1310  subdetMEs.SubDetClusterDBxCycleProf = ibooker.bookProfile(HistoName,HistoName,
1311  Parameters.getParameter<int32_t>("Nbins"),
1312  Parameters.getParameter<double>("xmin"),
1313  Parameters.getParameter<double>("xmax"),
1314  200, //that parameter should not be there !?
1315  Parameters.getParameter<double>("ymin"),
1316  Parameters.getParameter<double>("ymax"),
1317  "" );
1318  subdetMEs.SubDetClusterDBxCycleProf->setAxisTitle("Delta Bunch Crossing Cycle",1);
1319  }
1320  // DeltaBx vs ApvCycle - 2DProfile
1321  if(subdetswitchapvcycledbxprof2on){
1322  edm::ParameterSet Parameters = conf_.getParameter<edm::ParameterSet>("TProf2ApvCycleVsDBx");
1323  HistoName = "DeltaBx_vs_ApvCycle__" + label;
1324  subdetMEs.SubDetApvDBxProf2 = ibooker.bookProfile2D(HistoName,HistoName,
1325  Parameters.getParameter<int32_t>("Nbinsx"),
1326  Parameters.getParameter<double>("xmin"),
1327  Parameters.getParameter<double>("xmax"),
1328  Parameters.getParameter<int32_t>("Nbinsy"),
1329  Parameters.getParameter<double>("ymin"),
1330  Parameters.getParameter<double>("ymax"),
1331  Parameters.getParameter<double>("zmin"),
1332  Parameters.getParameter<double>("zmax"),
1333  "" );
1334  subdetMEs.SubDetApvDBxProf2->setAxisTitle("APV Cycle (Corrected Absolute Bx % 70)",1);
1335  subdetMEs.SubDetApvDBxProf2->setAxisTitle("Delta Bunch Crossing Cycle",2);
1336  }
1337  SubDetMEsMap[label]=subdetMEs;
1338 }
1339 
1340 //
1341 // -- Fill Module Level Histograms
1342 //
1343 void SiStripMonitorCluster::fillModuleMEs(ModMEs& mod_mes, ClusterProperties& cluster) {
1344 
1345  if(moduleswitchclusposon && (mod_mes.ClusterPosition)) // position of cluster
1346  (mod_mes.ClusterPosition)->Fill(cluster.position);
1347 
1348  // position of digis in cluster
1349  if(moduleswitchclusdigiposon && (mod_mes.ClusterDigiPosition)) {
1350  for(int ipos=cluster.start+1; ipos<=cluster.start+cluster.width; ipos++){
1351  (mod_mes.ClusterDigiPosition)->Fill(ipos);
1352  }
1353  }
1354 
1355  if(moduleswitchcluswidthon && (mod_mes.ClusterWidth)) // width of cluster
1356  (mod_mes.ClusterWidth)->Fill(static_cast<float>(cluster.width));
1357 
1358  if(moduleswitchclusstonon && (mod_mes.ClusterSignalOverNoise)) {// SignalToNoise
1359  if (cluster.noise > 0)
1360  (mod_mes.ClusterSignalOverNoise)->Fill(cluster.charge/cluster.noise);
1361  }
1362 
1363  if(moduleswitchclusstonVsposon && (mod_mes.ClusterSignalOverNoiseVsPos)) {// SignalToNoise
1364  if (cluster.noise > 0)
1365  (mod_mes.ClusterSignalOverNoiseVsPos)->Fill(cluster.position,cluster.charge/cluster.noise);
1366  }
1367 
1368  if(moduleswitchclusnoiseon && (mod_mes.ClusterNoise)) // Noise
1369  (mod_mes.ClusterNoise)->Fill(cluster.noise);
1370 
1371  if(moduleswitchcluschargeon && (mod_mes.ClusterCharge)) // charge of cluster
1372  (mod_mes.ClusterCharge)->Fill(cluster.charge);
1373 
1374  if(module_clusterWidth_vs_amplitude_on) (mod_mes.Module_ClusWidthVsAmpTH2)->Fill(cluster.charge, cluster.width);
1375 }
1376 //
1377 // -- Fill Layer Level MEs
1378 //
1379 void SiStripMonitorCluster::fillLayerMEs(LayerMEs& layerMEs, ClusterProperties& cluster) {
1380 
1381  if(layerswitchclusstonon) {
1382  fillME(layerMEs.LayerClusterStoN ,cluster.charge/cluster.noise);
1383  if (createTrendMEs) {
1384  fillME(layerMEs.LayerClusterStoNTrend,trendVar,cluster.charge/cluster.noise);
1385  }
1386  }
1387 
1388  if(layerswitchcluschargeon) {
1389  fillME(layerMEs.LayerClusterCharge,cluster.charge);
1390  if (createTrendMEs) {
1391  fillME(layerMEs.LayerClusterChargeTrend,trendVar,cluster.charge);
1392  }
1393  }
1394 
1395  if(layerswitchclusnoiseon) {
1396  fillME(layerMEs.LayerClusterNoise ,cluster.noise);
1397  if (createTrendMEs) {
1398  fillME(layerMEs.LayerClusterNoiseTrend,trendVar,cluster.noise);
1399  }
1400  }
1401 
1402  if(layerswitchcluswidthon) {
1403  fillME(layerMEs.LayerClusterWidth ,cluster.width);
1404  if (createTrendMEs) {
1405  fillME(layerMEs.LayerClusterWidthTrend,trendVar,cluster.width);
1406  }
1407  }
1408 
1409  if( layer_clusterWidth_vs_amplitude_on ) {
1410  fillME(layerMEs.LayerClusWidthVsAmpTH2, cluster.charge, cluster.width);
1411  }
1412 
1413  if (layerswitchclusposon) {
1414  fillME(layerMEs.LayerClusterPosition,cluster.position);
1415  }
1416 
1417 }
1418 //------------------------------------------------------------------------------------------
1419 MonitorElement* SiStripMonitorCluster::bookMETrend(const char* HistoName , DQMStore::IBooker & ibooker)
1420 {
1421  edm::ParameterSet ParametersTrend = trendVsLs_ ? conf_.getParameter<edm::ParameterSet>("TrendingLS") : conf_.getParameter<edm::ParameterSet>("Trending");
1422  MonitorElement* me = ibooker.bookProfile(HistoName,HistoName,
1423  ParametersTrend.getParameter<int32_t>("Nbins"),
1424  ParametersTrend.getParameter<double>("xmin"),
1425  ParametersTrend.getParameter<double>("xmax"),
1426  0 , 0 , "" );
1427  if(!me) return me;
1428  me->setAxisTitle(ParametersTrend.getParameter<std::string>("xaxis"),1);
1429  if (me->kind() == MonitorElement::DQM_KIND_TPROFILE) me->getTH1()->SetCanExtend(TH1::kAllAxes);
1430  return me;
1431 }
1432 
1433 //------------------------------------------------------------------------------------------
1434 MonitorElement* SiStripMonitorCluster::bookME1D(const char* ParameterSetLabel, const char* HistoName , DQMStore::IBooker & ibooker)
1435 {
1436  Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
1437  return ibooker.book1D(HistoName,HistoName,
1438  Parameters.getParameter<int32_t>("Nbinx"),
1439  Parameters.getParameter<double>("xmin"),
1440  Parameters.getParameter<double>("xmax")
1441  );
1442 }
1443 
1444 //------------------------------------------------------------------------------------------
1445 MonitorElement* SiStripMonitorCluster::bookME2D(const char* ParameterSetLabel, const char* HistoName , DQMStore::IBooker & ibooker)
1446 {
1447  Parameters = conf_.getParameter<edm::ParameterSet>(ParameterSetLabel);
1448  return ibooker.book2D(HistoName,HistoName,
1449  Parameters.getParameter<int32_t>("Nbinsx"),
1450  Parameters.getParameter<double>("xmin"),
1451  Parameters.getParameter<double>("xmax"),
1452  Parameters.getParameter<int32_t>("Nbinsy"),
1453  Parameters.getParameter<double>("ymin"),
1454  Parameters.getParameter<double>("ymax")
1455  );
1456 }
1457 
1458 
1459 int SiStripMonitorCluster::FindRegion(int nstrip,int npix){
1460 
1461  double kplus= k0*(1+dk0/100);
1462  double kminus=k0*(1-dk0/100);
1463  int region=0;
1464 
1465  if (nstrip!=0 && npix >= (nstrip*kminus-q0) && npix <=(nstrip*kplus+q0)) region=1;
1466  else if (nstrip!=0 && npix < (nstrip*kminus-q0) && nstrip <= maxClus) region=2;
1467  else if (nstrip!=0 && npix < (nstrip*kminus-q0) && nstrip > maxClus) region=3;
1468  else if (nstrip!=0 && npix > (nstrip*kplus+q0)) region=4;
1469  else if (npix > minPix && nstrip==0) region=5;
1470  return region;
1471 
1472 }
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
SiStripMonitorCluster::GlobalABXTH1_CSCP
MonitorElement * GlobalABXTH1_CSCP
Definition: SiStripMonitorCluster.h:108
SiStripMonitorCluster::ResetModuleMEs
void ResetModuleMEs(uint32_t idet)
Definition: SiStripMonitorCluster.cc:960
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
SiStripMonitorCluster::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &)
Definition: SiStripMonitorCluster.cc:585
SiStripMonitorCluster::bookME1D
MonitorElement * bookME1D(const char *ParameterSetLabel, const char *HistoName, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:1434
SiStripMonitorCluster::FindRegion
int FindRegion(int nstrip, int npixel)
Definition: SiStripMonitorCluster.cc:1459
edmNew::DetSetVector::const_iterator
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
Definition: DetSetVectorNew.h:178
SiStripMonitorCluster::NumberOfPixelClus_vs_BX
MonitorElement * NumberOfPixelClus_vs_BX
Definition: SiStripMonitorCluster.h:121
SiStripFolderOrganizer::setSiStripFolderName
void setSiStripFolderName(std::string name)
Definition: SiStripFolderOrganizer.cc:45
SiStripMonitorCluster::PixVsStripMultiplicityRegions
MonitorElement * PixVsStripMultiplicityRegions
Definition: SiStripMonitorCluster.h:109
SiStripFolderOrganizer::getSubDetFolderAndTag
std::pair< const std::string, const char * > getSubDetFolderAndTag(const uint32_t &detid, const TrackerTopology *tTopo)
Definition: SiStripFolderOrganizer.cc:456
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventIDconst &, edm::Timestampconst & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
DQMStore::IBooker::bookProfile
MonitorElement * bookProfile(Args &&...args)
Definition: DQMStore.h:157
SiStripMonitorCluster::GlobalTH2MainDiagonalPosition_vs_BX
MonitorElement * GlobalTH2MainDiagonalPosition_vs_BX
Definition: SiStripMonitorCluster.h:112
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
SiStripMonitorCluster::GlobalCStripVsCpix
MonitorElement * GlobalCStripVsCpix
Definition: SiStripMonitorCluster.h:107
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
SiStripMonitorCluster::SiStripDetCabling_
edm::ESHandle< SiStripDetCabling > SiStripDetCabling_
Definition: SiStripMonitorCluster.h:154
SiStripFolderOrganizer::GetSubDetAndLayer
std::pair< std::string, int32_t > GetSubDetAndLayer(const uint32_t &detid, const TrackerTopology *tTopo, bool ring_flag=0)
Definition: SiStripFolderOrganizer.cc:119
NULL
#define NULL
Definition: scimark2.h:8
SiStripMonitorCluster::ModMEs::ClusterSignalOverNoiseVsPos
MonitorElement * ClusterSignalOverNoiseVsPos
Definition: SiStripMonitorCluster.h:56
GenericTriggerEventFlag
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Definition: GenericTriggerEventFlag.h:41
SiStripMonitorCluster::StripNoise2Cycle
MonitorElement * StripNoise2Cycle
Definition: SiStripMonitorCluster.h:113
edmNew::DetSet::const_iterator
data_type const * const_iterator
Definition: DetSetNew.h:30
SiStripMonitorCluster::createLayerMEs
void createLayerMEs(std::string label, int ndets, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:1098
SiStripMonitorCluster::LayerMEsMap
std::map< std::string, LayerMEs > LayerMEsMap
Definition: SiStripMonitorCluster.h:145
edm::EventBase::isRealData
bool isRealData() const
Definition: EventBase.h:63
hlt_dqm_clientPB-live_cfg.me
tuple me
Definition: hlt_dqm_clientPB-live_cfg.py:43
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
SiStripMonitorCluster::LayerDetMap
std::map< std::string, std::vector< uint32_t > > LayerDetMap
Definition: SiStripMonitorCluster.h:146
SiStripMonitorClusterAlca_cfi.Nbinsx
tuple Nbinsx
Definition: SiStripMonitorClusterAlca_cfi.py:147
SiStripMonitorCluster::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: SiStripMonitorCluster.cc:563
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
SiStripMonitorCluster::NumberOfPixelClus
MonitorElement * NumberOfPixelClus
Definition: SiStripMonitorCluster.h:115
SiStripFolderOrganizer::setDetectorFolder
void setDetectorFolder(uint32_t rawdetid, const TrackerTopology *tTopo)
Definition: SiStripFolderOrganizer.cc:211
SiStripMonitorCluster::m_cacheID_
unsigned long long m_cacheID_
Definition: SiStripMonitorCluster.h:152
SiStripMonitorCluster::createSubDetMEs
void createSubDetMEs(std::string label, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:1192
SiStripMonitorCluster::ModuleMEsMap
std::map< uint32_t, ModMEs > ModuleMEsMap
Definition: SiStripMonitorCluster.h:144
SiStripDCSStatus::getStatus
bool getStatus(edm::Event const &e, edm::EventSetup const &eSetup)
Definition: SiStripDCSStatus.cc:45
iEvent
int iEvent
Definition: GenABIO.cc:230
SiStripMonitorCluster::createMEs
void createMEs(const edm::EventSetup &es, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:232
SiStripMonitorCluster::SubDetPhasePartMap
std::map< std::string, std::string > SubDetPhasePartMap
Definition: SiStripMonitorCluster.h:148
SiStripMonitorCluster::GlobalApvCycleDBxTH2
MonitorElement * GlobalApvCycleDBxTH2
Definition: SiStripMonitorCluster.h:104
Parameters
vector< ParameterSet > Parameters
Definition: HLTMuonPlotter.cc:31
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
DQMStore::IBooker::bookProfile2D
MonitorElement * bookProfile2D(Args &&...args)
Definition: DQMStore.h:163
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
SiStripMonitorCluster::genTriggerEventFlagPixelDCSfilter_
GenericTriggerEventFlag * genTriggerEventFlagPixelDCSfilter_
Definition: SiStripMonitorCluster.h:249
edm::EDConsumerBase::consumesCollector
ConsumesCollector consumesCollector()
Use a ConsumesCollector to gather consumes information from helper functions.
Definition: EDConsumerBase.cc:72
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:18
TkHistoMap::fill
void fill(uint32_t &detid, float value)
Definition: TkHistoMap.cc:180
GenericTriggerEventFlag::accept
bool accept(const edm::Event &event, const edm::EventSetup &setup)
To be called from analyze/filter() methods.
Definition: GenericTriggerEventFlag.cc:239
edmNew::DetSetVector::data
data_type const * data(size_t cell) const
Definition: DetSetVectorNew.h:561
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
DQMStore::IBooker::tag
void tag(MonitorElement *, unsigned int)
Definition: DQMStore.cc:288
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
SiStripApvGain::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripApvGain.h:27
SiStripMonitorCluster::bookMETrend
MonitorElement * bookMETrend(const char *, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:1419
MonitorElement::getTH1
TH1 * getTH1(void) const
Definition: MonitorElement.cc:1450
HLT_25ns10e33_v2_cff.noise
tuple noise
Definition: HLT_25ns10e33_v2_cff.py:4362
edm::EventBase::orbitNumber
int orbitNumber() const
Definition: EventBase.h:66
SiStripMonitorCluster::NumberOfStripClus
MonitorElement * NumberOfStripClus
Definition: SiStripMonitorCluster.h:116
MonitorElement::kind
Kind kind(void) const
Get the type of the monitor element.
Definition: MonitorElement.h:109
SiStripHistoId::getSubdetid
std::string getSubdetid(uint32_t id, const TrackerTopology *tTopo, bool flag_ring)
Definition: SiStripHistoId.cc:84
SiStripMonitorCluster::dqmBeginRun
void dqmBeginRun(const edm::Run &, const edm::EventSetup &)
Definition: SiStripMonitorCluster.cc:220
SiStripMonitorCluster::GlobalMainDiagonalPosition_vs_BX
MonitorElement * GlobalMainDiagonalPosition_vs_BX
Definition: SiStripMonitorCluster.h:111
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
SiStripMonitorCluster::NumberOfStripClus_vs_BX
MonitorElement * NumberOfStripClus_vs_BX
Definition: SiStripMonitorCluster.h:120
SiStripMonitorCluster::historyProducerToken_
edm::EDGetTokenT< EventWithHistory > historyProducerToken_
Definition: SiStripMonitorCluster.h:230
SiStripMonitorCluster::ClusWidthVsAmpTH2
MonitorElement * ClusWidthVsAmpTH2
Definition: SiStripMonitorCluster.h:119
SiStripMonitorCluster::GlobalMainDiagonalPosition
MonitorElement * GlobalMainDiagonalPosition
Definition: SiStripMonitorCluster.h:110
SiStripMonitorCluster::bookME2D
MonitorElement * bookME2D(const char *ParameterSetLabel, const char *HistoName, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:1445
SiStripMonitorCluster::dcsStatus_
SiStripDCSStatus * dcsStatus_
Definition: SiStripMonitorCluster.h:245
edm::HandleBase::failedToGet
bool failedToGet() const
Definition: HandleBase.h:79
SiStripMonitorCluster::genTriggerEventFlagBPTXfilter_
GenericTriggerEventFlag * genTriggerEventFlagBPTXfilter_
Definition: SiStripMonitorCluster.h:248
SiStripMonitorCluster::LayerMEs::LayerNumberOfClusterPerRingTrend
MonitorElement * LayerNumberOfClusterPerRingTrend
Definition: SiStripMonitorCluster.h:74
HLT_25ns10e33_v2_cff.region
tuple region
Definition: HLT_25ns10e33_v2_cff.py:61599
SiStripFolderOrganizer::setLayerFolder
void setLayerFolder(uint32_t rawdetid, const TrackerTopology *tTopo, int32_t layer=0, bool ring_flag=0)
Definition: SiStripFolderOrganizer.cc:300
SiStripMonitorCluster::fillME
void fillME(MonitorElement *ME, float value1)
Definition: SiStripMonitorCluster.h:135
SiStripMonitorCluster::SiStripMonitorCluster
SiStripMonitorCluster(const edm::ParameterSet &)
Definition: SiStripMonitorCluster.cc:45
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
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
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
cscdqm::HistoName
std::string HistoName
Definition: CSCDQM_HistoDef.h:32
SiStripMonitorCluster::StripNoise3Cycle
MonitorElement * StripNoise3Cycle
Definition: SiStripMonitorCluster.h:114
SiStripMonitorCluster::createModuleMEs
void createModuleMEs(ModMEs &mod_single, uint32_t detid, DQMStore::IBooker &ibooker)
Definition: SiStripMonitorCluster.cc:979
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
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
edmNew::DetSet::end
iterator end()
Definition: DetSetNew.h:70
SiStripMonitorCluster::clusterProducerPixToken_
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > > clusterProducerPixToken_
Definition: SiStripMonitorCluster.h:229
SiStripMonitorCluster::SubDetMEsMap
std::map< std::string, SubDetMEs > SubDetMEsMap
Definition: SiStripMonitorCluster.h:147
SiStripMonitorCluster::fillLayerMEs
void fillLayerMEs(LayerMEs &, ClusterProperties &cluster)
Definition: SiStripMonitorCluster.cc:1379
SiStripBadStrip::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripBadStrip.h:54
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
SiStripMonitorCluster::clusterProducerStripToken_
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > clusterProducerStripToken_
Definition: SiStripMonitorCluster.h:228
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
edm::false
volatile std::atomic< bool > shutdown_flag false
Definition: UnixSignalHandlers.cc:22
GenericTriggerEventFlag::initRun
void initRun(const edm::Run &run, const edm::EventSetup &setup)
To be called from beginRun() methods.
Definition: GenericTriggerEventFlag.cc:134
SiStripHistoId::createHistoLayer
std::string createHistoLayer(std::string description, std::string id_type, std::string path, std::string flag)
Definition: SiStripHistoId.cc:57
SiStripNoises::Range
std::pair< ContainerIterator, ContainerIterator > Range
Definition: SiStripNoises.h:48
edmNew::DetSet::size
size_type size() const
Definition: DetSetNew.h:87
SiStripMonitorCluster::GlobalDBxCycleTH1
MonitorElement * GlobalDBxCycleTH1
Definition: SiStripMonitorCluster.h:106
SiStripMonitorCluster::apvPhaseProducerToken_
edm::EDGetTokenT< APVCyclePhaseCollection > apvPhaseProducerToken_
Definition: SiStripMonitorCluster.h:231
SiStripMonitorCluster::genTriggerEventFlagStripDCSfilter_
GenericTriggerEventFlag * genTriggerEventFlagStripDCSfilter_
Definition: SiStripMonitorCluster.h:250
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
SiStripMonitorCluster::NclusVsCycleTimeProf2D
MonitorElement * NclusVsCycleTimeProf2D
Definition: SiStripMonitorCluster.h:118
SiStripMonitorCluster::fillModuleMEs
void fillModuleMEs(ModMEs &mod_mes, ClusterProperties &cluster)
Definition: SiStripMonitorCluster.cc:1343
SiStripMonitorCluster::SubDetMEs::SubDetNumberOfClusterPerLayerTrend
MonitorElement * SubDetNumberOfClusterPerLayerTrend
Definition: SiStripMonitorCluster.h:93
edm::Run
Definition: Run.h:43
SiPixelCluster::size
int size() const
Definition: SiPixelCluster.h:136
edmNew::DetSet::begin
iterator begin()
Definition: DetSetNew.h:67
Generated for CMSSW Reference Manual by doxygen 1.8.5