SeedClusterRemover Class Reference

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
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final

Private Types
typedef edm::ContainerMask< edmNew::DetSetVector< SiPixelCluster > >	PixelMaskContainer
typedef edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > >	StripMaskContainer

Private Member Functions
void	process (const OmniClusterRef &cluRef, uint32_t subdet)
void	process (const TrackingRecHit *hit, float chi2, const TrackerGeometry *tg)
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_
const edm::ESGetToken< TrackerGeometry, TrackerDigiGeometryRecord >	tTrackerGeom_

Static Private Attributes
static const unsigned int	NumberOfParamBlocks = 6

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache

Detailed Description

Definition at line 35 of file SeedClusterRemover.cc.

Member Typedef Documentation

PixelMaskContainer

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

private

Definition at line 60 of file SeedClusterRemover.cc.

StripMaskContainer

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

private

Definition at line 61 of file SeedClusterRemover.cc.

Constructor & Destructor Documentation

SeedClusterRemover()

SeedClusterRemover::SeedClusterRemover

(

const edm::ParameterSet &

iConfig

)

Definition at line 113 of file SeedClusterRemover.cc.

    : 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> > >();
}

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

Member Function Documentation

process() [1/2]

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

inlineprivate

Definition at line 160 of file SeedClusterRemover.cc.

                                                                                {
  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;
}

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.

process() [2/2]

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

inlineprivate

Definition at line 176 of file SeedClusterRemover.cc.

                                                                                                 {
  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";
  }
}

References cms::cuda::assert(), hltPixelTracks_cff::chi2, SiPixelRecHit::cluster(), collectedPixels_, doPixel_, doStrip_, Exception, 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().

produce()

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

override

Definition at line 238 of file SeedClusterRemover.cc.

                                                                        {
  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();
}

References cms::cuda::assert(), collectedPixels_, collectedStrips_, edm::ContainerMask< T >::copyMaskTo(), submitPVResolutionJobs::count, doPixel_, doStrip_, ntuplemaker::fill, edm::EventSetup::getData(), iEvent, LogDebug, mergeOld_, oldPxlMaskToken_, oldStrMaskToken_, LumiMonitor_cff::pixelClusters, pixelClusters_, pixels, pixelSourceProdID, process(), fileCollector::seed, InitialStep_cff::seeds, edm::ContainerMask< T >::size(), HLT_FULL_cff::stripClusters, stripClusters_, strips, stripSourceProdID, trajectories_, and tTrackerGeom_.

readPSet()

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.

                                                                                                                 {
  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;
    }
  }
}

References edm::ParameterSet::exists(), ntuplemaker::fill, edm::ParameterSet::getParameter(), globals_cff::id1, globals_cff::id2, and Skims_PA_cff::name.

Referenced by SeedClusterRemover().

Member Data Documentation

collectedPixels_

std::vector<bool> SeedClusterRemover::collectedPixels_

private

Definition at line 85 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

collectedStrips_

std::vector<bool> SeedClusterRemover::collectedStrips_

private

Definition at line 84 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

doPixel_

bool SeedClusterRemover::doPixel_

private

Definition at line 57 of file SeedClusterRemover.cc.

Referenced by process(), produce(), and SeedClusterRemover().

doStrip_

bool SeedClusterRemover::doStrip_

private

Definition at line 57 of file SeedClusterRemover.cc.

Referenced by process(), produce(), and SeedClusterRemover().

mergeOld_

bool SeedClusterRemover::mergeOld_

private

Definition at line 58 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

NumberOfParamBlocks

const unsigned int SeedClusterRemover::NumberOfParamBlocks = 6

staticprivate

Definition at line 55 of file SeedClusterRemover.cc.

Referenced by process(), and SeedClusterRemover().

oldPxlMaskToken_

edm::EDGetTokenT<PixelMaskContainer> SeedClusterRemover::oldPxlMaskToken_

private

Definition at line 64 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

oldStrMaskToken_

edm::EDGetTokenT<StripMaskContainer> SeedClusterRemover::oldStrMaskToken_

private

Definition at line 65 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

pblocks_

ParamBlock SeedClusterRemover::pblocks_[NumberOfParamBlocks]

private

Definition at line 68 of file SeedClusterRemover.cc.

Referenced by process(), and SeedClusterRemover().

pixelClusters_

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

private

Definition at line 62 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

pixels

std::vector<uint8_t> SeedClusterRemover::pixels

private

Definition at line 78 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

pixelSourceProdID

edm::ProductID SeedClusterRemover::pixelSourceProdID

private

Definition at line 79 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

stripClusters_

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

private

Definition at line 63 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

strips

std::vector<uint8_t> SeedClusterRemover::strips

private

Definition at line 78 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

stripSourceProdID

edm::ProductID SeedClusterRemover::stripSourceProdID

private

Definition at line 79 of file SeedClusterRemover.cc.

Referenced by process(), and produce().

trajectories_

edm::EDGetTokenT<TrajectorySeedCollection> SeedClusterRemover::trajectories_

private

Definition at line 66 of file SeedClusterRemover.cc.

Referenced by produce(), and SeedClusterRemover().

tTrackerGeom_

const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> SeedClusterRemover::tTrackerGeom_

private

Definition at line 41 of file SeedClusterRemover.cc.

Referenced by produce().