#include <SiPixelPhase1TrackEfficiency.h>

Inheritance diagram for SiPixelPhase1TrackEfficiency:

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

	SiPixelPhase1TrackEfficiency (const edm::ParameterSet &conf)

Public Member Functions inherited from SiPixelPhase1Base
void	analyze (edm::Event const &e, edm::EventSetup const &eSetup) override=0

void	bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &iSetup) override

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

	~SiPixelPhase1Base () override

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

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

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

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

	~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)

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

Private Types
enum	{ VALID, MISSING, INACTIVE, EFFICIENCY, VERTICES }

Private Attributes
bool	applyVertexCut_

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

edm::EDGetTokenT< reco::TrackCollection >	tracksToken_

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 Types inherited from SiPixelPhase1Base
enum	{ DCS }

Protected Member Functions inherited from SiPixelPhase1Base
bool	checktrigger (const edm::Event &iEvent, const edm::EventSetup &iSetup, const unsigned trgidx) const

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
VALID
MISSING
INACTIVE
EFFICIENCY
VERTICES

Definition at line 17 of file SiPixelPhase1TrackEfficiency.h.

        {
     VALID,
     MISSING,
     INACTIVE,
     EFFICIENCY,
     VERTICES
   };

Constructor & Destructor Documentation

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

explicit

Definition at line 24 of file SiPixelPhase1TrackEfficiency.cc.

References applyVertexCut_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), tracksToken_, and vtxToken_.

                                                                                          :
   SiPixelPhase1Base(iConfig) 
 { 
   tracksToken_ = consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracks"));
   vtxToken_ = consumes<reco::VertexCollection>(iConfig.getParameter<edm::InputTag>("primaryvertices"));
   applyVertexCut_=iConfig.getUntrackedParameter<bool>("VertexCut",true);
 
 }

Member Function Documentation

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

override

Definition at line 33 of file SiPixelPhase1TrackEfficiency.cc.

References funct::abs(), applyVertexCut_, SiPixelPhase1Base::checktrigger(), SiPixelPhase1Base::DCS, DEFINE_FWK_MODULE, EFFICIENCY, edm::EventSetup::get(), edm::Event::getByToken(), h, hcalSimParameters_cfi::hb, HistogramManagerHolder::histo, iEvent, INACTIVE, TrackingRecHit::inactive, edm::ESHandleBase::isValid(), edm::HandleBase::isValid(), MISSING, TrackingRecHit::missing, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, StripSubdetector::TEC, StripSubdetector::TIB, StripSubdetector::TID, StripSubdetector::TOB, HiIsolationCommonParameters_cff::track, mixOne_simraw_on_sim_cfi::tracker, l1t::tracks, tracksToken_, VALID, TrackingRecHit::valid, electrons_cff::vertices, VERTICES, and vtxToken_.

                                                                                               {
   if( !checktrigger(iEvent,iSetup,DCS) ) return;
 
   // 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()) return;
 
   histo[VERTICES].fill(vertices->size(),DetId(0),&iEvent);
 
   if (applyVertexCut_ &&  vertices->empty()) return;
 
   // should be used for weird cuts
   //const auto primaryVertex = vertices->at(0); 
 
   // get the map
   edm::Handle<reco::TrackCollection> tracks;
   iEvent.getByToken( tracksToken_, tracks);
   if (!tracks.isValid()) return;
 
   for (auto const & track : *tracks) {
 
     //this cut is needed to be consisten with residuals calculation
     if (applyVertexCut_ && (track.pt() < 0.75 || std::abs( track.dxy(vertices->at(0).position()) ) > 5*track.dxyError())) continue; 
 
     bool isBpixtrack = false, isFpixtrack = false;
     int nStripHits = 0;
 
     // first, look at the full track to see whether it is good
     // auto const & trajParams = track.extra()->trajParams();
     auto hb = track.recHitsBegin();
     for(unsigned int h=0;h<track.recHitsSize();h++){
       
       auto hit = *(hb+h);
       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(unsigned int h=0;h<track.recHitsSize();h++){
       auto hit = *(hb+h);
 
       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;
       bool isHitInactive = hit->getType()==TrackingRecHit::inactive;
 
       /*
       const SiPixelRecHit* pixhit = dynamic_cast<const SiPixelRecHit*>(hit);
       const PixelGeomDetUnit* geomdetunit = dynamic_cast<const PixelGeomDetUnit*> ( tracker->idToDet(id) );
       const PixelTopology& topol = geomdetunit->specificTopology();
 
       // this commented part is useful if one wants ROC level maps of hits, however the local position may fall out of a ROC and the ROC maps will look very strange (with no white cross)
       LocalPoint lp;
 
       if (pixhit) {
         lp = pixhit->localPosition();
       } else {
         lp = trajParams[h].position();
       }
 
       MeasurementPoint mp = topol.measurementPosition(lp);
       int row = (int) mp.x();
       int col = (int) mp.y();
       */
 
       if (isHitValid)   {
         histo[VALID].fill(id, &iEvent);
         histo[EFFICIENCY].fill(1, id, &iEvent);
       }
       if (isHitMissing) {
         histo[MISSING].fill(id, &iEvent);
         histo[EFFICIENCY].fill(0, id, &iEvent);
       }
       if (isHitInactive)   {
         histo[INACTIVE].fill(id, &iEvent);
       }
     }
   }
   histo[VALID   ].executePerEventHarvesting(&iEvent);
   histo[MISSING ].executePerEventHarvesting(&iEvent);
   histo[INACTIVE].executePerEventHarvesting(&iEvent);
 }