#include <SiPixelPhase1TrackingParticleV.h>

Inheritance diagram for SiPixelPhase1TrackingParticleV:

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

	SiPixelPhase1TrackingParticleV (const edm::ParameterSet &conf)

Public Member Functions inherited from SiPixelPhase1Base
void	bookHistograms (DQMStore::IBooker &iBooker, edm::Run const &run, edm::EventSetup const &) override

	SiPixelPhase1Base (const edm::ParameterSet &iConfig)

	~SiPixelPhase1Base () override

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

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

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

Public Member Functions inherited from HistogramManagerHolder
	HistogramManagerHolder (const edm::ParameterSet &iConfig, edm::ConsumesCollector &&iC, edm::Transition transition=edm::Transition::BeginRun)

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_

edm::EDGetTokenT< TrackingParticleCollection >	vec_TrackingParticle_Token_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
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

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Types inherited from SiPixelPhase1Base
enum	{ DCS }

Protected Member Functions inherited from SiPixelPhase1Base
bool	checktrigger (const edm::Event &iEvent, const edm::EventSetup &, const unsigned trgidx) const

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

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

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::SiPixelPhase1TrackingParticleV ( const edm::ParameterSet & conf )

explicit

Definition at line 45 of file SiPixelPhase1TrackingParticleV.cc.

     : 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));
   }
 }

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

Member Function Documentation

◆ analyze()

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

overridevirtual

Implements SiPixelPhase1Base.

Definition at line 53 of file SiPixelPhase1TrackingParticleV.cc.

                                                                                                 {
   edm::Handle<TrackingParticleCollection> TruthTrackContainer;
   iEvent.getByToken(vec_TrackingParticle_Token_, TruthTrackContainer);
   const TrackingParticleCollection *tPC = TruthTrackContainer.product();
   std::vector<std::pair<unsigned int, const PSimHit *>> trackIdToHitPtr;
  
   // 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);
     }
   }
 }