DTSegmentsTask Class Reference

#include <DTSegmentsTask.h>

Inheritance diagram for DTSegmentsTask:

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	accumulate (edm::Event const &event, edm::EventSetup const &setup) final
void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final
void	beginStream (edm::StreamID id) final
virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)
	DQMEDAnalyzer ()
void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final
void	endRun (edm::Run const &run, edm::EventSetup const &setup) final
virtual bool	getCanSaveByLumi ()
Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

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 DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore
typedef dqm::reco::MonitorElement	MonitorElement
Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
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
Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)
static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)
static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)
static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)
Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const
Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_
edm::EDPutTokenT< DQMToken >	runToken_
unsigned int	streamId_

Detailed Description

DQM Analysis of 4D DT segments

Author

G. Mila - INFN Torino

Definition at line 25 of file DTSegmentsTask.h.

Constructor & Destructor Documentation

DTSegmentsTask()

DTSegmentsTask::DTSegmentsTask

(

const edm::ParameterSet &

pset

)

Constructor.

Definition at line 27 of file DTSegmentsTask.cc.

                                                          {
  debug = pset.getUntrackedParameter<bool>("debug", false);
  parameters = pset;
 
  // should be init from pset, but it was commented out...
  // better false than undefined
  checkNoisyChannels = false;
 
  // the name of the 4D rec hits collection
  theRecHits4DLabel_ = consumes<DTRecSegment4DCollection>(parameters.getParameter<std::string>("recHits4DLabel"));
}

References dtDigiTask_cfi::checkNoisyChannels, debug, muonDTDigis_cfi::pset, and AlCaHLTBitMon_QueryRunRegistry::string.

~DTSegmentsTask()

DTSegmentsTask::~DTSegmentsTask ( )

override

Destructor.

Definition at line 39 of file DTSegmentsTask.cc.

{}

Member Function Documentation

analyze()

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

overridevirtual

book the histos

Reimplemented from DQMEDAnalyzer.

Definition at line 110 of file DTSegmentsTask.cc.

                                                                              {
  // Get the map of noisy channels
  //  bool checkNoisyChannels = parameters.getUntrackedParameter<bool>("checkNoisyChannels",false);
  ESHandle<DTStatusFlag> statusMap;
  if (checkNoisyChannels) {
    setup.get<DTStatusFlagRcd>().get(statusMap);
  }
 
  // Get the 4D segment collection from the event
  edm::Handle<DTRecSegment4DCollection> all4DSegments;
  event.getByToken(theRecHits4DLabel_, all4DSegments);
 
  // Loop over all chambers containing a segment
  DTRecSegment4DCollection::id_iterator chamberId;
  for (chamberId = all4DSegments->id_begin(); chamberId != all4DSegments->id_end(); ++chamberId) {
    // Get the range for the corresponding ChamerId
    DTRecSegment4DCollection::range range = all4DSegments->get(*chamberId);
    int nsegm = distance(range.first, range.second);
    if (debug)
      cout << "   Chamber: " << *chamberId << " has " << nsegm << " 4D segments" << endl;
 
    // Loop over the rechits of this ChamerId
    for (DTRecSegment4DCollection::const_iterator segment4D = range.first; segment4D != range.second; ++segment4D) {
      //FOR NOISY CHANNELS////////////////////////////////
      bool segmNoisy = false;
      if ((*segment4D).hasPhi()) {
        const DTChamberRecSegment2D* phiSeg = (*segment4D).phiSegment();
        vector<DTRecHit1D> phiHits = phiSeg->specificRecHits();
        map<DTSuperLayerId, vector<DTRecHit1D> > hitsBySLMap;
        for (vector<DTRecHit1D>::const_iterator hit = phiHits.begin(); hit != phiHits.end(); ++hit) {
          DTWireId wireId = (*hit).wireId();
 
          // Check for noisy channels to skip them
          if (checkNoisyChannels) {
            bool isNoisy = false;
            bool isFEMasked = false;
            bool isTDCMasked = false;
            bool isTrigMask = false;
            bool isDead = false;
            bool isNohv = false;
            statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
            if (isNoisy) {
              if (debug)
                cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
              segmNoisy = true;
            }
          }
        }
      }
 
      if ((*segment4D).hasZed()) {
        const DTSLRecSegment2D* zSeg = (*segment4D).zSegment();  // zSeg lives in the SL RF
        // Check for noisy channels to skip them
        vector<DTRecHit1D> zHits = zSeg->specificRecHits();
        for (vector<DTRecHit1D>::const_iterator hit = zHits.begin(); hit != zHits.end(); ++hit) {
          DTWireId wireId = (*hit).wireId();
          if (checkNoisyChannels) {
            bool isNoisy = false;
            bool isFEMasked = false;
            bool isTDCMasked = false;
            bool isTrigMask = false;
            bool isDead = false;
            bool isNohv = false;
            //cout<<"wire id "<<wireId<<endl;
            statusMap->cellStatus(wireId, isNoisy, isFEMasked, isTDCMasked, isTrigMask, isDead, isNohv);
            if (isNoisy) {
              if (debug)
                cout << "Wire: " << wireId << " is noisy, skipping!" << endl;
              segmNoisy = true;
            }
          }
        }
      }
 
      if (segmNoisy) {
        if (debug)
          cout << "skipping the segment: it contains noisy cells" << endl;
        continue;
      }
      //END FOR NOISY CHANNELS////////////////////////////////
 
      // Fill the histos
      int nHits = 0;
      if ((*segment4D).hasPhi()) {
        nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
        if (debug)
          cout << "Phi segment with number of hits: " << nHits << endl;
        phiHistos[0]->Fill((*chamberId).wheel(), nHits);
        phiHistos[1]->Fill((*chamberId).sector(), nHits);
        phiHistos[2]->Fill((*chamberId).station(), nHits);
      }
      if ((*segment4D).hasZed()) {
        nHits = (((*segment4D).zSegment())->specificRecHits()).size();
        if (debug)
          cout << "Zed segment with number of hits: " << nHits << endl;
        thetaHistos[0]->Fill((*chamberId).wheel(), nHits);
        thetaHistos[1]->Fill((*chamberId).sector(), nHits);
        thetaHistos[2]->Fill((*chamberId).station(), nHits);
      }
 
    }  //loop over segments
  }    // loop over chambers
}

References DTStatusFlag::cellStatus(), dtDigiTask_cfi::checkNoisyChannels, gather_cfg::cout, debug, HLT_FULL_cff::distance, edm::get(), nHits, FastTimerService_cff::range, singleTopDQM_cfi::setup, findQualityFiles::size, and DTRecSegment2D::specificRecHits().

bookHistograms()

