Inheritance diagram for ClusterCompatibilityProducer:

Classes
struct	ContainedHits

struct	VertexHit

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

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

	~ClusterCompatibilityProducer () override

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

Private Member Functions
ContainedHits	getContainedHits (const std::vector< VertexHit > &hits, double z0) const

Private Attributes
edm::InputTag	inputTag_

edm::EDGetTokenT< SiPixelRecHitCollection >	inputToken_

double	maxZ_

double	minZ_

edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	trackerToken_

double	zStep_

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

Definition at line 37 of file ClusterCompatibilityProducer.cc.

Constructor & Destructor Documentation

◆ ClusterCompatibilityProducer()

ClusterCompatibilityProducer::ClusterCompatibilityProducer ( const edm::ParameterSet & config )

explicit

Definition at line 67 of file ClusterCompatibilityProducer.cc.

     : inputTag_(config.getParameter<edm::InputTag>("inputTag")),
       minZ_(config.getParameter<double>("minZ")),
       maxZ_(config.getParameter<double>("maxZ")),
       zStep_(config.getParameter<double>("zStep")) {
   inputToken_ = consumes<SiPixelRecHitCollection>(inputTag_);
   trackerToken_ = esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>();
   LogDebug("") << "Using the " << inputTag_ << " input collection";
   produces<reco::ClusterCompatibility>();
 }

References inputTag_, inputToken_, LogDebug, and trackerToken_.

◆ ~ClusterCompatibilityProducer()

ClusterCompatibilityProducer::~ClusterCompatibilityProducer ( )

override

Definition at line 78 of file ClusterCompatibilityProducer.cc.

78 {}

Member Function Documentation

◆ getContainedHits()

ClusterCompatibilityProducer::ContainedHits ClusterCompatibilityProducer::getContainedHits	(	const std::vector< VertexHit > &	hits,
		double	z0
	)		const

private

Definition at line 140 of file ClusterCompatibilityProducer.cc.

                                                        {
   // Calculate number of hits contained in v-shaped window in cluster y-width vs. z-position.
   int n = 0;
   double chi = 0.;
  
   for (std::vector<VertexHit>::const_iterator hit = hits.begin(); hit != hits.end(); hit++) {
     // the calculation of the predicted cluster width p was
     // marked 'FIXME' in the HLTPixelClusterShapeFilter. It should
     // be revisited but is retained as it was for compatibility with the
     // older filter.
     double p = 2 * fabs(hit->z - z0) / hit->r + 0.5;
     if (fabs(p - hit->w) <= 1.) {
       chi += fabs(p - hit->w);
       n++;
     }
   }
   ClusterCompatibilityProducer::ContainedHits output;
   output.z0 = z0;
   output.nHit = n;
   output.chi = chi;
   return output;
 }

References hfClusterShapes_cfi::hits, dqmiodumpmetadata::n, convertSQLitetoXML_cfg::output, AlCaHLTBitMon_ParallelJobs::p, hit::z, and HLTMuonOfflineAnalyzer_cfi::z0.

Referenced by produce().

◆ produce()

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

override

Definition at line 80 of file ClusterCompatibilityProducer.cc.

                                                                                         {
   auto creco = std::make_unique<reco::ClusterCompatibility>();
  
   // get hold of products from Event
   edm::Handle<SiPixelRecHitCollection> hRecHits;
   iEvent.getByToken(inputToken_, hRecHits);
  
   // get tracker geometry
   if (hRecHits.isValid()) {
     edm::ESHandle<TrackerGeometry> trackerHandle = iSetup.getHandle(trackerToken_);
     const TrackerGeometry *tgeo = trackerHandle.product();
     const SiPixelRecHitCollection *hits = hRecHits.product();
  
     // loop over pixel rechits
     int nPxlHits = 0;
     std::vector<VertexHit> vhits;
     for (SiPixelRecHitCollection::DataContainer::const_iterator hit = hits->data().begin(), end = hits->data().end();
          hit != end;
          ++hit) {
       if (!hit->isValid())
         continue;
       ++nPxlHits;
       DetId id(hit->geographicalId());
       if (id.subdetId() != int(PixelSubdetector::PixelBarrel))
         continue;
       const PixelGeomDetUnit *pgdu = static_cast<const PixelGeomDetUnit *>(tgeo->idToDet(id));
       const PixelTopology *pixTopo = &(pgdu->specificTopology());
       std::vector<SiPixelCluster::Pixel> pixels(hit->cluster()->pixels());
       bool pixelOnEdge = false;
       for (std::vector<SiPixelCluster::Pixel>::const_iterator pixel = pixels.begin(); pixel != pixels.end(); ++pixel) {
         int pixelX = pixel->x;
         int pixelY = pixel->y;
         if (pixTopo->isItEdgePixelInX(pixelX) || pixTopo->isItEdgePixelInY(pixelY)) {
           pixelOnEdge = true;
           break;
         }
       }
       if (pixelOnEdge)
         continue;
  
       LocalPoint lpos = LocalPoint(hit->localPosition().x(), hit->localPosition().y(), hit->localPosition().z());
       GlobalPoint gpos = pgdu->toGlobal(lpos);
       VertexHit vh;
       vh.z = gpos.z();
       vh.r = gpos.perp();
       vh.w = hit->cluster()->sizeY();
       vhits.push_back(vh);
     }
  
     creco->setNValidPixelHits(nPxlHits);
  
     // append cluster compatibility  for each z-position
     for (double z0 = minZ_; z0 <= maxZ_; z0 += zStep_) {
       ContainedHits c = getContainedHits(vhits, z0);
       creco->append(c.z0, c.nHit, c.chi);
     }
   }
   iEvent.put(std::move(creco));
 }