#include <DTRecHitQuality.h>

Inheritance diagram for DTRecHitQuality:

Public Member Functions
	DTRecHitQuality (const edm::ParameterSet &pset)
	Constructor. More...

Public Member Functions inherited from DQMGlobalEDAnalyzer< dtrechit::Histograms >
virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &, dtrechit::Histograms const &) const

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

Public Member Functions inherited from DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >
void	accumulate (edm::StreamID id, edm::Event const &event, edm::EventSetup const &setup) const final

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

	DQMGlobalEDAnalyzerBase ()

std::shared_ptr< dtrechit::Histograms >	globalBeginRun (edm::Run const &run, edm::EventSetup const &setup) const final

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

Public Member Functions inherited from edm::global::EDProducer< edm::RunCache< dtrechit::Histograms >, edm::EndRunProducer, edm::Accumulator, Args... >
	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

EDProducer &	operator= (const EDProducer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ModuleDescription const &	moduleDescription () 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 const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &, dtrechit::Histograms &) const override
	Book the DQM plots. More...

template<typename type >
void	compute (const DTGeometry &dtGeom, const std::map< DTWireId, std::vector< PSimHit >> &simHitsPerWire, const std::map< DTWireId, std::vector< type >> &recHitsPerWire, dtrechit::Histograms const &histograms, int step) const

void	dqmAnalyze (edm::Event const &, edm::EventSetup const &, dtrechit::Histograms const &) const override
	Perform the real analysis. More...

template<typename type >
const type *	findBestRecHit (const DTLayer *layer, const DTWireId &wireId, const std::vector< type > &recHits, float simHitDist) const

std::map< DTWireId, std::vector< DTRecHit1DPair > >	map1DRecHitsPerWire (const DTRecHitCollection *dt1DRecHitPairs) const

std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment2DCollection *segment2Ds) const

std::map< DTWireId, std::vector< DTRecHit1D > >	map1DRecHitsPerWire (const DTRecSegment4DCollection *segment4Ds) const

float	recHitDistFromWire (const DTRecHit1DPair &hitPair, const DTLayer *layer) const

float	recHitDistFromWire (const DTRecHit1D &recHit, const DTLayer *layer) const

float	recHitPositionError (const DTRecHit1DPair &recHit) const

float	recHitPositionError (const DTRecHit1D &recHit) const

