#include <DTSegment2DSLPhiQuality.h>

Inheritance diagram for DTSegment2DSLPhiQuality:

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

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

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

Private Attributes
bool	debug_

bool	doall_

bool	local_

edm::InputTag	segment4DLabel_

edm::EDGetTokenT< DTRecSegment4DCollection >	segment4DToken_

double	sigmaResAngle_

double	sigmaResPos_

edm::InputTag	simHitLabel_

edm::EDGetTokenT< edm::PSimHitContainer >	simHitToken_

Additional Inherited Members
Public Types inherited from DQMGlobalEDAnalyzer< dtsegment2dsl::Histograms >
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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 reconstructed with both SL Phi and plot resolution comparing reconstructed and simulated quantities

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

Definition at line 32 of file DTSegment2DSLPhiQuality.h.

Constructor & Destructor Documentation

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

Constructor.

Definition at line 38 of file DTSegment2DSLPhiQuality.cc.

References 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
   segment4DLabel_ = pset.getUntrackedParameter<InputTag>("segment4DLabel");
   segment4DToken_ = consumes<DTRecSegment4DCollection>(pset.getUntrackedParameter<InputTag>("segment4DLabel"));
 
   // sigma resolution on position
   sigmaResPos_ = pset.getParameter<double>("sigmaResPos");
   // sigma resolution on angle
   sigmaResAngle_ = pset.getParameter<double>("sigmaResAngle");
   doall_ = pset.getUntrackedParameter<bool>("doall", false);
   local_ = pset.getUntrackedParameter<bool>("local", false);
 }

Member Function Documentation

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

overrideprivatevirtual

Book the DQM plots.

Implements DQMGlobalEDAnalyzer< dtsegment2dsl::Histograms >.

Definition at line 58 of file DTSegment2DSLPhiQuality.cc.

References dtsegment2dsl::Histograms::h2DHitEff_SuperPhi, and dtsegment2dsl::Histograms::h2DHitSuperPhi.

                                                                            {
   // Book the histos
   histograms.h2DHitSuperPhi = std::make_unique<HRes2DHit>("SuperPhi", booker, doall_, local_);
   if (doall_) {
     histograms.h2DHitEff_SuperPhi = std::make_unique<HEff2DHit>("SuperPhi", booker);
   }
 }

void DTSegment2DSLPhiQuality::dqmAnalyze	(	edm::Event const &	event,
		edm::EventSetup const &	setup,
		dtsegment2dsl::Histograms const &	histograms
	)		const

overrideprivate

Perform the real analysis.

Definition at line 70 of file DTSegment2DSLPhiQuality.cc.

