#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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	bookHistograms (DQMStore::ConcurrentBooker &, 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 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 ESProduct , Transition Tr = Transition::Event>
auto	esConsumes (eventsetup::EventSetupRecordKey const &, 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 33 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::ConcurrentBooker &	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 66 of file DTSegment2DSLPhiQuality.cc.

References relativeConstraints::chamber, DTGeometry::chamber(), gather_cfg::cout, DEFINE_FWK_MODULE, SoftLeptonByDistance_cfi::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(), rpcPointValidation_cfi::simHit, trackerHits::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 62 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::doall_

private

Definition at line 58 of file DTSegment2DSLPhiQuality.h.

bool DTSegment2DSLPhiQuality::local_

private

Definition at line 59 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::segment4DLabel_

private

Definition at line 48 of file DTSegment2DSLPhiQuality.h.

edm::EDGetTokenT<DTRecSegment4DCollection> DTSegment2DSLPhiQuality::segment4DToken_

private

Definition at line 50 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResAngle_

private

Definition at line 56 of file DTSegment2DSLPhiQuality.h.

double DTSegment2DSLPhiQuality::sigmaResPos_

private

Definition at line 53 of file DTSegment2DSLPhiQuality.h.

edm::InputTag DTSegment2DSLPhiQuality::simHitLabel_

private

Definition at line 47 of file DTSegment2DSLPhiQuality.h.

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

private

Definition at line 49 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