#include <DTSegment2DQuality.h>

Inheritance diagram for DTSegment2DQuality:

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

Public Member Functions inherited from edm::global::EDAnalyzer< edm::RunCache< dtsegment2d::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 &, dtsegment2d::Histograms &) const override
	Book the DQM plots. More...

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

Private Attributes
bool	debug_

edm::InputTag	segment2DLabel_

edm::EDGetTokenT< DTRecSegment2DCollection >	segment2DToken_

double	sigmaResAngle_

double	sigmaResPos_

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 2D DTSegments and plot resolution comparing reconstructed and simulated quantities

Author: S. Bolognesi and G. Cerminara - INFN Torino

Definition at line 33 of file DTSegment2DQuality.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 47 of file DTSegment2DQuality.cc.

References gather_cfg::cout, DTHitQualityUtils::debug, edm::ParameterSet::getParameter(), and edm::ParameterSet::getUntrackedParameter().

                                                                {
   // get the debug parameter for verbose output
   debug_ = pset.getUntrackedParameter<bool>("debug");
   DTHitQualityUtils::debug = debug_;
   // the name of the simhit collection
   simHitLabel_ = pset.getUntrackedParameter<InputTag>("simHitLabel");
   simHitToken_ = consumes<PSimHitContainer>(pset.getUntrackedParameter<InputTag>("simHitLabel"));
   // the name of the 2D rec hit collection
   segment2DLabel_ = pset.getUntrackedParameter<InputTag>("segment2DLabel");
   segment2DToken_ = consumes<DTRecSegment2DCollection>(pset.getUntrackedParameter<InputTag>("segment2DLabel"));
 
   // sigma resolution on position
   sigmaResPos_ = pset.getParameter<double>("sigmaResPos");
   // sigma resolution on angle
   sigmaResAngle_ = pset.getParameter<double>("sigmaResAngle");
 
   if (debug_) {
     cout << "[DTSegment2DQuality] Constructor called " << endl;
   }
 }

Member Function Documentation

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

overrideprivatevirtual

Book the DQM plots.

Implements DQMGlobalEDAnalyzer< dtsegment2d::Histograms >.

Definition at line 68 of file DTSegment2DQuality.cc.

References gather_cfg::cout, dtsegment2d::Histograms::h2DHitEff_RPhi, dtsegment2d::Histograms::h2DHitEff_RZ, dtsegment2d::Histograms::h2DHitEff_RZ_W0, dtsegment2d::Histograms::h2DHitEff_RZ_W1, dtsegment2d::Histograms::h2DHitEff_RZ_W2, dtsegment2d::Histograms::h2DHitRPhi, dtsegment2d::Histograms::h2DHitRZ, dtsegment2d::Histograms::h2DHitRZ_W0, dtsegment2d::Histograms::h2DHitRZ_W1, and dtsegment2d::Histograms::h2DHitRZ_W2.

                                                                       {
   histograms.h2DHitRPhi = std::make_unique<HRes2DHit>("RPhi", booker, true, true);
   histograms.h2DHitRZ = std::make_unique<HRes2DHit>("RZ", booker, true, true);
   histograms.h2DHitRZ_W0 = std::make_unique<HRes2DHit>("RZ_W0", booker, true, true);
   histograms.h2DHitRZ_W1 = std::make_unique<HRes2DHit>("RZ_W1", booker, true, true);
   histograms.h2DHitRZ_W2 = std::make_unique<HRes2DHit>("RZ_W2", booker, true, true);
 
   histograms.h2DHitEff_RPhi = std::make_unique<HEff2DHit>("RPhi", booker);
   histograms.h2DHitEff_RZ = std::make_unique<HEff2DHit>("RZ", booker);
   histograms.h2DHitEff_RZ_W0 = std::make_unique<HEff2DHit>("RZ_W0", booker);
   histograms.h2DHitEff_RZ_W1 = std::make_unique<HEff2DHit>("RZ_W1", booker);
   histograms.h2DHitEff_RZ_W2 = std::make_unique<HEff2DHit>("RZ_W2", booker);
   if (debug_) {
     cout << "[DTSegment2DQuality] hitsos created " << endl;
   }
 }

void DTSegment2DQuality::dqmAnalyze	(	edm::Event const &	event,
		edm::EventSetup const &	setup,
		dtsegment2d::Histograms const &	histograms
	)		const

overrideprivate

Perform the real analysis.