References relativeConstraints::chamber, DTGeometry::chamber(), gather_cfg::cout, DEFINE_FWK_MODULE, HLT_2018_cff::distance, PV3DBase< T, PVType, FrameType >::eta(), DTHitQualityUtils::findMuSimSegment(), DTHitQualityUtils::findMuSimSegmentDirAndPos(), DTHitQualityUtils::findSegmentAlphaAndBeta(), edm::EventSetup::get(), dtsegment2dsl::Histograms::h2DHitEff_SuperPhi, dtsegment2dsl::Histograms::h2DHitSuperPhi, edm::HandleBase::isValid(), DTRecSegment2D::localDirection(), DTRecSegment2D::localDirectionError(), DTRecSegment2D::localPosition(), DTRecSegment2D::localPositionError(), DTHitQualityUtils::mapMuSimHitsPerWire(), DTHitQualityUtils::mapSimHitsPerWire(), PV3DBase< T, PVType, FrameType >::phi(), FastTimerService_cff::range, rpcPointValidation_cfi::simHit, FastTrackerRecHitCombiner_cfi::simHits, mathSSE::sqrt(), GeomDet::toGlobal(), PV3DBase< T, PVType, FrameType >::x(), and LocalError::xx().

                                                                              {
   // 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 chamber
   map<DTChamberId, PSimHitContainer> simHitsPerCh;
   for (const auto &simHit : *simHits) {
     // Create the id of the chamber (the simHits in the DT known their wireId)
     DTChamberId chamberId = (((DTWireId(simHit.detUnitId())).layerId()).superlayerId()).chamberId();
     // Fill the map
     simHitsPerCh[chamberId].push_back(simHit);
   }
 
   // Get the 4D rechits from the event
   Handle<DTRecSegment4DCollection> segment4Ds;
   event.getByToken(segment4DToken_, segment4Ds);
 
   if (!segment4Ds.isValid()) {
     if (debug_) {
       cout << "[DTSegment2DSLPhiQuality]**Warning: no 4D Segments with label: " << segment4DLabel_
            << " in this event, skipping!" << endl;
     }
     return;
   }
 
   // Loop over all chambers containing a segment
   DTRecSegment4DCollection::id_iterator chamberId;
   for (chamberId = segment4Ds->id_begin(); chamberId != segment4Ds->id_end(); ++chamberId) {
     //------------------------- simHits ---------------------------//
     // Get simHits of each chamber
     PSimHitContainer simHits = simHitsPerCh[(*chamberId)];
 
     // 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 || nMuSimHit == 1) {
       if (debug_ && nMuSimHit == 1) {
         cout << "[DTSegment2DSLPhiQuality] Only " << nMuSimHit << " mu SimHit in this chamber, skipping!" << endl;
       }
       continue;  // If no or only one mu SimHit is found skip this chamber
     }
     if (debug_) {
       cout << "=== Chamber " << (*chamberId) << " has " << nMuSimHit << " SimHits" << endl;
     }
 
     // Find outer and inner mu SimHit to build a segment
     pair<const PSimHit *, const PSimHit *> inAndOutSimHit = DTHitQualityUtils::findMuSimSegment(muSimHitPerWire);
 
     // Find direction and position of the sim Segment in Chamber RF
     pair<LocalVector, LocalPoint> dirAndPosSimSegm =
         DTHitQualityUtils::findMuSimSegmentDirAndPos(inAndOutSimHit, (*chamberId), &(*dtGeom));
 
     LocalVector simSegmLocalDir = dirAndPosSimSegm.first;
     LocalPoint simSegmLocalPos = dirAndPosSimSegm.second;
     const DTChamber *chamber = dtGeom->chamber(*chamberId);
     GlobalPoint simSegmGlobalPos = chamber->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 lobal RF
     float etaSimSeg = simSegmGlobalPos.eta();
     float phiSimSeg = simSegmGlobalPos.phi();
 
     if (debug_) {
       cout << "  Simulated segment:  local direction " << simSegmLocalDir << endl
            << "                      local position  " << simSegmLocalPos << endl
            << "                      angle           " << angleSimSeg << endl;
     }
 
     //---------------------------- recHits --------------------------//
     // Get the range of rechit for the corresponding chamberId
     bool recHitFound = false;
     DTRecSegment4DCollection::range range = segment4Ds->get(*chamberId);
     int nsegm = distance(range.first, range.second);
     if (debug_) {
       cout << "   Chamber: " << *chamberId << " has " << nsegm << " 4D segments" << endl;
     }
 
     if (nsegm != 0) {
       // Find the best RecHit: look for the 4D RecHit with the phi 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 chamberId
       for (DTRecSegment4DCollection::const_iterator segment4D = range.first; segment4D != range.second; ++segment4D) {
         // Check the dimension
         if ((*segment4D).dimension() != 4) {
           if (debug_) {
             cout << "[DTSegment2DSLPhiQuality]***Error: This is not 4D "
                     "segment!!!"
                  << endl;
           }
           continue;
         }
 
         // Get 2D superPhi segments from 4D segments
         const DTChamberRecSegment2D *phiSegment2D = (*segment4D).phiSegment();
         if ((*phiSegment2D).dimension() != 2) {
           if (debug_) {
             cout << "[DTSegment2DQuality]***Error: This is not 2D segment!!!" << endl;
           }
           abort();
         }
 
         // Segment Local Direction and position (in Chamber RF)
         LocalVector recSegDirection = (*phiSegment2D).localDirection();
 
         float recSegAlpha = DTHitQualityUtils::findSegmentAlphaAndBeta(recSegDirection).first;
         if (debug_) {
           cout << "  RecSegment direction: " << recSegDirection << endl
                << "             position : " << (*phiSegment2D).localPosition() << endl
                << "             alpha    : " << recSegAlpha << endl;
         }
 
         if (fabs(recSegAlpha - angleSimSeg) < deltaAlpha) {
           deltaAlpha = fabs(recSegAlpha - angleSimSeg);
           bestRecHit = &(*phiSegment2D);
           bestRecHitFound = true;
         }
       }  // End of Loop over all 4D RecHits of this chambers
 
       if (bestRecHitFound) {
         // Best rechit direction and position in Chamber 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_ &&
             fabs(bestRecHitLocalPos.x() - posSimSeg) < 5 * sigmaResPos_) {
           recHitFound = true;
         }
 
         // Fill Residual histos
         histograms.h2DHitSuperPhi->fill(angleSimSeg,
                                         angleBestRHit,
                                         posSimSeg,
                                         bestRecHitLocalPos.x(),
                                         etaSimSeg,
                                         phiSimSeg,
                                         sqrt(bestRecHitLocalPosErr.xx()),
                                         sqrt(bestRecHitLocalDirErr.xx()));
       }
     }  // end of if (nsegm!= 0)
 
     // Fill Efficiency plot
     if (doall_) {
       histograms.h2DHitEff_SuperPhi->fill(etaSimSeg, phiSimSeg, posSimSeg, angleSimSeg, recHitFound);
     }
   }  // End of loop over chambers
 }

Member Data Documentation

bool DTSegment2DSLPhiQuality::debug_

private

Definition at line 64 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::doall_

private

Definition at line 60 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::local_

private

Definition at line 61 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::segment4DLabel_

private

Definition at line 50 of file DTSegment2DSLPhiQuality.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegment2DSLPhiQuality::segment4DToken_

private

Definition at line 52 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResAngle_

private

Definition at line 58 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResPos_

private

Definition at line 55 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::simHitLabel_

private

Definition at line 49 of file DTSegment2DSLPhiQuality.h.

edm::EDGetTokenT<edm::PSimHitContainer> DTSegment2DSLPhiQuality::simHitToken_

private

Definition at line 51 of file DTSegment2DSLPhiQuality.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation