#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

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

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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 &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)=0

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

Protected Member Functions inherited from edm::ProducerBase
template<class ProductType >
BranchAliasSetterT< ProductType >	produces ()
	declare what type of product will make and with which optional label More...

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

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

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)

template<Transition B>
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 (std::string instanceName)

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

ProducesCollector	producesCollector ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (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)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

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

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

const DTRecSegment4D &	getBestSegment (const DTRecSegment4DCollection::range &segs) 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

edm::ESHandle< DTGeometry >	dtGeom

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

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
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 26 of file DTChamberEfficiencyTask.cc.

                                                                          {
   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.getParameter<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");
 }

References debug, dtNoiseAnalysis_cfi::detailedAnalysis, dtChamberEfficiencyTask_cfi::onlineMonitor, and muonDTDigis_cfi::pset.

◆ ~DTChamberEfficiencyTask()

DTChamberEfficiencyTask::~DTChamberEfficiencyTask ( )

override

Destructor.

Definition at line 50 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 147 of file DTChamberEfficiencyTask.cc.

                                                                                        {
   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
 }

References dtNoiseAnalysis_cfi::detailedAnalysis, DTRecSegment4D::hasPhi(), DTRecSegment4D::hasZed(), combine::histos, HiCaloJetParameters_cff::interpolate, DTRecSegment4D::localPosition(), mag(), DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), makeMuonMisalignmentScenario::wheel, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

◆ beginLuminosityBlock()

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

override

To reset the MEs.

Definition at line 54 of file DTChamberEfficiencyTask.cc.

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

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

◆ bookHistograms()

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

overrideprotected

Definition at line 75 of file DTChamberEfficiencyTask.cc.

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

References bookHistos().

◆ bookHistos()

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

private

Definition at line 90 of file DTChamberEfficiencyTask.cc.

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

References dtNoiseAnalysis_cfi::detailedAnalysis, combine::histos, DTChamberId::sector(), DTChamberId::station(), relativeConstraints::station, DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

◆ dqmBeginRun()

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

overridevirtual

BeginRun.

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

Definition at line 70 of file DTChamberEfficiencyTask.cc.

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

References edm::get(), and singleTopDQM_cfi::setup.

◆ 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 DTRecSegment4D *	s1,
		const DTRecSegment4D *	s2
	)		const

private

Definition at line 310 of file DTChamberEfficiencyTask.cc.

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

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

◆ getBestSegment() [2/2]

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

private

Definition at line 289 of file DTChamberEfficiencyTask.cc.

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

References DTRecSegment4D::chi2().

◆ interpolate()

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

private

Definition at line 332 of file DTChamberEfficiencyTask.cc.

                                                                               {
   // 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;
 }

References DTRecSegment4D::chamberId(), DeadROC_duringRun::dir, 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().

◆ isGoodSegment()

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

private

Definition at line 366 of file DTChamberEfficiencyTask.cc.

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

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

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(), pkg.AbstractPkg::generate(), util.rrapi.RRApi::get(), rrapi.RRApi::get(), pkg.AbstractPkg::get_kwds(), runTauIdMVA.TauIDEmbedder::loadMVA_WPs_run2_2017(), runTauIdMVA.TauIDEmbedder::runTauID(), and pkg.AbstractPkg::write().

◆ detailedAnalysis

bool DTChamberEfficiencyTask::detailedAnalysis

private

Definition at line 70 of file DTChamberEfficiencyTask.h.

◆ dtGeom

edm::ESHandle<DTGeometry> DTChamberEfficiencyTask::dtGeom

private

Definition at line 83 of file DTChamberEfficiencyTask.h.

◆ histosPerCh

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

private

Definition at line 77 of file DTChamberEfficiencyTask.h.

◆ onlineMonitor

bool DTChamberEfficiencyTask::onlineMonitor

private

Definition at line 68 of file DTChamberEfficiencyTask.h.

◆ parameters

edm::ParameterSet DTChamberEfficiencyTask::parameters

private

Definition at line 75 of file DTChamberEfficiencyTask.h.

◆ recHits4DToken_

edm::EDGetTokenT<DTRecSegment4DCollection> DTChamberEfficiencyTask::recHits4DToken_

private

Definition at line 73 of file DTChamberEfficiencyTask.h.

◆ segs

edm::Handle<DTRecSegment4DCollection> DTChamberEfficiencyTask::segs

private

Definition at line 84 of file DTChamberEfficiencyTask.h.

◆ theMinChi2NormSegment

double DTChamberEfficiencyTask::theMinChi2NormSegment

private

Definition at line 80 of file DTChamberEfficiencyTask.h.

◆ theMinCloseDist

double DTChamberEfficiencyTask::theMinCloseDist

private

Definition at line 81 of file DTChamberEfficiencyTask.h.

◆ theMinHitsSegment

unsigned int DTChamberEfficiencyTask::theMinHitsSegment

private

Definition at line 79 of file DTChamberEfficiencyTask.h.

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

◆ onlineMonitor

◆ parameters

◆ recHits4DToken_

◆ segs

◆ theMinChi2NormSegment

◆ theMinCloseDist

◆ theMinHitsSegment