void DTSegmentsTask::bookHistograms ( DQMStore::IBooker &  ibooker, edm::Run const &  , edm::EventSetup const &    )

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 41 of file DTSegmentsTask.cc.

                                                            {
  ibooker.cd();
  ibooker.setCurrentFolder("Muons/DTSegmentsMonitor");
 
  // histos for phi segments
  phiHistos.push_back(
      ibooker.book2D("phiSegments_numHitsVsWheel", "phiSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
  phiHistos[0]->setBinLabel(1, "W-2", 1);
  phiHistos[0]->setBinLabel(2, "W-1", 1);
  phiHistos[0]->setBinLabel(3, "W0", 1);
  phiHistos[0]->setBinLabel(4, "W1", 1);
  phiHistos[0]->setBinLabel(5, "W2", 1);
  phiHistos.push_back(
      ibooker.book2D("phiSegments_numHitsVsSector", "phiSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
  phiHistos[1]->setBinLabel(1, "Sec1", 1);
  phiHistos[1]->setBinLabel(2, "Sec2", 1);
  phiHistos[1]->setBinLabel(3, "Sec3", 1);
  phiHistos[1]->setBinLabel(4, "Sec4", 1);
  phiHistos[1]->setBinLabel(5, "Sec5", 1);
  phiHistos[1]->setBinLabel(6, "Sec6", 1);
  phiHistos[1]->setBinLabel(7, "Sec7", 1);
  phiHistos[1]->setBinLabel(8, "Sec8", 1);
  phiHistos[1]->setBinLabel(9, "Sec9", 1);
  phiHistos[1]->setBinLabel(10, "Sec10", 1);
  phiHistos[1]->setBinLabel(11, "Sec11", 1);
  phiHistos[1]->setBinLabel(12, "Sec12", 1);
  phiHistos[1]->setBinLabel(13, "Sec13", 1);
  phiHistos[1]->setBinLabel(14, "Sec14", 1);
  phiHistos.push_back(
      ibooker.book2D("phiSegments_numHitsVsStation", "phiSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
  phiHistos[2]->setBinLabel(1, "St1", 1);
  phiHistos[2]->setBinLabel(2, "St2", 1);
  phiHistos[2]->setBinLabel(3, "St3", 1);
  phiHistos[2]->setBinLabel(4, "St4", 1);
 
  // histos for theta segments
  thetaHistos.push_back(
      ibooker.book2D("thetaSegments_numHitsVsWheel", "thetaSegments_numHitsVsWheel", 5, -2.5, 2.5, 20, 0, 20));
  thetaHistos[0]->setBinLabel(1, "W-2", 1);
  thetaHistos[0]->setBinLabel(2, "W-1", 1);
  thetaHistos[0]->setBinLabel(3, "W0", 1);
  thetaHistos[0]->setBinLabel(4, "W1", 1);
  thetaHistos[0]->setBinLabel(5, "W2", 1);
  thetaHistos.push_back(
      ibooker.book2D("thetaSegments_numHitsVsSector", "thetaSegments_numHitsVsSector", 14, 0.5, 14.5, 20, 0, 20));
  thetaHistos[1]->setBinLabel(1, "Sec1", 1);
  thetaHistos[1]->setBinLabel(2, "Sec2", 1);
  thetaHistos[1]->setBinLabel(3, "Sec3", 1);
  thetaHistos[1]->setBinLabel(4, "Sec4", 1);
  thetaHistos[1]->setBinLabel(5, "Sec5", 1);
  thetaHistos[1]->setBinLabel(6, "Sec6", 1);
  thetaHistos[1]->setBinLabel(7, "Sec7", 1);
  thetaHistos[1]->setBinLabel(8, "Sec8", 1);
  thetaHistos[1]->setBinLabel(9, "Sec9", 1);
  thetaHistos[1]->setBinLabel(10, "Sec10", 1);
  thetaHistos[1]->setBinLabel(11, "Sec11", 1);
  thetaHistos[1]->setBinLabel(12, "Sec12", 1);
  thetaHistos[1]->setBinLabel(13, "Sec13", 1);
  thetaHistos[1]->setBinLabel(14, "Sec14", 1);
  thetaHistos.push_back(
      ibooker.book2D("thetaSegments_numHitsVsStation", "thetaSegments_numHitsVsStation", 4, 0.5, 4.5, 20, 0, 20));
  thetaHistos[2]->setBinLabel(1, "St1", 1);
  thetaHistos[2]->setBinLabel(2, "St2", 1);
  thetaHistos[2]->setBinLabel(3, "St3", 1);
  thetaHistos[2]->setBinLabel(4, "St4", 1);
}

References dqm::implementation::IBooker::book2D(), dqm::implementation::NavigatorBase::cd(), and dqm::implementation::NavigatorBase::setCurrentFolder().

Member Data Documentation

checkNoisyChannels

bool DTSegmentsTask::checkNoisyChannels

private

Definition at line 41 of file DTSegmentsTask.h.

debug

bool DTSegmentsTask::debug

private

Definition at line 40 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(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder::runTauID(), and pkg.AbstractPkg::write().

parameters

edm::ParameterSet DTSegmentsTask::parameters

private

Definition at line 42 of file DTSegmentsTask.h.

phiHistos

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

private

Definition at line 45 of file DTSegmentsTask.h.

theRecHits4DLabel_

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegmentsTask::theRecHits4DLabel_

private

Definition at line 49 of file DTSegmentsTask.h.

thetaHistos

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

private

Definition at line 46 of file DTSegmentsTask.h.