#include <Calibration/TkAlCaRecoProducers/plugins/NearbyPixelClustersProducer.cc>

Inheritance diagram for NearbyPixelClustersProducer:

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

	~NearbyPixelClustersProducer () override=default

Public Member Functions inherited from edm::stream::EDProducer<>
	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)

Private Member Functions
bool	detidIsOnPixel (const DetId &detid)

const std::vector< edmNew::DetSet< SiPixelCluster >::const_iterator >	findAllNearbyClusters (const SiPixelClusterCollectionNew::const_iterator &clusterSet, const uint32_t rawId, const std::vector< LocalPoint > &vLocalPos)

const std::vector< edmNew::DetSet< SiPixelCluster >::const_iterator >	findAllNearbyClusters (const SiPixelClusterCollectionNew &clusterSet, const uint32_t rawId, const std::vector< LocalPoint > &vLocalPos)

const trajCrossings_t	findAllTrajectoriesCrossings (const edm::Handle< TrajTrackAssociationCollection > &trajTrackCollectionHandle)

TrajectoryStateOnSurface	getTrajectoryStateOnSurface (const TrajectoryMeasurement &measurement)

void	produce (edm::Event &, const edm::EventSetup &) override

Private Attributes
const edm::ESGetToken< SiPixelQuality, SiPixelQualityFromDbRcd >	badModuleToken_

edm::EDPutTokenT< SiPixelClusterCollectionNew >	clusterPutToken_

edm::EDGetTokenT< SiPixelClusterCollectionNew >	clustersToken_

const bool	dumpWholeDetId_

const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	geomEsToken_

const PixelClusterParameterEstimator *	pixelCPE_

const edm::ESGetToken< PixelClusterParameterEstimator, TkPixelCPERecord >	pixelCPEEsToken_

const bool	throwBadComponents_

const TrackerGeometry *	trackerGeometry_

edm::EDGetTokenT< TrajTrackAssociationCollection >	trajTrackCollectionToken_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
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

Detailed Description

Description: Class to produce the collection of SiPixelClusters closest, hit by hit, to the trajectory measurements of a given track

Implementation: Implementation of this class is heavily endebted to https://github.com/jkarancs/PhaseIPixelNtuplizer/blob/master/plugins/PhaseIPixelNtuplizer.h

Definition at line 63 of file NearbyPixelClustersProducer.cc.

Constructor & Destructor Documentation

◆ NearbyPixelClustersProducer()

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

explicit

Definition at line 112 of file NearbyPixelClustersProducer.cc.

References dumpWholeDetId_.

     : throwBadComponents_(iConfig.getParameter<bool>("throwBadComponents")),
       dumpWholeDetId_(iConfig.getParameter<bool>("dumpWholeDetIds")),
       geomEsToken_(esConsumes()),
       pixelCPEEsToken_(esConsumes(edm::ESInputTag("", "PixelCPEGeneric"))),
       badModuleToken_(esConsumes()),
       clustersToken_(consumes<SiPixelClusterCollectionNew>(iConfig.getParameter<edm::InputTag>("clusterCollection"))),
       trajTrackCollectionToken_(
           consumes<TrajTrackAssociationCollection>(iConfig.getParameter<edm::InputTag>("trajectoryInput"))),
       clusterPutToken_(produces<SiPixelClusterCollectionNew>()) {
   if (dumpWholeDetId_) {
     edm::LogPrint("NearbyPixelClustersProducer") << "WARNING: selected to dump the whole DetId's worth of clusters.\n "
                                                     "This will have consequences on the event size!"
                                                  << std::endl;
   }
 }

◆ ~NearbyPixelClustersProducer()

NearbyPixelClustersProducer::~NearbyPixelClustersProducer ( )

overridedefault

Member Function Documentation

◆ detidIsOnPixel()

bool NearbyPixelClustersProducer::detidIsOnPixel ( const DetId & detid )

private

Definition at line 383 of file NearbyPixelClustersProducer.cc.

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

