#include <DTSegmentAnalysisTask.h>

Inheritance diagram for DTSegmentAnalysisTask:

Public Member Functions
void	analyze (const edm::Event &event, const edm::EventSetup &setup) override

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override
	BeginRun. More...

	DTSegmentAnalysisTask (const edm::ParameterSet &pset)
	Constructor. More...

	~DTSegmentAnalysisTask () 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

	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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Protected Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Private Member Functions
void	bookHistos (DQMStore::IBooker &ibooker, DTChamberId chamberId)

void	fillHistos (DTChamberId chamberId, int nHits, float chi2)

Private Attributes
bool	checkNoisyChannels

bool	detailedAnalysis

const DTGeometry *	dtGeom

std::map< DTChamberId, std::vector< MonitorElement * > >	histosPerCh

bool	hltDQMMode

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	muonGeomToken_

MonitorElement *	nEventMonitor

int	nevents

int	nhitsCut

double	phiSegmCut

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

const DTStatusFlag *	statusMap

edm::ESGetToken< DTStatusFlag, DTStatusFlagRcd >	statusMapToken_

std::map< int, MonitorElement * >	summaryHistos

std::string	topHistoFolder

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::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 Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

DQM Analysis of 4D DT segments, it produces plots about:

number of segments per event
number of hits per segment
position of the segments in chamber RF
direction of the segments (theta and phi projections)
reduced chi-square
All histos are produce per Chamber

Author: G. Cerminara - INFN Torino

Definition at line 41 of file DTSegmentAnalysisTask.h.

Constructor & Destructor Documentation

◆ DTSegmentAnalysisTask()

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

Constructor.

Definition at line 33 of file DTSegmentAnalysisTask.cc.

References checkNoisyChannels, detailedAnalysis, hltDQMMode, ProducerED_cfi::InputTag, nhitsCut, phiSegmCut, muonDTDigis_cfi::pset, recHits4DToken_, and topHistoFolder.

     : muonGeomToken_(esConsumes<edm::Transition::BeginRun>()), statusMapToken_(esConsumes()), nevents(0) {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Constructor called!";
 
   // switch for detailed analysis
   detailedAnalysis = pset.getUntrackedParameter<bool>("detailedAnalysis", false);
   // the name of the 4D rec hits collection
   recHits4DToken_ =
       consumes<DTRecSegment4DCollection>(edm::InputTag(pset.getUntrackedParameter<string>("recHits4DLabel")));
   // Get the map of noisy channels
   checkNoisyChannels = pset.getUntrackedParameter<bool>("checkNoisyChannels", false);
   phiSegmCut = pset.getUntrackedParameter<double>("phiSegmCut", 30.);
   nhitsCut = pset.getUntrackedParameter<int>("nhitsCut", 12);
 
   // top folder for the histograms in DQMStore
   topHistoFolder = pset.getUntrackedParameter<string>("topHistoFolder", "DT/02-Segments");
   // hlt DQM mode
   hltDQMMode = pset.getUntrackedParameter<bool>("hltDQMMode", false);
 }

◆ ~DTSegmentAnalysisTask()

DTSegmentAnalysisTask::~DTSegmentAnalysisTask ( )

override

Destructor.

Definition at line 53 of file DTSegmentAnalysisTask.cc.

                                               {
   //FR moved fron endjob
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "[DTSegmentAnalysisTask] Destructor called!";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 106 of file DTSegmentAnalysisTask.cc.

References DTStatusFlag::cellStatus(), checkNoisyChannels, HLT_2023v12_cff::distance, dqm::impl::MonitorElement::Fill(), fillHistos(), edm::HandleBase::isValid(), nEventMonitor, nevents, nHits, nhitsCut, phiSegmCut, Pi, FastTimerService_cff::range, recHits4DToken_, singleTopDQM_cfi::setup, findQualityFiles::size, DTRecSegment2D::specificRecHits(), statusMap, statusMapToken_, and xdir.

                                                                                      {
   nevents++;
   nEventMonitor->Fill(nevents);
 
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
       << "[DTSegmentAnalysisTask] Analyze #Run: " << event.id().run() << " #Event: " << event.id().event();
   if (!(event.id().event() % 1000))
     edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
         << "[DTSegmentAnalysisTask] Analyze #Run: " << event.id().run() << " #Event: " << event.id().event();
 
   if (checkNoisyChannels) {
     statusMap = &setup.getData(statusMapToken_);
   }
 
   // -- 4D segment analysis  -----------------------------------------------------
 
   // Get the 4D segment collection from the event
   edm::Handle<DTRecSegment4DCollection> all4DSegments;
   event.getByToken(recHits4DToken_, all4DSegments);
 
   if (!all4DSegments.isValid())
     return;
 
   // 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);
 
     edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
         << "   Chamber: " << *chamberId << " has " << distance(range.first, range.second) << " 4D segments";
 
     // 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 (checkNoisyChannels) {
         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
             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) {
               edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
                   << "Wire: " << wireId << " is noisy, skipping!";
               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();
             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) {
               edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
                   << "Wire: " << wireId << " is noisy, skipping!";
               segmNoisy = true;
             }
           }
         }
 
       }  // end of switch on noisy channels
       if (segmNoisy) {
         edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask")
             << "skipping the segment: it contains noisy cells";
         continue;
       }
       //END FOR NOISY CHANNELS////////////////////////////////
 
       int nHits = 0;
       if ((*segment4D).hasPhi())
         nHits = (((*segment4D).phiSegment())->specificRecHits()).size();
       if ((*segment4D).hasZed())
         nHits = nHits + ((((*segment4D).zSegment())->specificRecHits()).size());
 
       double anglePhiSegm(0.);
       if ((*segment4D).hasPhi()) {
         double xdir = (*segment4D).phiSegment()->localDirection().x();
         double zdir = (*segment4D).phiSegment()->localDirection().z();
 
         anglePhiSegm = atan(xdir / zdir) * 180. / TMath::Pi();
       }
       if (fabs(anglePhiSegm) > phiSegmCut)
         continue;
       // If the segment is in Wh+-2/SecX/MB1, get the DT chambers just above and check if there is a segment
       // to validate the segment present in MB1
       if (fabs((*chamberId).wheel()) == 2 && (*chamberId).station() == 1) {
         bool segmOk = false;
         int mb(2);
         while (mb < 4) {
           DTChamberId checkMB((*chamberId).wheel(), mb, (*chamberId).sector());
           DTRecSegment4DCollection::range ckrange = all4DSegments->get(checkMB);
 
           for (DTRecSegment4DCollection::const_iterator cksegment4D = ckrange.first; cksegment4D != ckrange.second;
                ++cksegment4D) {
             int nHits = 0;
             if ((*cksegment4D).hasPhi())
               nHits = (((*cksegment4D).phiSegment())->specificRecHits()).size();
             if ((*cksegment4D).hasZed())
               nHits = nHits + ((((*cksegment4D).zSegment())->specificRecHits()).size());
 
             if (nHits >= nhitsCut)
               segmOk = true;
           }
           mb++;
         }
 
         if (!segmOk)
           continue;
       }
       fillHistos(*chamberId, nHits, (*segment4D).chi2() / (*segment4D).degreesOfFreedom());
     }
   }
 
   // -----------------------------------------------------------------------------
 }

◆ bookHistograms()

void DTSegmentAnalysisTask::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	context
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 63 of file DTSegmentAnalysisTask.cc.

