#include <SiPixelPhase1TrackingParticleV.h>

Inheritance diagram for SiPixelPhase1TrackingParticleV:

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

	SiPixelPhase1TrackingParticleV (const edm::ParameterSet &conf)

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

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

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

virtual	~SiPixelPhase1Base ()

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

virtual void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer (void)

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

virtual 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

virtual	~EDAnalyzerBase ()

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

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	{ MASS, CHARGE, ID, NHITS, MATCHED, PT, PHI, ETA, VTX, VTY, VYZ, TIP, LIP }

Private Attributes
std::vector< edm::EDGetTokenT< std::vector< PSimHit > > >	simHitTokens_

std::vector< std::pair< unsigned int, const PSimHit * > >	trackIdToHitPtr_

edm::EDGetTokenT< TrackingParticleCollection >	vec_TrackingParticle_Token_

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 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 21 of file SiPixelPhase1TrackingParticleV.h.

Member Enumeration Documentation

anonymous enum

private

Enumerator
MASS
CHARGE
ID
NHITS
MATCHED
PT
PHI
ETA
VTX
VTY
VYZ
TIP
LIP

Definition at line 22 of file SiPixelPhase1TrackingParticleV.h.

        {
     MASS,
     CHARGE,
     ID,
     NHITS,
     MATCHED,
     PT,
     PHI,
     ETA,
     VTX,
     VTY,
     VYZ,
     TIP,
     LIP,
   };

Constructor & Destructor Documentation

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

explicit

Definition at line 40 of file SiPixelPhase1TrackingParticleV.cc.

References edm::EDConsumerBase::consumes(), edm::ParameterSet::getParameter(), simHitTokens_, and GlobalPosition_Frontier_DevDB_cff::tag.

                                                                                              :
   SiPixelPhase1Base(iConfig),
   vec_TrackingParticle_Token_( consumes<TrackingParticleCollection>( iConfig.getParameter<edm::InputTag>( "src" ) ) )
 {
   for(const auto& tag: iConfig.getParameter<std::vector<edm::InputTag>>("simHitToken")) {
     simHitTokens_.push_back(consumes<std::vector<PSimHit>>(tag));
   }
 }

Member Function Documentation

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

Definition at line 49 of file SiPixelPhase1TrackingParticleV.cc.

References begin, CHARGE, DEFINE_FWK_MODULE, PSimHit::detUnitId(), end, ETA, PSimHit::eventId(), edm::Event::getByToken(), HistogramManagerHolder::histo, hcalTTPDigis_cfi::id, ID, iEvent, LIP, MASS, MATCHED, NHITS, PHI, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, edm::Handle< T >::product(), PT, rpcPointValidation_cfi::simHit, SiPixelPhase1TrackingParticleV_cfi::simHitToken, simHitTokens_, simTrackMatching_cfi::simTrack, mathSSE::sqrt(), lumiQTWidget::t, TIP, trackIdToHitPtr_, vec_TrackingParticle_Token_, VTX, VTY, and VYZ.

                                                                                                 {
 
   edm::Handle<TrackingParticleCollection>  TruthTrackContainer;
   iEvent.getByToken( vec_TrackingParticle_Token_, TruthTrackContainer );
   const TrackingParticleCollection *tPC   = TruthTrackContainer.product();
 
   // A multimap linking SimTrack::trackId() to a pointer to PSimHit
   // Similar to TrackingTruthAccumulator
   for(const auto& simHitToken: simHitTokens_) {
     edm::Handle<std::vector<PSimHit> > hsimhits;
     iEvent.getByToken(simHitToken, hsimhits);
     trackIdToHitPtr_.reserve(trackIdToHitPtr_.size()+hsimhits->size());
     for(const auto& simHit: *hsimhits) {
       trackIdToHitPtr_.emplace_back(simHit.trackId(), &simHit);
     }
   }
   std::stable_sort(trackIdToHitPtr_.begin(), trackIdToHitPtr_.end(), trackIdHitPairLessSort);
 
 
   // Loop over TrackingParticle's
   for (TrackingParticleCollection::const_iterator t = tPC -> begin(); t != tPC -> end(); ++t) {
 
     // histo manager requires a det ID, use first tracker hit
 
     bool isBpixtrack = false, isFpixtrack = false;
     DetId id;
 
     for(const SimTrack& simTrack: t->g4Tracks()) {
       // Logic is from TrackingTruthAccumulator
       auto range = std::equal_range(trackIdToHitPtr_.begin(), trackIdToHitPtr_.end(), std::pair<unsigned int, const PSimHit *>(simTrack.trackId(), nullptr), trackIdHitPairLess);
       if(range.first == range.second) continue;
 
       auto iHitPtr = range.first;
       for(; iHitPtr != range.second; ++iHitPtr) {
         const PSimHit& simHit = *(iHitPtr->second);
         if(simHit.eventId() != t->eventId())
           continue;
         id = DetId( simHit.detUnitId() );
 
         // check we are in pixel
         uint32_t subdetid = (id.subdetId());
     if (subdetid == PixelSubdetector::PixelBarrel) isBpixtrack = true;
     if (subdetid == PixelSubdetector::PixelEndcap) isFpixtrack = true;
     if (subdetid != PixelSubdetector::PixelBarrel && subdetid != PixelSubdetector::PixelEndcap) continue;
       }
     }
 
     if ( isBpixtrack || isFpixtrack ) {
       histo[MASS].fill(t->mass(), id, &iEvent);
       histo[CHARGE].fill(t->charge(), id, &iEvent);
       histo[ID].fill(t->pdgId(), id, &iEvent);
       histo[NHITS].fill(t->numberOfTrackerHits(), id, &iEvent);
       histo[MATCHED].fill(t->numberOfTrackerLayers(), id, &iEvent);
       histo[PT].fill(sqrt(t->momentum().perp2()), id, &iEvent);
       histo[PHI].fill(t->momentum().Phi(), id, &iEvent);
       histo[ETA].fill(t->momentum().eta(), id, &iEvent);
       histo[VTX].fill(t->vx(), id, &iEvent);
       histo[VTY].fill(t->vy(), id, &iEvent);
       histo[VYZ].fill(t->vz(), id, &iEvent);
       histo[TIP].fill(sqrt(t->vertex().perp2()), id, &iEvent);
       histo[LIP].fill(t->vz(), id, &iEvent);
     }
   }
 }