Definition at line 89 of file DTSegment2DQuality.cc.

                                                                         {
   // Get the DT Geometry
   ESHandle<DTGeometry> dtGeom;
   setup.get<MuonGeometryRecord>().get(dtGeom);
 
   // Get the SimHit collection from the event
   edm::Handle<PSimHitContainer> simHits;
   event.getByToken(simHitToken_, simHits);  // FIXME: second string to be removed
 
   // Map simHits by sl
   map<DTSuperLayerId, PSimHitContainer> simHitsPerSl;
   for (const auto &simHit : *simHits) {
     // Create the id of the sl (the simHits in the DT known their wireId)
     DTSuperLayerId slId = ((DTWireId(simHit.detUnitId())).layerId()).superlayerId();
     // Fill the map
     simHitsPerSl[slId].push_back(simHit);
   }
 
   // Get the 2D rechits from the event
   Handle<DTRecSegment2DCollection> segment2Ds;
   event.getByToken(segment2DToken_, segment2Ds);
 
   if (not segment2Ds.isValid()) {
     if (debug_) {
       cout << "[DTSegment2DQuality]**Warning: no 2DSegments with label: " << segment2DLabel_
            << " in this event, skipping !" << endl;
     }
     return;
   }
 
   // Loop over all superlayers containing a segment
   DTRecSegment2DCollection::id_iterator slId;
   for (slId = segment2Ds->id_begin(); slId != segment2Ds->id_end(); ++slId) {
     //------------------------- simHits ---------------------------//
     // Get simHits of each superlayer
     PSimHitContainer simHits = simHitsPerSl[(*slId)];
 
     // Map simhits per wire
     map<DTWireId, PSimHitContainer> simHitsPerWire = DTHitQualityUtils::mapSimHitsPerWire(simHits);
     map<DTWireId, const PSimHit *> muSimHitPerWire = DTHitQualityUtils::mapMuSimHitsPerWire(simHitsPerWire);
     int nMuSimHit = muSimHitPerWire.size();
     if (nMuSimHit == 0 or nMuSimHit == 1) {
       if (debug_ and nMuSimHit == 1) {
         cout << "[DTSegment2DQuality] Only " << nMuSimHit << " mu SimHit in this SL, skipping !" << endl;
       }
       continue;  // If no or only one mu SimHit is found skip this SL
     }
     if (debug_) {
       cout << "=== SL " << (*slId) << " has " << nMuSimHit << " SimHits" << endl;
     }
 
     // Find outer and inner mu SimHit to build a segment
     pair<const PSimHit *, const PSimHit *> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire);
     // Check that outermost and innermost SimHit are not the same
     if (inAndOutSimHit.first == inAndOutSimHit.second) {
       cout << "[DTHitQualityUtils]***Warning: outermost and innermost SimHit "
               "are the same !"
            << endl;
       continue;
     }
 
     // Find direction and position of the sim Segment in SL RF
     pair<LocalVector, LocalPoint> dirAndPosSimSegm =
         DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit, (*slId), &(*dtGeom));
 
     LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
     LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
     if (debug_) {
       cout << "  Simulated segment:  local direction " << simSegmLocalDir << endl
            << "                      local position  " << simSegmLocalPos << endl;
     }
     const DTSuperLayer *superLayer = dtGeom->superLayer(*slId);
     GlobalPoint simSegmGlobalPos = superLayer->toGlobal(simSegmLocalPos);
 
     // Atan(x/z) angle and x position in SL RF
     float angleSimSeg = DTHitQualityUtils::findSegmentAlphaAndBeta(simSegmLocalDir).first;
     float posSimSeg = simSegmLocalPos.x();
     // Position (in eta, phi coordinates) in the global RF
     float etaSimSeg = simSegmGlobalPos.eta();
     float phiSimSeg = simSegmGlobalPos.phi();
 
     //---------------------------- recHits --------------------------//
     // Get the range of rechit for the corresponding slId
     bool recHitFound = false;
     DTRecSegment2DCollection::range range = segment2Ds->get(*slId);
     int nsegm = distance(range.first, range.second);
     if (debug_) {
       cout << "   Sl: " << *slId << " has " << nsegm << " 2D segments" << endl;
     }
 
     if (nsegm != 0) {
       // Find the best RecHit: look for the 2D RecHit with the angle closest
       // to that of segment made of SimHits.
       // RecHits must have delta alpha and delta position within 5 sigma of
       // the residual distribution (we are looking for residuals of segments
       // usefull to the track fit) for efficency purpose
       const DTRecSegment2D *bestRecHit = nullptr;
       bool bestRecHitFound = false;
       double deltaAlpha = 99999;
 
       // Loop over the recHits of this slId
       for (DTRecSegment2DCollection::const_iterator segment2D = range.first; segment2D != range.second; ++segment2D) {
         // Check the dimension
         if ((*segment2D).dimension() != 2) {
           if (debug_) {
             cout << "[DTSegment2DQuality]***Error: This is not 2D segment !!!" << endl;
           }
           abort();
         }
         // Segment Local Direction and position (in SL RF)
         LocalVector recSegDirection = (*segment2D).localDirection();
         LocalPoint recSegPosition = (*segment2D).localPosition();
 
         float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
         if (debug_) {
           cout << "  RecSegment direction: " << recSegDirection << endl
                << "             position : " << recSegPosition << endl
                << "             alpha    : " << recSegAlpha << endl;
         }
 
         if (fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
           deltaAlpha = fabs(recSegAlpha - angleSimSeg);
           bestRecHit = &(*segment2D);
           bestRecHitFound = true;
         }
       }  // End of Loop over all 2D RecHits
 
       if (bestRecHitFound) {
         // Best rechit direction and position in SL RF
         LocalPoint bestRecHitLocalPos = bestRecHit->localPosition();
         LocalVector bestRecHitLocalDir = bestRecHit->localDirection();
 
         LocalError bestRecHitLocalPosErr = bestRecHit->localPositionError();
         LocalError bestRecHitLocalDirErr = bestRecHit->localDirectionError();
 
         float angleBestRHit = DTHitQualityUtils::findSegmentAlphaAndBeta(bestRecHitLocalDir).first;
 
         if (fabs(angleBestRHit - angleSimSeg) < 5 * sigmaResAngle_ and
             fabs(bestRecHitLocalPos.x() - posSimSeg) < 5 * sigmaResPos_) {
           recHitFound = true;
         }
 
         // Fill Residual histos
         HRes2DHit *hRes = nullptr;
         if ((*slId).superlayer() == 1 or (*slId).superlayer() == 3) {  // RPhi SL
           hRes = histograms.h2DHitRPhi.get();
         } else if ((*slId).superlayer() == 2) {  // RZ SL
           histograms.h2DHitRZ->fill(angleSimSeg,
                                     angleBestRHit,
                                     posSimSeg,
                                     bestRecHitLocalPos.x(),
                                     etaSimSeg,
                                     phiSimSeg,
                                     sqrt(bestRecHitLocalPosErr.xx()),
                                     sqrt(bestRecHitLocalDirErr.xx()));
           if (abs((*slId).wheel()) == 0) {
             hRes = histograms.h2DHitRZ_W0.get();
           } else if (abs((*slId).wheel()) == 1) {
             hRes = histograms.h2DHitRZ_W1.get();
           } else if (abs((*slId).wheel()) == 2) {
             hRes = histograms.h2DHitRZ_W2.get();
           }
         }
         hRes->fill(angleSimSeg,
                    angleBestRHit,
                    posSimSeg,
                    bestRecHitLocalPos.x(),
                    etaSimSeg,
                    phiSimSeg,
                    sqrt(bestRecHitLocalPosErr.xx()),
                    sqrt(bestRecHitLocalDirErr.xx()));
       }
     }  // end of if (nsegm != 0)
 
     // Fill Efficiency plot
     HEff2DHit *hEff = nullptr;
     if ((*slId).superlayer() == 1 or (*slId).superlayer() == 3) {  // RPhi SL
       hEff = histograms.h2DHitEff_RPhi.get();
     } else if ((*slId).superlayer() == 2) {  // RZ SL
       histograms.h2DHitEff_RZ->fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
       if (abs((*slId).wheel()) == 0) {
         hEff = histograms.h2DHitEff_RZ_W0.get();
       } else if (abs((*slId).wheel()) == 1) {
         hEff = histograms.h2DHitEff_RZ_W1.get();
       } else if (abs((*slId).wheel()) == 2) {
         hEff = histograms.h2DHitEff_RZ_W2.get();
       }
     }
     hEff->fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
   }  // End of loop over superlayers
 }

Member Data Documentation

bool DTSegment2DQuality::debug_

private

Definition at line 62 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::segment2DLabel_

private

Definition at line 51 of file DTSegment2DQuality.h.

edm::EDGetTokenT<DTRecSegment2DCollection> DTSegment2DQuality::segment2DToken_

private

Definition at line 53 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResAngle_

private

Definition at line 59 of file DTSegment2DQuality.h.

double DTSegment2DQuality::sigmaResPos_

private

Definition at line 56 of file DTSegment2DQuality.h.

edm::InputTag DTSegment2DQuality::simHitLabel_

private

Definition at line 50 of file DTSegment2DQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTSegment2DQuality::simHitToken_

private

Definition at line 52 of file DTSegment2DQuality.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation