#include <DTSegmentAnalysisTask.h>

Inheritance diagram for DTSegmentAnalysisTask:

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &eSetup) override
	Summary. More...

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

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

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &eSetup) override

	~DTSegmentAnalysisTask () override
	Destructor. More...

Public Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

	DQMOneEDAnalyzer ()

void	endRun (edm::Run const &, edm::EventSetup const &) final

void	endRunProduce (edm::Run &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::Accumulator, Args... >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

Protected Member Functions inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
virtual void	analyze (edm::Event const &, edm::EventSetup const &)

virtual void	dqmBeginRun (edm::Run const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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)

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

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

Private Attributes
bool	checkNoisyChannels

bool	detailedAnalysis

edm::ESHandle< DTGeometry >	dtGeom

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

std::map< int, std::map< int, DTTimeEvolutionHisto * > >	histoTimeEvol

bool	hltDQMMode

DTTimeEvolutionHisto *	hNevtPerLS

MonitorElement *	nEventMonitor

int	nevents

int	nEventsInLS

int	nhitsCut

int	nLSTimeBin

int	nTimeBins

edm::ParameterSet	parameters

double	phiSegmCut

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

bool	slideTimeBins

std::map< int, MonitorElement * >	summaryHistos

std::string	topHistoFolder

Additional Inherited Members
Public Types inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Attributes inherited from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >
edm::EDPutTokenT< DQMToken >	runToken_

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

Constructor & Destructor Documentation

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

Constructor.

Definition at line 38 of file DTSegmentAnalysisTask.cc.

References checkNoisyChannels, detailedAnalysis, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), hltDQMMode, HLT_2018_cff::InputTag, nhitsCut, nLSTimeBin, nTimeBins, phiSegmCut, recHits4DToken_, slideTimeBins, and topHistoFolder.

     : nevents(0), nEventsInLS(0), hNevtPerLS(nullptr) {
   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);
   // # of bins in the time histos
   nTimeBins = pset.getUntrackedParameter<int>("nTimeBins", 100);
   // # of LS per bin in the time histos
   nLSTimeBin = pset.getUntrackedParameter<int>("nLSTimeBin", 2);
   // switch on/off sliding bins in time histos
   slideTimeBins = pset.getUntrackedParameter<bool>("slideTimeBins", true);
   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 63 of file DTSegmentAnalysisTask.cc.

References hNevtPerLS.

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

Member Function Documentation

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

override

Definition at line 145 of file DTSegmentAnalysisTask.cc.

References DTStatusFlag::cellStatus(), checkNoisyChannels, HLT_2018_cff::distance, edm::EventID::event(), dqm::impl::MonitorElement::Fill(), fillHistos(), edm::EventSetup::get(), edm::EventBase::id(), edm::HandleBase::isValid(), nEventMonitor, nevents, nEventsInLS, nhitsCut, phiSegmCut, Pi, FastTimerService_cff::range, recHits4DToken_, findQualityFiles::size, DTRecSegment2D::specificRecHits(), and xdir.

                                                                                      {
   nevents++;
   nEventMonitor->Fill(nevents);
 
   nEventsInLS++;
   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();
 
   ESHandle<DTStatusFlag> statusMap;
   if (checkNoisyChannels) {
     setup.get<DTStatusFlagRcd>().get(statusMap);
   }
 
   // -- 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::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	eSetup
	)

override

Summary.

Definition at line 339 of file DTSegmentAnalysisTask.cc.

References nEventsInLS.

                                                                                                          {
   nEventsInLS = 0;
 }

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

overrideprotectedvirtual

Implements DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 74 of file DTSegmentAnalysisTask.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), dqm::dqmstoreimpl::DQMStore::IBooker::bookFloat(), bookHistos(), DTGeometry::chambers(), chambers, dtGeom, HltBtagPostValidation_cff::histoName, histoTimeEvol, hltDQMMode, hNevtPerLS, nEventMonitor, nLSTimeBin, nTimeBins, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), slideTimeBins, 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());
   }
 
   // book sector time-evolution histos
   int modeTimeHisto = 0;
   if (!slideTimeBins)
     modeTimeHisto = 1;
   for (int wheel = -2; wheel != 3; ++wheel) {       // loop over wheels
     for (int sector = 1; sector <= 12; ++sector) {  // loop over sectors
 
       stringstream wheelstr;
       wheelstr << wheel;
       stringstream sectorstr;
       sectorstr << sector;
       string sectorHistoName = "NSegmPerEvent_W" + wheelstr.str() + "_Sec" + sectorstr.str();
       string sectorHistoTitle = "# segm. W" + wheelstr.str() + " Sect." + sectorstr.str();
 
       ibooker.setCurrentFolder(topHistoFolder + "/Wheel" + wheelstr.str() + "/Sector" + sectorstr.str());
 
       histoTimeEvol[wheel][sector] = new DTTimeEvolutionHisto(
           ibooker, sectorHistoName, sectorHistoTitle, nTimeBins, nLSTimeBin, slideTimeBins, modeTimeHisto);
     }
   }
 
   if (hltDQMMode)
     ibooker.setCurrentFolder(topHistoFolder);
   else
     ibooker.setCurrentFolder("DT/EventInfo/");
 
   hNevtPerLS = new DTTimeEvolutionHisto(ibooker, "NevtPerLS", "# evt.", nTimeBins, nLSTimeBin, slideTimeBins, 2);
 }

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

private

Definition at line 285 of file DTSegmentAnalysisTask.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), detailedAnalysis, combine::histos, histosPerCh, DTChamberId::sector(), dqm::dqmstoreimpl::DQMStore::IBooker::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;
 }

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

override

BeginRun.

Definition at line 69 of file DTSegmentAnalysisTask.cc.

References dtGeom, and edm::EventSetup::get().

                                                                                     {
   // Get the DT Geometry
   context.get<MuonGeometryRecord>().get(dtGeom);
 }

void DTSegmentAnalysisTask::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	eSetup
	)

override

Definition at line 329 of file DTSegmentAnalysisTask.cc.

References histoTimeEvol, hNevtPerLS, edm::LuminosityBlockBase::luminosityBlock(), nEventsInLS, DTTimeEvolutionHisto::updateTimeSlot(), and makeMuonMisalignmentScenario::wheel.

                                                                                                        {
   hNevtPerLS->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
   // book sector time-evolution histos
   for (int wheel = -2; wheel != 3; ++wheel) {
     for (int sector = 1; sector <= 12; ++sector) {
       histoTimeEvol[wheel][sector]->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
     }
   }
 }

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

private

Definition at line 311 of file DTSegmentAnalysisTask.cc.

References detailedAnalysis, combine::histos, histosPerCh, histoTimeEvol, 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());
   histoTimeEvol[chamberId.wheel()][sector]->accumulateValueTimeSlot(1);
 
   vector<MonitorElement*> histos = histosPerCh[chamberId];
   histos[0]->Fill(nHits);
   if (detailedAnalysis) {
     histos[1]->Fill(chi2);
   }
 }