#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 42 of file DTSegmentAnalysisTask.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 33 of file DTSegmentAnalysisTask.cc.

References checkNoisyChannels, detailedAnalysis, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltDQMMode, HLT_FULL_cff::InputTag, nhitsCut, phiSegmCut, 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.getParameter<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 ( )

override

Destructor.

Definition at line 52 of file DTSegmentAnalysisTask.cc.

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

Member Function Documentation

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 105 of file DTSegmentAnalysisTask.cc.

References DTStatusFlag::cellStatus(), checkNoisyChannels, HLT_FULL_cff::distance, edm::EventID::event(), dqm::impl::MonitorElement::Fill(), fillHistos(), edm::EventSetup::getData(), edm::EventBase::id(), edm::HandleBase::isValid(), nEventMonitor, nevents, nHits, nhitsCut, phiSegmCut, Pi, sistrip::SpyUtilities::range(), recHits4DToken_, 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());
     }
   }
 
   // -----------------------------------------------------------------------------
 }

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 62 of file DTSegmentAnalysisTask.cc.

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

                                                                          {
   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());
   }
 }

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

private

Definition at line 243 of file DTSegmentAnalysisTask.cc.

References dqm::implementation::IBooker::book1D(), detailedAnalysis, mergeVDriftHistosByStation::histos, histosPerCh, DTChamberId::sector(), dqm::implementation::NavigatorBase::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, topHistoFolder, and DTChamberId::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;
 }

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

overridevirtual

BeginRun.

Reimplemented from DQMEDAnalyzer.

Definition at line 57 of file DTSegmentAnalysisTask.cc.

References dtGeom, edm::EventSetup::getData(), and muonGeomToken_.

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

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

private

Definition at line 269 of file DTSegmentAnalysisTask.cc.

References detailedAnalysis, mergeVDriftHistosByStation::histos, histosPerCh, 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);
   }
 }