#include <SiPixelPhase1TrackEfficiency.h>

Inheritance diagram for SiPixelPhase1TrackEfficiency:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &)

	SiPixelPhase1TrackEfficiency (const edm::ParameterSet &conf)

Public Member Functions inherited from SiPixelPhase1Base
void	bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &iSetup)

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

virtual	~SiPixelPhase1Base ()

Public Member Functions inherited from DQMEDAnalyzer
virtual void	beginRun (edm::Run const &, edm::EventSetup const &) final

virtual void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

virtual void	endLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

virtual void	endRunSummary (edm::Run const &, edm::EventSetup const &, dqmDetails::NoCache *) const final

uint32_t	streamId () const

Public Member Functions inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
	EDAnalyzer ()=default

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDAnalyzerBase ()

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

	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 &	itemsToGetFromEvent () 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

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

virtual	~EDConsumerBase () noexcept(false)

Public Member Functions inherited from HistogramManagerHolder
	HistogramManagerHolder (const edm::ParameterSet &iConfig)

Private Types
enum	{ CLUSTER_PROB, VALID, MISSING, EFFICIENCY }

Private Attributes
edm::EDGetTokenT < edmNew::DetSetVector < SiPixelCluster > >	clustersToken_

edm::EDGetTokenT < TrajTrackAssociationCollection >	trackAssociationToken_

edm::EDGetTokenT < reco::VertexCollection >	vtxToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDAnalyzer< edm::RunSummaryCache< dqmDetails::NoCache >, edm::LuminosityBlockSummaryCache< dqmDetails::NoCache > >
typedef CacheContexts< T...>	CacheTypes

typedef CacheTypes::GlobalCache	GlobalCache

typedef AbilityChecker< T...>	HasAbility

typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache

typedef LuminosityBlockContextT < LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext

typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache

typedef CacheTypes::RunCache	RunCache

typedef RunContextT< RunCache, GlobalCache >	RunContext

typedef CacheTypes::RunSummaryCache	RunSummaryCache

Public Types inherited from edm::stream::EDAnalyzerBase
typedef EDAnalyzerAdaptorBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from DQMEDAnalyzer
static std::shared_ptr < dqmDetails::NoCache >	globalBeginLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const *)

static std::shared_ptr < dqmDetails::NoCache >	globalBeginRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const *)

static void	globalEndLuminosityBlockSummary (edm::LuminosityBlock const &, edm::EventSetup const &, LuminosityBlockContext const , dqmDetails::NoCache )

static void	globalEndRunSummary (edm::Run const &, edm::EventSetup const &, RunContext const , dqmDetails::NoCache )

Static Public Member Functions inherited from edm::stream::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 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)

Protected Attributes inherited from HistogramManagerHolder
GeometryInterface	geometryInterface

std::vector< HistogramManager >	histo

Detailed Description

Definition at line 16 of file SiPixelPhase1TrackEfficiency.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
CLUSTER_PROB
VALID
MISSING
EFFICIENCY

Definition at line 17 of file SiPixelPhase1TrackEfficiency.h.

        {
     CLUSTER_PROB,
     VALID,
     MISSING,
     EFFICIENCY
   };

Constructor & Destructor Documentation

SiPixelPhase1TrackEfficiency::SiPixelPhase1TrackEfficiency ( const edm::ParameterSet & conf )

explicit

Definition at line 24 of file SiPixelPhase1TrackEfficiency.cc.

References edm::ParameterSet::getParameter(), trackAssociationToken_, and vtxToken_.

                                                                                          :
   SiPixelPhase1Base(iConfig) 
 {
   trackAssociationToken_ = consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("trajectories"));
   vtxToken_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryvertices"));
 }

Member Function Documentation

void SiPixelPhase1TrackEfficiency::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::stream::EDAnalyzerBase.

Definition at line 31 of file SiPixelPhase1TrackEfficiency.cc.

