#include <DebugTools/OverlapProblem/plugins/OverlapProblemTPAnalyzer.cc>

Inheritance diagram for OverlapProblemTPAnalyzer:

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

	~OverlapProblemTPAnalyzer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::WatchRuns, edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const noexcept final

bool	wantsGlobalRuns () const noexcept final

bool	wantsInputProcessBlocks () const noexcept final

bool	wantsProcessBlocks () const noexcept final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const noexcept

bool	wantsStreamRuns () const noexcept

	~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

ESResolverIndex const *	esGetTokenIndices (edm::Transition iTrans) const

std::vector< ESResolverIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

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

void	updateLookup (eventsetup::ESRecordsToProductResolverIndices const &)

virtual	~EDConsumerBase () noexcept(false)

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

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

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

Private Attributes
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	m_associatorToken

std::vector< TH1F * >	m_assosimhitytecr

TH1F *	m_etatp

TH1F *	m_llbit

TH1F *	m_nassotk

TH1F *	m_nhits

TH1F *	m_nrechits

TH2F *	m_nrecvssimhits

TH1F *	m_pdgid

TH1F *	m_ptp

std::vector< TH1F * >	m_simhitytecr

TH1F *	m_statustp

edm::EDGetTokenT< TrackingParticleCollection >	m_tpcollToken

edm::EDGetTokenT< edm::View< reco::Track > >	m_trkcollToken

Additional Inherited Members
Public Types inherited from edm::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

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

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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

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

Definition at line 58 of file OverlapProblemTPAnalyzer.cc.

Constructor & Destructor Documentation

◆ OverlapProblemTPAnalyzer()

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

explicit

Definition at line 99 of file OverlapProblemTPAnalyzer.cc.

References TFileService::kSharedResource, m_assosimhitytecr, m_etatp, m_llbit, m_nassotk, m_nhits, m_nrechits, m_nrecvssimhits, m_pdgid, m_ptp, m_simhitytecr, m_statustp, TFileService::make(), Skims_PA_cff::name, relativeConstraints::ring, and runGCPTkAlMap::title.

     : m_simhitytecr(),
       m_assosimhitytecr(),
       m_tpcollToken(
           consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("trackingParticlesCollection"))),
       m_trkcollToken(consumes<edm::View<reco::Track> >(iConfig.getParameter<edm::InputTag>("trackCollection"))),
       m_associatorToken(consumes<reco::TrackToTrackingParticleAssociator>(edm::InputTag("trackAssociatorByHits"))) {
   //now do what ever initialization is needed
 
   usesResource(TFileService::kSharedResource);
 
   edm::Service<TFileService> tfserv;
 
   m_ptp = tfserv->make<TH1F>("tpmomentum", "Tracking Particle momentum", 100, 0., 200.);
   m_etatp = tfserv->make<TH1F>("tpeta", "Tracking Particle pseudorapidity", 100, -4., 4.);
   m_nhits = tfserv->make<TH1F>("nhits", "Tracking Particle associated hits", 100, -0.5, 99.5);
   m_nrechits = tfserv->make<TH1F>("nrechits", "Tracking Particle associated rec hits", 100, -0.5, 99.5);
   m_nrecvssimhits = tfserv->make<TH2F>(
       "nrecvssimhits", "Tracking Particle associated rec hits vs sim hits", 100, -0.5, 99.5, 100, -0.5, 99.5);
   m_nassotk = tfserv->make<TH1F>("nassotk", "Number of assocated reco tracks", 10, -0.5, 9.5);
 
   m_pdgid = tfserv->make<TH1F>("pdgid", "Tracking Particle PDG id", 1000, -500.5, 499.5);
   m_llbit = tfserv->make<TH1F>("llbit", "Tracking Particle LongLived bit", 2, -0.5, 1.5);
   m_statustp = tfserv->make<TH1F>("statustp", "Tracking Particle status", 2000, -1000.5, 999.5);
 
   for (unsigned int ring = 0; ring < 7; ++ring) {
     char name[100];
     char title[100];
 
     sprintf(name, "simytecr_%d", ring + 1);
     sprintf(title, "SimHit local Y TEC ring %d", ring + 1);
 
     m_simhitytecr.push_back(tfserv->make<TH1F>(name, title, 200, -20., 20.));
 
     sprintf(name, "assosimytecr_%d", ring + 1);
     sprintf(title, "SimHit local Y TEC ring %d with associated RecHit", ring + 1);
 
     m_assosimhitytecr.push_back(tfserv->make<TH1F>(name, title, 200, -20., 20.));
   }
 }

