#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 final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const 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

bool	wantsStreamRuns () 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

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

std::vector< ESProxyIndex > 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

bool	registeredToConsumeMany (TypeID const &, 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::ESRecordsToProxyIndices 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 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<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, HI_PhotonSkim_cff::rechits, 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();
   }
 }