#include <DTChamberEfficiencyTask.h>

Inheritance diagram for DTChamberEfficiencyTask:

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

void	beginLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context) override
	To reset the MEs. More...

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

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

void	endLuminosityBlock (edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context) final

	~DTChamberEfficiencyTask () override
	Destructor. More...

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

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

	EDProducer (const EDProducer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

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

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~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

std::vector< bool > const &	recordProvenanceList () 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)

TypeLabelList const &	typeLabelList () const
	used by the fwk to register the list of products of this module More...

	~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

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

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices 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	dqmEndRun (edm::Run const &, edm::EventSetup const &)

Protected Member Functions inherited from edm::ProducerBase
template<Transition Tr = Transition::Event>
auto	produces (std::string instanceName) noexcept
	declare what type of product will make and with which optional label More...

template<Transition B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<BranchType B>
BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

BranchAliasSetter	produces (const TypeID &id, std::string instanceName=std::string(), bool recordProvenance=true)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , BranchType B>
BranchAliasSetterT< ProductType >	produces ()

template<class ProductType >
BranchAliasSetterT< ProductType >	produces (std::string instanceName)

template<typename ProductType , Transition B>
BranchAliasSetterT< ProductType >	produces ()

template<Transition Tr = Transition::Event>
auto	produces () noexcept

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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

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

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

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

const DTRecSegment4D &	getBestSegment (const DTRecSegment4DCollection::range &segs) const

const DTRecSegment4D *	getBestSegment (const DTRecSegment4D s1, const DTRecSegment4D s2) const

LocalPoint	interpolate (const DTRecSegment4D &seg1, const DTRecSegment4D &seg3, const DTChamberId &MB2) const

bool	isGoodSegment (const DTRecSegment4D &seg) const

Private Attributes
bool	debug

bool	detailedAnalysis

const DTGeometry *	dtGeom

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

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	muonGeomToken_

bool	onlineMonitor

edm::ParameterSet	parameters

edm::EDGetTokenT< DTRecSegment4DCollection >	recHits4DToken_

edm::Handle< DTRecSegment4DCollection >	segs

double	theMinChi2NormSegment

double	theMinCloseDist

unsigned int	theMinHitsSegment

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
template<typename T >
using	BranchAliasSetterT = ProductRegistryHelper::BranchAliasSetterT< T >

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:

single chamber efficiency All histos are produced per Chamber

Class based on the code written by S. Lacaprara : RecoLocalMuon / DTSegment / test / DTEffAnalyzer.h

Author: G. Mila - INFN Torino

Definition at line 35 of file DTChamberEfficiencyTask.h.

Constructor & Destructor Documentation

◆ DTChamberEfficiencyTask()

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

Constructor.

Definition at line 23 of file DTChamberEfficiencyTask.cc.

References debug, detailedAnalysis, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), ProducerED_cfi::InputTag, onlineMonitor, parameters, muonDTDigis_cfi::pset, recHits4DToken_, theMinChi2NormSegment, theMinCloseDist, and theMinHitsSegment.

     : muonGeomToken_(esConsumes<edm::Transition::BeginRun>()) {
   debug = pset.getUntrackedParameter<bool>("debug", false);
 
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "[DTChamberEfficiencyTask] Constructor called!";
 
   parameters = pset;
 
   // the name of the 4D rec hits collection
   recHits4DToken_ =
       consumes<DTRecSegment4DCollection>(edm::InputTag(parameters.getUntrackedParameter<string>("recHits4DLabel")));
 
   // parameters to use for the segment quality check
   theMinHitsSegment = static_cast<unsigned int>(parameters.getParameter<int>("minHitsSegment"));
   theMinChi2NormSegment = parameters.getParameter<double>("minChi2NormSegment");
   // parameter to use for the exstrapolated segment check
   theMinCloseDist = parameters.getParameter<double>("minCloseDist");
 
   // the running modality
   onlineMonitor = parameters.getUntrackedParameter<bool>("onlineMonitor");
 
   // the analysis mode
   detailedAnalysis = parameters.getUntrackedParameter<bool>("detailedAnalysis");
 }

◆ ~DTChamberEfficiencyTask()

DTChamberEfficiencyTask::~DTChamberEfficiencyTask ( )

