#include <DQMOffline/CalibTracker/plugins/SiPixelCalSingleMuonAnalyzer.cc>

Inheritance diagram for SiPixelCalSingleMuonAnalyzer:

Public Member Functions
	SiPixelCalSingleMuonAnalyzer (const edm::ParameterSet &)

	~SiPixelCalSingleMuonAnalyzer () override=default

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

	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 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 Types
typedef dqm::reco::MonitorElement	MonitorElement

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

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	countClusters (const edm::Handle< SiPixelClusterCollectionNew > &handle, unsigned int &nClusGlobal)

const bool	detidIsOnPixel (const DetId &detid)

std::pair< float, float >	findClosestCluster (const edm::Handle< SiPixelClusterCollectionNew > &handle, const PixelClusterParameterEstimator pixelCPE_, const TrackerGeometry trackerGeometry_, uint32_t rawId, float traj_lx, float traj_ly)

TrajectoryStateOnSurface	getTrajectoryStateOnSurface (const TrajectoryMeasurement &measurement)

Private Attributes
const edm::EDGetTokenT< SiPixelClusterCollectionNew >	clustersToken_

const edm::EDGetTokenT< edm::ValueMap< std::vector< float > > >	distanceToken_

const std::string	dqm_path_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomEsToken_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomEsTokenBR_

MonitorElement *	h_distClosestInactive

MonitorElement *	h_distClosestMissing

MonitorElement *	h_distClosestTrack

MonitorElement *	h_distClosestValid

MonitorElement *	h_nALCARECOClusters

MonitorElement *	h_nCloseByClusters

const edm::EDGetTokenT< edm::View< reco::Track > >	muonTracksToken_

const edm::EDGetTokenT< SiPixelClusterCollectionNew >	nearByClustersToken_

bool	phase_

const edm::ESGetToken< PixelClusterParameterEstimator, TkPixelCPERecord >	pixelCPEEsToken_

SiPixelDetInfoFileReader	reader_

const edm::ESGetToken< TrackerTopology, TrackerTopologyRcd >	topoEsTokenBR_

const edm::EDGetTokenT< TrajTrackAssociationCollection >	trajTrackCollectionToken_

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 DQMEDAnalyzer
uint64_t	meId () const

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Description: Analysis of the closebyPixelClusters collections for Pixel Hit Efficiency mearurements purposes

Definition at line 50 of file SiPixelCalSingleMuonAnalyzer.cc.

Member Typedef Documentation

◆ MonitorElement

typedef dqm::reco::MonitorElement SiPixelCalSingleMuonAnalyzer::MonitorElement

private

Definition at line 85 of file SiPixelCalSingleMuonAnalyzer.cc.

Constructor & Destructor Documentation

◆ SiPixelCalSingleMuonAnalyzer()

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

explicit

Definition at line 98 of file SiPixelCalSingleMuonAnalyzer.cc.

     : geomEsToken_(esConsumes()),
       pixelCPEEsToken_(esConsumes(edm::ESInputTag("", "PixelCPEGeneric"))),
       geomEsTokenBR_(esConsumes<edm::Transition::BeginRun>()),
       topoEsTokenBR_(esConsumes<edm::Transition::BeginRun>()),
       clustersToken_(consumes<SiPixelClusterCollectionNew>(iConfig.getParameter<edm::InputTag>("clusterCollection"))),
       nearByClustersToken_(
           consumes<SiPixelClusterCollectionNew>(iConfig.getParameter<edm::InputTag>("nearByClusterCollection"))),
       trajTrackCollectionToken_(
           consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("trajectoryInput"))),
       distanceToken_(consumes<edm::ValueMap<std::vector<float>>>(iConfig.getParameter<edm::InputTag>("distToTrack"))),
       muonTracksToken_(consumes<edm::View<reco::Track>>(iConfig.getParameter<edm::InputTag>("muonTracks"))),
       dqm_path_(iConfig.getParameter<std::string>("dqmPath")),
       reader_(edm::FileInPath(iConfig.getParameter<std::string>("skimmedGeometryPath")).fullPath()) {}

◆ ~SiPixelCalSingleMuonAnalyzer()

SiPixelCalSingleMuonAnalyzer::~SiPixelCalSingleMuonAnalyzer ( )

overridedefault

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 118 of file SiPixelCalSingleMuonAnalyzer.cc.