Referenced by findAllTrajectoriesCrossings().

 {
   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 NearbyPixelClustersProducer::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 421 of file NearbyPixelClustersProducer.cc.

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

                                                                                              {
   edm::ParameterSetDescription desc;
   desc.setComment(
       "Produces the collection of SiPixelClusters closest, hit by hit, to the trajectory measurements of a given "
       "track");
   desc.add<bool>("throwBadComponents", false)
       ->setComment(
           "do not consider modules flagged as bad components. Careful, it changes the efficiency denominator!");
   desc.add<bool>("dumpWholeDetIds", false)
       ->setComment("put in the event all the pixel cluster on the impacted module, by default only the closest");
   ;
   desc.add<edm::InputTag>("clusterCollection", edm::InputTag("siPixelClusters"));
   desc.add<edm::InputTag>("trajectoryInput", edm::InputTag("myRefitter"));
   descriptions.addWithDefaultLabel(desc);
 }

◆ findAllNearbyClusters() [1/2]

const std::vector< edmNew::DetSet< SiPixelCluster >::const_iterator > NearbyPixelClustersProducer::findAllNearbyClusters	(	const SiPixelClusterCollectionNew::const_iterator &	clusterSet,
		const uint32_t	rawId,
		const std::vector< LocalPoint > &	vLocalPos
	)

private

Definition at line 247 of file NearbyPixelClustersProducer.cc.

References pv::closest(), ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), submitPVResolutionJobs::count, dumpWholeDetId_, edmNew::DetSet< T >::end(), PixelClusterParameterEstimator::getParameters(), TrackerGeometry::idToDetUnit(), LogDebug, slimmedTrackExtras_cff::outputClusters, submitPVValidationJobs::params, pixelCPE_, nano_mu_digi_cff::rawId, mathSSE::sqrt(), trackerGeometry_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

Referenced by produce().

 {
   std::vector<edmNew::DetSet<SiPixelCluster>::const_iterator> outputClusters;
 
   static constexpr unsigned int k_maxClustersInDet = 1024;
 
   // something funny is going on here ...
   if ((*clusterSet).size() > k_maxClustersInDet) {
     edm::LogWarning("NearbyPixelClustersProducer")
         << __func__ << "() number of clusters in det " << rawId /*(*clusterSet).id()*/ << " is " << (*clusterSet).size()
         << ", which is larger than maximum (1024).\n Something funny with the data is going on!" << std::endl;
     return outputClusters;
   }
 
   const PixelGeomDetUnit* pixdet = (const PixelGeomDetUnit*)trackerGeometry_->idToDetUnit(rawId);
   edmNew::DetSet<SiPixelCluster>::const_iterator itCluster = clusterSet->begin();
 
   // just copy straight the whole set of clusters in the detid
   if (dumpWholeDetId_) {
     for (; itCluster != clusterSet->end(); ++itCluster) {
       outputClusters.push_back(itCluster);
     }
     return outputClusters;
   }
 
   int count = 0;
   for (const auto& localPos : vLocalPos) {
     count++;
     //trajectory crossing local coordinates
     const auto& traj_lx = localPos.x();
     const auto& traj_ly = localPos.y();
 
     float minD = 10000.;
     edmNew::DetSet<SiPixelCluster>::const_iterator closest = nullptr;
 
     //std::cout<< rawId << " count: " << count << " n. clusters: " << (*clusterSet).size() << std::endl;
     LogDebug("NearbyPixelClustersProducer")
         << __func__ << rawId << " count: " << count << " n. clusters: " << (*clusterSet).size() << std::endl;
 
     for (; itCluster != clusterSet->end(); ++itCluster) {
       LocalPoint lp(itCluster->x(), itCluster->y(), 0.);
       const auto& params = pixelCPE_->getParameters(*itCluster, *pixdet);
       lp = std::get<0>(params);
 
       float D = sqrt((lp.x() - traj_lx) * (lp.x() - traj_lx) + (lp.y() - traj_ly) * (lp.y() - traj_ly));
       if (D < minD) {
         closest = itCluster;
         minD = D;
       }
     }  // loop on cluster sets
 
     if (closest) {
       outputClusters.push_back(closest);
     }
   }  // loop on the trajectory crossings
 
   return outputClusters;
 }

◆ findAllNearbyClusters() [2/2]

const std::vector< edmNew::DetSet< SiPixelCluster >::const_iterator > NearbyPixelClustersProducer::findAllNearbyClusters	(	const SiPixelClusterCollectionNew &	clusterSet,
		const uint32_t	rawId,
		const std::vector< LocalPoint > &	vLocalPos
	)

private

Definition at line 312 of file NearbyPixelClustersProducer.cc.