override

Destructor.

Definition at line 48 of file DTChamberEfficiencyTask.cc.

                                                   {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "[DTChamberEfficiencyTask] Destructor called!";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 145 of file DTChamberEfficiencyTask.cc.

References DTGeometry::chamber(), DTGeometry::chambers(), detailedAnalysis, dtGeom, getBestSegment(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), mergeVDriftHistosByStation::histos, histosPerCh, interpolate(), isGoodSegment(), DTRecSegment4D::localPosition(), mag(), recHits4DToken_, nano_mu_digi_cff::sector, DTChamberId::sector(), segs, DTChamberId::station(), relativeConstraints::station, theMinCloseDist, DTChamberId::wheel(), makeMuonMisalignmentScenario::wheel, x, PV3DBase< T, PVType, FrameType >::x(), y, and PV3DBase< T, PVType, FrameType >::y().

                                                                                        {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiencyTask")
       << "[DTChamberEfficiencyTask] Analyze #Run: " << event.id().run() << " #Event: " << event.id().event();
 
   // Get the 4D rechit collection from the event
   event.getByToken(recHits4DToken_, segs);
 
   int bottom = 0, top = 0;
 
   // Loop over all the chambers
   vector<const DTChamber*>::const_iterator ch_it = dtGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator ch_end = dtGeom->chambers().end();
   for (; ch_it != ch_end; ++ch_it) {
     DTChamberId ch = (*ch_it)->id();
     int wheel = ch.wheel();
     int sector = ch.sector();
     int station = ch.station();
 
     DTChamberId MidId(wheel, station, sector);
 
     // get efficiency for MB1 using MB2 and MB3
     if (station == 1) {
       bottom = 2;
       top = 3;
     }
 
     // get efficiency for MB2 using MB1 and MB3
     if (station == 2) {
       bottom = 1;
       top = 3;
     }
 
     // get efficiency for MB2 using MB2 and MB4
     if (station == 3) {
       bottom = 2;
       top = 4;
     }
 
     // get efficiency for MB4 using MB2 and MB3
     if (station == 4) {
       bottom = 2;
       top = 3;
     }
 
     // Select events with (good) segments in Bot and Top
     DTChamberId BotId(wheel, bottom, sector);
     DTChamberId TopId(wheel, top, sector);
 
     // Get segments in the bottom chambers (if any)
     DTRecSegment4DCollection::range segsBot = segs->get(BotId);
     int nSegsBot = segsBot.second - segsBot.first;
     // check if any segments is there
     if (nSegsBot == 0)
       continue;
 
     vector<MonitorElement*> histos = histosPerCh[MidId];
 
     // Get segments in the top chambers (if any)
     DTRecSegment4DCollection::range segsTop = segs->get(TopId);
     int nSegsTop = segsTop.second - segsTop.first;
 
     // Select one segment for the bottom chamber
     const DTRecSegment4D& bestBotSeg = getBestSegment(segsBot);
 
     // Select one segment for the top chamber
     DTRecSegment4D* pBestTopSeg = nullptr;
     if (nSegsTop > 0)
       pBestTopSeg = const_cast<DTRecSegment4D*>(&getBestSegment(segsTop));
     //if top chamber is MB4 sector 10, consider also sector 14
     if (TopId.station() == 4 && TopId.sector() == 10) {
       DTChamberId TopId14(wheel, top, 14);
       DTRecSegment4DCollection::range segsTop14 = segs->get(TopId14);
       int nSegsTop14 = segsTop14.second - segsTop14.first;
       nSegsTop += nSegsTop14;
       if (nSegsTop14) {
         DTRecSegment4D* pBestTopSeg14 = const_cast<DTRecSegment4D*>(&getBestSegment(segsTop14));
         // get best between sector 10 and 14
         pBestTopSeg = const_cast<DTRecSegment4D*>(getBestSegment(pBestTopSeg, pBestTopSeg14));
       }
     }
     if (!pBestTopSeg)
       continue;
     const DTRecSegment4D& bestTopSeg = *pBestTopSeg;
 
     DTRecSegment4DCollection::range segsMid = segs->get(MidId);
     int nSegsMid = segsMid.second - segsMid.first;
 
     if (detailedAnalysis) {
       // very trivial efficiency, just count segments
       histos[3]->Fill(0);
       if (nSegsMid > 0)
         histos[3]->Fill(1);
     }
 
     // get position at Mid by interpolating the position (not direction) of best
     // segment in Bot and Top to Mid surface
     LocalPoint posAtMid = interpolate(bestBotSeg, bestTopSeg, MidId);
 
     // is best segment good enough?
     if (isGoodSegment(bestBotSeg) && isGoodSegment(bestTopSeg)) {
       if (detailedAnalysis)
         histos[4]->Fill(posAtMid.x(), posAtMid.y());
       //check if interpolation fall inside middle chamber
       if ((dtGeom->chamber(MidId))->surface().bounds().inside(posAtMid)) {
         histos[0]->Fill(posAtMid.x(), posAtMid.y());
         if (nSegsMid > 0) {
           if (detailedAnalysis) {
             histos[3]->Fill(2);
             histos[5]->Fill(posAtMid.x(), posAtMid.y());
           }
 
           const DTRecSegment4D& bestMidSeg = getBestSegment(segsMid);
           // check if middle segments is good enough
           if (isGoodSegment(bestMidSeg)) {
             if (detailedAnalysis)
               histos[6]->Fill(posAtMid.x(), posAtMid.y());
             LocalPoint midSegPos = bestMidSeg.localPosition();
 
             // check if middle segments is also close enough
             double dist;
             if (bestMidSeg.hasPhi()) {
               if (bestTopSeg.hasZed() && bestBotSeg.hasZed() && bestMidSeg.hasZed()) {
                 dist = (midSegPos - posAtMid).mag();
               } else {
                 dist = fabs((midSegPos - posAtMid).x());
               }
             } else {
               dist = fabs((midSegPos - posAtMid).y());
             }
             if (dist < theMinCloseDist) {
               histos[1]->Fill(posAtMid.x(), posAtMid.y());
             }
             if (detailedAnalysis)
               histos[2]->Fill(dist);
           }
         }
       }
     }
   }  // loop over stations
 }

◆ beginLuminosityBlock()

void DTChamberEfficiencyTask::beginLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

override

To reset the MEs.

Definition at line 52 of file DTChamberEfficiencyTask.cc.

References edm::ParameterSet::getUntrackedParameter(), timingPdfMaker::histo, histosPerCh, mps_fire::i, edm::LuminosityBlockBase::id(), edm::LuminosityBlockID::luminosityBlock(), genParticles_cff::map, onlineMonitor, parameters, and findQualityFiles::size.

                                                                                                             {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiencyTask")
       << "[DTChamberEfficiencyTask]: Begin of LS transition";
 
   if (lumiSeg.id().luminosityBlock() % parameters.getUntrackedParameter<int>("ResetCycle", 3) == 0 && onlineMonitor) {
     for (map<DTChamberId, vector<MonitorElement*> >::const_iterator histo = histosPerCh.begin();
          histo != histosPerCh.end();
          histo++) {
       int size = (*histo).second.size();
       for (int i = 0; i < size; i++) {
         (*histo).second[i]->Reset();
       }
     }
   }
 }

◆ bookHistograms()

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

overrideprotectedvirtual

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

Definition at line 73 of file DTChamberEfficiencyTask.cc.

References bookHistos(), DTGeometry::chambers(), dtGeom, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                            {
   ibooker.setCurrentFolder("DT/DTChamberEfficiencyTask");
 
   // Loop over all the chambers
   vector<const DTChamber*>::const_iterator ch_it = dtGeom->chambers().begin();
   vector<const DTChamber*>::const_iterator ch_end = dtGeom->chambers().end();
   for (; ch_it != ch_end; ++ch_it) {
     // histo booking
     bookHistos(ibooker, (*ch_it)->id());
   }
 }

◆ bookHistos()

void DTChamberEfficiencyTask::bookHistos	(	DQMStore::IBooker &	ibooker,
		DTChamberId	chId
	)

private

Definition at line 88 of file DTChamberEfficiencyTask.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book2D(), detailedAnalysis, mergeVDriftHistosByStation::histos, histosPerCh, nano_mu_digi_cff::sector, DTChamberId::sector(), dqm::implementation::NavigatorBase::setCurrentFolder(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by bookHistograms().

                                                                                    {
   edm::LogVerbatim("DTDQM|DTMonitorModule|DTChamberEfficiencyTask") << "   Booking histos for CH : " << chId;
 
   // Compose the chamber name
   stringstream wheel;
   wheel << chId.wheel();
   stringstream station;
   station << chId.station();
   stringstream sector;
   sector << chId.sector();
 
   string HistoName = "_W" + wheel.str() + "_St" + station.str() + "_Sec" + sector.str();
 
   ibooker.setCurrentFolder("DT/01-DTChamberEfficiency/Task/Wheel" + wheel.str() + "/Sector" + sector.str() +
                            "/Station" + station.str());
 
   // Create the monitor elements
   vector<MonitorElement*> histos;
 
   //efficiency selection cuts
   // a- number of segments of the top chamber > 0 && number of segments of the bottom chamber > 0
   // b- number of segments of the middle chamber > 0
   // c- check of the top and bottom segment quality
   // d- check if interpolation falls inside the middle chamber
   // e- check of the middle segment quality
   // f- check if the distance between the reconstructed and the exstrapolated segments is ok
 
   // histo for efficiency with cuts a-/c-/d-
   histos.push_back(ibooker.book2D(
       "hEffGoodSegVsPosDen" + HistoName, "Eff vs local position (good) ", 25, -250., 250., 25, -250., 250.));
   // histo for efficiency with cuts a-/b-/c-/d-/e-/f-
   histos.push_back(ibooker.book2D("hEffGoodCloseSegVsPosNum" + HistoName,
                                   "Eff vs local position (good and close segs) ",
                                   25,
                                   -250.,
                                   250.,
                                   25,
                                   -250.,
                                   250.));
   if (detailedAnalysis) {
     histos.push_back(
         ibooker.book1D("hDistSegFromExtrap" + HistoName, "Distance segments from extrap position ", 200, 0., 200.));
     // histo for efficiency from segment counting
     histos.push_back(ibooker.book1D("hNaiveEffSeg" + HistoName, "Naive eff ", 10, 0., 10.));
     // histo for efficiency with cuts a-/c-
     histos.push_back(ibooker.book2D(
         "hEffSegVsPosDen" + HistoName, "Eff vs local position (all) ", 25, -250., 250., 25, -250., 250.));
     // histo for efficiency with cuts a-/b-/c-/d-
     histos.push_back(
         ibooker.book2D("hEffSegVsPosNum" + HistoName, "Eff vs local position ", 25, -250., 250., 25, -250., 250.));
     // histo for efficiency with cuts a-/b-/c-/d-/e-
     histos.push_back(ibooker.book2D(
         "hEffGoodSegVsPosNum" + HistoName, "Eff vs local position (good segs) ", 25, -250., 250., 25, -250., 250.));
   }
   histosPerCh[chId] = histos;
 }

◆ dqmBeginRun()

void DTChamberEfficiencyTask::dqmBeginRun	(	const edm::Run &	run,
		const edm::EventSetup &	setup
	)

overridevirtual

BeginRun.

Reimplemented from DQMOneEDAnalyzer< edm::one::WatchLuminosityBlocks >.

Definition at line 68 of file DTChamberEfficiencyTask.cc.

References dtGeom, muonGeomToken_, and singleTopDQM_cfi::setup.

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

◆ endLuminosityBlock()

void DTChamberEfficiencyTask::endLuminosityBlock	(	edm::LuminosityBlock const &	lumiSeg,
		edm::EventSetup const &	context
	)

inlinefinal

Definition at line 48 of file DTChamberEfficiencyTask.h.

48 {}

◆ getBestSegment() [1/2]

const DTRecSegment4D & DTChamberEfficiencyTask::getBestSegment ( const DTRecSegment4DCollection::range & segs ) const

private

Definition at line 287 of file DTChamberEfficiencyTask.cc.

References nHits, and segs.

Referenced by analyze().

                                                                                                              {
   DTRecSegment4DCollection::const_iterator bestIter;
   unsigned int nHitBest = 0;
   double chi2Best = 99999.;
   for (DTRecSegment4DCollection::const_iterator seg = segs.first; seg != segs.second; ++seg) {
     unsigned int nHits = ((*seg).hasPhi() ? (*seg).phiSegment()->recHits().size() : 0);
     nHits += ((*seg).hasZed() ? (*seg).zSegment()->recHits().size() : 0);
 
     if (nHits == nHitBest) {
       if ((*seg).chi2() / (*seg).degreesOfFreedom() < chi2Best) {
         chi2Best = (*seg).chi2() / (*seg).degreesOfFreedom();
         bestIter = seg;
       }
     } else if (nHits > nHitBest) {
       nHitBest = nHits;
       bestIter = seg;
     }
   }
   return *bestIter;
 }

◆ getBestSegment() [2/2]

const DTRecSegment4D * DTChamberEfficiencyTask::getBestSegment	(	const DTRecSegment4D *	s1,
		const DTRecSegment4D *	s2
	)		const

private

Definition at line 308 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::chi2(), DTRecSegment4D::degreesOfFreedom(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), DTRecSegment4D::phiSegment(), DTRecSegment2D::recHits(), and DTRecSegment4D::zSegment().

                                                                                               {
   if (!s1)
     return s2;
   if (!s2)
     return s1;
   unsigned int nHits1 = (s1->hasPhi() ? s1->phiSegment()->recHits().size() : 0);
   nHits1 += (s1->hasZed() ? s1->zSegment()->recHits().size() : 0);
 
   unsigned int nHits2 = (s2->hasPhi() ? s2->phiSegment()->recHits().size() : 0);
   nHits2 += (s2->hasZed() ? s2->zSegment()->recHits().size() : 0);
 
   if (nHits1 == nHits2) {
     if (s1->chi2() / s1->degreesOfFreedom() < s2->chi2() / s2->degreesOfFreedom())
       return s1;
     else
       return s2;
   } else if (nHits1 > nHits2)
     return s1;
   return s2;
 }

◆ interpolate()

LocalPoint DTChamberEfficiencyTask::interpolate	(	const DTRecSegment4D &	seg1,
		const DTRecSegment4D &	seg3,
		const DTChamberId &	MB2
	)		const

private

Definition at line 330 of file DTChamberEfficiencyTask.cc.

References DTGeometry::chamber(), DTRecSegment4D::chamberId(), DeadROC_duringRun::dir, dtGeom, DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), globals_cff::id2, DTRecSegment4D::localPosition(), toLocal(), unit(), PV3DBase< T, PVType, FrameType >::x(), PV3DBase< T, PVType, FrameType >::y(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by analyze().

                                                                               {
   // Get GlobalPoition of Seg in MB1
   GlobalPoint gpos1 = (dtGeom->chamber(seg1.chamberId()))->toGlobal(seg1.localPosition());
 
   // Get GlobalPoition of Seg in MB3
   GlobalPoint gpos3 = (dtGeom->chamber(seg3.chamberId()))->toGlobal(seg3.localPosition());
 
   // interpolate
   // get all in MB2 frame
   LocalPoint pos1 = (dtGeom->chamber(id2))->toLocal(gpos1);
   LocalPoint pos3 = (dtGeom->chamber(id2))->toLocal(gpos3);
 
   // case 1: 1 and 3 has both projection. No problem
 
   // case 2: one projection is missing for one of the segments. Keep the other's segment position
   if (!seg1.hasZed())
     pos1 = LocalPoint(pos1.x(), pos3.y(), pos1.z());
   if (!seg3.hasZed())
     pos3 = LocalPoint(pos3.x(), pos1.y(), pos3.z());
 
   if (!seg1.hasPhi())
     pos1 = LocalPoint(pos3.x(), pos1.y(), pos1.z());
   if (!seg3.hasPhi())
     pos3 = LocalPoint(pos1.x(), pos3.y(), pos3.z());
 
   // direction
   LocalVector dir = (pos3 - pos1).unit();  // z points inward!
   LocalPoint pos2 = pos1 + dir * pos1.z() / (-dir.z());
 
   return pos2;
 }

◆ isGoodSegment()

bool DTChamberEfficiencyTask::isGoodSegment ( const DTRecSegment4D & seg ) const

private

Definition at line 364 of file DTChamberEfficiencyTask.cc.

References DTRecSegment4D::chamberId(), DTRecSegment4D::chi2(), DTRecSegment4D::degreesOfFreedom(), DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), nHits, DTRecSegment4D::phiSegment(), DTRecSegment2D::recHits(), DTChamberId::station(), theMinChi2NormSegment, theMinHitsSegment, and DTRecSegment4D::zSegment().

Referenced by analyze().

                                                                            {
   if (seg.chamberId().station() != 4 && !seg.hasZed())
     return false;
   unsigned int nHits = (seg.hasPhi() ? seg.phiSegment()->recHits().size() : 0);
   nHits += (seg.hasZed() ? seg.zSegment()->recHits().size() : 0);
   return (nHits >= theMinHitsSegment && seg.chi2() / seg.degreesOfFreedom() < theMinChi2NormSegment);
 }

Member Data Documentation

◆ debug

bool DTChamberEfficiencyTask::debug

private

Definition at line 66 of file DTChamberEfficiencyTask.h.

Referenced by util.rrapi.RRApi::dprint(), rrapi.RRApi::dprint(), DTChamberEfficiencyTask(), pkg.AbstractPkg::generate(), rrapi.RRApi::get(), util.rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::load_againstElectronMVA6(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), and pkg.AbstractPkg::write().

◆ detailedAnalysis

bool DTChamberEfficiencyTask::detailedAnalysis

private

Definition at line 70 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), bookHistos(), and DTChamberEfficiencyTask().