References clustersToken_, countClusters(), detidIsOnPixel(), HLT_2024v14_cff::distance, distanceToken_, dqm::impl::MonitorElement::Fill(), findClosestCluster(), geomEsToken_, edm::EventSetup::getData(), getTrajectoryStateOnSurface(), h_distClosestInactive, h_distClosestMissing, h_distClosestTrack, h_distClosestValid, h_nALCARECOClusters, h_nCloseByClusters, iEvent, TrackingRecHit::inactive, TrajectoryStateOnSurface::isValid(), edm::Ref< C, T, F >::key(), TrajectoryStateOnSurface::localPosition(), LogDebug, LogTrace, TrackingRecHit::missing, DiMuonV_cfg::muonTracks, muonTracksToken_, nearByClustersToken_, custom_jme_cff::nMuons, pixelCPEEsToken_, DetId::rawId(), rpcPointValidation_cfi::recHit, mathSSE::sqrt(), HLT_2024v14_cff::track, trajTrackCollectionToken_, TrackingRecHit::valid, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

                                                                                               {
   using namespace edm;
 
   // get the Tracker geometry from event setup
   const TrackerGeometry* trackerGeometry_ = &iSetup.getData(geomEsToken_);
 
   // get the Pixel CPE from event setup
   const PixelClusterParameterEstimator* pixelCPE_ = &iSetup.getData(pixelCPEEsToken_);
 
   // get the muon track collection
   const auto& muonTrackCollectionHandle = iEvent.getHandle(muonTracksToken_);
   auto const& muonTracks = *muonTrackCollectionHandle;
 
   // get the track distances
   const auto& distancesToTrack = iEvent.getHandle(distanceToken_);
 
   unsigned int nMuons = muonTracks.size();
   for (unsigned int ij = 0; ij < nMuons; ij++) {
     auto muon = muonTrackCollectionHandle->ptrAt(ij);
     edm::RefToBase<reco::Track> trackRef = muonTrackCollectionHandle->refAt(ij);
     const auto& distances = (*distancesToTrack)[trackRef];
 
     LogDebug("SiPixelCalSingleMuonAnalyzer") << "distances size: " << distances.size() << std::endl;
 
     unsigned counter = 0;
     double closestDR = 999.;
     for (const auto& distance : distances) {
       counter++;
       LogDebug("SiPixelCalSingleMuonAnalyzer")
           << "track: " << counter << " distance:" << std::sqrt(distance) << std::endl;
       if (distance < closestDR && distance > 0) {
         closestDR = distance;
       }
     }
 
     h_distClosestTrack->Fill(std::sqrt(closestDR));
   }
 
   // Get cluster collection
   const auto& clusterCollectionHandle = iEvent.getHandle(clustersToken_);
 
   unsigned int nClusGlobal = 0;
   countClusters(clusterCollectionHandle, nClusGlobal);
 
   h_nALCARECOClusters->Fill(nClusGlobal);
   LogTrace("SiPixelCalSingleMuonAnalyzer") << "total ALCARECO clusters: " << nClusGlobal << std::endl;
 
   // Get nearby cluster collection
   const auto& nearByClusterCollectionHandle = iEvent.getHandle(nearByClustersToken_);
 
   unsigned int nNearByClusGlobal = 0;
   countClusters(nearByClusterCollectionHandle, nNearByClusGlobal);
 
   h_nCloseByClusters->Fill(nNearByClusGlobal);
   LogTrace("SiPixelCalSingleMuonAnalyzer") << "total close-by clusters: " << nNearByClusGlobal << std::endl;
 
   // Get Traj-Track Collection
   const auto& trajTrackCollectionHandle = iEvent.getHandle(trajTrackCollectionToken_);
 
   if (!trajTrackCollectionHandle.isValid())
     return;
 
   for (const auto& pair : *trajTrackCollectionHandle) {
     const edm::Ref<std::vector<Trajectory>> traj = pair.key;
     const reco::TrackRef track = pair.val;
 
     for (const TrajectoryMeasurement& measurement : traj->measurements()) {
       if (!measurement.updatedState().isValid())
         return;
 
       const TransientTrackingRecHit::ConstRecHitPointer& recHit = measurement.recHit();
 
       // Only looking for pixel hits
       DetId r_rawId = recHit->geographicalId();
 
       if (!this->detidIsOnPixel(r_rawId))
         continue;
 
       // Skipping hits with undeterminable positions
       TrajectoryStateOnSurface trajStateOnSurface = this->getTrajectoryStateOnSurface(measurement);
 
       if (!(trajStateOnSurface.isValid()))
         continue;
 
       // Position measurements
       // Looking for valid and missing hits
       LocalPoint localPosition = trajStateOnSurface.localPosition();
 
       const auto& traj_lx = localPosition.x();
       const auto& traj_ly = localPosition.y();
 
       const auto loc = this->findClosestCluster(
           nearByClusterCollectionHandle, pixelCPE_, trackerGeometry_, r_rawId.rawId(), traj_lx, traj_ly);
 
       float dist = (loc.first != -999.) ? std::sqrt(loc.first * loc.first + loc.second * loc.second) : -0.1;
 
       if (recHit->getType() == TrackingRecHit::valid) {
         LogTrace("SiPixelCalSingleMuonAnalyzer")
             << "RawID:" << r_rawId.rawId() << " (valid hit), distance: " << dist << std::endl;
         h_distClosestValid->Fill(dist);
       }
 
       if (recHit->getType() == TrackingRecHit::missing) {
         LogTrace("SiPixelCalSingleMuonAnalyzer")
             << "RawID:" << r_rawId.rawId() << " (missing hit), distance: " << dist << std::endl;
         h_distClosestMissing->Fill(dist);
       }
 
       if (recHit->getType() == TrackingRecHit::inactive) {
         LogTrace("SiPixelCalSingleMuonAnalyzer")
             << "RawID:" << r_rawId.rawId() << " (inactive hit), distance: " << dist << std::endl;
         h_distClosestInactive->Fill(dist);
       }
     }
   }
 }