float	simHitDistFromFE (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const

float	simHitDistFromWire (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const

float	simHitImpactAngle (const DTLayer *layer, const DTWireId &wireId, const PSimHit &hit) const

Private Attributes
bool	debug_

bool	doall_

bool	doStep1_

bool	doStep2_

bool	doStep3_

bool	local_

edm::ESGetToken< DTGeometry, MuonGeometryRecord >	muonGeomToken_

edm::InputTag	recHitLabel_

edm::EDGetTokenT< DTRecHitCollection >	recHitToken_

edm::InputTag	segment2DLabel_

edm::EDGetTokenT< DTRecSegment2DCollection >	segment2DToken_

edm::InputTag	segment4DLabel_

edm::EDGetTokenT< DTRecSegment4DCollection >	segment4DToken_

edm::InputTag	simHitLabel_

edm::EDGetTokenT< edm::PSimHitContainer >	simHitToken_

Additional Inherited Members
Public Types inherited from DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::global::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::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >
uint64_t	meId (edm::Run const &run) 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 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<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)

Protected Attributes inherited from DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >
DQMStore *	dqmstore_

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Basic analyzer class which accesses 1D DTRecHits and plot resolution comparing reconstructed and simulated quantities

Residual/pull plots are filled for the rechit with distance from wire closer to that of the muon simhit.

Efficiencies are defined as the fraction of muon simhits with a rechit in the same cell, for the given reconstruction step. Hence, for S2 and S3 the definition incorporate the segment reconstruction efficiency.

Author: G. Cerminara - INFN Torino

Definition at line 47 of file DTRecHitQuality.h.

Constructor & Destructor Documentation

◆ DTRecHitQuality()

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

Constructor.

Definition at line 87 of file DTRecHitQuality.cc.

References debug_, doall_, doStep1_, doStep2_, doStep3_, local_, muonDTDigis_cfi::pset, recHitLabel_, recHitToken_, segment2DLabel_, segment2DToken_, segment4DLabel_, segment4DToken_, simHitLabel_, and simHitToken_.

                                                          : muonGeomToken_(esConsumes()) {
   // Get the debug parameter for verbose output
   debug_ = pset.getUntrackedParameter<bool>("debug");
   // the name of the simhit collection
   simHitLabel_ = pset.getUntrackedParameter<InputTag>("simHitLabel");
   simHitToken_ = consumes<PSimHitContainer>(pset.getUntrackedParameter<InputTag>("simHitLabel"));
   // the name of the 1D rec hit collection
   recHitLabel_ = pset.getUntrackedParameter<InputTag>("recHitLabel");
   recHitToken_ = consumes<DTRecHitCollection>(pset.getUntrackedParameter<InputTag>("recHitLabel"));
   // the name of the 2D rec hit collection
   segment2DLabel_ = pset.getUntrackedParameter<InputTag>("segment2DLabel");
   segment2DToken_ = consumes<DTRecSegment2DCollection>(pset.getUntrackedParameter<InputTag>("segment2DLabel"));
   // the name of the 4D rec hit collection
   segment4DLabel_ = pset.getUntrackedParameter<InputTag>("segment4DLabel");
   segment4DToken_ = consumes<DTRecSegment4DCollection>(pset.getUntrackedParameter<InputTag>("segment4DLabel"));
   // Switches for analysis at various steps
   doStep1_ = pset.getUntrackedParameter<bool>("doStep1", false);
   doStep2_ = pset.getUntrackedParameter<bool>("doStep2", false);
   doStep3_ = pset.getUntrackedParameter<bool>("doStep3", false);
   doall_ = pset.getUntrackedParameter<bool>("doall", false);
   local_ = pset.getUntrackedParameter<bool>("local", true);
 }

Member Function Documentation

◆ bookHistograms()

void DTRecHitQuality::bookHistograms	(	DQMStore::IBooker &	booker,
		edm::Run const &	run,
		edm::EventSetup const &	setup,
		dtrechit::Histograms &	histograms
	)		const

overrideprivatevirtual

Book the DQM plots.

Implements DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >.

Definition at line 110 of file DTRecHitQuality.cc.

References doall_, doStep1_, doStep2_, doStep3_, local_, hfnoseParametersInitialization_cfi::name2, alignCSCRings::s, and w().

                                                                    {
   if (doall_ && doStep1_) {
     histograms.hRes_S1RPhi = std::make_unique<HRes1DHit>("S1RPhi", booker, true, local_);  // RecHits, 1. step, RPhi
     histograms.hRes_S1RPhi_W0 =
         std::make_unique<HRes1DHit>("S1RPhi_W0", booker, true, local_);  // RecHits, 1. step, RZ, wheel 0
     histograms.hRes_S1RPhi_W1 =
         std::make_unique<HRes1DHit>("S1RPhi_W1", booker, true, local_);  // RecHits, 1. step, RZ, wheel +-1
     histograms.hRes_S1RPhi_W2 =
         std::make_unique<HRes1DHit>("S1RPhi_W2", booker, true, local_);  // RecHits, 1. step, RZ, wheel +-2
     histograms.hRes_S1RZ = std::make_unique<HRes1DHit>("S1RZ", booker, true, local_);  // RecHits, 1. step, RZ
     histograms.hRes_S1RZ_W0 =
         std::make_unique<HRes1DHit>("S1RZ_W0", booker, true, local_);  // RecHits, 1. step, RZ, wheel 0
     histograms.hRes_S1RZ_W1 =
         std::make_unique<HRes1DHit>("S1RZ_W1", booker, true, local_);  // RecHits, 1. step, RZ, wheel +-1
     histograms.hRes_S1RZ_W2 =
         std::make_unique<HRes1DHit>("S1RZ_W2", booker, true, local_);          // RecHits, 1. step, RZ, wheel +-2
     histograms.hEff_S1RPhi = std::make_unique<HEff1DHit>("S1RPhi", booker);    // RecHits, 1. step, RPhi
     histograms.hEff_S1RZ = std::make_unique<HEff1DHit>("S1RZ", booker);        // RecHits, 1. step, RZ
     histograms.hEff_S1RZ_W0 = std::make_unique<HEff1DHit>("S1RZ_W0", booker);  // RecHits, 1. step, RZ, wheel 0
     histograms.hEff_S1RZ_W1 = std::make_unique<HEff1DHit>("S1RZ_W1", booker);  // RecHits, 1. step, RZ, wheel +-1
     histograms.hEff_S1RZ_W2 = std::make_unique<HEff1DHit>("S1RZ_W2", booker);  // RecHits, 1. step, RZ, wheel +-2
   }
   if (doall_ && doStep2_) {
     histograms.hRes_S2RPhi = std::make_unique<HRes1DHit>("S2RPhi", booker, true, local_);  // RecHits, 2. step, RPhi
     histograms.hRes_S2RPhi_W0 =
         std::make_unique<HRes1DHit>("S2RPhi_W0", booker, true, local_);  // RecHits, 2. step, RPhi, wheel 0
     histograms.hRes_S2RPhi_W1 =
         std::make_unique<HRes1DHit>("S2RPhi_W1", booker, true, local_);  // RecHits, 2. step, RPhi, wheel +-1
     histograms.hRes_S2RPhi_W2 =
         std::make_unique<HRes1DHit>("S2RPhi_W2", booker, true, local_);  // RecHits, 2. step, RPhi, wheel +-2
     histograms.hRes_S2RZ = std::make_unique<HRes1DHit>("S2RZ", booker, true, local_);  // RecHits, 2. step, RZ
     histograms.hRes_S2RZ_W0 =
         std::make_unique<HRes1DHit>("S2RZ_W0", booker, true, local_);  // RecHits, 2. step, RZ, wheel 0
     histograms.hRes_S2RZ_W1 =
         std::make_unique<HRes1DHit>("S2RZ_W1", booker, true, local_);  // RecHits, 2. step, RZ, wheel +-1
     histograms.hRes_S2RZ_W2 =
         std::make_unique<HRes1DHit>("S2RZ_W2", booker, true, local_);          // RecHits, 2. step, RZ, wheel +-2
     histograms.hEff_S2RPhi = std::make_unique<HEff1DHit>("S2RPhi", booker);    // RecHits, 2. step, RPhi
     histograms.hEff_S2RZ_W0 = std::make_unique<HEff1DHit>("S2RZ_W0", booker);  // RecHits, 2. step, RZ, wheel 0
     histograms.hEff_S2RZ_W1 = std::make_unique<HEff1DHit>("S2RZ_W1", booker);  // RecHits, 2. step, RZ, wheel +-1
     histograms.hEff_S2RZ_W2 = std::make_unique<HEff1DHit>("S2RZ_W2", booker);  // RecHits, 2. step, RZ, wheel +-2
     histograms.hEff_S2RZ = std::make_unique<HEff1DHit>("S2RZ", booker);        // RecHits, 2. step, RZ
   }
   if (doStep3_) {
     histograms.hRes_S3RPhi = std::make_unique<HRes1DHit>("S3RPhi", booker, doall_, local_);  // RecHits, 3. step, RPhi
     histograms.hRes_S3RPhi_W0 =
         std::make_unique<HRes1DHit>("S3RPhi_W0", booker, doall_, local_);  // RecHits, 3. step, RPhi, wheel 0
     histograms.hRes_S3RPhi_W1 = std::make_unique<HRes1DHit>("S3RPhi_W1",
                                                             booker,
                                                             doall_,
                                                             local_);  // RecHits, 3. step, RPhi, wheel +-1
     histograms.hRes_S3RPhi_W2 = std::make_unique<HRes1DHit>("S3RPhi_W2",
                                                             booker,
                                                             doall_,
                                                             local_);  // RecHits, 3. step, RPhi, wheel +-2
     histograms.hRes_S3RZ = std::make_unique<HRes1DHit>("S3RZ", booker, doall_, local_);  // RecHits, 3. step, RZ
     histograms.hRes_S3RZ_W0 =
         std::make_unique<HRes1DHit>("S3RZ_W0", booker, doall_, local_);  // RecHits, 3. step, RZ, wheel 0
     histograms.hRes_S3RZ_W1 =
         std::make_unique<HRes1DHit>("S3RZ_W1", booker, doall_, local_);  // RecHits, 3. step, RZ, wheel +-1
     histograms.hRes_S3RZ_W2 =
         std::make_unique<HRes1DHit>("S3RZ_W2", booker, doall_, local_);  // RecHits, 3. step, RZ, wheel +-2
 
     if (local_) {
       // Plots with finer granularity, not to be included in DQM
       TString name1 = "RPhi_W";
       TString name2 = "RZ_W";
       for (long w = 0; w <= 2; ++w) {
         for (long s = 1; s <= 4; ++s) {
           histograms.hRes_S3RPhiWS[w][s - 1] =
               std::make_unique<HRes1DHit>(("S3" + name1 + w + "_St" + s).Data(), booker, doall_, local_);
           histograms.hEff_S1RPhiWS[w][s - 1] =
               std::make_unique<HEff1DHit>(("S1" + name1 + w + "_St" + s).Data(), booker);
           histograms.hEff_S3RPhiWS[w][s - 1] =
               std::make_unique<HEff1DHit>(("S3" + name1 + w + "_St" + s).Data(), booker);
           if (s != 4) {
             histograms.hRes_S3RZWS[w][s - 1] =
                 std::make_unique<HRes1DHit>(("S3" + name2 + w + "_St" + s).Data(), booker, doall_, local_);
             histograms.hEff_S1RZWS[w][s - 1] =
                 std::make_unique<HEff1DHit>(("S1" + name2 + w + "_St" + s).Data(), booker);
             histograms.hEff_S3RZWS[w][s - 1] =
                 std::make_unique<HEff1DHit>(("S3" + name2 + w + "_St" + s).Data(), booker);
           }
         }
       }
     }
 
     if (doall_) {
       histograms.hEff_S3RPhi = std::make_unique<HEff1DHit>("S3RPhi", booker);    // RecHits, 3. step, RPhi
       histograms.hEff_S3RZ = std::make_unique<HEff1DHit>("S3RZ", booker);        // RecHits, 3. step, RZ
       histograms.hEff_S3RZ_W0 = std::make_unique<HEff1DHit>("S3RZ_W0", booker);  // RecHits, 3. step, RZ, wheel 0
       histograms.hEff_S3RZ_W1 = std::make_unique<HEff1DHit>("S3RZ_W1", booker);  // RecHits, 3. step, RZ, wheel +-1
       histograms.hEff_S3RZ_W2 = std::make_unique<HEff1DHit>("S3RZ_W2", booker);  // RecHits, 3. step, RZ, wheel +-2
     }
   }
 }

◆ compute()

template<typename type >

void DTRecHitQuality::compute	(	const DTGeometry &	dtGeom,
		const std::map< DTWireId, std::vector< PSimHit >> &	simHitsPerWire,
		const std::map< DTWireId, std::vector< type >> &	recHitsPerWire,
		dtrechit::Histograms const &	histograms,
		int	step
	)		const

private

Definition at line 414 of file DTRecHitQuality.cc.

References funct::abs(), gather_cfg::cout, debug_, doall_, PV3DBase< T, PVType, FrameType >::eta(), HRes1DHit::fill(), HEff1DHit::fill(), findBestRecHit(), DTHitQualityUtils::findMuSimHit(), PixelMapPlotter::hRes, DTGeometry::layer(), pixelTopology::layer, local_, PSimHit::localPosition(), PV3DBase< T, PVType, FrameType >::phi(), recHitDistFromWire(), recHitPositionError(), FastTrackerRecHitMaskProducer_cfi::recHits, simHitDistFromFE(), simHitDistFromWire(), simHitImpactAngle(), DTChamberId::station(), DTSuperLayerId::superLayer(), DTChamberId::wheel(), and makeMuonMisalignmentScenario::wheel.

Referenced by dqmAnalyze().

                                               {
   // Loop over cells with a muon SimHit
   for (const auto &wireAndSHits : simHitsPerWire) {
     DTWireId wireId = wireAndSHits.first;
     int wheel = wireId.wheel();
     int sl = wireId.superLayer();
 
     vector<PSimHit> simHitsInCell = wireAndSHits.second;
 
     // Get the layer
     const DTLayer *layer = dtGeom.layer(wireId);
 
     // Look for a mu hit in the cell
     const PSimHit *muSimHit = DTHitQualityUtils::findMuSimHit(simHitsInCell);
     if (muSimHit == nullptr) {
       if (debug_) {
         cout << "   No mu SimHit in channel: " << wireId << ", skipping! " << endl;
       }
       continue;  // Skip this cell
     }
 
     // Find the distance of the simhit from the wire
     float simHitWireDist = simHitDistFromWire(layer, wireId, *muSimHit);
     // Skip simhits out of the cell
     if (simHitWireDist > 2.1) {
       if (debug_) {
         cout << "  [DTRecHitQuality]###Warning: The mu SimHit in out of the "
                 "cell, skipping!"
              << endl;
       }
       continue;  // Skip this cell
     }
     GlobalPoint simHitGlobalPos = layer->toGlobal(muSimHit->localPosition());
 
     // find SH impact angle
     float simHitTheta = simHitImpactAngle(layer, wireId, *muSimHit);
 
     // find SH distance from FE
     float simHitFEDist = simHitDistFromFE(layer, wireId, *muSimHit);
 
     bool recHitReconstructed = false;
 
     // Look for RecHits in the same cell
     if (recHitsPerWire.find(wireId) == recHitsPerWire.end()) {
       // No RecHit found in this cell
       if (debug_) {
         cout << "   No RecHit found at Step: " << step << " in cell: " << wireId << endl;
       }
     } else {
       recHitReconstructed = true;
       // vector<type> recHits = (*wireAndRecHits).second;
       const vector<type> &recHits = recHitsPerWire.at(wireId);
       if (debug_) {
         cout << "   " << recHits.size() << " RecHits, Step " << step << " in channel: " << wireId << endl;
       }
 
       // Find the best RecHit
       const type *theBestRecHit = findBestRecHit(layer, wireId, recHits, simHitWireDist);
 
       float recHitWireDist = recHitDistFromWire(*theBestRecHit, layer);
       if (debug_) {
         cout << "    SimHit distance from wire: " << simHitWireDist << endl
              << "    SimHit distance from FE:   " << simHitFEDist << endl
              << "    SimHit angle in layer RF:  " << simHitTheta << endl
              << "    RecHit distance from wire: " << recHitWireDist << endl;
       }
       float recHitErr = recHitPositionError(*theBestRecHit);
       HRes1DHit *hRes = nullptr;
       HRes1DHit *hResTot = nullptr;
 
       // Mirror angle in phi so that + and - wheels can be plotted together
       if (mirrorMinusWheels && wheel < 0 && sl != 2) {
         simHitTheta *= -1.;
         // Note: local X, if used, would have to be mirrored as well
       }
 
       // Fill residuals and pulls
       // Select the histo to be filled
       if (step == 1) {
         // Step 1
         if (sl != 2) {
           hResTot = histograms.hRes_S1RPhi.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S1RPhi_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S1RPhi_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S1RPhi_W2.get();
           }
         } else {
           hResTot = histograms.hRes_S1RZ.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S1RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S1RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S1RZ_W2.get();
           }
         }
 
       } else if (step == 2) {
         // Step 2
         if (sl != 2) {
           hRes = histograms.hRes_S2RPhi.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S2RPhi_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S2RPhi_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S2RPhi_W2.get();
           }
         } else {
           hResTot = histograms.hRes_S2RZ.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S2RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S2RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S2RZ_W2.get();
           }
         }
 
       } else if (step == 3) {
         // Step 3
         if (sl != 2) {
           hResTot = histograms.hRes_S3RPhi.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S3RPhi_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S3RPhi_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S3RPhi_W2.get();
           }
           if (local_) {
             histograms.hRes_S3RPhiWS[abs(wheel)][wireId.station() - 1]->fill(simHitWireDist,
                                                                              simHitTheta,
                                                                              simHitFEDist,
                                                                              recHitWireDist,
                                                                              simHitGlobalPos.eta(),
                                                                              simHitGlobalPos.phi(),
                                                                              recHitErr,
                                                                              wireId.station());
           }
         } else {
           hResTot = histograms.hRes_S3RZ.get();
           if (wheel == 0) {
             hRes = histograms.hRes_S3RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hRes = histograms.hRes_S3RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hRes = histograms.hRes_S3RZ_W2.get();
           }
           if (local_) {
             histograms.hRes_S3RZWS[abs(wheel)][wireId.station() - 1]->fill(simHitWireDist,
                                                                            simHitTheta,
                                                                            simHitFEDist,
                                                                            recHitWireDist,
                                                                            simHitGlobalPos.eta(),
                                                                            simHitGlobalPos.phi(),
                                                                            recHitErr,
                                                                            wireId.station());
           }
         }
       }
 
       // Fill
       hRes->fill(simHitWireDist,
                  simHitTheta,
                  simHitFEDist,
                  recHitWireDist,
                  simHitGlobalPos.eta(),
                  simHitGlobalPos.phi(),
                  recHitErr,
                  wireId.station());
       if (hResTot != nullptr) {
         hResTot->fill(simHitWireDist,
                       simHitTheta,
                       simHitFEDist,
                       recHitWireDist,
                       simHitGlobalPos.eta(),
                       simHitGlobalPos.phi(),
                       recHitErr,
                       wireId.station());
       }
     }
 
     // Fill Efficiencies
     if (doall_) {
       HEff1DHit *hEff = nullptr;
       HEff1DHit *hEffTot = nullptr;
       if (step == 1) {
         // Step 1
         if (sl != 2) {
           hEff = histograms.hEff_S1RPhi.get();
           if (local_) {
             histograms.hEff_S1RPhiWS[abs(wheel)][wireId.station() - 1]->fill(
                 simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
           }
         } else {
           hEffTot = histograms.hEff_S1RZ.get();
           if (wheel == 0) {
             hEff = histograms.hEff_S1RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hEff = histograms.hEff_S1RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hEff = histograms.hEff_S1RZ_W2.get();
           }
           if (local_) {
             histograms.hEff_S1RZWS[abs(wheel)][wireId.station() - 1]->fill(
                 simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
           }
         }
 
       } else if (step == 2) {
         // Step 2
         if (sl != 2) {
           hEff = histograms.hEff_S2RPhi.get();
         } else {
           hEffTot = histograms.hEff_S2RZ.get();
           if (wheel == 0) {
             hEff = histograms.hEff_S2RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hEff = histograms.hEff_S2RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hEff = histograms.hEff_S2RZ_W2.get();
           }
         }
 
       } else if (step == 3) {
         // Step 3
         if (sl != 2) {
           hEff = histograms.hEff_S3RPhi.get();
           if (local_) {
             histograms.hEff_S3RPhiWS[abs(wheel)][wireId.station() - 1]->fill(
                 simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
           }
         } else {
           hEffTot = histograms.hEff_S3RZ.get();
           if (wheel == 0) {
             hEff = histograms.hEff_S3RZ_W0.get();
           }
           if (abs(wheel) == 1) {
             hEff = histograms.hEff_S3RZ_W1.get();
           }
           if (abs(wheel) == 2) {
             hEff = histograms.hEff_S3RZ_W2.get();
           }
           if (local_) {
             histograms.hEff_S3RZWS[abs(wheel)][wireId.station() - 1]->fill(
                 simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
           }
         }
       }
       // Fill
       hEff->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
       if (hEffTot != nullptr) {
         hEffTot->fill(simHitWireDist, simHitGlobalPos.eta(), simHitGlobalPos.phi(), recHitReconstructed);
       }
     }
   }
 }

◆ dqmAnalyze()

void DTRecHitQuality::dqmAnalyze	(	edm::Event const &	event,
		edm::EventSetup const &	setup,
		dtrechit::Histograms const &	histograms
	)		const

overrideprivatevirtual

Perform the real analysis.

Implements DQMGlobalEDAnalyzerBase< dtrechit::Histograms, Args... >.

Definition at line 211 of file DTRecHitQuality.cc.

References compute(), gather_cfg::cout, debug_, doall_, doStep1_, doStep2_, doStep3_, edm::HandleBase::isValid(), map1DRecHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), muonGeomToken_, edm::Handle< T >::product(), recHitLabel_, recHitToken_, segment2DLabel_, segment2DToken_, segment4DLabel_, segment4DToken_, singleTopDQM_cfi::setup, FastTrackerRecHitCombiner_cfi::simHits, and simHitToken_.

                                                                      {
   if (debug_) {
     cout << "--- [DTRecHitQuality] Analysing Event: #Run: " << event.id().run() << " #Event: " << event.id().event()
          << endl;
   }
 
   // Get the DT Geometry
   const DTGeometry &dtGeom = setup.getData(muonGeomToken_);
 
   // Get the SimHit collection from the event
   Handle<PSimHitContainer> simHits;
   event.getByToken(simHitToken_, simHits);
 
   // Map simhits per wire
   map<DTWireId, PSimHitContainer> simHitsPerWire = DTHitQualityUtils::mapSimHitsPerWire(*(simHits.product()));
 
   //=======================================================================================
   // RecHit analysis at Step 1
   if (doStep1_ && doall_) {
     if (debug_) {
       cout << "  -- DTRecHit S1: begin analysis:" << endl;
     }
     // Get the rechit collection from the event
     Handle<DTRecHitCollection> dtRecHits;
     event.getByToken(recHitToken_, dtRecHits);
 
     if (!dtRecHits.isValid()) {
       if (debug_) {
         cout << "[DTRecHitQuality]**Warning: no 1DRechits with label: " << recHitLabel_ << " in this event, skipping!"
              << endl;
       }
       return;
     }
 
     // Map rechits per wire
     auto const &recHitsPerWire = map1DRecHitsPerWire(dtRecHits.product());
     compute(dtGeom, simHitsPerWire, recHitsPerWire, histograms, 1);
   }
 
   //=======================================================================================
   // RecHit analysis at Step 2
   if (doStep2_ && doall_) {
     if (debug_) {
       cout << "  -- DTRecHit S2: begin analysis:" << endl;
     }
 
     // Get the 2D rechits from the event
     Handle<DTRecSegment2DCollection> segment2Ds;
     event.getByToken(segment2DToken_, segment2Ds);
 
     if (!segment2Ds.isValid()) {
       if (debug_) {
         cout << "[DTRecHitQuality]**Warning: no 2DSegments with label: " << segment2DLabel_
              << " in this event, skipping!" << endl;
       }
 
     } else {
       // Map rechits per wire
       auto const &recHitsPerWire = map1DRecHitsPerWire(segment2Ds.product());
       compute(dtGeom, simHitsPerWire, recHitsPerWire, histograms, 2);
     }
   }
 
   //=======================================================================================
   // RecHit analysis at Step 3
   if (doStep3_) {
     if (debug_) {
       cout << "  -- DTRecHit S3: begin analysis:" << endl;
     }
 
     // Get the 4D rechits from the event
     Handle<DTRecSegment4DCollection> segment4Ds;
     event.getByToken(segment4DToken_, segment4Ds);
 
     if (!segment4Ds.isValid()) {
       if (debug_) {
         cout << "[DTRecHitQuality]**Warning: no 4D Segments with label: " << segment4DLabel_
              << " in this event, skipping!" << endl;
       }
       return;
     }
 
     // Map rechits per wire
     auto const &recHitsPerWire = map1DRecHitsPerWire(segment4Ds.product());
     compute(dtGeom, simHitsPerWire, recHitsPerWire, histograms, 3);
   }
 }

◆ findBestRecHit()

template<typename type >

const type * DTRecHitQuality::findBestRecHit	(	const DTLayer *	layer,
		const DTWireId &	wireId,
		const std::vector< type > &	recHits,
		float	simHitDist
	)		const

private

Definition at line 384 of file DTRecHitQuality.cc.

References pixelTopology::layer, rpcPointValidation_cfi::recHit, recHitDistFromWire(), and FastTrackerRecHitMaskProducer_cfi::recHits.

Referenced by compute().

                                                                           {
   float res = 99999;
   const type *theBestRecHit = nullptr;
   // Loop over RecHits within the cell
   for (auto recHit = recHits.begin(); recHit != recHits.end(); recHit++) {
     float distTmp = recHitDistFromWire(*recHit, layer);
     if (fabs(distTmp - simHitDist) < res) {
       res = fabs(distTmp - simHitDist);
       theBestRecHit = &(*recHit);
     }
   }  // End of loop over RecHits within the cell
 
   return theBestRecHit;
 }

◆ map1DRecHitsPerWire() [1/3]

map< DTWireId, vector< DTRecHit1DPair > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecHitCollection * dt1DRecHitPairs ) const

