Inheritance diagram for SeedClusterRemover:

Classes
struct	ParamBlock

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

	SeedClusterRemover (const edm::ParameterSet &iConfig)

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 Types
typedef edm::ContainerMask < edmNew::DetSetVector < SiPixelCluster > >	PixelMaskContainer

typedef edm::ContainerMask < edmNew::DetSetVector < SiStripCluster > >	StripMaskContainer

Private Member Functions
void	process (const TrackingRecHit hit, float chi2, const TrackerGeometry tg)

void	process (const OmniClusterRef &cluRef, uint32_t subdet)

void	readPSet (const edm::ParameterSet &iConfig, const std::string &name, int id1=-1, int id2=-1, int id3=-1, int id4=-1, int id5=-1, int id6=-1)

Private Attributes
std::vector< bool >	collectedPixels_

std::vector< bool >	collectedStrips_

bool	doPixel_

bool	doStrip_

bool	mergeOld_

edm::EDGetTokenT < PixelMaskContainer >	oldPxlMaskToken_

edm::EDGetTokenT < StripMaskContainer >	oldStrMaskToken_

ParamBlock	pblocks_ [NumberOfParamBlocks]

edm::EDGetTokenT < edmNew::DetSetVector < SiPixelCluster > >	pixelClusters_

std::vector< uint8_t >	pixels

edm::ProductID	pixelSourceProdID

edm::EDGetTokenT < edmNew::DetSetVector < SiStripCluster > >	stripClusters_

std::vector< uint8_t >	strips

edm::ProductID	stripSourceProdID

edm::EDGetTokenT < TrajectorySeedCollection >	trajectories_

edm::ESGetToken < TrackerGeometry, TrackerDigiGeometryRecord > const	tTrackerGeom_

Static Private Attributes
static const unsigned int	NumberOfParamBlocks = 6

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 35 of file SeedClusterRemover.cc.

Member Typedef Documentation

typedef edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > SeedClusterRemover::PixelMaskContainer

private

Definition at line 60 of file SeedClusterRemover.cc.

typedef edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > SeedClusterRemover::StripMaskContainer

private

Definition at line 61 of file SeedClusterRemover.cc.

Constructor & Destructor Documentation

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

Definition at line 113 of file SeedClusterRemover.cc.

References doPixel_, doStrip_, Exception, PVValHelper::fill(), edm::ParameterSet::getParameter(), mps_fire::i, mergeOld_, NumberOfParamBlocks, oldPxlMaskToken_, oldStrMaskToken_, pblocks_, pixelClusters_, readPSet(), stripClusters_, and trajectories_.

     : tTrackerGeom_(esConsumes<TrackerGeometry, TrackerDigiGeometryRecord>()),
       doStrip_(iConfig.existsAs<bool>("doStrip") ? iConfig.getParameter<bool>("doStrip") : true),
       doPixel_(iConfig.existsAs<bool>("doPixel") ? iConfig.getParameter<bool>("doPixel") : true),
       mergeOld_(iConfig.exists("oldClusterRemovalInfo")) {
   fill(pblocks_, pblocks_ + NumberOfParamBlocks, ParamBlock());
   readPSet(iConfig, "Common", -1);
   if (doPixel_) {
     readPSet(iConfig, "Pixel", 0, 1);
     readPSet(iConfig, "PXB", 0);
     readPSet(iConfig, "PXE", 1);
   }
   if (doStrip_) {
     readPSet(iConfig, "Strip", 2, 3, 4, 5);
     readPSet(iConfig, "StripInner", 2, 3);
     readPSet(iConfig, "StripOuter", 4, 5);
     readPSet(iConfig, "TIB", 2);
     readPSet(iConfig, "TID", 3);
     readPSet(iConfig, "TOB", 4);
     readPSet(iConfig, "TEC", 5);
   }
 
   bool usingCharge = false;
   for (size_t i = 0; i < NumberOfParamBlocks; ++i) {
     if (!pblocks_[i].isSet_)
       throw cms::Exception("Configuration Error")
           << "SeedClusterRemover: Missing configuration for detector with subDetID = " << (i + 1);
     if (pblocks_[i].usesCharge_ && !usingCharge) {
       throw cms::Exception("Configuration Error") << "SeedClusterRemover: Configuration for subDetID = " << (i + 1)
                                                   << " uses cluster charge, which is not enabled.";
     }
   }
 
   trajectories_ = consumes<TrajectorySeedCollection>(iConfig.getParameter<InputTag>("trajectories"));
   if (doPixel_)
     pixelClusters_ = consumes<edmNew::DetSetVector<SiPixelCluster> >(iConfig.getParameter<InputTag>("pixelClusters"));
   if (doStrip_)
     stripClusters_ = consumes<edmNew::DetSetVector<SiStripCluster> >(iConfig.getParameter<InputTag>("stripClusters"));
   if (mergeOld_) {
     oldPxlMaskToken_ = consumes<PixelMaskContainer>(iConfig.getParameter<InputTag>("oldClusterRemovalInfo"));
     oldStrMaskToken_ = consumes<StripMaskContainer>(iConfig.getParameter<InputTag>("oldClusterRemovalInfo"));
   }
 
   produces<edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > >();
   produces<edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > >();
 }

