#include <DTRecHitQuality.h>

Inheritance diagram for DTRecHitQuality:

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

Public Member Functions inherited from edm::global::EDAnalyzer< edm::RunCache< dtrechit::Histograms >, Args... >
	EDAnalyzer ()=default

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

	~EDAnalyzerBase () override

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

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

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

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

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHistograms (DQMStore::ConcurrentBooker &, 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::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 edm::global::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

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

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

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

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

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

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

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

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

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

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 ( const edm::ParameterSet & pset )

Constructor.

Definition at line 88 of file DTRecHitQuality.cc.

References edm::ParameterSet::getUntrackedParameter().

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

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

overrideprivatevirtual

Book the DQM plots.

Implements DQMGlobalEDAnalyzer< dtrechit::Histograms >.

Definition at line 111 of file DTRecHitQuality.cc.

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

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 416 of file DTRecHitQuality.cc.

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

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

overrideprivate

Perform the real analysis.

Definition at line 212 of file DTRecHitQuality.cc.

References bookConverter::compute(), gather_cfg::cout, edm::EventSetup::get(), ecalSelectiveReadoutValidation_cfi::histograms, edm::HandleBase::isValid(), DTHitQualityUtils::mapSimHitsPerWire(), edm::Handle< T >::product(), edm::ESHandle< T >::product(), and trackerHits::simHits.

                                                                      {
   if (debug_) {
     cout << "--- [DTRecHitQuality] Analysing Event: #Run: " << event.id().run() << " #Event: " << event.id().event()
          << endl;
   }
 
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
   // 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.product(), 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.product(), 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.product(), simHitsPerWire, recHitsPerWire, histograms, 3);
   }
 }

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 386 of file DTRecHitQuality.cc.

References rpcPointValidation_cfi::recHit.

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

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

private

Definition at line 304 of file DTRecHitQuality.cc.

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

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

private

Definition at line 316 of file DTRecHitQuality.cc.

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

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

private

Definition at line 332 of file DTRecHitQuality.cc.

References hfClusterShapes_cfi::hits.

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

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

private

Definition at line 405 of file DTRecHitQuality.cc.

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

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

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

private

Definition at line 411 of file DTRecHitQuality.cc.

References DTRecHit1D::localPosition(), DTLayer::specificTopology(), DTWireId::wire(), DTRecHit1D::wireId(), DTTopology::wirePosition(), and PV3DBase< T, PVType, FrameType >::x().

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

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

private

Definition at line 699 of file DTRecHitQuality.cc.

References DTEnums::Left, DTRecHit1DPair::localPositionError(), mathSSE::sqrt(), and LocalError::xx().

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

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

private

Definition at line 704 of file DTRecHitQuality.cc.

References DEFINE_FWK_MODULE, DTRecHit1D::localPositionError(), mathSSE::sqrt(), and LocalError::xx().

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

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

private

Definition at line 374 of file DTRecHitQuality.cc.

References DTTopology::cellLenght(), PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), and PV3DBase< T, PVType, FrameType >::y().

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

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

private

Definition at line 355 of file DTRecHitQuality.cc.

References PSimHit::entryPoint(), PSimHit::exitPoint(), DTLayer::specificTopology(), DTWireId::wire(), DTTopology::wirePosition(), PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::z().

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

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

private

Definition at line 366 of file DTRecHitQuality.cc.

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

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

Member Data Documentation

bool DTRecHitQuality::debug_

private

Definition at line 64 of file DTRecHitQuality.h.

bool DTRecHitQuality::doall_

private

Definition at line 83 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep1_

private

Definition at line 79 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep2_

private

Definition at line 80 of file DTRecHitQuality.h.

bool DTRecHitQuality::doStep3_

private

Definition at line 81 of file DTRecHitQuality.h.

bool DTRecHitQuality::local_

private

Definition at line 82 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::recHitLabel_

private

Definition at line 74 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecHitCollection> DTRecHitQuality::recHitToken_

private

Definition at line 68 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::segment2DLabel_

private

Definition at line 75 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTRecHitQuality::segment2DToken_

private

Definition at line 69 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::segment4DLabel_

private

Definition at line 76 of file DTRecHitQuality.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTRecHitQuality::segment4DToken_

private

Definition at line 70 of file DTRecHitQuality.h.

edm::InputTag DTRecHitQuality::simHitLabel_

private

Definition at line 71 of file DTRecHitQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTRecHitQuality::simHitToken_

private

Definition at line 67 of file DTRecHitQuality.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation