#include <Phase2ITValidateTrackingRecHit.cc>

Inheritance diagram for Phase2ITValidateTrackingRecHit:

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

	Phase2ITValidateTrackingRecHit (const edm::ParameterSet &)

	~Phase2ITValidateTrackingRecHit () override

Public Member Functions inherited from Phase2ITValidateRecHitBase
void	bookHistograms (DQMStore::IBooker &ibooker, edm::Run const &iRun, edm::EventSetup const &iSetup) override

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

	Phase2ITValidateRecHitBase (const edm::ParameterSet &)

	~Phase2ITValidateRecHitBase () 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

	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

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from Phase2ITValidateRecHitBase
static void	fillPSetDescription (edm::ParameterSetDescription &desc)

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 &)

Private Member Functions
void	fillITHistos (const edm::Event &iEvent, const TrackerHitAssociator &associateRecHit, const std::vector< edm::Handle< edm::PSimHitContainer >> &simHits, const std::map< unsigned int, SimTrack > &selectedSimTrackMap)

Private Attributes
edm::ParameterSet	config_

std::vector< edm::EDGetTokenT< edm::PSimHitContainer > >	simHitTokens_

const double	simtrackminpt_

const edm::EDGetTokenT< edm::SimTrackContainer >	simTracksToken_

const edm::EDGetTokenT< reco::TrackCollection >	tokenTracks_

TrackerHitAssociator::Config	trackerHitAssociatorConfig_

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 >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Protected Member Functions inherited from Phase2ITValidateRecHitBase
void	bookLayerHistos (DQMStore::IBooker &ibooker, unsigned int det_id, std::string &subdir)

void	fillRechitHistos (const PSimHit simhitClosest, const SiPixelRecHit rechit, const std::map< unsigned int, SimTrack > &selectedSimTrackMap, std::map< std::string, unsigned int > &nrechitLayerMap_primary)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from Phase2ITValidateRecHitBase
edm::ParameterSet	config_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomToken_

std::map< std::string, RecHitME >	layerMEs_

const TrackerGeometry *	tkGeom_ = nullptr

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoToken_

const TrackerTopology *	tTopo_ = nullptr

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: Plugin for Phase2 TrackingRecHit validation

Definition at line 52 of file Phase2ITValidateTrackingRecHit.cc.

Constructor & Destructor Documentation

◆ Phase2ITValidateTrackingRecHit()

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

explicit

Definition at line 72 of file Phase2ITValidateTrackingRecHit.cc.

References config_, edm::ParameterSet::getParameter(), and simHitTokens_.

     : Phase2ITValidateRecHitBase(iConfig),
       config_(iConfig),
       trackerHitAssociatorConfig_(iConfig, consumesCollector()),
       simtrackminpt_(iConfig.getParameter<double>("SimTrackMinPt")),
       tokenTracks_(consumes<reco::TrackCollection>(iConfig.getParameter<edm::InputTag>("tracksSrc"))),
       simTracksToken_(consumes<edm::SimTrackContainer>(iConfig.getParameter<edm::InputTag>("simTracksSrc"))) {
   edm::LogInfo("Phase2ITValidateRecHit") << ">>> Construct Phase2ITValidateRecHit ";
   for (const auto& itName : config_.getParameter<std::vector<std::string>>("ROUList")) {
     simHitTokens_.push_back(consumes<std::vector<PSimHit>>(edm::InputTag("g4SimHits", itName)));
   }
 }

◆ ~Phase2ITValidateTrackingRecHit()

Phase2ITValidateTrackingRecHit::~Phase2ITValidateTrackingRecHit ( )

override

Definition at line 85 of file Phase2ITValidateTrackingRecHit.cc.

                                                                 {
   edm::LogInfo("Phase2ITValidateTrackingRecHit") << ">>> Destroy Phase2ITValidateTrackingRecHit ";
 }

Member Function Documentation

◆ analyze()

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 89 of file Phase2ITValidateTrackingRecHit.cc.

References fillITHistos(), iEvent, FastTrackerRecHitCombiner_cfi::simHits, simHitTokens_, simtrackminpt_, TrackCandidateProducer_cfi::simTracks, simTracksToken_, and trackerHitAssociatorConfig_.

                                                                                                 {
   std::vector<edm::Handle<edm::PSimHitContainer>> simHits;
   for (const auto& itoken : simHitTokens_) {
     const auto& simHitHandle = iEvent.getHandle(itoken);
     if (!simHitHandle.isValid())
       continue;
     simHits.emplace_back(simHitHandle);
   }
   // Get the SimTracks and push them in a map of id, SimTrack
   const auto& simTracks = iEvent.getHandle(simTracksToken_);
 
   std::map<unsigned int, SimTrack> selectedSimTrackMap;
   for (const auto& simTrackIt : *simTracks) {
     if (simTrackIt.momentum().pt() > simtrackminpt_) {
       selectedSimTrackMap.emplace(simTrackIt.trackId(), simTrackIt);
     }
   }
   TrackerHitAssociator associateRecHit(iEvent, trackerHitAssociatorConfig_);
   fillITHistos(iEvent, associateRecHit, simHits, selectedSimTrackMap);
 }

◆ fillDescriptions()

void Phase2ITValidateTrackingRecHit::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 182 of file Phase2ITValidateTrackingRecHit.cc.