Member Function Documentation

void SeedClusterRemover::process	(	const TrackingRecHit *	hit,
		float	chi2,
		const TrackerGeometry *	tg
	)

inlineprivate

Definition at line 176 of file SeedClusterRemover.cc.

References cms::cuda::assert(), SiPixelRecHit::cluster(), collectedPixels_, doPixel_, doStrip_, Exception, TrackingRecHit::geographicalId(), TrackerGeometry::geomDetSubDetector(), edm::Ref< C, T, F >::id(), GeomDetEnumerators::isTrackerPixel(), edm::Ref< C, T, F >::key(), SeedClusterRemover::ParamBlock::maxSize_, SiStripMatchedRecHit2D::monoClusterRef(), NumberOfParamBlocks, TrackerSingleRecHit::omniClusterRef(), ProjectedSiStripRecHit2D::originalHit(), pblocks_, pixels, pixelSourceProdID, DetId::rawId(), SiStripMatchedRecHit2D::stereoClusterRef(), and DetId::subdetId().

Referenced by produce().

                                                                                                  {
   DetId detid = hit->geographicalId();
   uint32_t subdet = detid.subdetId();
 
   assert((subdet > 0) && (subdet <= NumberOfParamBlocks));
 
   // chi2 cut
   if (chi2 > pblocks_[subdet - 1].maxChi2_)
     return;
 
   if (GeomDetEnumerators::isTrackerPixel(tg->geomDetSubDetector(subdet))) {
     if (!doPixel_)
       return;
     // this is a pixel, and i *know* it is
     const SiPixelRecHit *pixelHit = static_cast<const SiPixelRecHit *>(hit);
 
     SiPixelRecHit::ClusterRef cluster = pixelHit->cluster();
 
     if (cluster.id() != pixelSourceProdID)
       throw cms::Exception("Inconsistent Data")
           << "SeedClusterRemover: pixel cluster ref from Product ID = " << cluster.id()
           << " does not match with source cluster collection (ID = " << pixelSourceProdID << ")\n.";
 
     assert(cluster.id() == pixelSourceProdID);
     //DBG// cout << "HIT NEW PIXEL DETID = " << detid.rawId() << ", Cluster [ " << cluster.key().first << " / " <<  cluster.key().second << " ] " << endl;
 
     // if requested, cut on cluster size
     if (pblocks_[subdet - 1].usesSize_ && (cluster->pixels().size() > pblocks_[subdet - 1].maxSize_))
       return;
 
     // mark as used
     pixels[cluster.key()] = false;
 
     assert(collectedPixels_.size() > cluster.key());
     collectedPixels_[cluster.key()] = true;
 
   } else {  // aka Strip
     if (!doStrip_)
       return;
     const type_info &hitType = typeid(*hit);
     if (hitType == typeid(SiStripRecHit2D)) {
       const SiStripRecHit2D *stripHit = static_cast<const SiStripRecHit2D *>(hit);
       //DBG//     cout << "Plain RecHit 2D: " << endl;
       process(stripHit->omniClusterRef(), subdet);
     } else if (hitType == typeid(SiStripRecHit1D)) {
       const SiStripRecHit1D *hit1D = static_cast<const SiStripRecHit1D *>(hit);
       process(hit1D->omniClusterRef(), subdet);
     } else if (hitType == typeid(SiStripMatchedRecHit2D)) {
       const SiStripMatchedRecHit2D *matchHit = static_cast<const SiStripMatchedRecHit2D *>(hit);
       //DBG//     cout << "Matched RecHit 2D: " << endl;
       process(matchHit->monoClusterRef(), subdet);
       process(matchHit->stereoClusterRef(), subdet);
     } else if (hitType == typeid(ProjectedSiStripRecHit2D)) {
       const ProjectedSiStripRecHit2D *projHit = static_cast<const ProjectedSiStripRecHit2D *>(hit);
       //DBG//     cout << "Projected RecHit 2D: " << endl;
       process(projHit->originalHit().omniClusterRef(), subdet);
     } else
       throw cms::Exception("NOT IMPLEMENTED")
           << "Don't know how to handle " << hitType.name() << " on detid " << detid.rawId() << "\n";
   }
 }