◆ ~OverlapProblemTPAnalyzer()

OverlapProblemTPAnalyzer::~OverlapProblemTPAnalyzer ( )

override

Definition at line 140 of file OverlapProblemTPAnalyzer.cc.

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

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 150 of file OverlapProblemTPAnalyzer.cc.

References funct::abs(), reco::TrackToTrackingParticleAssociator::associateSimToReco(), edm::AssociationMap< Tag >::end(), edm::AssociationMap< Tag >::find(), mps_fire::i, iEvent, LogDebug, LogTrace, m_associatorToken, m_etatp, m_llbit, m_nassotk, m_nhits, m_nrechits, m_nrecvssimhits, m_pdgid, m_ptp, m_statustp, m_tpcollToken, m_trkcollToken, cmsswSequenceInfo::tp, and hltEgammaHLTExtra_cfi::trks.

                                                                                           {
   using namespace edm;
 
   // reco track Handle
 
   //   Handle<reco::TrackCollection> trkcoll;
   Handle<edm::View<reco::Track> > trkcoll;
   iEvent.getByToken(m_trkcollToken, trkcoll);
 
   Handle<TrackingParticleCollection> tpcoll;
   iEvent.getByToken(m_tpcollToken, tpcoll);
 
   Handle<reco::TrackToTrackingParticleAssociator> tahandle;
   iEvent.getByToken(m_associatorToken, tahandle);
 
   // associate reco to sim tracks
 
   reco::SimToRecoCollection srcoll = tahandle->associateSimToReco(trkcoll, tpcoll);
 
   // loop on Handle with index and use find
 
   for (unsigned int index = 0; index != tpcoll->size(); ++index) {
     // get TrackingParticleRef
 
     const TrackingParticleRef tp(tpcoll, index);
 
     if (std::abs(tp->pdgId()) != 13)
       continue;
 
     // get the SimHIt from tracker only
 
     // Commented since the new TP's do not have this method
     //    std::vector<PSimHit> tksimhits = tp->trackPSimHit(DetId::Tracker);
 
     m_ptp->Fill(tp->p());
     m_etatp->Fill(tp->eta());
     //     m_nhits->Fill(tp->matchedHit());
     // With the new Tracking Particles I have to use a different method
     //    m_nhits->Fill(tksimhits.size());
     m_nhits->Fill(tp->numberOfTrackerHits());
 
     m_pdgid->Fill(tp->pdgId());
     m_llbit->Fill(tp->longLived());
     m_statustp->Fill(tp->status());
 
     // prepare a vector of TrackingRecHit from the associated reco tracks
 
     std::vector<DetId> rechits;
 
     // look at associated tracks
 
     if (srcoll.find(tp) != srcoll.end()) {
       reco::SimToRecoCollection::result_type trks = srcoll[tp];
       m_nassotk->Fill(trks.size());
 
       // loop on associated tracks and fill TrackingRecHit vector
       for (reco::SimToRecoCollection::result_type::const_iterator trk = trks.begin(); trk != trks.end(); ++trk) {
         for (trackingRecHit_iterator rh = trk->first->recHitsBegin(); rh != trk->first->recHitsEnd(); ++rh) {
           rechits.push_back((*rh)->geographicalId());
         }
       }
 
     } else {
       m_nassotk->Fill(0.);
       edm::LogInfo("NoAssociatedRecoTrack")
           << "No associated reco track for TP with p = " << tp->p() << " and eta = " << tp->eta();
     }
 
     m_nrechits->Fill(rechits.size());
     // new method used to be compliant with the new TP's
     m_nrecvssimhits->Fill(tp->numberOfTrackerHits(), rechits.size());
 
     LogDebug("RecHitDetId") << "List of " << rechits.size() << " rechits detid from muon with p = " << tp->p()
                             << "and eta = " << tp->eta();
     for (unsigned int i = 0; i < rechits.size(); ++i) {
       LogTrace("RecHitDetId") << rechits[i].rawId();
     }
 
     // loop on sim hits
 
     LogDebug("SimHitDetId") << "List of " << tp->numberOfTrackerHits()
                             << " simhits detid from muon with p = " << tp->p() << "and eta = " << tp->eta();
   }
 }