References edm::ConfigurationDescriptions::add(), submitPVResolutionJobs::desc, Phase2ITValidateRecHitBase::fillPSetDescription(), HLT_2022v12_cff::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                                 {
   // rechitValidIT
   edm::ParameterSetDescription desc;
 
   // call the base fillPsetDescription for the plots bookings
   Phase2ITValidateRecHitBase::fillPSetDescription(desc);
 
   //to be used in TrackerHitAssociator
   desc.add<bool>("associatePixel", true);
   desc.add<bool>("associateStrip", false);
   desc.add<bool>("usePhase2Tracker", true);
   desc.add<bool>("associateRecoTracks", false);
   desc.add<bool>("associateHitbySimTrack", true);
   desc.add<edm::InputTag>("pixelSimLinkSrc", edm::InputTag("simSiPixelDigis", "Pixel"));
   desc.add<std::vector<std::string>>("ROUList",
                                      {
                                          "TrackerHitsPixelBarrelLowTof",
                                          "TrackerHitsPixelBarrelHighTof",
                                          "TrackerHitsPixelEndcapLowTof",
                                          "TrackerHitsPixelEndcapHighTof",
                                      });
   //
   desc.add<edm::InputTag>("simTracksSrc", edm::InputTag("g4SimHits"));
   desc.add<edm::InputTag>("SimVertexSource", edm::InputTag("g4SimHits"));
   desc.add<double>("SimTrackMinPt", 2.0);
   desc.add<edm::InputTag>("tracksSrc", edm::InputTag("generalTracks"));
   desc.add<std::string>("TopFolderName", "TrackerPhase2ITTrackingRecHitV");
   desc.add<bool>("Verbosity", false);
   descriptions.add("Phase2ITValidateTrackingRecHit", desc);
 }

◆ fillITHistos()

void Phase2ITValidateTrackingRecHit::fillITHistos	(	const edm::Event &	iEvent,
		const TrackerHitAssociator &	associateRecHit,
		const std::vector< edm::Handle< edm::PSimHitContainer >> &	simHits,
		const std::map< unsigned int, SimTrack > &	selectedSimTrackMap
	)

private

Definition at line 110 of file Phase2ITValidateTrackingRecHit.cc.

References funct::abs(), TrackerHitAssociator::associateHitId(), Phase2ITValidateRecHitBase::fillRechitHistos(), phase2tkutil::getITHistoId(), TrackerGeometry::idToDetUnit(), iEvent, crabWrapper::key, Phase2ITValidateRecHitBase::layerMEs_, BaseTrackerRecHit::localPosition(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, FastTrackerRecHitCombiner_cfi::simHits, AlCaHLTBitMon_QueryRunRegistry::string, Phase2ITValidateRecHitBase::tkGeom_, tokenTracks_, HLT_2022v12_cff::track, tracks, Phase2ITValidateRecHitBase::tTopo_, and PV3DBase< T, PVType, FrameType >::x().

Referenced by analyze().

                                                                                                              {
   const auto& tracks = iEvent.getHandle(tokenTracks_);
   if (!tracks.isValid())
     return;
 
   std::map<std::string, unsigned int> nrechitLayerMap_primary;
 
   // loop over tracks
   for (const auto& track : *tracks) {
     // loop over hits
     for (auto const& hit : track.recHits()) {
       if (!hit->isValid())
         continue;
 
       auto id = hit->geographicalId();
       // check that we are in the pixel
       auto subdetid = (id.subdetId());
       if (!(subdetid == PixelSubdetector::PixelBarrel) && !(subdetid == PixelSubdetector::PixelEndcap))
         continue;
 
       const GeomDetUnit* geomDetunit(tkGeom_->idToDetUnit(id));
       if (!geomDetunit)
         continue;
       // determine the detector we are in
       std::string key = phase2tkutil::getITHistoId(id.rawId(), tTopo_);
       if (nrechitLayerMap_primary.find(key) == nrechitLayerMap_primary.end()) {
         nrechitLayerMap_primary.emplace(key, 1);
       } else {
         nrechitLayerMap_primary[key] += 1;
       }
 
       const SiPixelRecHit* rechit = dynamic_cast<const SiPixelRecHit*>(hit);
       if (!rechit)
         continue;
 
       const std::vector<SimHitIdpr>& matchedId = associateRecHit.associateHitId(*rechit);
       const PSimHit* simhitClosest = nullptr;
       float minx = 10000;
       LocalPoint lp = rechit->localPosition();
       for (const auto& simHitCol : simHits) {
         for (const auto& simhitIt : *simHitCol) {
           if (id.rawId() != simhitIt.detUnitId())
             continue;
           for (const auto& mId : matchedId) {
             if (simhitIt.trackId() == mId.first) {
               if (!simhitClosest || std::abs(simhitIt.localPosition().x() - lp.x()) < minx) {
                 minx = std::abs(simhitIt.localPosition().x() - lp.x());
                 simhitClosest = &simhitIt;
               }
             }
           }
         }  //end loop over PSimhitcontainers
       }    //end loop over simHits
 
       if (!simhitClosest)
         continue;
 
       // call the base class method to fill the plots
       fillRechitHistos(simhitClosest, rechit, selectedSimTrackMap, nrechitLayerMap_primary);
 
     }  // loop over tracking rechits
   }    // loop over tracks
 
   //fill nRecHit counter per layer
   for (const auto& lme : nrechitLayerMap_primary) {
     layerMEs_[lme.first].numberRecHitsprimary->Fill(nrechitLayerMap_primary[lme.first]);
   }
 }