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	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () 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_FULL_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().