private

Definition at line 302 of file DTRecHitQuality.cc.

References runTheMatrix::ret.

Referenced by dqmAnalyze().

                                                      {
   map<DTWireId, vector<DTRecHit1DPair>> ret;
 
   for (const auto &dt1DRecHitPair : *dt1DRecHitPairs) {
     ret[dt1DRecHitPair.wireId()].push_back(dt1DRecHitPair);
   }
 
   return ret;
 }

◆ map1DRecHitsPerWire() [2/3]

map< DTWireId, vector< DTRecHit1D > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecSegment2DCollection * segment2Ds ) const

private

Definition at line 314 of file DTRecHitQuality.cc.

References runTheMatrix::ret.

                                                       {
   map<DTWireId, vector<DTRecHit1D>> ret;
 
   // Loop over all 2D segments
   for (const auto &segment2D : *segment2Ds) {
     vector<DTRecHit1D> component1DHits = segment2D.specificRecHits();
     // Loop over all component 1D hits
     for (auto &component1DHit : component1DHits) {
       ret[component1DHit.wireId()].push_back(component1DHit);
     }
   }
   return ret;
 }

◆ map1DRecHitsPerWire() [3/3]

map< DTWireId, std::vector< DTRecHit1D > > DTRecHitQuality::map1DRecHitsPerWire ( const DTRecSegment4DCollection * segment4Ds ) const