References dqm::implementation::IBooker::book2D(), dqm::implementation::IBooker::bookFloat(), bookHistos(), DTGeometry::chambers(), chambers, dtGeom, HltBtagPostValidation_cff::histoName, hltDQMMode, nEventMonitor, dqm::implementation::NavigatorBase::setCurrentFolder(), summaryHistos, topHistoFolder, and makeMuonMisalignmentScenario::wheel.

                                                                          {
   if (!hltDQMMode) {
     ibooker.setCurrentFolder("DT/EventInfo/Counters");
     nEventMonitor = ibooker.bookFloat("nProcessedEventsSegment");
   }
 
   for (int wh = -2; wh <= 2; wh++) {
     stringstream wheel;
     wheel << wh;
     ibooker.setCurrentFolder(topHistoFolder + "/Wheel" + wheel.str());
     string histoName = "numberOfSegments_W" + wheel.str();
 
     summaryHistos[wh] = ibooker.book2D(histoName.c_str(), histoName.c_str(), 12, 1, 13, 4, 1, 5);
     summaryHistos[wh]->setAxisTitle("Sector", 1);
     summaryHistos[wh]->setBinLabel(1, "1", 1);
     summaryHistos[wh]->setBinLabel(2, "2", 1);
     summaryHistos[wh]->setBinLabel(3, "3", 1);
     summaryHistos[wh]->setBinLabel(4, "4", 1);
     summaryHistos[wh]->setBinLabel(5, "5", 1);
     summaryHistos[wh]->setBinLabel(6, "6", 1);
     summaryHistos[wh]->setBinLabel(7, "7", 1);
     summaryHistos[wh]->setBinLabel(8, "8", 1);
     summaryHistos[wh]->setBinLabel(9, "9", 1);
     summaryHistos[wh]->setBinLabel(10, "10", 1);
     summaryHistos[wh]->setBinLabel(11, "11", 1);
     summaryHistos[wh]->setBinLabel(12, "12", 1);
     summaryHistos[wh]->setBinLabel(1, "MB1", 2);
     summaryHistos[wh]->setBinLabel(2, "MB2", 2);
     summaryHistos[wh]->setBinLabel(3, "MB3", 2);
     summaryHistos[wh]->setBinLabel(4, "MB4", 2);
   }
 
   // loop over all the DT chambers & book the histos
   const vector<const DTChamber*>& chambers = dtGeom->chambers();
   vector<const DTChamber*>::const_iterator ch_it = chambers.begin();
   vector<const DTChamber*>::const_iterator ch_end = chambers.end();
   for (; ch_it != ch_end; ++ch_it) {
     bookHistos(ibooker, (*ch_it)->id());
   }
 }

◆ bookHistos()

void DTSegmentAnalysisTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		DTChamberId	chamberId
	)

private

Definition at line 244 of file DTSegmentAnalysisTask.cc.

References dqm::implementation::IBooker::book1D(), detailedAnalysis, combine::histos, histosPerCh, hgcalTBTopologyTester_cfi::sector, DTChamberId::sector(), dqm::implementation::NavigatorBase::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, topHistoFolder, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                                       {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTSegmentAnalysisTask") << "   Booking histos for chamber: " << chamberId;
 
   // Compose the chamber name
   stringstream wheel;
   wheel << chamberId.wheel();
   stringstream station;
   station << chamberId.station();
   stringstream sector;
   sector << chamberId.sector();
 
   string chamberHistoName = "_W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str();
 
   ibooker.setCurrentFolder(topHistoFolder + "/Wheel" + wheel.str() + "/Sector" + sector.str() + "/Station" +
                            station.str());
 
   // Create the monitor elements
   vector<MonitorElement*> histos;
   histos.push_back(ibooker.book1D("h4DSegmNHits" + chamberHistoName, "# of hits per segment", 16, 0.5, 16.5));
   if (detailedAnalysis) {
     histos.push_back(ibooker.book1D("h4DChi2" + chamberHistoName, "4D Segment reduced Chi2", 20, 0, 20));
   }
   histosPerCh[chamberId] = histos;
 }

◆ dqmBeginRun()

void DTSegmentAnalysisTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	context
	)

overridevirtual

BeginRun.

Reimplemented from DQMEDAnalyzer.

Definition at line 58 of file DTSegmentAnalysisTask.cc.

References visDQMUpload::context, dtGeom, and muonGeomToken_.

                                                                                     {
   // Get the DT Geometry
   dtGeom = &context.getData(muonGeomToken_);
 }

◆ fillHistos()

void DTSegmentAnalysisTask::fillHistos	(	DTChamberId	chamberId,
		int	nHits,
		float	chi2
	)

private

Definition at line 270 of file DTSegmentAnalysisTask.cc.

References hltPixelTracks_cff::chi2, detailedAnalysis, combine::histos, histosPerCh, nHits, hgcalTBTopologyTester_cfi::sector, DTChamberId::sector(), DTChamberId::station(), summaryHistos, and DTChamberId::wheel().

Referenced by analyze().

                                                                                    {
   int sector = chamberId.sector();
   if (chamberId.sector() == 13) {
     sector = 4;
   } else if (chamberId.sector() == 14) {
     sector = 10;
   }
 
   summaryHistos[chamberId.wheel()]->Fill(sector, chamberId.station());
 
   vector<MonitorElement*> histos = histosPerCh[chamberId];
   histos[0]->Fill(nHits);
   if (detailedAnalysis) {
     histos[1]->Fill(chi2);
   }
 }