◆ dtGeom

const DTGeometry* DTChamberEfficiencyTask::dtGeom

private

Definition at line 85 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), bookHistograms(), dqmBeginRun(), and interpolate().

◆ histosPerCh

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

private

Definition at line 77 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), beginLuminosityBlock(), and bookHistos().

◆ muonGeomToken_

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

private

Definition at line 84 of file DTChamberEfficiencyTask.h.

Referenced by dqmBeginRun().

◆ onlineMonitor

bool DTChamberEfficiencyTask::onlineMonitor

private

Definition at line 68 of file DTChamberEfficiencyTask.h.

Referenced by beginLuminosityBlock(), and DTChamberEfficiencyTask().

◆ parameters

edm::ParameterSet DTChamberEfficiencyTask::parameters

private

Definition at line 75 of file DTChamberEfficiencyTask.h.

Referenced by beginLuminosityBlock(), and DTChamberEfficiencyTask().

◆ recHits4DToken_

edm::EDGetTokenT<DTRecSegment4DCollection> DTChamberEfficiencyTask::recHits4DToken_

private

Definition at line 73 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), and DTChamberEfficiencyTask().

◆ segs

edm::Handle<DTRecSegment4DCollection> DTChamberEfficiencyTask::segs

private

Definition at line 86 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), and getBestSegment().