private

Definition at line 330 of file DTRecHitQuality.cc.

References hfClusterShapes_cfi::hits, and runTheMatrix::ret.

                                                       {
   map<DTWireId, vector<DTRecHit1D>> ret;
   // Loop over all 4D segments
   for (const auto &segment4D : *segment4Ds) {
     // Get component 2D segments
     vector<const TrackingRecHit *> segment2Ds = segment4D.recHits();
     // Loop over 2D segments:
     for (auto &segment2D : segment2Ds) {
       // Get 1D component rechits
       vector<const TrackingRecHit *> hits = segment2D->recHits();
       // Loop over them
       for (auto &hit : hits) {
         const auto *hit1D = dynamic_cast<const DTRecHit1D *>(hit);
         ret[hit1D->wireId()].push_back(*hit1D);
       }
     }
   }
 
   return ret;
 }

◆ recHitDistFromWire() [1/2]

float DTRecHitQuality::recHitDistFromWire	(	const DTRecHit1DPair &	hitPair,
		const DTLayer *	layer
	)		const

private

Definition at line 403 of file DTRecHitQuality.cc.

References DTEnums::Left, DTRecHit1DPair::localPosition(), DTEnums::Right, and PV3DBase< T, PVType, FrameType >::x().

Referenced by compute(), and findBestRecHit().

                                                                                                    {
   // Compute the rechit distance from wire
   return fabs(hitPair.localPosition(DTEnums::Left).x() - hitPair.localPosition(DTEnums::Right).x()) / 2.;
 }

