|
|
|
template<typename T > |
std::unique_ptr< edmNew::DetSetVector< T > > | cleanup (const edmNew::DetSetVector< T > &oldClusters, const std::vector< uint8_t > &isGood, reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices *oldRefs) |
|
void | mergeOld (reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices &oldRefs) |
|
void | process (const OmniClusterRef &cluRef, SiStripDetId &detid, bool fromTrack) |
|
void | process (const TrackingRecHit *hit, unsigned char chi2, const TrackerGeometry *tg) |
|
void | readPSet (const edm::ParameterSet &iConfig, const std::string &name, int id1=-1, int id2=-1, int id3=-1, int id4=-1, int id5=-1, int id6=-1) |
|
float | sensorThickness (const SiStripDetId &detid) const |
|
Definition at line 42 of file HITrackClusterRemover.cc.
◆ PixelMaskContainer
◆ StripMaskContainer
◆ HITrackClusterRemover()
Definition at line 145 of file HITrackClusterRemover.cc.
152 iConfig.
existsAs<
bool>(
"doStripChargeCheck") ? iConfig.
getParameter<
bool>(
"doStripChargeCheck") :
false),
154 iConfig.
existsAs<
bool>(
"doPixelChargeCheck") ? iConfig.
getParameter<
bool>(
"doPixelChargeCheck") :
false),
160 if (iConfig.
exists(
"overrideTrkQuals"))
162 if (iConfig.
exists(
"clusterLessSolution"))
166 <<
"HITrackClusterRemover: Charge check asked without cluster collection ";
169 <<
"HITrackClusterRemover: Pixel cluster charge check not yet implemented";
172 produces<edmNew::DetSetVector<SiPixelCluster> >();
174 produces<edmNew::DetSetVector<SiStripCluster> >();
176 produces<ClusterRemovalInfo>();
188 readPSet(iConfig,
"Strip", 2, 3, 4, 5);
189 readPSet(iConfig,
"StripInner", 2, 3);
190 readPSet(iConfig,
"StripOuter", 4, 5);
197 bool usingCharge =
false;
201 <<
"HITrackClusterRemover: Missing configuration for detector with subDetID = " << (
i + 1);
202 if (
pblocks_[
i].usesCharge_ && !usingCharge) {
203 throw cms::Exception(
"Configuration Error") <<
"HITrackClusterRemover: Configuration for subDetID = " << (
i + 1)
204 <<
" uses cluster charge, which is not enabled.";
209 produces<edm::ContainerMask<edmNew::DetSetVector<SiPixelCluster> > >();
210 produces<edm::ContainerMask<edmNew::DetSetVector<SiStripCluster> > >();
214 if (iConfig.
exists(
"TrackQuality")) {
218 iConfig.
existsAs<
int>(
"minNumberOfLayersWithMeasBeforeFiltering")
219 ? iConfig.
getParameter<
int>(
"minNumberOfLayersWithMeasBeforeFiltering")
References cms::cuda::assert(), clusterWasteSolution_, doPixel_, doPixelChargeCheck_, doStrip_, doStripChargeCheck_, doTracks_, Exception, edm::ParameterSet::exists(), edm::ParameterSet::existsAs(), ntuplemaker::fill, filterTracks_, edm::ParameterSet::getParameter(), mps_fire::i, HLT_2018_cff::InputTag, label, mergeOld_, minNumberOfLayersWithMeasBeforeFiltering_, NumberOfParamBlocks, oldPxlMaskToken_, oldRemovalInfo_, oldStrMaskToken_, overrideTrkQuals_, pblocks_, pixelClusters_, reco::TrackBase::qualityByName(), readPSet(), rphiRecHitToken_, stereoRecHitToken_, AlCaHLTBitMon_QueryRunRegistry::string, stripClusters_, stripRecHits_, trackQuality_, tracks_, and reco::TrackBase::undefQuality.
◆ ~HITrackClusterRemover()
HITrackClusterRemover::~HITrackClusterRemover |
( |
| ) |
|
|
override |
◆ cleanup()
◆ mergeOld()
◆ process() [1/2]
Definition at line 309 of file HITrackClusterRemover.cc.
313 <<
"HITrackClusterRemover: strip cluster ref from Product ID = " << cluster.
id()
314 <<
" does not match with source cluster collection (ID = " <<
stripSourceProdID <<
")\n.";
322 for (std::vector<uint8_t>::const_iterator iAmp = cluster->amplitudes().begin(); iAmp != cluster->amplitudes().end();
326 if (
pblocks_[subdet - 1].cutOnStripCharge_ &&
References cms::cuda::assert(), OmniClusterRef::cluster_strip(), clusterWasteSolution_, collectedStrips_, edm::Ref< C, T, F >::id(), edm::Ref< C, T, F >::key(), HITrackClusterRemover::ParamBlock::maxSize_, pblocks_, sensorThickness(), strips, stripSourceProdID, and DetId::subdetId().
◆ process() [2/2]
Definition at line 340 of file HITrackClusterRemover.cc.
360 <<
"HITrackClusterRemover: pixel cluster ref from Product ID = " << cluster.
id()
361 <<
" does not match with source cluster collection (ID = " <<
pixelSourceProdID <<
")\n.";
382 const type_info &hitType =
typeid(*hit);
401 <<
"Don't know how to handle " << hitType.name() <<
" on detid " << detid.
rawId() <<
"\n";
References cms::cuda::assert(), hltPixelTracks_cff::chi2, SiPixelRecHit::cluster(), clusterWasteSolution_, collectedPixels_, doPixel_, doStrip_, Exception, TrackerGeometry::geomDetSubDetector(), edm::Ref< C, T, F >::id(), GeomDetEnumerators::isTrackerPixel(), edm::Ref< C, T, F >::key(), HITrackClusterRemover::ParamBlock::maxSize_, SiStripMatchedRecHit2D::monoClusterRef(), NumberOfParamBlocks, TrackerSingleRecHit::omniClusterRef(), ProjectedSiStripRecHit2D::originalHit(), pblocks_, pixels, pixelSourceProdID, DetId::rawId(), SiStripMatchedRecHit2D::stereoClusterRef(), DetId::subdetId(), and Traj2TrackHits::toChi2x5().
Referenced by produce().
◆ produce()
Definition at line 415 of file HITrackClusterRemover.cc.
416 ProductID pixelOldProdID, stripOldProdID;
435 std::unique_ptr<ClusterRemovalInfo> cri;
451 cri->getOldClustersFrom(*oldRemovalInfo);
460 <<
"HITrackClusterRemover: "
463 <<
"\t but the *old* ClusterRemovalInfo " <<
labels.productInstance <<
" refers as 'new product ids' to "
466 <<
"NOTA BENE: when running HITrackClusterRemover with an old ClusterRemovalInfo the hits in the trajectory "
467 "MUST be already re-keyed.\n";
487 LogDebug(
"HITrackClusterRemover") <<
"to merge in, " << oldStrMask->
size() <<
" strp and " << oldPxlMask->
size()
502 std::vector<Handle<edm::ValueMap<int> > > quals;
511 if (!quals.empty()) {
512 int qual = (*(quals[0])).get(it++);
526 auto const &chi2sX5 =
track.extra()->chi2sX5();
528 auto hb =
track.recHitsBegin();
529 for (
unsigned int h = 0;
h <
track.recHitsSize();
h++) {
567 cri->setNewPixelClusters(newPixels);
573 cri->setNewStripClusters(newStrips);
590 LogDebug(
"HITrackClusterRemover") <<
"total strip to skip: "
596 LogDebug(
"HITrackClusterRemover") <<
"total pxl to skip: "
References cms::cuda::assert(), cleanup(), clusterWasteSolution_, collectedPixels_, collectedStrips_, edm::ContainerMask< T >::copyMaskTo(), KineDebug3::count(), doPixel_, doStrip_, doStripChargeCheck_, doTracks_, Exception, ntuplemaker::fill, filterTracks_, edm::EventSetup::get(), get, h, es_hardcode_cfi::hb, edm::RefProd< C >::id(), iEvent, SummaryClient_cfi::labels, LogDebug, mergeOld_, minNumberOfLayersWithMeasBeforeFiltering_, eostools::move(), oldPxlMaskToken_, oldRemovalInfo_, oldStrMaskToken_, overrideTrkQuals_, LumiMonitor_cff::pixelClusters, pixelClusters_, reco::ClusterRemovalInfo::pixelIndices(), reco::ClusterRemovalInfo::pixelNewRefProd(), reco::ClusterRemovalInfo::pixelRefProd(), pixels, pixelSourceProdID, process(), edm::ESHandle< T >::product(), rpcPointValidation_cfi::recHit, LaserSeedGenerator_cfi::rphiRecHits, rphiRecHitToken_, edm::ContainerMask< T >::size(), LaserSeedGenerator_cfi::stereoRecHits, stereoRecHitToken_, HLT_2018_cff::stripClusters, stripClusters_, reco::ClusterRemovalInfo::stripIndices(), reco::ClusterRemovalInfo::stripNewRefProd(), reco::ClusterRemovalInfo::stripRefProd(), strips, stripSourceProdID, HLT_2018_cff::track, trackQuality_, PDWG_EXOHSCP_cff::tracks, and tracks_.
◆ readPSet()
void HITrackClusterRemover::readPSet |
( |
const edm::ParameterSet & |
iConfig, |
|
|
const std::string & |
name, |
|
|
int |
id1 = -1 , |
|
|
int |
id2 = -1 , |
|
|
int |
id3 = -1 , |
|
|
int |
id4 = -1 , |
|
|
int |
id5 = -1 , |
|
|
int |
id6 = -1 |
|
) |
| |
|
private |
◆ sensorThickness()
float HITrackClusterRemover::sensorThickness |
( |
const SiStripDetId & |
detid | ) |
const |
|
private |
◆ clusterWasteSolution_
bool HITrackClusterRemover::clusterWasteSolution_ |
|
private |
◆ collectedPixels_
std::vector<bool> HITrackClusterRemover::collectedPixels_ |
|
private |
◆ collectedStrips_
std::vector<bool> HITrackClusterRemover::collectedStrips_ |
|
private |
◆ doPixel_
bool HITrackClusterRemover::doPixel_ |
|
private |
◆ doPixelChargeCheck_
bool HITrackClusterRemover::doPixelChargeCheck_ |
|
private |
◆ doStrip_
bool HITrackClusterRemover::doStrip_ |
|
private |
◆ doStripChargeCheck_
bool HITrackClusterRemover::doStripChargeCheck_ |
|
private |
◆ doTracks_
bool HITrackClusterRemover::doTracks_ |
|
private |
◆ filterTracks_
bool HITrackClusterRemover::filterTracks_ |
|
private |
◆ mergeOld_
bool HITrackClusterRemover::mergeOld_ |
|
private |
◆ minNumberOfLayersWithMeasBeforeFiltering_
int HITrackClusterRemover::minNumberOfLayersWithMeasBeforeFiltering_ |
|
private |
◆ NumberOfParamBlocks
const unsigned int HITrackClusterRemover::NumberOfParamBlocks = 6 |
|
staticprivate |
◆ oldPxlMaskToken_
◆ oldRemovalInfo_
◆ oldStrMaskToken_
◆ overrideTrkQuals_
◆ pblocks_
◆ pixelClusters_
◆ pixelRecHits_
std::string HITrackClusterRemover::pixelRecHits_ |
|
private |
◆ pixels
std::vector<uint8_t> HITrackClusterRemover::pixels |
|
private |
◆ pixelSourceProdID
◆ rphiRecHitToken_
◆ stereoRecHitToken_
◆ stripClusters_
◆ stripRecHits_
std::string HITrackClusterRemover::stripRecHits_ |
|
private |
◆ strips
std::vector<uint8_t> HITrackClusterRemover::strips |
|
private |
◆ stripSourceProdID
◆ trackQuality_
◆ tracks_
T const * product() const
SiStripModuleGeometry moduleGeometry() const
edm::EDGetTokenT< edmNew::DetSetVector< SiPixelCluster > > pixelClusters_
ProductID id() const
Accessor for product ID.
edm::EDGetTokenT< PixelMaskContainer > oldPxlMaskToken_
static unsigned char toChi2x5(float chi2)
ParamBlock pblocks_[NumberOfParamBlocks]
FWCore Framework interface EventSetupRecordImplementation h
Helper function to determine trigger accepts.
const SiStripClusterRefProd & stripRefProd() const
boost::transform_iterator< IterHelp, const_IdIter > const_iterator
void process(const TrackingRecHit *hit, unsigned char chi2, const TrackerGeometry *tg)
const_iterator begin(bool update=false) const
reco::TrackBase::TrackQuality trackQuality_
bool existsAs(std::string const ¶meterName, bool trackiness=true) const
checks if a parameter exists as a given type
ClusterRef cluster() const
std::string pixelRecHits_
std::vector< uint8_t > strips
edm::ProductID pixelSourceProdID
void readPSet(const edm::ParameterSet &iConfig, const std::string &name, int id1=-1, int id2=-1, int id3=-1, int id4=-1, int id5=-1, int id6=-1)
edm::ProductID stripSourceProdID
std::vector< bool > collectedPixels_
std::vector< bool > collectedStrips_
float sensorThickness(const SiStripDetId &detid) const
SiStripRecHit2D originalHit() const
int minNumberOfLayersWithMeasBeforeFiltering_
edm::EDGetTokenT< SiStripRecHit2DCollection > rphiRecHitToken_
edm::EDGetTokenT< edmNew::DetSetVector< SiStripCluster > > stripClusters_
size_type dataSize() const
constexpr int subdetId() const
get the contents of the subdetector field (not cast into any detector's numbering enum)
static const unsigned int NumberOfParamBlocks
edm::EDGetTokenT< reco::TrackCollection > tracks_
bool exists(std::string const ¶meterName) const
checks if a parameter exists
std::vector< edm::EDGetTokenT< edm::ValueMap< int > > > overrideTrkQuals_
const GeomDetEnumerators::SubDetector geomDetSubDetector(int subdet) const
OmniClusterRef const & stereoClusterRef() const
edm::EDGetTokenT< reco::ClusterRemovalInfo > oldRemovalInfo_
std::string stripRecHits_
bool isTrackerPixel(GeomDetEnumerators::SubDetector m)
static constexpr auto TOB
std::vector< data_type > DataContainer
const SiStripClusterRefProd & stripNewRefProd() const
OmniClusterRef const & monoClusterRef() const
const SiPixelClusterRefProd & pixelRefProd() const
edm::EDGetTokenT< StripMaskContainer > oldStrMaskToken_
static TrackQuality qualityByName(const std::string &name)
T getParameter(std::string const &) const
void mergeOld(reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices &oldRefs)
ProductID id() const
Accessor for product ID.
constexpr uint32_t rawId() const
get the raw id
data_type const * data(size_t cell) const
void copyMaskTo(std::vector< bool > &) const
const SiPixelClusterRefProd & pixelNewRefProd() const
const_iterator end(bool update=false) const
std::unique_ptr< edmNew::DetSetVector< T > > cleanup(const edmNew::DetSetVector< T > &oldClusters, const std::vector< uint8_t > &isGood, reco::ClusterRemovalInfo::Indices &refs, const reco::ClusterRemovalInfo::Indices *oldRefs)
static constexpr auto TIB
std::vector< uint8_t > pixels
edm::EDGetTokenT< SiStripRecHit2DCollection > stereoRecHitToken_
key_type key() const
Accessor for product key.
Detector identifier class for the strip tracker.
bool clusterWasteSolution_
OmniClusterRef const & omniClusterRef() const