◆ bookHistograms()

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

overrideprivatevirtual

Implements DQMEDAnalyzer.

Definition at line 236 of file SiPixelCalSingleMuonAnalyzer.cc.

References dqm::implementation::IBooker::book1D(), dqm::implementation::IBooker::book1I(), dqm_path_, h_distClosestInactive, h_distClosestMissing, h_distClosestTrack, h_distClosestValid, h_nALCARECOClusters, h_nCloseByClusters, and dqm::implementation::NavigatorBase::setCurrentFolder().

                                                                                                                {
   // book the overall event count and event types histograms
   iBooker.setCurrentFolder(dqm_path_ + "/ClusterCounts");
   h_nALCARECOClusters = iBooker.book1I(
       "h_nALCARECOClusters", "Number of Pixel clusters per event (ALCARECO) ;N_{clusters};events", 20, 0, 20);
   h_nCloseByClusters = iBooker.book1I(
       "h_nCloseByClusters", "Number of Pixel clusters per event (close-by) ;N_{clusters};events", 20, 0, 20);
 
   iBooker.setCurrentFolder(dqm_path_ + "/TrajDistance");
   h_distClosestValid =
       iBooker.book1D("h_distClosestValid",
                      "Distance of Closest cluster to trajectory (valid);distance (cm); valid trajectory measurements",
                      110,
                      -0.105,
                      0.995);
   h_distClosestMissing = iBooker.book1D(
       "h_distClosestMissing",
       "Distance of Closest cluster to trajectory (missing);distance (cm);missing trajectory measurements",
       110,
       -0.105,
       0.995);
   h_distClosestInactive = iBooker.book1D(
       "h_distClosestInactive",
       "Distance of Closest cluster to trajectory (inactive);distance (cm);inactive trajectory measurements",
       110,
       -0.105,
       0.995);
 
   iBooker.setCurrentFolder(dqm_path_ + "/OtherTrackDistance");
   h_distClosestTrack =
       iBooker.book1D("h_distClosestTrack",
                      "#DeltaR Distance of Closest track to the muon trajectory;#DeltaR distance; muon trajectories",
                      100,
                      0.,
                      5.);
 }

◆ countClusters()

void SiPixelCalSingleMuonAnalyzer::countClusters	(	const edm::Handle< SiPixelClusterCollectionNew > &	handle,
		unsigned int &	nClusGlobal
	)

private

Definition at line 312 of file SiPixelCalSingleMuonAnalyzer.cc.

References hcalRecHitTable_cff::detId, patZpeak::handle, and LogDebug.

Referenced by analyze().

 {
   for (const auto& DSVItr : *handle) {
     uint32_t rawid(DSVItr.detId());
     DetId detId(rawid);
     LogDebug("SiPixelCalSingleMuonAnalyzer") << "DetId: " << detId.rawId() << " size: " << DSVItr.size() << std::endl;
     nClusGlobal += DSVItr.size();
   }
 }

◆ detidIsOnPixel()

const bool SiPixelCalSingleMuonAnalyzer::detidIsOnPixel ( const DetId & detid )

private

Definition at line 274 of file SiPixelCalSingleMuonAnalyzer.cc.

References DetId::det(), PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, DetId::subdetId(), and DetId::Tracker.

Referenced by analyze().

 {
   if (detid.det() != DetId::Tracker)
     return false;
   if (detid.subdetId() == PixelSubdetector::PixelBarrel)
     return true;
   if (detid.subdetId() == PixelSubdetector::PixelEndcap)
     return true;
   return false;
 }

◆ fillDescriptions()

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

static

Definition at line 373 of file SiPixelCalSingleMuonAnalyzer.cc.