◆ recHitDistFromWire() [2/2]

float DTRecHitQuality::recHitDistFromWire	(	const DTRecHit1D &	recHit,
		const DTLayer *	layer
	)		const

private

Definition at line 409 of file DTRecHitQuality.cc.

References pixelTopology::layer, and rpcPointValidation_cfi::recHit.

                                                                                               {
   return fabs(recHit.localPosition().x() - layer->specificTopology().wirePosition(recHit.wireId().wire()));
 }

◆ recHitPositionError() [1/2]

float DTRecHitQuality::recHitPositionError ( const DTRecHit1DPair & recHit ) const

private

Definition at line 697 of file DTRecHitQuality.cc.

References DTEnums::Left, rpcPointValidation_cfi::recHit, and mathSSE::sqrt().

Referenced by compute().

                                                                              {
   return sqrt(recHit.localPositionError(DTEnums::Left).xx());
 }

◆ recHitPositionError() [2/2]

float DTRecHitQuality::recHitPositionError ( const DTRecHit1D & recHit ) const

private

Definition at line 702 of file DTRecHitQuality.cc.

References rpcPointValidation_cfi::recHit, and mathSSE::sqrt().

                                                                          {
   return sqrt(recHit.localPositionError().xx());
 }

◆ simHitDistFromFE()

