CMS 3D CMS Logo
List of all members | Public Member Functions | Private Attributes
DTSegmentsTask Class Reference

#include <DTSegmentsTask.h>

Inheritance diagram for DTSegmentsTask:
DQMEDAnalyzer edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > > edm::stream::EDAnalyzerBase edm::EDConsumerBase

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &) override
 book the histos More...
 
void bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
 
 DTSegmentsTask (const edm::ParameterSet &pset)
 Constructor. More...
 
 ~DTSegmentsTask () override
 Destructor. More...
 
- Public Member Functions inherited from DQMEDAnalyzer
void beginRun (edm::Run const &, edm::EventSetup const &) final
 
void beginStream (edm::StreamID id) final
 
virtual void dqmBeginRun (edm::Run const &, edm::EventSetup const &)
 
 DQMEDAnalyzer ()
 
void endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
void endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final
 
uint32_t streamId () const
 
- Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
 EDAnalyzer ()=default
 
- Public Member Functions inherited from edm::stream::EDAnalyzerBase
void callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
 
 EDAnalyzerBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDAnalyzerBase () override
 
- Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo > consumesInfo () const
 
void convertCurrentProcessAlias (std::string const &processName)
 Convert "@currentProcess" in InputTag process names to the actual current process name. More...
 
 EDConsumerBase ()
 
 EDConsumerBase (EDConsumerBase const &)=delete
 
 EDConsumerBase (EDConsumerBase &&)=default
 
ProductResolverIndexAndSkipBit indexFrom (EDGetToken, BranchType, TypeID const &) const
 
void itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
void itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
 
std::vector< ProductResolverIndexAndSkipBit > const & itemsToGetFrom (BranchType iType) const
 
void labelsForToken (EDGetToken iToken, Labels &oLabels) const
 
void modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
 
EDConsumerBase const & operator= (EDConsumerBase const &)=delete
 
EDConsumerBase & operator= (EDConsumerBase &&)=default
 
bool registeredToConsume (ProductResolverIndex, bool, BranchType) const
 
bool registeredToConsumeMany (TypeID const &, BranchType) const
 
ProductResolverIndexAndSkipBit uncheckedIndexFrom (EDGetToken) const
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Private Attributes

bool checkNoisyChannels
 
bool debug
 
edm::ParameterSet parameters
 
std::vector< MonitorElement * > phiHistos
 
edm::EDGetTokenT< DTRecSegment4DCollection > theRecHits4DLabel_
 
std::vector< MonitorElement * > thetaHistos
 

Additional Inherited Members

- Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T... > CacheTypes
 
typedef CacheTypes::GlobalCache GlobalCache
 
typedef AbilityChecker< T... > HasAbility
 
typedef CacheTypes::LuminosityBlockCache LuminosityBlockCache
 
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache > LuminosityBlockContext
 
typedef CacheTypes::LuminosityBlockSummaryCache LuminosityBlockSummaryCache
 
typedef CacheTypes::RunCache RunCache
 
typedef RunContextT< RunCache, GlobalCache > RunContext
 
typedef CacheTypes::RunSummaryCache RunSummaryCache
 
- Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase ModuleType
 
- Public Types inherited from edm::EDConsumerBase
typedef ProductLabels Labels
 
- Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr< dqmDetails::NoCache > globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)
 
static std::shared_ptr< dqmDetails::NoCache > globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)
 
static void globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *, dqmDetails::NoCache *)
 
static void globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *, dqmDetails::NoCache *)
 
- Static Public Member Functions inherited from edm::stream::EDAnalyzerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > consumes (edm::InputTag const &tag)
 
