SiPixelPhase1HitsV Class Reference

#include <SiPixelPhase1HitsV.h>

Inheritance diagram for SiPixelPhase1HitsV:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override
	SiPixelPhase1HitsV (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	{   TOF_R, ELOSS, ENTRY_EXIT_X, ENTRY_EXIT_Y,   ENTRY_EXIT_Z, LOCAL_X, LOCAL_Y, LOCAL_Z,   LOCAL_PHI, LOCAL_ETA, EFFICIENCY_TRACK }

Private Attributes
edm::EDGetTokenT< edm::PSimHitContainer >	pixelBarrelHighToken_
edm::EDGetTokenT< edm::PSimHitContainer >	pixelBarrelLowToken_
edm::EDGetTokenT< edm::PSimHitContainer >	pixelForwardHighToken_
edm::EDGetTokenT< edm::PSimHitContainer >	pixelForwardLowToken_
edm::EDGetTokenT< TrackingParticleCollection >	tpToken_
edm::EDGetTokenT< reco::TrackToTrackingParticleAssociator >	trackAssociatorByHitsToken_
edm::EDGetTokenT< edm::View< reco::Track > >	tracksToken_

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 SiPixelPhase1HitsV.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
TOF_R
ELOSS
ENTRY_EXIT_X
ENTRY_EXIT_Y
ENTRY_EXIT_Z
LOCAL_X
LOCAL_Y
LOCAL_Z
LOCAL_PHI
LOCAL_ETA
EFFICIENCY_TRACK

Definition at line 22 of file SiPixelPhase1HitsV.h.

       {
    TOF_R,
    ELOSS,
    ENTRY_EXIT_X,
    ENTRY_EXIT_Y,
    ENTRY_EXIT_Z,
    LOCAL_X,
    LOCAL_Y,
    LOCAL_Z,
    LOCAL_PHI,
    LOCAL_ETA,
    EFFICIENCY_TRACK,
  };

Constructor & Destructor Documentation

SiPixelPhase1HitsV()

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

explicit

Definition at line 26 of file SiPixelPhase1HitsV.cc.

    : SiPixelPhase1Base(iConfig),
      pixelBarrelLowToken_(consumes<edm::PSimHitContainer>(iConfig.getParameter<edm::InputTag>("pixBarrelLowSrc"))),
      pixelBarrelHighToken_(consumes<edm::PSimHitContainer>(iConfig.getParameter<edm::InputTag>("pixBarrelHighSrc"))),
      pixelForwardLowToken_(consumes<edm::PSimHitContainer>(iConfig.getParameter<edm::InputTag>("pixForwardLowSrc"))),
      pixelForwardHighToken_(consumes<edm::PSimHitContainer>(iConfig.getParameter<edm::InputTag>("pixForwardHighSrc"))),
 
      tracksToken_(consumes<edm::View<reco::Track>>(iConfig.getParameter<edm::InputTag>("tracksTag"))),
      tpToken_(consumes<TrackingParticleCollection>(iConfig.getParameter<edm::InputTag>("tpTag"))),
      trackAssociatorByHitsToken_(consumes<reco::TrackToTrackingParticleAssociator>(
          iConfig.getParameter<edm::InputTag>("trackAssociatorByHitsTag"))) {}

Member Function Documentation

analyze()

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

overridevirtual

Implements SiPixelPhase1Base.

Definition at line 38 of file SiPixelPhase1HitsV.cc.

                                                                                    {
  edm::Handle<edm::PSimHitContainer> barrelLowInput;
  iEvent.getByToken(pixelBarrelLowToken_, barrelLowInput);
  if (!barrelLowInput.isValid())
    return;
 
  edm::Handle<edm::PSimHitContainer> barrelHighInput;
  iEvent.getByToken(pixelBarrelHighToken_, barrelHighInput);
  if (!barrelHighInput.isValid())
    return;
 
  edm::Handle<edm::PSimHitContainer> forwardLowInput;
  iEvent.getByToken(pixelForwardLowToken_, forwardLowInput);
  if (!forwardLowInput.isValid())
    return;
 
  edm::Handle<edm::PSimHitContainer> forwardHighInput;
  iEvent.getByToken(pixelForwardHighToken_, forwardHighInput);
  if (!forwardHighInput.isValid())
    return;
 
  edm::PSimHitContainer::const_iterator it;
 
  // Get geometry information
 
  edm::ESHandle<TrackerGeometry> tracker;
  iSetup.get<TrackerDigiGeometryRecord>().get(tracker);
 
  // get low barrel info
  for (it = barrelLowInput->begin(); it != barrelLowInput->end(); ++it) {
    auto id = DetId(it->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(id);
    GlobalPoint gpos = det->toGlobal(it->localPosition());
 
    float tof = it->timeOfFlight();
    float globalR = gpos.mag();
 
    float energyLoss = it->energyLoss();
 
    float entryExitX = (it->entryPoint().x() - it->exitPoint().x());
    float entryExitY = (it->entryPoint().y() - it->exitPoint().y());
    float entryExitZ = std::abs(it->entryPoint().z() - it->exitPoint().z());
 
    float localX = it->localPosition().x();
    float localY = it->localPosition().y();
    float localZ = it->localPosition().z();
    float localPhi = it->localPosition().phi();
    float localEta = it->localPosition().eta();
 
    histo[TOF_R].fill(globalR, tof, id, &iEvent);
    histo[ELOSS].fill(energyLoss, id, &iEvent);
    histo[ENTRY_EXIT_X].fill(entryExitX, id, &iEvent);
    histo[ENTRY_EXIT_Y].fill(entryExitY, id, &iEvent);
    histo[ENTRY_EXIT_Z].fill(entryExitZ, id, &iEvent);
    histo[LOCAL_X].fill(localX, id, &iEvent);
    histo[LOCAL_Y].fill(localY, id, &iEvent);
    histo[LOCAL_Z].fill(localZ, id, &iEvent);
    histo[LOCAL_PHI].fill(localPhi, id, &iEvent);
    histo[LOCAL_ETA].fill(localEta, id, &iEvent);
  }
  // get high barrel info
  for (it = barrelHighInput->begin(); it != barrelHighInput->end(); ++it) {
    auto id = DetId(it->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(id);
    GlobalPoint gpos = det->toGlobal(it->localPosition());
 
    float tof = it->timeOfFlight();
    float globalR = gpos.mag();
 
    float energyLoss = it->energyLoss();
 
    float entryExitX = (it->entryPoint().x() - it->exitPoint().x());
    float entryExitY = (it->entryPoint().y() - it->exitPoint().y());
    float entryExitZ = std::abs(it->entryPoint().z() - it->exitPoint().z());
 
    float localX = it->localPosition().x();
    float localY = it->localPosition().y();
    float localZ = it->localPosition().z();
    float localPhi = it->localPosition().phi();
    float localEta = it->localPosition().eta();
 
    histo[TOF_R].fill(globalR, tof, id, &iEvent);
    histo[ELOSS].fill(energyLoss, id, &iEvent);
    histo[ENTRY_EXIT_X].fill(entryExitX, id, &iEvent);
    histo[ENTRY_EXIT_Y].fill(entryExitY, id, &iEvent);
    histo[ENTRY_EXIT_Z].fill(entryExitZ, id, &iEvent);
    histo[LOCAL_X].fill(localX, id, &iEvent);
    histo[LOCAL_Y].fill(localY, id, &iEvent);
    histo[LOCAL_Z].fill(localZ, id, &iEvent);
    histo[LOCAL_PHI].fill(localPhi, id, &iEvent);
    histo[LOCAL_ETA].fill(localEta, id, &iEvent);
  }
 
  // get low forward info
  for (it = forwardLowInput->begin(); it != forwardLowInput->end(); ++it) {
    auto id = DetId(it->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(id);
    GlobalPoint gpos = det->toGlobal(it->localPosition());
 
    float tof = it->timeOfFlight();
    float globalR = gpos.mag();
 
    float energyLoss = it->energyLoss();
 
    float entryExitX = (it->entryPoint().x() - it->exitPoint().x());
    float entryExitY = (it->entryPoint().y() - it->exitPoint().y());
    float entryExitZ = std::abs(it->entryPoint().z() - it->exitPoint().z());
 
    float localX = it->localPosition().x();
    float localY = it->localPosition().y();
    float localZ = it->localPosition().z();
    float localPhi = it->localPosition().phi();
    float localEta = it->localPosition().eta();
 
    histo[TOF_R].fill(globalR, tof, id, &iEvent);
    histo[ELOSS].fill(energyLoss, id, &iEvent);
    histo[ENTRY_EXIT_X].fill(entryExitX, id, &iEvent);
    histo[ENTRY_EXIT_Y].fill(entryExitY, id, &iEvent);
    histo[ENTRY_EXIT_Z].fill(entryExitZ, id, &iEvent);
    histo[LOCAL_X].fill(localX, id, &iEvent);
    histo[LOCAL_Y].fill(localY, id, &iEvent);
    histo[LOCAL_Z].fill(localZ, id, &iEvent);
    histo[LOCAL_PHI].fill(localPhi, id, &iEvent);
    histo[LOCAL_ETA].fill(localEta, id, &iEvent);
  }
 
  // get high forward info
  for (it = forwardHighInput->begin(); it != forwardHighInput->end(); ++it) {
    auto id = DetId(it->detUnitId());
    const GeomDetUnit *det = (const GeomDetUnit *)tracker->idToDetUnit(id);
    GlobalPoint gpos = det->toGlobal(it->localPosition());
 
    float tof = it->timeOfFlight();
    float globalR = gpos.mag();
 
    float energyLoss = it->energyLoss();
 
    float entryExitX = (it->entryPoint().x() - it->exitPoint().x());
    float entryExitY = (it->entryPoint().y() - it->exitPoint().y());
    float entryExitZ = std::abs(it->entryPoint().z() - it->exitPoint().z());
 
    float localX = it->localPosition().x();
    float localY = it->localPosition().y();
    float localZ = it->localPosition().z();
    float localPhi = it->localPosition().phi();
    float localEta = it->localPosition().eta();
 
    histo[TOF_R].fill(globalR, tof, id, &iEvent);
    histo[ELOSS].fill(energyLoss, id, &iEvent);
    histo[ENTRY_EXIT_X].fill(entryExitX, id, &iEvent);
    histo[ENTRY_EXIT_Y].fill(entryExitY, id, &iEvent);
    histo[ENTRY_EXIT_Z].fill(entryExitZ, id, &iEvent);
    histo[LOCAL_X].fill(localX, id, &iEvent);
    histo[LOCAL_Y].fill(localY, id, &iEvent);
    histo[LOCAL_Z].fill(localZ, id, &iEvent);
    histo[LOCAL_PHI].fill(localPhi, id, &iEvent);
    histo[LOCAL_ETA].fill(localEta, id, &iEvent);
  }
 
  // Sim Hit efficiency info
  edm::Handle<edm::View<reco::Track>> trackCollectionH;
  iEvent.getByToken(tracksToken_, trackCollectionH);
  const edm::View<reco::Track> &tC = *(trackCollectionH.product());
 
  edm::Handle<TrackingParticleCollection> TPCollectionH;
  iEvent.getByToken(tpToken_, TPCollectionH);
 
  edm::Handle<reco::TrackToTrackingParticleAssociator> theHitsAssociator;
  iEvent.getByToken(trackAssociatorByHitsToken_, theHitsAssociator);
  if (!theHitsAssociator.isValid()) {
    throw cms::Exception("NO VALID HIT ASSOCIATOR");
  }
  reco::TrackToTrackingParticleAssociator const *associatorByHits = theHitsAssociator.product();
 
  if (TPCollectionH.isValid() && trackCollectionH.isValid()) {
    reco::RecoToSimCollection const &p = associatorByHits->associateRecoToSim(trackCollectionH, TPCollectionH);
 
    for (edm::View<reco::Track>::size_type i = 0; i < tC.size(); ++i) {
      edm::RefToBase<reco::Track> track(trackCollectionH, i);
      //      const reco::Track& t = *track;
      auto id = DetId(track->innerDetId());  // histo manager requires a det ID,
                                             // use innermost ID for ease
 
      auto iter = p.find(track);
      histo[EFFICIENCY_TRACK].fill(iter != p.end() ? 1 : 0, id, &iEvent);
    }
  }
}

References funct::abs(), reco::TrackToTrackingParticleAssociator::associateRecoToSim(), EFFICIENCY_TRACK, ELOSS, fastSimProducer_cff::energyLoss, ENTRY_EXIT_X, ENTRY_EXIT_Y, ENTRY_EXIT_Z, Exception, edm::EventSetup::get(), get, HistogramManagerHolder::histo, mps_fire::i, triggerObjects_cff::id, iEvent, edm::HandleBase::isValid(), LOCAL_ETA, LOCAL_PHI, LOCAL_X, LOCAL_Y, LOCAL_Z, phase1PixelTopology::localX(), phase1PixelTopology::localY(), PV3DBase< T, PVType, FrameType >::mag(), AlCaHLTBitMon_ParallelJobs::p, pixelBarrelHighToken_, pixelBarrelLowToken_, pixelForwardHighToken_, pixelForwardLowToken_, edm::Handle< T >::product(), edm::View< T >::size(), TOF_R, tpToken_, HLT_2018_cff::track, trackAssociatorByHitsToken_, PbPb_ZMuSkimMuonDPG_cff::tracker, and tracksToken_.

Member Data Documentation

pixelBarrelHighToken_

edm::EDGetTokenT<edm::PSimHitContainer> SiPixelPhase1HitsV::pixelBarrelHighToken_

private

Definition at line 42 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

pixelBarrelLowToken_

edm::EDGetTokenT<edm::PSimHitContainer> SiPixelPhase1HitsV::pixelBarrelLowToken_

private

Definition at line 41 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

pixelForwardHighToken_

edm::EDGetTokenT<edm::PSimHitContainer> SiPixelPhase1HitsV::pixelForwardHighToken_

private

Definition at line 44 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

pixelForwardLowToken_

edm::EDGetTokenT<edm::PSimHitContainer> SiPixelPhase1HitsV::pixelForwardLowToken_

private

Definition at line 43 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

tpToken_

edm::EDGetTokenT<TrackingParticleCollection> SiPixelPhase1HitsV::tpToken_

private

Definition at line 47 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

trackAssociatorByHitsToken_

edm::EDGetTokenT<reco::TrackToTrackingParticleAssociator> SiPixelPhase1HitsV::trackAssociatorByHitsToken_

private

Definition at line 48 of file SiPixelPhase1HitsV.h.

Referenced by analyze().

tracksToken_

edm::EDGetTokenT<edm::View<reco::Track> > SiPixelPhase1HitsV::tracksToken_

private

Definition at line 46 of file SiPixelPhase1HitsV.h.

Referenced by analyze().