References pv::closest(), ALCARECOSiPixelCalSingleMuonDQM_cff::clusterCollection, submitPVResolutionJobs::count, hcalRecHitTable_cff::detId, dumpWholeDetId_, edmNew::DetSet< T >::end(), PixelClusterParameterEstimator::getParameters(), TrackerGeometry::idToDetUnit(), LogDebug, slimmedTrackExtras_cff::outputClusters, submitPVValidationJobs::params, PixelSubdetector::PixelBarrel, pixelCPE_, PixelSubdetector::PixelEndcap, nano_mu_digi_cff::rawId, mathSSE::sqrt(), GeomDetEnumerators::subDetId, trackerGeometry_, PV3DBase< T, PVType, FrameType >::x(), and PV3DBase< T, PVType, FrameType >::y().

 {
   std::vector<edmNew::DetSet<SiPixelCluster>::const_iterator> outputClusters;
 
   int count = 0;
   for (const auto& localPos : vLocalPos) {
     count++;
 
     //trajectory crossing local coordinates
     const auto& traj_lx = localPos.x();
     const auto& traj_ly = localPos.y();
 
     float minD = 10000.;
 
     SiPixelClusterCollectionNew::const_iterator itClusterSet = clusterCollection.begin();
     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("NearByPixelClusterProducer")
             << "ERROR: not a pixel cluster!!!" << std::endl;  // should not happen
         continue;
       }
 
       edmNew::DetSet<SiPixelCluster>::const_iterator itCluster = itClusterSet->begin();
 
       // just copy straight the whole set of clusters in the detid
       if (dumpWholeDetId_ && count == 1) {
         for (; itCluster != itClusterSet->end(); ++itCluster) {
           outputClusters.push_back(itCluster);
         }
         return outputClusters;
       }
 
       //std::cout<< rawId << " count: " << count << " n. clusters: " << (*clusterSet).size() << std::endl;
       LogDebug("NearbyPixelClustersProducer")
           << __func__ << rawId << " count: " << count << " n. clusters: " << (*itClusterSet).size() << std::endl;
 
       const PixelGeomDetUnit* pixdet = (const PixelGeomDetUnit*)trackerGeometry_->idToDetUnit(rawId);
 
       edmNew::DetSet<SiPixelCluster>::const_iterator closest = nullptr;
 
       for (; itCluster != itClusterSet->end(); ++itCluster) {
         LocalPoint lp(itCluster->x(), itCluster->y(), 0.);
         const auto& params = pixelCPE_->getParameters(*itCluster, *pixdet);
         lp = std::get<0>(params);
 
         float D = sqrt((lp.x() - traj_lx) * (lp.x() - traj_lx) + (lp.y() - traj_ly) * (lp.y() - traj_ly));
         if (D < minD) {
           closest = itCluster;
           minD = D;
         }
       }  // loop on cluster sets
 
       if (closest) {
         outputClusters.push_back(closest);
       }
     }  // loop on all clusters
   }    // loop on the trajectory crossings
 
   return outputClusters;
 }

◆ findAllTrajectoriesCrossings()

const trajCrossings_t NearbyPixelClustersProducer::findAllTrajectoriesCrossings ( const edm::Handle< TrajTrackAssociationCollection > & trajTrackCollectionHandle )

private

Definition at line 199 of file NearbyPixelClustersProducer.cc.

References detidIsOnPixel(), spr::find(), getTrajectoryStateOnSurface(), TrajectoryStateOnSurface::isValid(), edm::Ref< C, T, F >::key(), TrajectoryStateOnSurface::localPosition(), nano_mu_digi_cff::rawId, DetId::rawId(), and rpcPointValidation_cfi::recHit.

Referenced by produce().

 {
   trajCrossings_t crossings;
 
   std::vector<uint32_t> treatedIds;
 
   for (const auto& pair : *trajTrackCollectionHandle) {
     const edm::Ref<std::vector<Trajectory>> traj = pair.key;
 
     for (const TrajectoryMeasurement& measurement : traj->measurements()) {
       //Check if the measurement infos can be read
       if (!measurement.updatedState().isValid())
         continue;
 
       const TransientTrackingRecHit::ConstRecHitPointer& recHit = measurement.recHit();
 
       // Only looking for pixel hits
       DetId recHitDetid = recHit->geographicalId();
       const auto& rawId = recHitDetid.rawId();
 
       if (!this->detidIsOnPixel(recHitDetid))
         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();
 
       if (std::find(treatedIds.begin(), treatedIds.end(), rawId) != treatedIds.end()) {
         crossings.at(rawId).push_back(localPosition);
       } else {
         crossings.insert(std::pair<uint32_t, std::vector<LocalPoint>>(rawId, {localPosition}));
         treatedIds.push_back(rawId);
       }
     }  // loop on measurements in trajectory
   }    // loop on trajectories
 
   return crossings;
 }