void SeedClusterRemover::process	(	const OmniClusterRef &	cluRef,
		uint32_t	subdet
	)

inlineprivate

Definition at line 160 of file SeedClusterRemover.cc.

References cms::cuda::assert(), OmniClusterRef::cluster_strip(), collectedStrips_, edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), SeedClusterRemover::ParamBlock::maxSize_, pblocks_, strips, and stripSourceProdID.

                                                                                 {
   SiStripRecHit2D::ClusterRef cluster = ocluster.cluster_strip();
   if (cluster.id() != stripSourceProdID)
     throw cms::Exception("Inconsistent Data")
         << "SeedClusterRemover: strip cluster ref from Product ID = " << cluster.id()
         << " does not match with source cluster collection (ID = " << stripSourceProdID << ")\n.";
 
   assert(cluster.id() == stripSourceProdID);
   if (pblocks_[subdet - 1].usesSize_ && (cluster->amplitudes().size() > pblocks_[subdet - 1].maxSize_))
     return;
 
   strips[cluster.key()] = false;
   assert(collectedStrips_.size() > cluster.key());
   collectedStrips_[cluster.key()] = true;
 }

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

override

Definition at line 238 of file SeedClusterRemover.cc.

References cms::cuda::assert(), collectedPixels_, collectedStrips_, submitPVResolutionJobs::count, doPixel_, doStrip_, PVValHelper::fill(), edm::Event::getByToken(), edm::EventSetup::getData(), edm::HandleBase::id(), LogDebug, mergeOld_, oldPxlMaskToken_, oldStrMaskToken_, HLT_FULL_cff::pixelClusters, pixelClusters_, pixels, pixelSourceProdID, process(), edm::Event::put(), fileCollector::seed, DetachedQuadStep_cff::seeds, HLT_FULL_cff::stripClusters, stripClusters_, strips, stripSourceProdID, trajectories_, and tTrackerGeom_.

                                                                         {
   const auto &tgh = &iSetup.getData(tTrackerGeom_);
 
   Handle<edmNew::DetSetVector<SiPixelCluster> > pixelClusters;
   if (doPixel_) {
     iEvent.getByToken(pixelClusters_, pixelClusters);
     pixelSourceProdID = pixelClusters.id();
   }
   //DBG// std::cout << "SeedClusterRemover: Read pixel " << pixelClusters_.encode() << " = ID " << pixelSourceProdID << std::endl;
 
   Handle<edmNew::DetSetVector<SiStripCluster> > stripClusters;
   if (doStrip_) {
     iEvent.getByToken(stripClusters_, stripClusters);
     stripSourceProdID = stripClusters.id();
   }
   //DBG// std::cout << "SeedClusterRemover: Read strip " << stripClusters_.encode() << " = ID " << stripSourceProdID << std::endl;
 
   if (doStrip_) {
     strips.resize(stripClusters->dataSize());
     fill(strips.begin(), strips.end(), true);
   }
   if (doPixel_) {
     pixels.resize(pixelClusters->dataSize());
     fill(pixels.begin(), pixels.end(), true);
   }
   if (mergeOld_) {
     edm::Handle<PixelMaskContainer> oldPxlMask;
     edm::Handle<StripMaskContainer> oldStrMask;
     iEvent.getByToken(oldPxlMaskToken_, oldPxlMask);
     iEvent.getByToken(oldStrMaskToken_, oldStrMask);
     LogDebug("SeedClusterRemover") << "to merge in, " << oldStrMask->size() << " strp and " << oldPxlMask->size()
                                    << " pxl";
     oldStrMask->copyMaskTo(collectedStrips_);
     oldPxlMask->copyMaskTo(collectedPixels_);
     assert(stripClusters->dataSize() >= collectedStrips_.size());
     collectedStrips_.resize(stripClusters->dataSize(), false);  // for ondemand
   } else {
     collectedStrips_.resize(stripClusters->dataSize(), false);
     collectedPixels_.resize(pixelClusters->dataSize(), false);
   }
 
   edm::Handle<TrajectorySeedCollection> seeds;
   iEvent.getByToken(trajectories_, seeds);
 
   for (auto const &seed : (*seeds)) {
     for (auto const &hit : seed.recHits()) {
       if (!hit.isValid())
         continue;
       process(&hit, 0., tgh);
     }
   }
 
   pixels.clear();
   strips.clear();
 
   LogDebug("SeedClusterRemover") << "total strip to skip: "
                                  << std::count(collectedStrips_.begin(), collectedStrips_.end(), true);
   iEvent.put(std::make_unique<StripMaskContainer>(edm::RefProd<edmNew::DetSetVector<SiStripCluster> >(stripClusters),
                                                   collectedStrips_));
 
   LogDebug("SeedClusterRemover") << "total pxl to skip: "
                                  << std::count(collectedPixels_.begin(), collectedPixels_.end(), true);
   iEvent.put(std::make_unique<PixelMaskContainer>(edm::RefProd<edmNew::DetSetVector<SiPixelCluster> >(pixelClusters),
                                                   collectedPixels_));
 
   collectedStrips_.clear();
   collectedPixels_.clear();
 }

void SeedClusterRemover::readPSet	(	const edm::ParameterSet &	iConfig,
		const std::string &	name,
		int	id1 = `-1`,
		int	id2 = `-1`,
		int	id3 = `-1`,
		int	id4 = `-1`,
		int	id5 = `-1`,
		int	id6 = `-1`
	)

private

Definition at line 91 of file SeedClusterRemover.cc.

References edm::ParameterSet::exists(), PVValHelper::fill(), and edm::ParameterSet::getParameter().

Referenced by SeedClusterRemover().

                                                                                                                  {
   if (iConfig.exists(name)) {
     ParamBlock pblock(iConfig.getParameter<ParameterSet>(name));
     if (id1 == -1) {
       fill(pblocks_, pblocks_ + NumberOfParamBlocks, pblock);
     } else {
       pblocks_[id1] = pblock;
       if (id2 != -1)
         pblocks_[id2] = pblock;
       if (id3 != -1)
         pblocks_[id3] = pblock;
       if (id4 != -1)
         pblocks_[id4] = pblock;
       if (id5 != -1)
         pblocks_[id5] = pblock;
       if (id6 != -1)
         pblocks_[id6] = pblock;
     }
   }
 }