#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

◆ SiStripClusterRef

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

Definition at line 28 of file HitExtractorSTRP.h.

◆ SkipClustersCollection

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

private

Definition at line 86 of file HitExtractorSTRP.h.

◆ SkipPhase2ClustersCollection

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

private

Definition at line 87 of file HitExtractorSTRP.h.

Constructor & Destructor Documentation

◆ HitExtractorSTRP()

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

◆ ~HitExtractorSTRP()

ctfseeding::HitExtractorSTRP::~HitExtractorSTRP ( )

inlineoverride

Definition at line 35 of file HitExtractorSTRP.h.

35 {}

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 112 of file HitExtractorSTRP.cc.

References cms::cuda::assert(), runTheMatrix::const, relativeConstraints::empty, makeMEIFBenchmarkPlots::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().

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

◆ clone()

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

inlineoverridevirtual

Implements ctfseeding::HitExtractor.

Definition at line 40 of file HitExtractorSTRP.h.

References HitExtractorSTRP().

40 { return new HitExtractorSTRP(*this); }

ctfseeding::HitExtractorSTRP::HitExtractorSTRP

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

Definition: HitExtractorSTRP.cc:23

◆ getMinMaxRing()

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

inline

Definition at line 81 of file HitExtractorSTRP.h.

References theMaxRing, and theMinRing.

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

ctfseeding::HitExtractorSTRP::theMaxRing

int theMaxRing

Definition: HitExtractorSTRP.h:95

ctfseeding::HitExtractorSTRP::theMinRing

int theMinRing

Definition: HitExtractorSTRP.h:95

◆ hits()

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(), edmNew::DetSetVector< T >::equal_range(), makeMEIFBenchmarkPlots::ev, edm::EventSetup::getData(), hasMatchedHits, hasRPhiHits, hasSimpleRphiHitsCleaner, hasStereoHits, hasVectorHits, LogDebug, LogTrace, minAbsZ, ctfseeding::HitExtractor::minGoodCharge, 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, theTtopo, theVectorHits, GeomDetEnumerators::TIB, TrackerTopology::tibDetIdLayerComparator(), GeomDetEnumerators::TID, TrackerTopology::tidDetIdWheelComparator(), TrackerTopology::tidRing(), GeomDetEnumerators::TOB, TrackerTopology::tobDetIdLayerComparator(), HLT_2023v12_cff::vectorHits, and hit::z.

Referenced by cleanedOfClusters().

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

◆ ringRange()

bool HitExtractorSTRP::ringRange ( int ring ) const

private

Definition at line 57 of file HitExtractorSTRP.cc.

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

Referenced by hits().

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

◆ setMinAbsZ()

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

inline

Definition at line 78 of file HitExtractorSTRP.h.

References minAbsZ.

78 { minAbsZ = minZToSet; }

ctfseeding::HitExtractorSTRP::minAbsZ

double minAbsZ

Definition: HitExtractorSTRP.h:94

◆ setNoProjection()

void ctfseeding::HitExtractorSTRP::setNoProjection ( )

inline

Definition at line 77 of file HitExtractorSTRP.h.

References failProjection.

77 { failProjection = true; }

ctfseeding::HitExtractorSTRP::failProjection

bool failProjection

Definition: HitExtractorSTRP.h:109

◆ 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 75 of file HitExtractorSTRP.cc.

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

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

◆ skipThis() [2/2]

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

◆ useMatchedHits()

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

◆ useRingSelector() [1/2]

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

Definition at line 51 of file HitExtractorSTRP.cc.

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

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

◆ useRingSelector() [2/2]

bool ctfseeding::HitExtractorSTRP::useRingSelector ( ) const

inline

Definition at line 80 of file HitExtractorSTRP.h.

References hasRingSelector.

80 { return hasRingSelector; }

ctfseeding::HitExtractorSTRP::hasRingSelector

bool hasRingSelector

Definition: HitExtractorSTRP.h:107

◆ useRPhiHits()

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

◆ useSimpleRphiHitsCleaner()

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

inline

Definition at line 60 of file HitExtractorSTRP.h.

References hasSimpleRphiHitsCleaner.

60 { hasSimpleRphiHitsCleaner = use; }

ctfseeding::HitExtractorSTRP::hasSimpleRphiHitsCleaner

bool hasSimpleRphiHitsCleaner

Definition: HitExtractorSTRP.h:108

◆ useSkipClusters_()

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

◆ useStereoHits()

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

◆ useVectorHits()

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

◆ failProjection

bool ctfseeding::HitExtractorSTRP::failProjection

private

Definition at line 109 of file HitExtractorSTRP.h.

Referenced by setNoProjection(), and skipThis().

◆ hasMatchedHits

bool ctfseeding::HitExtractorSTRP::hasMatchedHits

private

Definition at line 103 of file HitExtractorSTRP.h.

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

◆ hasRingSelector

bool ctfseeding::HitExtractorSTRP::hasRingSelector

private

Definition at line 107 of file HitExtractorSTRP.h.

Referenced by ringRange(), and useRingSelector().

◆ hasRPhiHits

bool ctfseeding::HitExtractorSTRP::hasRPhiHits

private

Definition at line 104 of file HitExtractorSTRP.h.

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

◆ hasSimpleRphiHitsCleaner

bool ctfseeding::HitExtractorSTRP::hasSimpleRphiHitsCleaner

private

Definition at line 108 of file HitExtractorSTRP.h.

Referenced by hits(), and useSimpleRphiHitsCleaner().

◆ hasStereoHits

bool ctfseeding::HitExtractorSTRP::hasStereoHits

private

Definition at line 105 of file HitExtractorSTRP.h.

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

◆ hasVectorHits

bool ctfseeding::HitExtractorSTRP::hasVectorHits

private

Definition at line 106 of file HitExtractorSTRP.h.

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

◆ minAbsZ

double ctfseeding::HitExtractorSTRP::minAbsZ

private

Definition at line 94 of file HitExtractorSTRP.h.

Referenced by hits(), and setMinAbsZ().

◆ theIdLayer

int ctfseeding::HitExtractorSTRP::theIdLayer

private

Definition at line 93 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theLayerSubDet

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

private

Definition at line 91 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theMatchedHits

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

private

Definition at line 98 of file HitExtractorSTRP.h.

Referenced by hits(), and useMatchedHits().

◆ theMaxRing

int ctfseeding::HitExtractorSTRP::theMaxRing

private

Definition at line 95 of file HitExtractorSTRP.h.

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

◆ theMinRing

int ctfseeding::HitExtractorSTRP::theMinRing

private

Definition at line 95 of file HitExtractorSTRP.h.

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

◆ theRPhiHits

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

private

Definition at line 99 of file HitExtractorSTRP.h.

Referenced by hits(), and useRPhiHits().

◆ theSide

TrackerDetSide ctfseeding::HitExtractorSTRP::theSide

private

Definition at line 92 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theSkipClusters

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

private

Definition at line 96 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

◆ theSkipPhase2Clusters

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

private

Definition at line 97 of file HitExtractorSTRP.h.

Referenced by cleanedOfClusters(), and useSkipClusters_().

◆ theStereoHits

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

private

Definition at line 100 of file HitExtractorSTRP.h.

Referenced by hits(), and useStereoHits().

◆ theTtopo

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

private

Definition at line 102 of file HitExtractorSTRP.h.

Referenced by hits().

◆ theVectorHits

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

private

Definition at line 101 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

◆ theTtopo

◆ theVectorHits