#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)

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)

bool	useRingSelector () const

void	useRingSelector (int minRing, int maxRing)

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::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 24 of file HitExtractorSTRP.h.

Member Typedef Documentation

◆ SiStripClusterRef

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

Definition at line 26 of file HitExtractorSTRP.h.

◆ SkipClustersCollection

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

private

Definition at line 80 of file HitExtractorSTRP.h.

◆ SkipPhase2ClustersCollection

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

private

Definition at line 81 of file HitExtractorSTRP.h.

Constructor & Destructor Documentation

◆ HitExtractorSTRP()

HitExtractorSTRP::HitExtractorSTRP	(	GeomDetEnumerators::SubDetector	subdet,
		TrackerDetSide	side,
		int	idLayer,
		float	iminGoodCharge
	)

Definition at line 29 of file HitExtractorSTRP.cc.

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

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

Referenced by clone().

◆ ~HitExtractorSTRP()

ctfseeding::HitExtractorSTRP::~HitExtractorSTRP ( )

inlineoverride

Definition at line 29 of file HitExtractorSTRP.h.

29 {}

Member Function Documentation

◆ cleanedOfClusters()

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

Definition at line 116 of file HitExtractorSTRP.cc.

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

References cms::cuda::assert(), relativeConstraints::empty, ev, hasMatchedHits, hasRPhiHits, hasStereoHits, hasVectorHits, hits(), sistrip::SpyUtilities::isValid(), OmniClusterRef::key(), dqmdumpme::last, LogTrace, VectorHit::lowerClusterRef(), edm::ContainerMask< T >::mask(), ctfseeding::HitExtractor::maskCluster, muonTagProbeFilters_cff::matched, AlCaHLTBitMon_ParallelJobs::p, python.rootplot.root2matplotlib::replace(), runEdmFileComparison::skipped, skipThis(), theSkipClusters, theSkipPhase2Clusters, and VectorHit::upperClusterRef().

Referenced by hits().

◆ clone()

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

inlineoverridevirtual

Implements ctfseeding::HitExtractor.

Definition at line 34 of file HitExtractorSTRP.h.

34 { return new HitExtractorSTRP(*this); }

References HitExtractorSTRP().

◆ getMinMaxRing()

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

inline

Definition at line 75 of file HitExtractorSTRP.h.

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

References theMaxRing, and theMinRing.

◆ hits()

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

overridevirtual

Implements ctfseeding::HitExtractor.

Definition at line 185 of file HitExtractorSTRP.cc.

                                                                          {
   LogDebug("HitExtractorSTRP") << "HitExtractorSTRP::hits";
   HitExtractor::Hits result;
   unsigned int cleanFrom = 0;
  
   //Retrieve tracker topology from geometry
   edm::ESHandle<TrackerTopology> tTopoHandle;
   es.get<TrackerTopologyRcd>().get(tTopoHandle);
   const TrackerTopology* const tTopo = tTopoHandle.product();
  
   //
   // 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;
 }

References cleanedOfClusters(), mps_fire::end, edmNew::DetSetVector< T >::equal_range(), ev, edm::EventSetup::get(), get, hasMatchedHits, hasRPhiHits, hasSimpleRphiHitsCleaner, hasStereoHits, hasVectorHits, LogDebug, LogTrace, minAbsZ, ctfseeding::HitExtractor::minGoodCharge, edm::ESHandle< T >::product(), FastTimerService_cff::range, ctfseeding::range2SeedingHits(), mps_fire::result, relativeConstraints::ring, ringRange(), ctfseeding::HitExtractor::skipClusters, GeomDetEnumerators::TEC, TrackerTopology::tecDetIdWheelComparator(), TrackerTopology::tecRing(), theIdLayer, theLayerSubDet, theMatchedHits, theRPhiHits, theSide, theStereoHits, theVectorHits, GeomDetEnumerators::TIB, TrackerTopology::tibDetIdLayerComparator(), GeomDetEnumerators::TID, TrackerTopology::tidDetIdWheelComparator(), TrackerTopology::tidRing(), GeomDetEnumerators::TOB, TrackerTopology::tobDetIdLayerComparator(), HLT_FULL_cff::vectorHits, and hit::z.

Referenced by cleanedOfClusters().

◆ ringRange()

bool HitExtractorSTRP::ringRange ( int ring ) const

private

Definition at line 61 of file HitExtractorSTRP.cc.

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

References hasRingSelector, relativeConstraints::ring, theMaxRing, and theMinRing.

Referenced by hits().

◆ setMinAbsZ()

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

inline

Definition at line 72 of file HitExtractorSTRP.h.

72 { minAbsZ = minZToSet; }

References minAbsZ.

◆ setNoProjection()

void ctfseeding::HitExtractorSTRP::setNoProjection ( )

inline

Definition at line 71 of file HitExtractorSTRP.h.

71 { failProjection = true; }

References failProjection.

◆ skipThis() [1/2]

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

Definition at line 79 of file HitExtractorSTRP.cc.

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

References cms::cuda::assert(), TkTransientTrackingRecHitBuilder::cloner(), failProjection, TrackingRecHit::geographicalId(), LogDebug, muonTagProbeFilters_cff::matched, and DetId::rawId().

Referenced by cleanedOfClusters().

◆ skipThis() [2/2]

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

Definition at line 67 of file HitExtractorSTRP.cc.

                                                                                                {
   if (maskCluster && (stripClusterMask->mask(clus.key())))
     return true;
  
   if UNLIKELY (minGoodCharge <= 0)
     return false;
   return siStripClusterTools::chargePerCM(id, *clus.cluster_strip()) <= minGoodCharge;
 }

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