◆ theMinChi2NormSegment

double DTChamberEfficiencyTask::theMinChi2NormSegment

private

Definition at line 80 of file DTChamberEfficiencyTask.h.

Referenced by DTChamberEfficiencyTask(), and isGoodSegment().

◆ theMinCloseDist

double DTChamberEfficiencyTask::theMinCloseDist

private

Definition at line 81 of file DTChamberEfficiencyTask.h.

Referenced by analyze(), and DTChamberEfficiencyTask().

◆ theMinHitsSegment

unsigned int DTChamberEfficiencyTask::theMinHitsSegment

private

Definition at line 79 of file DTChamberEfficiencyTask.h.

Referenced by DTChamberEfficiencyTask(), and isGoodSegment().

Public Member Functions

Protected Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ DTChamberEfficiencyTask()

◆ ~DTChamberEfficiencyTask()

Member Function Documentation

◆ analyze()

◆ beginLuminosityBlock()

◆ bookHistograms()

◆ bookHistos()

◆ dqmBeginRun()

◆ endLuminosityBlock()

◆ getBestSegment() [1/2]

◆ getBestSegment() [2/2]

◆ interpolate()

◆ isGoodSegment()

Member Data Documentation

◆ debug

◆ detailedAnalysis

◆ dtGeom

◆ histosPerCh

◆ muonGeomToken_

◆ onlineMonitor

◆ parameters

◆ recHits4DToken_

◆ segs

◆ theMinChi2NormSegment

◆ theMinCloseDist

◆ theMinHitsSegment