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

Inheritance diagram for OverlapProblemTPAnalyzer:

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

	~OverlapProblemTPAnalyzer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

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

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

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

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

virtual 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::EDAnalyzer
typedef EDAnalyzer	ModuleType

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

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)

Detailed Description

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

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

Definition at line 61 of file OverlapProblemTPAnalyzer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 103 of file OverlapProblemTPAnalyzer.cc.

References m_assosimhitytecr, m_etatp, m_llbit, m_nassotk, m_nhits, m_nrechits, m_nrecvssimhits, m_pdgid, m_ptp, m_simhitytecr, m_statustp, TFileService::make(), mergeVDriftHistosByStation::name, relativeConstraints::ring, and indexGen::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
 
 
 
   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 ( )

Definition at line 147 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

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 162 of file OverlapProblemTPAnalyzer.cc.

References funct::abs(), edm::Event::getByToken(), i, cmsHarvester::index, 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, and HI_PhotonSkim_cff::rechits.

 {
   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();
     
     // commented since with the new TP's I don't know how to loop on PSimHits
 
     /*
     for( std::vector<PSimHit>::const_iterator sh = tksimhits.begin(); sh!= tksimhits.end(); ++sh) {
       
       // check if the SimHit is Tracker and TEC
       
       LogTrace("SimHitDetId") << sh->detUnitId();
       
       TECDetId det(sh->detUnitId());
       if(det.subDetector() == SiStripDetId::TEC) {
     
     unsigned int iring = det.ring() - 1;
     m_simhitytecr[iring]->Fill(sh->entryPoint().y());
     
     // check if there is a TrackingRecHit in the same detid and if this is the case fill the histos
     
     for(std::vector<DetId>::const_iterator rhdet = rechits.begin(); rhdet!= rechits.end() ; ++rhdet) {
       if(det==*rhdet) {
         m_assosimhitytecr[iring]->Fill(sh->entryPoint().y());
         break;
       }
     }
     
       }
     }
     */    
   }
   
 }