References edm::ConfigurationDescriptions::addWithDefaultLabel(), submitPVResolutionJobs::desc, ProducerED_cfi::InputTag, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                                               {
   edm::ParameterSetDescription desc;
   desc.setComment("Analysis of the closebyPixelClusters collections");
   desc.add<edm::InputTag>("clusterCollection", edm::InputTag("ALCARECOSiPixelCalSingleMuonTight"));
   desc.add<edm::InputTag>("nearByClusterCollection", edm::InputTag("closebyPixelClusters"));
   desc.add<edm::InputTag>("trajectoryInput", edm::InputTag("refittedTracks"));
   desc.add<edm::InputTag>("muonTracks", edm::InputTag("ALCARECOSiPixelCalSingleMuonTight"));
   desc.add<edm::InputTag>("distToTrack", edm::InputTag("trackDistances"));
   desc.add<std::string>("dqmPath", "SiPixelCalSingleMuonTight");
   desc.add<std::string>("skimmedGeometryPath",
                         "SLHCUpgradeSimulations/Geometry/data/PhaseI/PixelSkimmedGeometry_phase1.txt");
   descriptions.addWithDefaultLabel(desc);
 }

◆ findClosestCluster()

std::pair< float, float > SiPixelCalSingleMuonAnalyzer::findClosestCluster	(	const edm::Handle< SiPixelClusterCollectionNew > &	handle,
		const PixelClusterParameterEstimator *	pixelCPE_,
		const TrackerGeometry *	trackerGeometry_,
		uint32_t	rawId,
		float	traj_lx,
		float	traj_ly
	)

private

Definition at line 326 of file SiPixelCalSingleMuonAnalyzer.cc.

References ALCARECOSiPixelCalSingleMuonDQM_cff::clusterCollection, hcalRecHitTable_cff::detId, edmNew::DetSet< T >::end(), PixelClusterParameterEstimator::getParameters(), patZpeak::handle, TrackerGeometry::idToDetUnit(), submitPVValidationJobs::params, PixelSubdetector::PixelBarrel, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, GeomDetEnumerators::subDetId, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by analyze().

 {
   const SiPixelClusterCollectionNew& clusterCollection = *handle;
   SiPixelClusterCollectionNew::const_iterator itClusterSet = clusterCollection.begin();
 
   float minD = 10e7;
 
   auto loc = std::make_pair(-999., -999.);
 
   for (; itClusterSet != clusterCollection.end(); itClusterSet++) {
     DetId detId(itClusterSet->id());
     if (detId.rawId() != rawId)
       continue;
 
     unsigned int subDetId = detId.subdetId();
     if (subDetId != PixelSubdetector::PixelBarrel && subDetId != PixelSubdetector::PixelEndcap) {
       edm::LogError("NearByPixelClustersAnalyzer")
           << "ERROR: not a pixel cluster!!!" << std::endl;  // should not happen
       continue;
     }
 
     const PixelGeomDetUnit* pixdet = (const PixelGeomDetUnit*)trackerGeometry_->idToDetUnit(detId);
     edmNew::DetSet<SiPixelCluster>::const_iterator itCluster = itClusterSet->begin();
     for (; itCluster != itClusterSet->end(); ++itCluster) {
       LocalPoint lp(itCluster->x(), itCluster->y(), 0.);
       PixelClusterParameterEstimator::ReturnType params = pixelCPE_->getParameters(*itCluster, *pixdet);
       lp = std::get<0>(params);
 
       float D = (lp.x() - traj_lx) * (lp.x() - traj_lx) + (lp.y() - traj_ly) * (lp.y() - traj_ly);
       if (D < minD) {
         minD = D;
         loc.first = (lp.x() - traj_lx);
         loc.second = (lp.y() - traj_ly);
       }
     }  // loop on cluster sets
   }
   return loc;
 }

◆ getTrajectoryStateOnSurface()

TrajectoryStateOnSurface SiPixelCalSingleMuonAnalyzer::getTrajectoryStateOnSurface ( const TrajectoryMeasurement & measurement )

private

Definition at line 287 of file SiPixelCalSingleMuonAnalyzer.cc.

References TrajectoryMeasurement::backwardPredictedState(), and TrajectoryMeasurement::forwardPredictedState().

Referenced by analyze().

 {
   const static TrajectoryStateCombiner trajStateCombiner;
 
   const auto& forwardPredictedState = measurement.forwardPredictedState();
   const auto& backwardPredictedState = measurement.backwardPredictedState();
 
   if (forwardPredictedState.isValid() && backwardPredictedState.isValid())
     return trajStateCombiner(forwardPredictedState, backwardPredictedState);
 
   else if (backwardPredictedState.isValid())
     return backwardPredictedState;
 
   else if (forwardPredictedState.isValid())
     return forwardPredictedState;
 
   edm::LogError("NearbyPixelClusterProducer") << "Error saving traj. measurement data."
                                               << " Trajectory state on surface cannot be determined." << std::endl;
 
   return TrajectoryStateOnSurface();
 }