Member Data Documentation

◆ checkNoisyChannels

bool DTSegmentAnalysisTask::checkNoisyChannels

private

Definition at line 75 of file DTSegmentAnalysisTask.h.

Referenced by analyze(), and DTSegmentAnalysisTask().

◆ detailedAnalysis

bool DTSegmentAnalysisTask::detailedAnalysis

private

Definition at line 61 of file DTSegmentAnalysisTask.h.

Referenced by bookHistos(), DTSegmentAnalysisTask(), and fillHistos().

◆ dtGeom

const DTGeometry* DTSegmentAnalysisTask::dtGeom

private

Definition at line 65 of file DTSegmentAnalysisTask.h.

Referenced by bookHistograms(), and dqmBeginRun().

◆ histosPerCh

std::map<DTChamberId, std::vector<MonitorElement*> > DTSegmentAnalysisTask::histosPerCh

private

Definition at line 83 of file DTSegmentAnalysisTask.h.

Referenced by bookHistos(), and fillHistos().

◆ hltDQMMode

bool DTSegmentAnalysisTask::hltDQMMode

private

Definition at line 90 of file DTSegmentAnalysisTask.h.

Referenced by bookHistograms(), and DTSegmentAnalysisTask().

◆ muonGeomToken_

edm::ESGetToken<DTGeometry, MuonGeometryRecord> DTSegmentAnalysisTask::muonGeomToken_

private

Definition at line 64 of file DTSegmentAnalysisTask.h.

Referenced by dqmBeginRun().

◆ nEventMonitor

MonitorElement* DTSegmentAnalysisTask::nEventMonitor

private

Definition at line 96 of file DTSegmentAnalysisTask.h.

Referenced by analyze(), and bookHistograms().

◆ nevents

int DTSegmentAnalysisTask::nevents

private

Definition at line 86 of file DTSegmentAnalysisTask.h.

Referenced by analyze().

◆ nhitsCut

int DTSegmentAnalysisTask::nhitsCut

private

Definition at line 94 of file DTSegmentAnalysisTask.h.

Referenced by analyze(), and DTSegmentAnalysisTask().

◆ phiSegmCut

double DTSegmentAnalysisTask::phiSegmCut

private

Definition at line 92 of file DTSegmentAnalysisTask.h.

Referenced by analyze(), and DTSegmentAnalysisTask().

◆ recHits4DToken_

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegmentAnalysisTask::recHits4DToken_

private

Definition at line 72 of file DTSegmentAnalysisTask.h.

Referenced by analyze(), and DTSegmentAnalysisTask().

◆ statusMap

const DTStatusFlag* DTSegmentAnalysisTask::statusMap

private

Definition at line 69 of file DTSegmentAnalysisTask.h.

Referenced by analyze().

◆ statusMapToken_

edm::ESGetToken<DTStatusFlag, DTStatusFlagRcd> DTSegmentAnalysisTask::statusMapToken_

private

Definition at line 68 of file DTSegmentAnalysisTask.h.

Referenced by analyze().

◆ summaryHistos

std::map<int, MonitorElement*> DTSegmentAnalysisTask::summaryHistos

private

Definition at line 84 of file DTSegmentAnalysisTask.h.

Referenced by bookHistograms(), and fillHistos().

◆ topHistoFolder

std::string DTSegmentAnalysisTask::topHistoFolder

private

Definition at line 88 of file DTSegmentAnalysisTask.h.

Referenced by bookHistograms(), bookHistos(), and DTSegmentAnalysisTask().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTSegmentAnalysisTask()

◆ ~DTSegmentAnalysisTask()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ bookHistos()

◆ dqmBeginRun()

◆ fillHistos()

Member Data Documentation

◆ checkNoisyChannels

◆ detailedAnalysis

◆ dtGeom

◆ histosPerCh

◆ hltDQMMode

◆ muonGeomToken_

◆ nEventMonitor

◆ nevents

◆ nhitsCut

◆ phiSegmCut

◆ recHits4DToken_

◆ statusMap

◆ statusMapToken_

◆ summaryHistos

◆ topHistoFolder