◆ getTrajectoryStateOnSurface()

TrajectoryStateOnSurface NearbyPixelClustersProducer::getTrajectoryStateOnSurface ( const TrajectoryMeasurement & measurement )

private

Definition at line 396 of file NearbyPixelClustersProducer.cc.

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

Referenced by findAllTrajectoriesCrossings().

 {
   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("NearbyPixelClustersProducer") << "Error saving traj. measurement data."
                                                << " Trajectory state on surface cannot be determined." << std::endl;
 
   return TrajectoryStateOnSurface();
 }

◆ produce()

void NearbyPixelClustersProducer::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 134 of file NearbyPixelClustersProducer.cc.

References edmNew::DetSetVector< T >::FastFiller::abort(), badModuleToken_, ALCARECOSiPixelCalSingleMuonDQM_cff::clusterCollection, clusterPutToken_, clustersToken_, edmNew::DetSetVector< T >::FastFiller::empty(), findAllNearbyClusters(), findAllTrajectoriesCrossings(), geomEsToken_, edm::EventSetup::getData(), iEvent, LogDebug, eostools::move(), slimmedTrackExtras_cff::outputClusters, pixelCPE_, pixelCPEEsToken_, edmNew::DetSetVector< T >::FastFiller::push_back(), throwBadComponents_, trackerGeometry_, and trajTrackCollectionToken_.

                                                                                        {
   using namespace edm;
 
   auto outputClusters = std::make_unique<SiPixelClusterCollectionNew>();
 
   // get the Tracker geometry from event setup
   trackerGeometry_ = &iSetup.getData(geomEsToken_);
 
   // get the Pixel CPE from event setup
   pixelCPE_ = &iSetup.getData(pixelCPEEsToken_);
 
   const auto& SiPixelBadModule_ = &iSetup.getData(badModuleToken_);
 
   // get cluster collection
   const auto& clusterCollectionHandle = iEvent.getHandle(clustersToken_);
   const SiPixelClusterCollectionNew& clusterCollection = *clusterCollectionHandle;
 
   // get Traj-Track Collection
   const auto& trajTrackCollectionHandle = iEvent.getHandle(trajTrackCollectionToken_);
   if (!trajTrackCollectionHandle.isValid())
     return;
 
   // find all trajectory crossings in the event
   const auto& allCrossings = this->findAllTrajectoriesCrossings(trajTrackCollectionHandle);
 
   LogDebug("NearbyPixelClustersProducer") << allCrossings.size() << std::endl;
 
   // now find all nearby clusters
   for (const auto& [id, vLocalPos] : allCrossings) {
     // prepare the filler
     SiPixelClusterCollectionNew::FastFiller spc(*outputClusters, id);
 
     // retrieve the clusters of the right detId
     const auto& clustersOnDet = clusterCollection.find(DetId(id));
 
     if (clustersOnDet == clusterCollection.end())
       continue;
 
     // if the cluster DetSet is not valid, move on
     if (!(*clustersOnDet).isValid())
       continue;
 
     // if the module is bad continue
     if (throwBadComponents_ && SiPixelBadModule_->IsModuleBad(id))
       continue;
 
     // find all the clusters to put into the event
     const auto& clustersToPut = this->findAllNearbyClusters(clustersOnDet, id, vLocalPos);
 
     // find all the clusters to put into the event (different interface)
     //const auto& clustersToPut = this->findAllNearbyClusters(clusterCollection, id, vLocalPos);
 
     for (const auto& cluster : clustersToPut) {
       spc.push_back(*cluster);
     }
 
     if (spc.empty())
       spc.abort();
 
   }  // loop on trajectory crossings
 
   iEvent.put(clusterPutToken_, std::move(outputClusters));
 }