float DTRecHitQuality::simHitDistFromFE	(	const DTLayer *	layer,
		const DTWireId &	wireId,
		const PSimHit &	hit
	)		const

private

Definition at line 372 of file DTRecHitQuality.cc.

References pixelTopology::layer, and PV3DBase< T, PVType, FrameType >::y().

Referenced by compute().

                                                                                                               {
   LocalPoint entryP = hit.entryPoint();
   LocalPoint exitP = hit.exitPoint();
   float wireLenght = layer->specificTopology().cellLenght();
   // FIXME: should take only wireLenght/2.;
   // moreover, pos+cellLenght/2. is shorter than the distance from FE.
   // In fact it would make more sense to make plots vs y.
   return (entryP.y() + exitP.y()) / 2. + wireLenght;
 }

◆ simHitDistFromWire()

float DTRecHitQuality::simHitDistFromWire	(	const DTLayer *	layer,
		const DTWireId &	wireId,
		const PSimHit &	hit
	)		const

private

Definition at line 353 of file DTRecHitQuality.cc.

References pixelTopology::layer, DTWireId::wire(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by compute().

                                                                                                                 {
   float xwire = layer->specificTopology().wirePosition(wireId.wire());
   LocalPoint entryP = hit.entryPoint();
   LocalPoint exitP = hit.exitPoint();
   float xEntry = entryP.x() - xwire;
   float xExit = exitP.x() - xwire;
 
   return fabs(xEntry - (entryP.z() * (xExit - xEntry)) / (exitP.z() - entryP.z()));  // FIXME: check...
 }

◆ simHitImpactAngle()

float DTRecHitQuality::simHitImpactAngle	(	const DTLayer *	layer,
		const DTWireId &	wireId,
		const PSimHit &	hit
	)		const

private

Definition at line 364 of file DTRecHitQuality.cc.

References theta(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

Referenced by compute().

                                                                                                                {
   LocalPoint entryP = hit.entryPoint();
   LocalPoint exitP = hit.exitPoint();
   float theta = (exitP.x() - entryP.x()) / (exitP.z() - entryP.z());
   return atan(theta);
 }