References assert(), CLUSTER_PROB, SiPixelRecHit::clusterProbability(), cuy::col, edm::EventSetup::get(), edm::Event::getByToken(), HistogramManagerHolder::histo, edm::ESHandleBase::isValid(), BaseTrackerRecHit::localPositionFast(), Topology::measurementPosition(), MISSING, TrackingRecHit::missing, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, PixelGeomDetUnit::specificTopology(), StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, trackAssociationToken_, patCandidatesForDimuonsSequences_cff::tracker, VALID, TrackingRecHit::valid, HLT_25ns10e33_v2_cff::vertices, vtxToken_, PV2DBase< T, PVType, FrameType >::x(), and PV2DBase< T, PVType, FrameType >::y().

                                                                                               {
 
   // get geometry
   edm::ESHandle<TrackerGeometry> tracker;
   iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
   assert(tracker.isValid());
 
   // get primary vertex
   edm::Handle<reco::VertexCollection> vertices;
   iEvent.getByToken( vtxToken_, vertices);
   if (!vertices.isValid() || vertices->size() == 0) return;
   // should be used for weird cuts
   //const auto primaryVertex = vertices->at(0); 
 
   // get the map
   edm::Handle<TrajTrackAssociationCollection> ttac;
   iEvent.getByToken( trackAssociationToken_, ttac);
 
   for (auto& item : *ttac) {
     auto trajectory_ref = item.key;
     reco::TrackRef track_ref = item.val;
 
     bool isBpixtrack = false, isFpixtrack = false;
     int nStripHits = 0;
 
     // first, look at the full track to see whether it is good
     for (auto& measurement : trajectory_ref->measurements()) {
       // check if things are all valid
       if (!measurement.updatedState().isValid()) continue;
       auto hit = measurement.recHit();
       if (!hit->isValid()) continue;
 
       DetId id = hit->geographicalId();
       uint32_t subdetid = (id.subdetId());
 
       // count strip hits
       if(subdetid==StripSubdetector::TIB) nStripHits++;
       if(subdetid==StripSubdetector::TOB) nStripHits++;
       if(subdetid==StripSubdetector::TID) nStripHits++;
       if(subdetid==StripSubdetector::TEC) nStripHits++;
 
       // check that we are in the pixel
       if (subdetid == PixelSubdetector::PixelBarrel) isBpixtrack = true;
       if (subdetid == PixelSubdetector::PixelEndcap) isFpixtrack = true;
     }
 
     if (!isBpixtrack && !isFpixtrack) continue;
 
     // then, look at each hit
     for (auto& measurement : trajectory_ref->measurements()) {
       if (!measurement.updatedState().isValid()) continue;
       auto hit = measurement.recHit();
 
       DetId id = hit->geographicalId();
       uint32_t subdetid = (id.subdetId());
       if (   subdetid != PixelSubdetector::PixelBarrel 
           && subdetid != PixelSubdetector::PixelEndcap) continue;
 
       bool isHitValid   = hit->getType()==TrackingRecHit::valid;
       bool isHitMissing = hit->getType()==TrackingRecHit::missing;
 
       const SiPixelRecHit* pixhit = dynamic_cast<const SiPixelRecHit*>(hit->hit());
       int row = 0, col = 0;
       if (pixhit) {
         double clusterProbability= pixhit->clusterProbability(0);
         if (clusterProbability > 0) 
           histo[CLUSTER_PROB].fill(log10(clusterProbability), id, &iEvent);
 
         const PixelGeomDetUnit* geomdetunit = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(id) );
         const PixelTopology& topol = geomdetunit->specificTopology();
         LocalPoint const& lp = pixhit->localPositionFast();
         MeasurementPoint mp = topol.measurementPosition(lp);
         row = (int) mp.x();
         col = (int) mp.y();
       }
 
       if (isHitValid)   histo[VALID  ].fill(id, &iEvent, col, row);
       if (isHitMissing) histo[MISSING].fill(id, &iEvent, col, row);
     }
   }
 }