◆ useMatchedHits()

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

inline

Definition at line 36 of file HitExtractorSTRP.h.

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

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

◆ useRingSelector() [1/2]

bool ctfseeding::HitExtractorSTRP::useRingSelector ( ) const

inline

Definition at line 74 of file HitExtractorSTRP.h.

74 { return hasRingSelector; }

References hasRingSelector.

◆ useRingSelector() [2/2]

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

Definition at line 55 of file HitExtractorSTRP.cc.

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

References hasRingSelector, HLT_FULL_cff::maxRing, HLT_FULL_cff::minRing, theMaxRing, and theMinRing.

◆ useRPhiHits()

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

inline

Definition at line 40 of file HitExtractorSTRP.h.

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

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

◆ useSimpleRphiHitsCleaner()

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

inline

Definition at line 54 of file HitExtractorSTRP.h.

54 { hasSimpleRphiHitsCleaner = use; }

References hasSimpleRphiHitsCleaner.

◆ useSkipClusters_()

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

overrideprivatevirtual

Implements ctfseeding::HitExtractor.

Definition at line 50 of file HitExtractorSTRP.cc.

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

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

◆ useStereoHits()

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

inline

Definition at line 44 of file HitExtractorSTRP.h.

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

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

◆ useVectorHits()

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

inline

Definition at line 49 of file HitExtractorSTRP.h.

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

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

Member Data Documentation

◆ failProjection

bool ctfseeding::HitExtractorSTRP::failProjection

private

Definition at line 102 of file HitExtractorSTRP.h.

Referenced by setNoProjection(), and skipThis().

◆ hasMatchedHits

bool ctfseeding::HitExtractorSTRP::hasMatchedHits

private

Definition at line 96 of file HitExtractorSTRP.h.

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

◆ hasRingSelector

bool ctfseeding::HitExtractorSTRP::hasRingSelector

private

Definition at line 100 of file HitExtractorSTRP.h.

Referenced by ringRange(), and useRingSelector().

◆ hasRPhiHits

bool ctfseeding::HitExtractorSTRP::hasRPhiHits

private

Definition at line 97 of file HitExtractorSTRP.h.

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

◆ hasSimpleRphiHitsCleaner

bool ctfseeding::HitExtractorSTRP::hasSimpleRphiHitsCleaner

private

Definition at line 101 of file HitExtractorSTRP.h.

Referenced by hits(), and useSimpleRphiHitsCleaner().

◆ hasStereoHits

bool ctfseeding::HitExtractorSTRP::hasStereoHits

private

Definition at line 98 of file HitExtractorSTRP.h.

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

◆ hasVectorHits

bool ctfseeding::HitExtractorSTRP::hasVectorHits

private

Definition at line 99 of file HitExtractorSTRP.h.

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

◆ minAbsZ

double ctfseeding::HitExtractorSTRP::minAbsZ

private

Definition at line 88 of file HitExtractorSTRP.h.

Referenced by hits(), and setMinAbsZ().

◆ theIdLayer

int ctfseeding::HitExtractorSTRP::theIdLayer

private

Definition at line 87 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theLayerSubDet

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

private

Definition at line 85 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theMatchedHits

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

private

Definition at line 92 of file HitExtractorSTRP.h.

Referenced by hits(), and useMatchedHits().

◆ theMaxRing

int ctfseeding::HitExtractorSTRP::theMaxRing

private

Definition at line 89 of file HitExtractorSTRP.h.

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

◆ theMinRing

int ctfseeding::HitExtractorSTRP::theMinRing

private

Definition at line 89 of file HitExtractorSTRP.h.

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

◆ theRPhiHits

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

private

Definition at line 93 of file HitExtractorSTRP.h.

Referenced by hits(), and useRPhiHits().

◆ theSide

TrackerDetSide ctfseeding::HitExtractorSTRP::theSide

private

Definition at line 86 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theSkipClusters

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

private

Definition at line 90 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

◆ theSkipPhase2Clusters

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

private

Definition at line 91 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

◆ theStereoHits

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

private

Definition at line 94 of file HitExtractorSTRP.h.

Referenced by hits(), and useStereoHits().

◆ theVectorHits

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

private

Definition at line 95 of file HitExtractorSTRP.h.

Referenced by hits(), and useVectorHits().

Public Types

Public Member Functions

Private Types

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ SiStripClusterRef

◆ SkipClustersCollection

◆ SkipPhase2ClustersCollection

Constructor & Destructor Documentation

◆ HitExtractorSTRP()

◆ ~HitExtractorSTRP()

Member Function Documentation

◆ cleanedOfClusters()

◆ clone()

◆ getMinMaxRing()

◆ hits()

◆ ringRange()

◆ setMinAbsZ()

◆ setNoProjection()

◆ skipThis() [1/2]

◆ skipThis() [2/2]

◆ useMatchedHits()

◆ useRingSelector() [1/2]

◆ useRingSelector() [2/2]

◆ useRPhiHits()

◆ useSimpleRphiHitsCleaner()

◆ useSkipClusters_()

◆ useStereoHits()

◆ useVectorHits()

Member Data Documentation

◆ failProjection

◆ hasMatchedHits

◆ hasRingSelector

◆ hasRPhiHits

◆ hasSimpleRphiHitsCleaner

◆ hasStereoHits

◆ hasVectorHits

◆ minAbsZ

◆ theIdLayer

◆ theLayerSubDet

◆ theMatchedHits

◆ theMaxRing

◆ theMinRing

◆ theRPhiHits

◆ theSide

◆ theSkipClusters

◆ theSkipPhase2Clusters

◆ theStereoHits

◆ theVectorHits