EDGetToken consumes (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken consumes (TypeToGet const &id, edm::InputTag const &tag)
 
ConsumesCollector consumesCollector ()
 Use a ConsumesCollector to gather consumes information from helper functions. More...
 
template<typename ProductType , BranchType B = InEvent>
void consumesMany ()
 
void consumesMany (const TypeToGet &id)
 
template<BranchType B>
void consumesMany (const TypeToGet &id)
 
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType > mayConsume (edm::InputTag const &tag)
 
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 
template<BranchType B>
EDGetToken mayConsume (const TypeToGet &id, edm::InputTag const &tag)
 

Detailed Description

DQM Analysis of 4D DT segments

Author
G. Mila - INFN Torino

Definition at line 26 of file DTSegmentsTask.h.

Constructor & Destructor Documentation

DTSegmentsTask::DTSegmentsTask ( const edm::ParameterSet &  pset)

Constructor.

Definition at line 28 of file DTSegmentsTask.cc.

References debug, edm::ParameterSet::getUntrackedParameter(), muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

28  {
29 
30  debug = pset.getUntrackedParameter<bool>("debug",false);
31  parameters = pset;
32 
33  // should be init from pset, but it was commented out...
34  // better false than undefined
35  checkNoisyChannels = false;
36 
37  // the name of the 4D rec hits collection
38  theRecHits4DLabel_ = consumes<DTRecSegment4DCollection>(parameters.getParameter<std::string>("recHits4DLabel"));
39 }
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
AlCaHLTBitMon_QueryRunRegistry.string
string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
muonDTDigis_cfi.pset
pset
Definition: muonDTDigis_cfi.py:27
DTSegmentsTask::theRecHits4DLabel_
edm::EDGetTokenT< DTRecSegment4DCollection > theRecHits4DLabel_
Definition: DTSegmentsTask.h:54
parameters
parameters
Definition: BeamSpot_PayloadInspector.cc:14
DTSegmentsTask::~DTSegmentsTask ( )
override

Destructor.

Definition at line 42 of file DTSegmentsTask.cc.

42  {
43 }

Member Function Documentation

void DTSegmentsTask::analyze ( const edm::Event &  event,
const edm::EventSetup &  setup 
)
override

book the histos

Definition at line 110 of file DTSegmentsTask.cc.

References DTStatusFlag::cellStatus(), gather_cfg::cout, debug, SoftLeptonByDistance_cfi::distance, edm::EventSetup::get(), findQualityFiles::size, and DTRecSegment2D::specificRecHits().

110  {
111 
112  // Get the map of noisy channels
113  // bool checkNoisyChannels = parameters.getUntrackedParameter<bool>("checkNoisyChannels",false);
114  ESHandle<DTStatusFlag> statusMap;
115  if(checkNoisyChannels) {
116  setup.get<DTStatusFlagRcd>().get(statusMap);
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
131  << " 4D segments" << endl;
132 
133  // Loop over the rechits of this ChamerId
134  for (DTRecSegment4DCollection::const_iterator segment4D = range.first; segment4D!=range.second; ++segment4D){
135 
136  //FOR NOISY CHANNELS////////////////////////////////
137  bool segmNoisy = false;
138  if((*segment4D).hasPhi()){
139  const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
140  vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
141  map<DTSuperLayerId,vector<DTRecHit1D> > hitsBySLMap;
142  for(vector<DTRecHit1D>::const_iterator hit = phiHits.begin(); hit != phiHits.end(); ++hit) {
143  DTWireId wireId = (*hit).wireId();
144 
145  // Check for noisy channels to skip them
146  if(checkNoisyChannels) {
147  bool isNoisy = false;
148  bool isFEMasked = false;
149  bool isTDCMasked = false;
150  bool isTrigMask = false;
151  bool isDead = false;
152  bool isNohv = false;
153  statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
154  if(isNoisy) {
155  if(debug)
156  cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
157  segmNoisy = true;
158  }
159  }
160  }
161  }
162 
163  if((*segment4D).hasZed()) {
164  const DTSLRecSegment2D* zSeg = (*segment4D).zSegment(); // zSeg lives in the SL RF
165  // Check for noisy channels to skip them
166  vector<DTRecHit1D> zHits = zSeg->specificRecHits();
167  for(vector<DTRecHit1D>::const_iterator hit = zHits.begin();
168  hit != zHits.end(); ++hit) {
169  DTWireId wireId = (*hit).wireId();
170  if(checkNoisyChannels) {
171  bool isNoisy = false;
172  bool isFEMasked = false;
173  bool isTDCMasked = false;
174  bool isTrigMask = false;
175  bool isDead = false;
176  bool isNohv = false;
177  //cout<<"wire id "<<wireId<<endl;
178  statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
179  if(isNoisy) {
180  if(debug)
181  cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
182  segmNoisy = true;
183  }
184  }
185  }
186  }
187 
188  if (segmNoisy) {
189  if(debug)
190  cout<<"skipping the segment: it contains noisy cells"<<endl;
191  continue;
192  }
193  //END FOR NOISY CHANNELS////////////////////////////////
194 
195  // Fill the histos
196  int nHits=0;
197  if((*segment4D).hasPhi()){
198  nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
199  if(debug)
200  cout<<"Phi segment with number of hits: "<<nHits<<endl;
201  phiHistos[0]->Fill((*chamberId).wheel(), nHits);
202  phiHistos[1]->Fill((*chamberId).sector(), nHits);
203  phiHistos[2]->Fill((*chamberId).station(), nHits);
204  }
205  if((*segment4D).hasZed()) {
206  nHits = (((*segment4D).zSegment())->specificRecHits()).size();
207  if(debug)
208  cout<<"Zed segment with number of hits: "<<nHits<<endl;
209  thetaHistos[0]->Fill((*chamberId).wheel(), nHits);
210  thetaHistos[1]->Fill((*chamberId).sector(), nHits);
211  thetaHistos[2]->Fill((*chamberId).station(), nHits);
212  }
213 
214  } //loop over segments
215  } // loop over chambers
216 }
findQualityFiles.size
size
Write out results.
Definition: findQualityFiles.py:442
edm::RangeMap< DTChamberId, edm::OwnVector< DTRecSegment4D > >::range
std::pair< const_iterator, const_iterator > range
iterator range
Definition: RangeMap.h:50
DTSegmentsTask::theRecHits4DLabel_
edm::EDGetTokenT< DTRecSegment4DCollection > theRecHits4DLabel_
Definition: DTSegmentsTask.h:54
SoftLeptonByDistance_cfi.distance
distance
Definition: SoftLeptonByDistance_cfi.py:6
DTSegmentsTask::thetaHistos
std::vector< MonitorElement * > thetaHistos
Definition: DTSegmentsTask.h:51
DTSegmentsTask::phiHistos
std::vector< MonitorElement * > phiHistos
Definition: DTSegmentsTask.h:50
DTRecSegment2D::specificRecHits
std::vector< DTRecHit1D > specificRecHits() const
Access to specific components.
Definition: DTRecSegment2D.cc:110
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:58
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:100
edm::OwnVector::const_iterator
Definition: OwnVector.h:41
gather_cfg.cout
cout
Definition: gather_cfg.py:145
void DTSegmentsTask::bookHistograms ( DQMStore::IBooker &  ibooker,
edm::Run const &  ,
edm::EventSetup const &   
)
overridevirtual

Implements DQMEDAnalyzer.

Definition at line 45 of file DTSegmentsTask.cc.

References DQMStore::IBooker::book2D(), DQMStore::IBooker::cd(), and DQMStore::IBooker::setCurrentFolder().

47  {
48 
49  ibooker.cd();
50  ibooker.setCurrentFolder("Muons/DTSegmentsMonitor");
51 
52  // histos for phi segments
53  phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsWheel", "phiSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
54  phiHistos[0]->setBinLabel(1,"W-2",1);
55  phiHistos[0]->setBinLabel(2,"W-1",1);
56  phiHistos[0]->setBinLabel(3,"W0",1);
57  phiHistos[0]->setBinLabel(4,"W1",1);
58  phiHistos[0]->setBinLabel(5,"W2",1);
59  phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsSector", "phiSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
60  phiHistos[1]->setBinLabel(1,"Sec1",1);
61  phiHistos[1]->setBinLabel(2,"Sec2",1);
62  phiHistos[1]->setBinLabel(3,"Sec3",1);
63  phiHistos[1]->setBinLabel(4,"Sec4",1);
64  phiHistos[1]->setBinLabel(5,"Sec5",1);
65  phiHistos[1]->setBinLabel(6,"Sec6",1);
66  phiHistos[1]->setBinLabel(7,"Sec7",1);
67  phiHistos[1]->setBinLabel(8,"Sec8",1);
68  phiHistos[1]->setBinLabel(9,"Sec9",1);
69  phiHistos[1]->setBinLabel(10,"Sec10",1);
70  phiHistos[1]->setBinLabel(11,"Sec11",1);
71  phiHistos[1]->setBinLabel(12,"Sec12",1);
72  phiHistos[1]->setBinLabel(13,"Sec13",1);
73  phiHistos[1]->setBinLabel(14,"Sec14",1);
74  phiHistos.push_back(ibooker.book2D("phiSegments_numHitsVsStation", "phiSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
75  phiHistos[2]->setBinLabel(1,"St1",1);
76  phiHistos[2]->setBinLabel(2,"St2",1);
77  phiHistos[2]->setBinLabel(3,"St3",1);
78  phiHistos[2]->setBinLabel(4,"St4",1);
79 
80  // histos for theta segments
81  thetaHistos.push_back(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(ibooker.book2D("thetaSegments_numHitsVsSector", "thetaSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
88  thetaHistos[1]->setBinLabel(1,"Sec1",1);
89  thetaHistos[1]->setBinLabel(2,"Sec2",1);
90  thetaHistos[1]->setBinLabel(3,"Sec3",1);
91  thetaHistos[1]->setBinLabel(4,"Sec4",1);
92  thetaHistos[1]->setBinLabel(5,"Sec5",1);
93  thetaHistos[1]->setBinLabel(6,"Sec6",1);
94  thetaHistos[1]->setBinLabel(7,"Sec7",1);
95  thetaHistos[1]->setBinLabel(8,"Sec8",1);
96  thetaHistos[1]->setBinLabel(9,"Sec9",1);
97  thetaHistos[1]->setBinLabel(10,"Sec10",1);
98  thetaHistos[1]->setBinLabel(11,"Sec11",1);
99  thetaHistos[1]->setBinLabel(12,"Sec12",1);
100  thetaHistos[1]->setBinLabel(13,"Sec13",1);
101  thetaHistos[1]->setBinLabel(14,"Sec14",1);
102  thetaHistos.push_back(ibooker.book2D("thetaSegments_numHitsVsStation", "thetaSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
103  thetaHistos[2]->setBinLabel(1,"St1",1);
104  thetaHistos[2]->setBinLabel(2,"St2",1);
105  thetaHistos[2]->setBinLabel(3,"St3",1);
106  thetaHistos[2]->setBinLabel(4,"St4",1);
107 
108 }
DTSegmentsTask::thetaHistos
std::vector< MonitorElement * > thetaHistos
Definition: DTSegmentsTask.h:51
DTSegmentsTask::phiHistos
std::vector< MonitorElement * > phiHistos
Definition: DTSegmentsTask.h:50
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:279
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:136

Member Data Documentation

bool DTSegmentsTask::checkNoisyChannels
private

Definition at line 46 of file DTSegmentsTask.h.

bool DTSegmentsTask::debug
private

Definition at line 45 of file DTSegmentsTask.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), and pkg.AbstractPkg::write().

edm::ParameterSet DTSegmentsTask::parameters
private

Definition at line 47 of file DTSegmentsTask.h.

Referenced by Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::inputTags(), Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::properties(), and Vispa.Plugins.ConfigEditor.ConfigDataAccessor.ConfigDataAccessor::recursePSetProperties().

std::vector<MonitorElement*> DTSegmentsTask::phiHistos
private

Definition at line 50 of file DTSegmentsTask.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegmentsTask::theRecHits4DLabel_
private

Definition at line 54 of file DTSegmentsTask.h.

std::vector<MonitorElement*> DTSegmentsTask::thetaHistos
private

Definition at line 51 of file DTSegmentsTask.h.

Generated for CMSSW Reference Manual by doxygen 1.8.11