#include <trackCount/DuplicateRecHits/src/DuplicateRecHits.cc>

Inheritance diagram for DuplicateRecHits:

Public Member Functions
	DuplicateRecHits (const edm::ParameterSet &)

	~DuplicateRecHits () override

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

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob () override

void	beginRun (const edm::Run &, const edm::EventSetup &) override

void	endJob () override

void	endRun (const edm::Run &, const edm::EventSetup &) override

Private Attributes
const TransientTrackingRecHitBuilder *	m_builder

std::string	m_buildername

TH1F *	m_nduplicate

TH1F *	m_nduplmod

edm::EDGetTokenT< reco::TrackCollection >	m_trkcollToken

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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

Description: <one line="" class="" summary>="">

Implementation: <Notes on="" implementation>="">

Definition at line 63 of file DuplicateRecHits.cc.

Constructor & Destructor Documentation

DuplicateRecHits::DuplicateRecHits ( const edm::ParameterSet & iConfig )

explicit

Definition at line 98 of file DuplicateRecHits.cc.

References edm::ParameterSet::getParameter(), label, m_nduplicate, m_nduplmod, and TFileService::make().

                                                                 :
   m_trkcollToken(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("trackCollection"))),
   m_buildername(iConfig.getParameter<std::string>("TTRHBuilder"))
   
 {
    //now do what ever initialization is needed
 
   // histogram parameters
 
   edm::LogInfo("TrackCollection") << "Using collection " << iConfig.getParameter<edm::InputTag>("trackCollection").label().c_str() ;
 
 
   edm::Service<TFileService> tfserv;
 
   m_nduplicate = tfserv->make<TH1F>("nduplicate","Number of duplicated clusters per track",10,-0.5,9.5);
   m_nduplmod = tfserv->make<TH1F>("nduplmod","Number of duplicated clusters per module",10,-0.5,9.5);
   m_nduplmod->SetCanExtend(TH1::kXaxis);
 }

DuplicateRecHits::~DuplicateRecHits ( )

override

Definition at line 118 of file DuplicateRecHits.cc.

 {
 
    // do anything here that needs to be done at desctruction time
    // (e.g. close files, deallocate resources etc.)
 
 }

Member Function Documentation

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

overrideprivate

Definition at line 133 of file DuplicateRecHits.cc.

References TransientTrackingRecHitBuilder::build(), SiPixelRecHit::cluster(), fastPrimaryVertexProducer_cfi::clusters, TrackingRecHit::det(), GeomDet::geographicalId(), edm::Event::getByToken(), edm::Ref< C, T, F >::index(), LogDebug, LogTrace, m_builder, m_nduplicate, m_nduplmod, m_trkcollToken, DetId::rawId(), and l1t::tracks.

 {
    using namespace edm;
    edm::Service<TFileService> tfserv;
 
    Handle<reco::TrackCollection> tracks;
    iEvent.getByToken(m_trkcollToken,tracks);
 
 
    for(reco::TrackCollection::const_iterator it = tracks->begin();it!=tracks->end();it++) {
 
      std::set<SiPixelRecHit::ClusterRef::key_type> clusters;
      int nduplicate=0;
      for(trackingRecHit_iterator rh= it->recHitsBegin(); rh!=it->recHitsEnd();++rh) {
        TransientTrackingRecHit::RecHitPointer ttrh = m_builder->build(&**rh);
         const SiPixelRecHit * pxrh = dynamic_cast<const SiPixelRecHit*>(ttrh->hit());
         if (pxrh) {
       //      LogTrace("DuplicateHitFinder") << ttrh->det()->geographicalId() << " " << pxrh->cluster().index();
       if(clusters.find(pxrh->cluster().index())!=clusters.end()) {
         nduplicate++;
         std::stringstream detidstr;
         detidstr << ttrh->det()->geographicalId().rawId();
         m_nduplmod->Fill(detidstr.str().c_str(),1.);
         LogDebug("DuplicateHitFinder") << "Track with " << it->recHitsSize() << " RecHits";
         LogTrace("DuplicateHitFinder") << "Duplicate found " << ttrh->det()->geographicalId().rawId() << " " << pxrh->cluster().index();
       }
       clusters.insert(pxrh->cluster().index());
     }
     m_nduplicate->Fill(nduplicate);
      }
 
    }
 
 }

void DuplicateRecHits::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 185 of file DuplicateRecHits.cc.

186 {

187 }

void DuplicateRecHits::beginRun	(	const edm::Run &	iRun,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 169 of file DuplicateRecHits.cc.

References edm::EventSetup::get(), m_builder, m_buildername, and edm::ESHandle< T >::product().

 {
   edm::ESHandle<TransientTrackingRecHitBuilder> theBuilder;
   iSetup.get<TransientRecHitRecord>().get(m_buildername, theBuilder);
   m_builder = theBuilder.product();
 }

void DuplicateRecHits::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 191 of file DuplicateRecHits.cc.

References DEFINE_FWK_MODULE.

192 {

193 }

void DuplicateRecHits::endRun	(	const edm::Run &	iRun,
		const edm::EventSetup &
	)

overrideprivate

Definition at line 177 of file DuplicateRecHits.cc.

178 {

179 }

Member Data Documentation

const TransientTrackingRecHitBuilder* DuplicateRecHits::m_builder

private

Definition at line 80 of file DuplicateRecHits.cc.

Referenced by analyze(), and beginRun().

std::string DuplicateRecHits::m_buildername

private

Definition at line 79 of file DuplicateRecHits.cc.

Referenced by beginRun().

TH1F* DuplicateRecHits::m_nduplicate

private

Definition at line 82 of file DuplicateRecHits.cc.

Referenced by analyze(), and DuplicateRecHits().

TH1F* DuplicateRecHits::m_nduplmod

private

Definition at line 83 of file DuplicateRecHits.cc.

Referenced by analyze(), and DuplicateRecHits().

edm::EDGetTokenT<reco::TrackCollection> DuplicateRecHits::m_trkcollToken

private

Definition at line 78 of file DuplicateRecHits.cc.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation