ctfseeding::HitExtractorSTRP Class Reference

#include <HitExtractorSTRP.h>

Inheritance diagram for ctfseeding::HitExtractorSTRP:

Public Types
typedef SiStripRecHit2D::ClusterRef	SiStripClusterRef
Public Types inherited from ctfseeding::HitExtractor
using	HitPointer = mayown_ptr< BaseTrackerRecHit >
using	Hits = std::vector< HitPointer >
using	TkHit = BaseTrackerRecHit
using	TkHitRef = BaseTrackerRecHit const &

Public Member Functions
void	cleanedOfClusters (const TkTransientTrackingRecHitBuilder &ttrhBuilder, const edm::Event &ev, HitExtractor::Hits &hits, bool matched, unsigned int cleanFrom=0) const
HitExtractorSTRP *	clone () const override
std::tuple< int, int >	getMinMaxRing () const
	HitExtractorSTRP (GeomDetEnumerators::SubDetector subdet, TrackerDetSide side, int idLayer, float iminGoodCharge, edm::ConsumesCollector &iC)
HitExtractor::Hits	hits (const TkTransientTrackingRecHitBuilder &ttrhBuilder, const edm::Event &, const edm::EventSetup &) const override
void	setMinAbsZ (double minZToSet)
void	setNoProjection ()
std::pair< bool, ProjectedSiStripRecHit2D * >	skipThis (const TkTransientTrackingRecHitBuilder &ttrhBuilder, TkHitRef matched, edm::Handle< edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > > > &stripClusterMask) const
bool	skipThis (DetId id, OmniClusterRef const &clus, edm::Handle< edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > > > &stripClusterMask) const
void	useMatchedHits (const edm::InputTag &m, edm::ConsumesCollector &iC)
void	useRingSelector (int minRing, int maxRing)
bool	useRingSelector () const
void	useRPhiHits (const edm::InputTag &m, edm::ConsumesCollector &iC)
void	useSimpleRphiHitsCleaner (bool use)
void	useStereoHits (const edm::InputTag &m, edm::ConsumesCollector &iC)
void	useVectorHits (const edm::InputTag &m, edm::ConsumesCollector &iC)
	~HitExtractorSTRP () override
Public Member Functions inherited from ctfseeding::HitExtractor
	HitExtractor ()
void	useSkipClusters (const edm::InputTag &m, edm::ConsumesCollector &iC)
virtual	~HitExtractor ()

Private Types
typedef edm::ContainerMask < edmNew::DetSetVector < SiStripCluster > >	SkipClustersCollection
typedef edm::ContainerMask < Phase2TrackerCluster1DCollectionNew >	SkipPhase2ClustersCollection

Private Member Functions
bool	ringRange (int ring) const
void	useSkipClusters_ (const edm::InputTag &m, edm::ConsumesCollector &iC) override

Private Attributes
bool	failProjection
bool	hasMatchedHits
bool	hasRingSelector
bool	hasRPhiHits
bool	hasSimpleRphiHitsCleaner
bool	hasStereoHits
bool	hasVectorHits
double	minAbsZ
int	theIdLayer
const GeomDetEnumerators::SubDetector	theLayerSubDet
edm::EDGetTokenT < SiStripMatchedRecHit2DCollection >	theMatchedHits
int	theMaxRing
int	theMinRing
edm::EDGetTokenT < SiStripRecHit2DCollection >	theRPhiHits
TrackerDetSide	theSide
edm::EDGetTokenT < SkipClustersCollection >	theSkipClusters
edm::EDGetTokenT < SkipPhase2ClustersCollection >	theSkipPhase2Clusters
edm::EDGetTokenT < SiStripRecHit2DCollection >	theStereoHits
edm::ESGetToken < TrackerTopology, TrackerTopologyRcd >	theTtopo
edm::EDGetTokenT < VectorHitCollection >	theVectorHits

Additional Inherited Members
Public Attributes inherited from ctfseeding::HitExtractor
bool	filterCluster = false
bool	maskCluster = false
float	minGoodCharge = 0
bool	skipClusters = false

Detailed Description

Definition at line 26 of file HitExtractorSTRP.h.

Member Typedef Documentation

typedef SiStripRecHit2D::ClusterRef ctfseeding::HitExtractorSTRP::SiStripClusterRef

Definition at line 28 of file HitExtractorSTRP.h.

typedef edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > ctfseeding::HitExtractorSTRP::SkipClustersCollection

private

Definition at line 86 of file HitExtractorSTRP.h.

typedef edm::ContainerMask<Phase2TrackerCluster1DCollectionNew> ctfseeding::HitExtractorSTRP::SkipPhase2ClustersCollection

private

Definition at line 87 of file HitExtractorSTRP.h.

Constructor & Destructor Documentation

HitExtractorSTRP::HitExtractorSTRP	(	GeomDetEnumerators::SubDetector	subdet,
		TrackerDetSide	side,
		int	idLayer,
		float	iminGoodCharge,
		edm::ConsumesCollector &	iC
	)

Definition at line 23 of file HitExtractorSTRP.cc.

References ctfseeding::HitExtractor::minGoodCharge, and ctfseeding::HitExtractor::skipClusters.

Referenced by clone().

    : theLayerSubDet(subdet),
      theSide(side),
      theIdLayer(idLayer),
      minAbsZ(0),
      theMinRing(1),
      theMaxRing(0),
      theTtopo(iC.esConsumes()),
      hasMatchedHits(false),
      hasRPhiHits(false),
      hasStereoHits(false),
      hasVectorHits(false),
      hasRingSelector(false),
      hasSimpleRphiHitsCleaner(true) {
  minGoodCharge = iminGoodCharge;
  if (minGoodCharge > 0)
    skipClusters = true;
}

ctfseeding::HitExtractorSTRP::~HitExtractorSTRP ( )

inlineoverride

Definition at line 35 of file HitExtractorSTRP.h.

{}

Member Function Documentation

void HitExtractorSTRP::cleanedOfClusters	(	const TkTransientTrackingRecHitBuilder &	ttrhBuilder,
		const edm::Event &	ev,
		HitExtractor::Hits &	hits,
		bool	matched,
		unsigned int	cleanFrom = 0
	)		const

Definition at line 112 of file HitExtractorSTRP.cc.

References cms::cuda::assert(), runTheMatrix::const, relativeConstraints::empty, edm::Event::getByToken(), hasMatchedHits, hasRPhiHits, hasStereoHits, hasVectorHits, sistrip::SpyUtilities::isValid(), OmniClusterRef::key(), dqmdumpme::last, LogTrace, VectorHit::lowerClusterRef(), ctfseeding::HitExtractor::maskCluster, AlCaHLTBitMon_ParallelJobs::p, python.rootplot.root2matplotlib::replace(), skipThis(), theSkipClusters, theSkipPhase2Clusters, and VectorHit::upperClusterRef().

Referenced by hits().

                                                                       {
  unsigned int skipped = 0;
  unsigned int projected = 0;
  if (hasMatchedHits || hasRPhiHits || hasStereoHits) {
    LogTrace("HitExtractorSTRP") << "getting " << hits.size() << " strip hit in input.";
    edm::Handle<SkipClustersCollection> stripClusterMask;
    if (maskCluster)
      ev.getByToken(theSkipClusters, stripClusterMask);
    for (unsigned int iH = cleanFrom; iH < hits.size(); ++iH) {
      assert(hits[iH]->isValid());
      auto id = hits[iH]->geographicalId();
      if (matched) {
        auto [replace, replaceMe] = skipThis(ttrhBuilder, *hits[iH], stripClusterMask);
        if (replace) {
          if (!replaceMe) {
            LogTrace("HitExtractorSTRP") << "skipping a matched hit on :" << hits[iH]->geographicalId().rawId();
            skipped++;
          } else {
            projected++;
          }
          hits[iH].reset(replaceMe);
          if (replaceMe == nullptr)
            assert(hits[iH].empty());
          else
            assert(hits[iH].isOwn());
        }
      } else if (skipThis(id, hits[iH]->firstClusterRef(), stripClusterMask)) {
        LogTrace("HitExtractorSTRP") << "skipping a hit on :" << hits[iH]->geographicalId().rawId() << " key: ";
        skipped++;
        hits[iH].reset();
      }
    }
  }
  if (hasVectorHits) {
    LogTrace("HitExtractorSTRP") << "getting " << hits.size() << " vector hit in input.";
    edm::Handle<SkipPhase2ClustersCollection> ph2ClusterMask;
    if (maskCluster)
      ev.getByToken(theSkipPhase2Clusters, ph2ClusterMask);
    for (unsigned int iH = cleanFrom; iH < hits.size(); ++iH) {
      LogTrace("HitExtractorSTRP") << "analizing hit on :" << hits[iH]->geographicalId().rawId();
      assert(hits[iH]->isValid());
      const VectorHit& vhit = dynamic_cast<VectorHit const&>(*hits[iH]);
      LogTrace("HitExtractorSTRP") << " key lower: " << vhit.lowerClusterRef().key()
                                   << " and key upper: " << vhit.upperClusterRef().key();
      LogTrace("HitExtractorSTRP") << " key lower: " << hits[iH]->firstClusterRef().key();

      //FIXME:: introduce a "projected" version later?
      if (maskCluster &&
          (ph2ClusterMask->mask(vhit.lowerClusterRef().key()) || ph2ClusterMask->mask(vhit.upperClusterRef().key()))) {
        LogTrace("HitExtractorSTRP") << "skipping a vector hit on :" << hits[iH]->geographicalId().rawId()
                                     << " key lower: " << vhit.lowerClusterRef().key()
                                     << " and key upper: " << vhit.upperClusterRef().key();
        skipped++;
        hits[iH].reset();
      }
    }
  }

  //  remove empty elements...
  auto last = std::remove_if(hits.begin() + cleanFrom, hits.end(), [](HitPointer const& p) { return p.empty(); });
  hits.resize(last - hits.begin());

  LogTrace("HitExtractorSTRP") << "skipped :" << skipped << " rechits because of clusters and projected: " << projected;
}

HitExtractorSTRP* ctfseeding::HitExtractorSTRP::clone ( void  ) const

inlineoverridevirtual

Implements ctfseeding::HitExtractor.

Definition at line 40 of file HitExtractorSTRP.h.

References HitExtractorSTRP().

{ return new HitExtractorSTRP(*this); }

std::tuple<int, int> ctfseeding::HitExtractorSTRP::getMinMaxRing ( ) const

inline

Definition at line 81 of file HitExtractorSTRP.h.

References theMaxRing, and theMinRing.

Referenced by SeedingLayerSetsBuilder::LayerSpec::print().

{ return std::make_tuple(theMinRing, theMaxRing); }

HitExtractor::Hits HitExtractorSTRP::hits ( const TkTransientTrackingRecHitBuilder &  ttrhBuilder, const edm::Event &  ev, const edm::EventSetup &  es  ) const

overridevirtual

Implements ctfseeding::HitExtractor.

Definition at line 181 of file HitExtractorSTRP.cc.

References cleanedOfClusters(), dataset::end, edm::Event::get(), edm::Event::getByToken(), edm::EventSetup::getData(), hasMatchedHits, hasRPhiHits, hasSimpleRphiHitsCleaner, hasStereoHits, hasVectorHits, LogDebug, LogTrace, minAbsZ, ctfseeding::HitExtractor::minGoodCharge, sistrip::SpyUtilities::range(), ctfseeding::range2SeedingHits(), mps_fire::result, relativeConstraints::ring, ringRange(), ctfseeding::HitExtractor::skipClusters, GeomDetEnumerators::TEC, TrackerTopology::tecDetIdWheelComparator(), TrackerTopology::tecRing(), theIdLayer, theLayerSubDet, theMatchedHits, theRPhiHits, theSide, theStereoHits, theTtopo, theVectorHits, GeomDetEnumerators::TIB, TrackerTopology::tibDetIdLayerComparator(), GeomDetEnumerators::TID, TrackerTopology::tidDetIdWheelComparator(), TrackerTopology::tidRing(), GeomDetEnumerators::TOB, TrackerTopology::tobDetIdLayerComparator(), vectorHits_cff::vectorHits, and hit::z.

                                                                         {
  LogDebug("HitExtractorSTRP") << "HitExtractorSTRP::hits";
  HitExtractor::Hits result;
  unsigned int cleanFrom = 0;

  //Retrieve tracker topology from geometry
  const TrackerTopology* const tTopo = &es.getData(theTtopo);

  //
  // TIB
  //
  if (theLayerSubDet == GeomDetEnumerators::TIB) {
    LogTrace("HitExtractorSTRP") << "Getting hits into the TIB";
    if (hasMatchedHits) {
      edm::Handle<SiStripMatchedRecHit2DCollection> matchedHits;
      ev.getByToken(theMatchedHits, matchedHits);
      if (skipClusters)
        cleanFrom = result.size();
      range2SeedingHits(*matchedHits, result, tTopo->tibDetIdLayerComparator(theIdLayer));
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, true, cleanFrom);
    }
    if (hasRPhiHits) {
      edm::Handle<SiStripRecHit2DCollection> rphiHits;
      ev.getByToken(theRPhiHits, rphiHits);
      if (hasMatchedHits) {
        if (!hasSimpleRphiHitsCleaner) {  // this is a brutal "cleaning". Add something smarter in the future
          if (skipClusters)
            cleanFrom = result.size();
          range2SeedingHits(*rphiHits, result, tTopo->tibDetIdLayerComparator(theIdLayer));
          if (skipClusters)
            cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
        }
      } else {
        if (skipClusters)
          cleanFrom = result.size();
        range2SeedingHits(*rphiHits, result, tTopo->tibDetIdLayerComparator(theIdLayer));
        if (skipClusters)
          cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
      }
    }
    if (hasStereoHits) {
      edm::Handle<SiStripRecHit2DCollection> stereoHits;
      ev.getByToken(theStereoHits, stereoHits);
      if (skipClusters)
        cleanFrom = result.size();
      range2SeedingHits(*stereoHits, result, tTopo->tibDetIdLayerComparator(theIdLayer));
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasVectorHits) {
      LogError("HitExtractorSTRP") << "TIB is not supposed to be in Phase2 TRK detector configuration. What follows "
                                      "have never been checked before! ";
      auto const& vectorHits = ev.get(theVectorHits);
      if (skipClusters)
        cleanFrom = result.size();
      range2SeedingHits(vectorHits, result, tTopo->tibDetIdLayerComparator(theIdLayer));
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }

  }

  //
  // TID
  //
  else if (theLayerSubDet == GeomDetEnumerators::TID) {
    LogTrace("HitExtractorSTRP") << "Getting hits into the TID";
    if (hasMatchedHits) {
      edm::Handle<SiStripMatchedRecHit2DCollection> matchedHits;
      ev.getByToken(theMatchedHits, matchedHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tidDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripMatchedRecHit2DCollection::Range range = matchedHits->equal_range(getter.first, getter.second);
      for (SiStripMatchedRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tidRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (SiStripMatchedRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end;
             ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, true, cleanFrom);
    }
    if (hasRPhiHits) {
      edm::Handle<SiStripRecHit2DCollection> rphiHits;
      ev.getByToken(theRPhiHits, rphiHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tidDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripRecHit2DCollection::Range range = rphiHits->equal_range(getter.first, getter.second);
      for (SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tidRing(it->detId());
        if (!ringRange(ring))
          continue;
        if ((SiStripDetId(it->detId()).partnerDetId() != 0) && hasSimpleRphiHitsCleaner)
          continue;  // this is a brutal "cleaning". Add something smarter in the future
        for (SiStripRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasStereoHits) {
      edm::Handle<SiStripRecHit2DCollection> stereoHits;
      ev.getByToken(theStereoHits, stereoHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tidDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripRecHit2DCollection::Range range = stereoHits->equal_range(getter.first, getter.second);
      for (SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tidRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (SiStripRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasVectorHits) {
      LogTrace("HitExtractorSTRP") << "Getting vector hits for IdLayer " << theIdLayer;
      auto const& vectorHits = ev.get(theVectorHits);
      //FIXME: check the skipClusters with VHits
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tidDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      VectorHitCollection::Range range = vectorHits.equal_range(getter.first, getter.second);
      for (VectorHitCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tidRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (VectorHitCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      LogTrace("HitExtractorSTRP") << "result size value:" << result.size();
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
  }
  //
  // TOB
  //
  else if (theLayerSubDet == GeomDetEnumerators::TOB) {
    LogTrace("HitExtractorSTRP") << "Getting hits into the TOB";
    if (hasMatchedHits) {
      edm::Handle<SiStripMatchedRecHit2DCollection> matchedHits;
      ev.getByToken(theMatchedHits, matchedHits);
      if (skipClusters)
        cleanFrom = result.size();
      if (minAbsZ > 0.) {
        auto getter = tTopo->tobDetIdLayerComparator(theIdLayer);
        SiStripMatchedRecHit2DCollection::Range range = matchedHits->equal_range(getter.first, getter.second);
        for (SiStripMatchedRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
          for (SiStripMatchedRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end;
               ++hit) {
            if (fabs(hit->globalPosition().z()) >= minAbsZ)
              result.emplace_back(*hit);
          }
        }
      } else {
        range2SeedingHits(*matchedHits, result, tTopo->tobDetIdLayerComparator(theIdLayer));
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, true, cleanFrom);
    }
    if (hasRPhiHits) {
      edm::Handle<SiStripRecHit2DCollection> rphiHits;
      ev.getByToken(theRPhiHits, rphiHits);
      if (hasMatchedHits) {
        if (!hasSimpleRphiHitsCleaner) {  // this is a brutal "cleaning". Add something smarter in the future
          if (skipClusters)
            cleanFrom = result.size();
          range2SeedingHits(*rphiHits, result, tTopo->tobDetIdLayerComparator(theIdLayer));
          if (skipClusters)
            cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
        }
      } else {
        if (skipClusters)
          cleanFrom = result.size();
        range2SeedingHits(*rphiHits, result, tTopo->tobDetIdLayerComparator(theIdLayer));
        if (skipClusters)
          cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
      }
    }
    if (hasStereoHits) {
      edm::Handle<SiStripRecHit2DCollection> stereoHits;
      ev.getByToken(theStereoHits, stereoHits);
      if (skipClusters)
        cleanFrom = result.size();
      range2SeedingHits(*stereoHits, result, tTopo->tobDetIdLayerComparator(theIdLayer));
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasVectorHits) {
      LogTrace("HitExtractorSTRP") << "Getting vector hits for IdLayer " << theIdLayer;
      edm::Handle<VectorHitCollection> vectorHits;
      ev.getByToken(theVectorHits, vectorHits);
      //FIXME: check the skipClusters with VHits
      if (skipClusters)
        cleanFrom = result.size();
      range2SeedingHits(*vectorHits, result, tTopo->tobDetIdLayerComparator(theIdLayer));
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }

  }

  //
  // TEC
  //
  else if (theLayerSubDet == GeomDetEnumerators::TEC) {
    LogTrace("HitExtractorSTRP") << "Getting hits into the TEC";
    if (hasMatchedHits) {
      edm::Handle<SiStripMatchedRecHit2DCollection> matchedHits;
      ev.getByToken(theMatchedHits, matchedHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tecDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripMatchedRecHit2DCollection::Range range = matchedHits->equal_range(getter.first, getter.second);
      for (SiStripMatchedRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tecRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (SiStripMatchedRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end;
             ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, true, cleanFrom);
    }
    if (hasRPhiHits) {
      edm::Handle<SiStripRecHit2DCollection> rphiHits;
      ev.getByToken(theRPhiHits, rphiHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tecDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripRecHit2DCollection::Range range = rphiHits->equal_range(getter.first, getter.second);
      for (SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tecRing(it->detId());
        if (!ringRange(ring))
          continue;
        if ((SiStripDetId(it->detId()).partnerDetId() != 0) && hasSimpleRphiHitsCleaner)
          continue;  // this is a brutal "cleaning". Add something smarter in the future
        for (SiStripRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasStereoHits) {
      edm::Handle<SiStripRecHit2DCollection> stereoHits;
      ev.getByToken(theStereoHits, stereoHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tecDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      SiStripRecHit2DCollection::Range range = stereoHits->equal_range(getter.first, getter.second);
      for (SiStripRecHit2DCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tecRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (SiStripRecHit2DCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
    if (hasVectorHits) {
      LogError("HitExtractorSTRP") << "TEC is not supposed to be in Phase2 TRK detector configuration. What follows "
                                      "have never been checked before! ";
      edm::Handle<VectorHitCollection> vectorHits;
      ev.getByToken(theVectorHits, vectorHits);
      if (skipClusters)
        cleanFrom = result.size();
      auto getter = tTopo->tidDetIdWheelComparator(static_cast<unsigned int>(theSide), theIdLayer);
      VectorHitCollection::Range range = vectorHits->equal_range(getter.first, getter.second);
      for (VectorHitCollection::const_iterator it = range.first; it != range.second; ++it) {
        int ring = tTopo->tidRing(it->detId());
        if (!ringRange(ring))
          continue;
        for (VectorHitCollection::DetSet::const_iterator hit = it->begin(), end = it->end(); hit != end; ++hit) {
          result.emplace_back(*hit);
        }
      }
      if (skipClusters)
        cleanedOfClusters(ttrhBuilder, ev, result, false, cleanFrom);
    }
  }

  LogDebug("HitExtractorSTRP") << " giving: " << result.size() << " out for charge cut " << minGoodCharge;
  return result;
}

bool HitExtractorSTRP::ringRange ( int  ring ) const

private

Definition at line 57 of file HitExtractorSTRP.cc.

References hasRingSelector, theMaxRing, and theMinRing.

Referenced by hits().

                                               {
  if (!hasRingSelector)
    return true;
  return (ring >= theMinRing) & (ring <= theMaxRing);
}

void ctfseeding::HitExtractorSTRP::setMinAbsZ ( double  minZToSet )

inline

Definition at line 78 of file HitExtractorSTRP.h.

References minAbsZ.

{ minAbsZ = minZToSet; }

void ctfseeding::HitExtractorSTRP::setNoProjection ( )

inline

Definition at line 77 of file HitExtractorSTRP.h.

References failProjection.

{ failProjection = true; }

std::pair< bool, ProjectedSiStripRecHit2D * > HitExtractorSTRP::skipThis	(	const TkTransientTrackingRecHitBuilder &	ttrhBuilder,
		TkHitRef	matched,
		edm::Handle< edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > > > &	stripClusterMask
	)		const

Definition at line 75 of file HitExtractorSTRP.cc.

References cms::cuda::assert(), TkTransientTrackingRecHitBuilder::cloner(), failProjection, TrackingRecHit::geographicalId(), OmniClusterRef::key(), LogDebug, SiStripMatchedRecHit2D::monoClusterRef(), DetId::rawId(), and SiStripMatchedRecHit2D::stereoClusterRef().

Referenced by cleanedOfClusters().

                                                                                               {
  const SiStripMatchedRecHit2D& hit = (SiStripMatchedRecHit2D const&)(matched);

  assert(dynamic_cast<SiStripMatchedRecHit2D const*>(&matched));

  auto id = hit.geographicalId();
  ProjectedSiStripRecHit2D* replaceMe = nullptr;
  bool rejectSt = skipThis(id, hit.stereoClusterRef(), stripClusterMask);
  bool rejectMono = skipThis(id, hit.monoClusterRef(), stripClusterMask);

  if ((!rejectSt) & (!rejectMono)) {
    // keepit
    return std::make_pair(false, replaceMe);
  }

  if (failProjection || (rejectSt & rejectMono)) {
    //only skip if both hits are done
    return std::make_pair(true, replaceMe);
  }

  // replace with one

  auto cloner = ttrhBuilder.cloner();
  replaceMe = cloner.project(hit, rejectSt, TrajectoryStateOnSurface()).release();
  if (rejectSt)
    LogDebug("HitExtractorSTRP") << "a matched hit is partially masked, and the mono hit got projected onto: "
                                 << replaceMe->geographicalId().rawId() << " key: " << hit.monoClusterRef().key();
  else
    LogDebug("HitExtractorSTRP") << "a matched hit is partially masked, and the stereo hit got projected onto: "
                                 << replaceMe->geographicalId().rawId() << " key: " << hit.stereoClusterRef().key();

  return std::make_pair(true, replaceMe);
}

bool HitExtractorSTRP::skipThis	(	DetId	id,
		OmniClusterRef const &	clus,
		edm::Handle< edm::ContainerMask< edmNew::DetSetVector< SiStripCluster > > > &	stripClusterMask
	)		const

Definition at line 63 of file HitExtractorSTRP.cc.

References siStripClusterTools::chargePerCM(), OmniClusterRef::cluster_strip(), OmniClusterRef::key(), ctfseeding::HitExtractor::maskCluster, ctfseeding::HitExtractor::minGoodCharge, and UNLIKELY.

                                                                                               {
  if (maskCluster && (stripClusterMask->mask(clus.key())))
    return true;

  if UNLIKELY (minGoodCharge <= 0)
    return false;
  return siStripClusterTools::chargePerCM(id, *clus.cluster_strip()) <= minGoodCharge;
}

void ctfseeding::HitExtractorSTRP::useMatchedHits ( const edm::InputTag &  m, edm::ConsumesCollector &  iC  )

inline

Definition at line 42 of file HitExtractorSTRP.h.

References edm::ConsumesCollector::consumes(), hasMatchedHits, visualization-live-secondInstance_cfg::m, and theMatchedHits.

                                                                        {
      hasMatchedHits = true;
      theMatchedHits = iC.consumes<SiStripMatchedRecHit2DCollection>(m);
    }

void HitExtractorSTRP::useRingSelector	(	int	minRing,
		int	maxRing
	)

Definition at line 51 of file HitExtractorSTRP.cc.

References hasRingSelector, theMaxRing, and theMinRing.

Referenced by SeedingLayerSetsBuilder::LayerSpec::print().

                                                               {
  hasRingSelector = true;
  theMinRing = minRing;
  theMaxRing = maxRing;
}

bool ctfseeding::HitExtractorSTRP::useRingSelector ( ) const

inline

Definition at line 80 of file HitExtractorSTRP.h.

References hasRingSelector.

{ return hasRingSelector; }

void ctfseeding::HitExtractorSTRP::useRPhiHits ( const edm::InputTag &  m, edm::ConsumesCollector &  iC  )

inline

Definition at line 46 of file HitExtractorSTRP.h.

References edm::ConsumesCollector::consumes(), hasRPhiHits, visualization-live-secondInstance_cfg::m, and theRPhiHits.

                                                                     {
      hasRPhiHits = true;
      theRPhiHits = iC.consumes<SiStripRecHit2DCollection>(m);
    }

void ctfseeding::HitExtractorSTRP::useSimpleRphiHitsCleaner ( bool  use )

inline

Definition at line 60 of file HitExtractorSTRP.h.

References hasSimpleRphiHitsCleaner.

{ hasSimpleRphiHitsCleaner = use; }

void HitExtractorSTRP::useSkipClusters_ ( const edm::InputTag &  m, edm::ConsumesCollector &  iC  )

overrideprivatevirtual

Implements ctfseeding::HitExtractor.

Definition at line 46 of file HitExtractorSTRP.cc.

References edm::ConsumesCollector::consumes(), visualization-live-secondInstance_cfg::m, theSkipClusters, and theSkipPhase2Clusters.

                                                                                      {
  theSkipClusters = iC.consumes<SkipClustersCollection>(m);
  theSkipPhase2Clusters = iC.consumes<SkipPhase2ClustersCollection>(m);
}

void ctfseeding::HitExtractorSTRP::useStereoHits ( const edm::InputTag &  m, edm::ConsumesCollector &  iC  )

inline

Definition at line 50 of file HitExtractorSTRP.h.

References edm::ConsumesCollector::consumes(), hasStereoHits, visualization-live-secondInstance_cfg::m, and theStereoHits.

                                                                       {
      hasStereoHits = true;
      theStereoHits = iC.consumes<SiStripRecHit2DCollection>(m);
    }

void ctfseeding::HitExtractorSTRP::useVectorHits ( const edm::InputTag &  m, edm::ConsumesCollector &  iC  )

inline

Definition at line 55 of file HitExtractorSTRP.h.

References edm::ConsumesCollector::consumes(), hasVectorHits, visualization-live-secondInstance_cfg::m, and theVectorHits.

                                                                       {
      hasVectorHits = true;
      theVectorHits = iC.consumes<VectorHitCollection>(m);
    }

Member Data Documentation

bool ctfseeding::HitExtractorSTRP::failProjection

private

Definition at line 109 of file HitExtractorSTRP.h.

Referenced by setNoProjection(), and skipThis().

bool ctfseeding::HitExtractorSTRP::hasMatchedHits

private

Definition at line 103 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), hits(), and useMatchedHits().

bool ctfseeding::HitExtractorSTRP::hasRingSelector

private

Definition at line 107 of file HitExtractorSTRP.h.

Referenced by ringRange(), and useRingSelector().

bool ctfseeding::HitExtractorSTRP::hasRPhiHits

private

Definition at line 104 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), hits(), and useRPhiHits().

bool ctfseeding::HitExtractorSTRP::hasSimpleRphiHitsCleaner

private

Definition at line 108 of file HitExtractorSTRP.h.

Referenced by hits(), and useSimpleRphiHitsCleaner().

bool ctfseeding::HitExtractorSTRP::hasStereoHits

private

Definition at line 105 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), hits(), and useStereoHits().

bool ctfseeding::HitExtractorSTRP::hasVectorHits

private

Definition at line 106 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), hits(), and useVectorHits().

double ctfseeding::HitExtractorSTRP::minAbsZ

private

Definition at line 94 of file HitExtractorSTRP.h.

Referenced by hits(), and setMinAbsZ().

int ctfseeding::HitExtractorSTRP::theIdLayer

private

Definition at line 93 of file HitExtractorSTRP.h.

Referenced by hits().

const GeomDetEnumerators::SubDetector ctfseeding::HitExtractorSTRP::theLayerSubDet

private

Definition at line 91 of file HitExtractorSTRP.h.

Referenced by hits().

edm::EDGetTokenT<SiStripMatchedRecHit2DCollection> ctfseeding::HitExtractorSTRP::theMatchedHits

private

Definition at line 98 of file HitExtractorSTRP.h.

Referenced by hits(), and useMatchedHits().

int ctfseeding::HitExtractorSTRP::theMaxRing

private

Definition at line 95 of file HitExtractorSTRP.h.

Referenced by getMinMaxRing(), ringRange(), and useRingSelector().

int ctfseeding::HitExtractorSTRP::theMinRing

private

Definition at line 95 of file HitExtractorSTRP.h.

Referenced by getMinMaxRing(), ringRange(), and useRingSelector().

edm::EDGetTokenT<SiStripRecHit2DCollection> ctfseeding::HitExtractorSTRP::theRPhiHits

private

Definition at line 99 of file HitExtractorSTRP.h.

Referenced by hits(), and useRPhiHits().

TrackerDetSide ctfseeding::HitExtractorSTRP::theSide

private

Definition at line 92 of file HitExtractorSTRP.h.

Referenced by hits().

edm::EDGetTokenT<SkipClustersCollection> ctfseeding::HitExtractorSTRP::theSkipClusters

private

Definition at line 96 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

edm::EDGetTokenT<SkipPhase2ClustersCollection> ctfseeding::HitExtractorSTRP::theSkipPhase2Clusters

private

Definition at line 97 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

edm::EDGetTokenT<SiStripRecHit2DCollection> ctfseeding::HitExtractorSTRP::theStereoHits

private

Definition at line 100 of file HitExtractorSTRP.h.

Referenced by hits(), and useStereoHits().

edm::ESGetToken<TrackerTopology, TrackerTopologyRcd> ctfseeding::HitExtractorSTRP::theTtopo

private

Definition at line 102 of file HitExtractorSTRP.h.

Referenced by hits().

edm::EDGetTokenT<VectorHitCollection> ctfseeding::HitExtractorSTRP::theVectorHits

private

Definition at line 101 of file HitExtractorSTRP.h.

Referenced by hits(), and useVectorHits().