CMS 3D CMS Logo
DTSegmentsTask.cc
Go to the documentation of this file.
1 /*
2  * See header file for a description of this class.
3  *
4  * \author G. Mila - INFN Torino
5  */
6 
7 #include "DQMOffline/Muon/interface/DTSegmentsTask.h"
8 
9 // Framework
10 #include "FWCore/Framework/interface/Event.h"
11 #include "FWCore/Framework/interface/EventSetup.h"
12 #include "FWCore/Framework/interface/ESHandle.h"
13 #include "FWCore/ParameterSet/interface/ParameterSet.h"
14 #include "FWCore/ServiceRegistry/interface/Service.h"
15 
16 #include "DQMServices/Core/interface/DQMStore.h"
17 
18 //Geometry
19 #include "DataFormats/GeometryVector/interface/Pi.h"
20 #include "DataFormats/MuonDetId/interface/DTChamberId.h"
21 
22 #include <iterator>
23 
24 using namespace edm;
25 using namespace std;
26 
27 DTSegmentsTask::DTSegmentsTask(const edm::ParameterSet& pset)
28  : statusMapToken_(esConsumes<DTStatusFlag, DTStatusFlagRcd>()) {
29  debug = pset.getUntrackedParameter<bool>("debug", false);
30  parameters = pset;
31 
32  // should be init from pset, but it was commented out...
33  // better false than undefined
34  checkNoisyChannels = false;
35 
36  // the name of the 4D rec hits collection
37  theRecHits4DLabel_ = consumes<DTRecSegment4DCollection>(parameters.getParameter<std::string>("recHits4DLabel"));
38 }
39 
40 DTSegmentsTask::~DTSegmentsTask() {}
41 
42 void DTSegmentsTask::bookHistograms(DQMStore::IBooker& ibooker,
43  edm::Run const& /*iRun*/,
44  edm::EventSetup const& /* iSetup */) {
45  ibooker.cd();
46  ibooker.setCurrentFolder("Muons/DTSegmentsMonitor");
47 
48  // histos for phi segments
49  phiHistos.push_back(
50  ibooker.book2D("phiSegments_numHitsVsWheel", "phiSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
51  phiHistos[0]->setBinLabel(1, "W-2", 1);
52  phiHistos[0]->setBinLabel(2, "W-1", 1);
53  phiHistos[0]->setBinLabel(3, "W0", 1);
54  phiHistos[0]->setBinLabel(4, "W1", 1);
55  phiHistos[0]->setBinLabel(5, "W2", 1);
56  phiHistos.push_back(
57  ibooker.book2D("phiSegments_numHitsVsSector", "phiSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
58  phiHistos[1]->setBinLabel(1, "Sec1", 1);
59  phiHistos[1]->setBinLabel(2, "Sec2", 1);
60  phiHistos[1]->setBinLabel(3, "Sec3", 1);
61  phiHistos[1]->setBinLabel(4, "Sec4", 1);
62  phiHistos[1]->setBinLabel(5, "Sec5", 1);
63  phiHistos[1]->setBinLabel(6, "Sec6", 1);
64  phiHistos[1]->setBinLabel(7, "Sec7", 1);
65  phiHistos[1]->setBinLabel(8, "Sec8", 1);
66  phiHistos[1]->setBinLabel(9, "Sec9", 1);
67  phiHistos[1]->setBinLabel(10, "Sec10", 1);
68  phiHistos[1]->setBinLabel(11, "Sec11", 1);
69  phiHistos[1]->setBinLabel(12, "Sec12", 1);
70  phiHistos[1]->setBinLabel(13, "Sec13", 1);
71  phiHistos[1]->setBinLabel(14, "Sec14", 1);
72  phiHistos.push_back(
73  ibooker.book2D("phiSegments_numHitsVsStation", "phiSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
74  phiHistos[2]->setBinLabel(1, "St1", 1);
75  phiHistos[2]->setBinLabel(2, "St2", 1);
76  phiHistos[2]->setBinLabel(3, "St3", 1);
77  phiHistos[2]->setBinLabel(4, "St4", 1);
78 
79  // histos for theta segments
80  thetaHistos.push_back(
81  ibooker.book2D("thetaSegments_numHitsVsWheel", "thetaSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
82  thetaHistos[0]->setBinLabel(1, "W-2", 1);
83  thetaHistos[0]->setBinLabel(2, "W-1", 1);
84  thetaHistos[0]->setBinLabel(3, "W0", 1);
85  thetaHistos[0]->setBinLabel(4, "W1", 1);
86  thetaHistos[0]->setBinLabel(5, "W2", 1);
87  thetaHistos.push_back(
88  ibooker.book2D("thetaSegments_numHitsVsSector", "thetaSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
89  thetaHistos[1]->setBinLabel(1, "Sec1", 1);
90  thetaHistos[1]->setBinLabel(2, "Sec2", 1);
91  thetaHistos[1]->setBinLabel(3, "Sec3", 1);
92  thetaHistos[1]->setBinLabel(4, "Sec4", 1);
93  thetaHistos[1]->setBinLabel(5, "Sec5", 1);
94  thetaHistos[1]->setBinLabel(6, "Sec6", 1);
95  thetaHistos[1]->setBinLabel(7, "Sec7", 1);
96  thetaHistos[1]->setBinLabel(8, "Sec8", 1);
97  thetaHistos[1]->setBinLabel(9, "Sec9", 1);
98  thetaHistos[1]->setBinLabel(10, "Sec10", 1);
99  thetaHistos[1]->setBinLabel(11, "Sec11", 1);
100  thetaHistos[1]->setBinLabel(12, "Sec12", 1);
101  thetaHistos[1]->setBinLabel(13, "Sec13", 1);
102  thetaHistos[1]->setBinLabel(14, "Sec14", 1);
103  thetaHistos.push_back(
104  ibooker.book2D("thetaSegments_numHitsVsStation", "thetaSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
105  thetaHistos[2]->setBinLabel(1, "St1", 1);
106  thetaHistos[2]->setBinLabel(2, "St2", 1);
107  thetaHistos[2]->setBinLabel(3, "St3", 1);
108  thetaHistos[2]->setBinLabel(4, "St4", 1);
109 }
110 
111 void DTSegmentsTask::analyze(const edm::Event& event, const edm::EventSetup& setup) {
112  // Get the map of noisy channels
113  // bool checkNoisyChannels = parameters.getUntrackedParameter<bool>("checkNoisyChannels",false);
114  ESHandle<DTStatusFlag> statusMap;
115  if (checkNoisyChannels) {
116  statusMap = setup.getHandle(statusMapToken_);
117  }
118 
119  // Get the 4D segment collection from the event
120  edm::Handle<DTRecSegment4DCollection> all4DSegments;
121  event.getByToken(theRecHits4DLabel_, all4DSegments);
122 
123  // Loop over all chambers containing a segment
124  DTRecSegment4DCollection::id_iterator chamberId;
125  for (chamberId = all4DSegments->id_begin(); chamberId != all4DSegments->id_end(); ++chamberId) {
126  // Get the range for the corresponding ChamerId
127  DTRecSegment4DCollection::range range = all4DSegments->get(*chamberId);
128  int nsegm = distance(range.first, range.second);
129  if (debug)
130  cout << " Chamber: " << *chamberId << " has " << nsegm << " 4D segments" << endl;
131 
132  // Loop over the rechits of this ChamerId
133  for (DTRecSegment4DCollection::const_iterator segment4D = range.first; segment4D != range.second; ++segment4D) {
134  //FOR NOISY CHANNELS////////////////////////////////
135  bool segmNoisy = false;
136  if ((*segment4D).hasPhi()) {
137  const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
138  vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
139  map<DTSuperLayerId, vector<DTRecHit1D> > hitsBySLMap;
140  for (vector<DTRecHit1D>::const_iterator hit = phiHits.begin(); hit != phiHits.end(); ++hit) {
141  DTWireId wireId = (*hit).wireId();
142 
143  // Check for noisy channels to skip them
144  if (checkNoisyChannels) {
145  bool isNoisy = false;
146  bool isFEMasked = false;
147  bool isTDCMasked = false;
148  bool isTrigMask = false;
149  bool isDead = false;
150  bool isNohv = false;
151  statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
152  if (isNoisy) {
153  if (debug)
154  cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
155  segmNoisy = true;
156  }
157  }
158  }
159  }
160 
161  if ((*segment4D).hasZed()) {
162  const DTSLRecSegment2D* zSeg = (*segment4D).zSegment(); // zSeg lives in the SL RF
163  // Check for noisy channels to skip them
164  vector<DTRecHit1D> zHits = zSeg->specificRecHits();
165  for (vector<DTRecHit1D>::const_iterator hit = zHits.begin(); hit != zHits.end(); ++hit) {
166  DTWireId wireId = (*hit).wireId();
167  if (checkNoisyChannels) {
168  bool isNoisy = false;
169  bool isFEMasked = false;
170  bool isTDCMasked = false;
171  bool isTrigMask = false;
172  bool isDead = false;
173  bool isNohv = false;
174  //cout<<"wire id "<<wireId<<endl;
175  statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
176  if (isNoisy) {
177  if (debug)
178  cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
179  segmNoisy = true;
180  }
181  }
182  }
183  }
184 
185  if (segmNoisy) {
186  if (debug)
187  cout << "skipping the segment: it contains noisy cells" << endl;
188  continue;
189  }
190  //END FOR NOISY CHANNELS////////////////////////////////
191 
192  // Fill the histos
193  int nHits = 0;
194  if ((*segment4D).hasPhi()) {
195  nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
196  if (debug)
197  cout << "Phi segment with number of hits: " << nHits << endl;
198  phiHistos[0]->Fill((*chamberId).wheel(), nHits);
199  phiHistos[1]->Fill((*chamberId).sector(), nHits);
200  phiHistos[2]->Fill((*chamberId).station(), nHits);
201  }
202  if ((*segment4D).hasZed()) {
203  nHits = (((*segment4D).zSegment())->specificRecHits()).size();
204  if (debug)
205  cout << "Zed segment with number of hits: " << nHits << endl;
206  thetaHistos[0]->Fill((*chamberId).wheel(), nHits);
207  thetaHistos[1]->Fill((*chamberId).sector(), nHits);
208  thetaHistos[2]->Fill((*chamberId).station(), nHits);
209  }
210 
211  } //loop over segments
212  } // loop over chambers
213 }
DTSLRecSegment2D
Definition: DTSLRecSegment2D.h:15
FastTimerService_cff.range
range
Definition: FastTimerService_cff.py:34
DTStatusFlag::cellStatus
int cellStatus(int wheelId, int stationId, int sectorId, int slId, int layerId, int cellId, bool &noiseFlag, bool &feMask, bool &tdcMask, bool &trigMask, bool &deadFlag, bool &nohvFlag) const
get content
Definition: DTStatusFlag.h:88
ESHandle.h
edm::Run
Definition: Run.h:45
edm
HLT enums.
Definition: AlignableModifier.h:19
gather_cfg.cout
cout
Definition: gather_cfg.py:144
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
DQMStore.h
edm::Handle< DTRecSegment4DCollection >
edm::RangeMap::id_iterator
identifier iterator
Definition: RangeMap.h:130
singleTopDQM_cfi.setup
setup
Definition: singleTopDQM_cfi.py:37
DTWireId
Definition: DTWireId.h:12
nHits
const caConstants::TupleMultiplicity *__restrict__ const HitsOnGPU *__restrict__ double *__restrict__ float *__restrict__ double *__restrict__ uint32_t nHits
Definition: BrokenLineFitOnGPU.h:27
Service.h
DTSegmentsTask.h
edm::ESHandle
Definition: DTSurvey.h:22
dqm::implementation::NavigatorBase::cd
virtual void cd()
Definition: DQMStore.cc:29
DTSegmentsTask::theRecHits4DLabel_
edm::EDGetTokenT< DTRecSegment4DCollection > theRecHits4DLabel_
Definition: DTSegmentsTask.h:51
DTChamberId.h
DTSegmentsTask::parameters
edm::ParameterSet parameters
Definition: DTSegmentsTask.h:42
DTSegmentsTask::checkNoisyChannels
bool checkNoisyChannels
Definition: DTSegmentsTask.h:41
DTSegmentsTask::debug
bool debug
Definition: DTSegmentsTask.h:40
DTStatusFlagRcd
Definition: DTStatusFlagRcd.h:5
DTSegmentsTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: DTSegmentsTask.cc:42
edm::ParameterSet
Definition: ParameterSet.h:47
Event.h
Pi.h
edm::RangeMap::const_iterator
C::const_iterator const_iterator
constant access iterator type
Definition: RangeMap.h:43
DTChamberRecSegment2D
Definition: DTChamberRecSegment2D.h:31
DTSegmentsTask::statusMapToken_
const edm::ESGetToken< DTStatusFlag, DTStatusFlagRcd > statusMapToken_
Definition: DTSegmentsTask.h:44
edm::EventSetup
Definition: EventSetup.h:58
DTSegmentsTask::thetaHistos
std::vector< MonitorElement * > thetaHistos
Definition: DTSegmentsTask.h:48
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
DTSegmentsTask::~DTSegmentsTask
~DTSegmentsTask() override
Destructor.
Definition: DTSegmentsTask.cc:40
edm::RangeMap::range
std::pair< const_iterator, const_iterator > range
iterator range
Definition: RangeMap.h:50
std
Definition: JetResolutionObject.h:76
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:177
DTRecSegment2D::specificRecHits
std::vector< DTRecHit1D > specificRecHits() const
Access to specific components.
Definition: DTRecSegment2D.cc:104
EventSetup.h
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
dqm::implementation::IBooker
Definition: DQMStore.h:43
ParameterSet.h
DTSegmentsTask::DTSegmentsTask
DTSegmentsTask(const edm::ParameterSet &pset)
Constructor.
Definition: DTSegmentsTask.cc:27
event
Definition: event.py:1
DeDxTools::esConsumes
ESGetTokenH3DDVariant esConsumes(std::string const &Reccord, edm::ConsumesCollector &)
Definition: DeDxTools.cc:283
edm::Event
Definition: Event.h:73
HLT_FULL_cff.distance
distance
Definition: HLT_FULL_cff.py:7746
hit
Definition: SiStripHitEffFromCalibTree.cc:88
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
DTSegmentsTask::phiHistos
std::vector< MonitorElement * > phiHistos
Definition: DTSegmentsTask.h:47
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:443
DTSegmentsTask::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
book the histos
Definition: DTSegmentsTask.cc:111
DTStatusFlag
Definition: DTStatusFlag.h:74
Generated for CMSSW Reference Manual by doxygen 1.8.16