#include <SiStripApproximateCluster.h>

Public Member Functions
cms_uint8_t	avgCharge () const

cms_uint16_t	barycenter () const

bool	filter () const

bool	isSaturated () const

bool	peakFilter () const

	SiStripApproximateCluster ()

	SiStripApproximateCluster (cms_uint16_t barycenter, cms_uint8_t width, cms_uint8_t avgCharge, bool filter, bool isSaturated, bool peakFilter=false)

	SiStripApproximateCluster (const SiStripCluster &cluster, unsigned int maxNSat, float hitPredPos, bool peakFilter)

cms_uint8_t	width () const

Private Attributes
cms_uint8_t	avgCharge_ = 0

cms_uint16_t	barycenter_ = 0

bool	filter_ = false

bool	isSaturated_ = false

bool	peakFilter_ = false

cms_uint8_t	width_ = 0

Static Private Attributes
static constexpr double	maxTrimmedSizeDiffNeg_ = .7

static constexpr double	maxTrimmedSizeDiffPos_ = 1.

static constexpr double	trimMaxADC_ = 30.

static constexpr double	trimMaxFracNeigh_ = .25

static constexpr double	trimMaxFracTotal_ = .15

Detailed Description

Definition at line 7 of file SiStripApproximateCluster.h.

Constructor & Destructor Documentation

◆ SiStripApproximateCluster() [1/3]

SiStripApproximateCluster::SiStripApproximateCluster ( )

inline

Definition at line 9 of file SiStripApproximateCluster.h.

9 {}

◆ SiStripApproximateCluster() [2/3]

SiStripApproximateCluster::SiStripApproximateCluster	(	cms_uint16_t	barycenter,
		cms_uint8_t	width,
		cms_uint8_t	avgCharge,
		bool	filter,
		bool	isSaturated,
		bool	peakFilter = `false`
	)

inlineexplicit

Definition at line 11 of file SiStripApproximateCluster.h.

       : barycenter_(barycenter),
         width_(width),
         avgCharge_(avgCharge),
         filter_(filter),
         isSaturated_(isSaturated),
         peakFilter_(peakFilter) {}

◆ SiStripApproximateCluster() [3/3]

SiStripApproximateCluster::SiStripApproximateCluster	(	const SiStripCluster &	cluster,
		unsigned int	maxNSat,
		float	hitPredPos,
		bool	peakFilter
	)

explicit

Definition at line 6 of file SiStripApproximateCluster.cc.

References funct::abs(), SiStripCluster::amplitudes(), avgCharge_, SiStripCluster::barycenter(), barycenter_, reco::ceil(), SiStripCluster::charge(), filter_, mps_fire::i, isSaturated_, dqmdumpme::last, HLT_2023v12_cff::maxNSat, maxTrimmedSizeDiffNeg_, maxTrimmedSizeDiffPos_, dqmiodumpmetadata::n, peakFilter(), peakFilter_, SiStripCluster::size(), trimMaxADC_, trimMaxFracNeigh_, trimMaxFracTotal_, and width_.

                                                                       {
   barycenter_ = std::round(cluster.barycenter() * 10);
   width_ = cluster.size();
   avgCharge_ = cluster.charge() / cluster.size();
   filter_ = false;
   isSaturated_ = false;
   peakFilter_ = peakFilter;
 
   //mimicing the algorithm used in StripSubClusterShapeTrajectoryFilter...
   //Looks for 3 adjacent saturated strips (ADC>=254)
   const auto& ampls = cluster.amplitudes();
   unsigned int thisSat = (ampls[0] >= 254), maxSat = thisSat;
   for (unsigned int i = 1, n = ampls.size(); i < n; ++i) {
     if (ampls[i] >= 254) {
       thisSat++;
     } else if (thisSat > 0) {
       maxSat = std::max<int>(maxSat, thisSat);
       thisSat = 0;
     }
   }
   if (thisSat > 0) {
     maxSat = std::max<int>(maxSat, thisSat);
   }
   if (maxSat >= maxNSat) {
     filter_ = true;
     isSaturated_ = true;
   }
 
   unsigned int hitStripsTrim = ampls.size();
   int sum = std::accumulate(ampls.begin(), ampls.end(), 0);
   uint8_t trimCut = std::min<uint8_t>(trimMaxADC_, std::floor(trimMaxFracTotal_ * sum));
   auto begin = ampls.begin();
   auto last = ampls.end() - 1;
   while (hitStripsTrim > 1 && (*begin < std::max<uint8_t>(trimCut, trimMaxFracNeigh_ * (*(begin + 1))))) {
     hitStripsTrim--;
     ++begin;
   }
   while (hitStripsTrim > 1 && (*last < std::max<uint8_t>(trimCut, trimMaxFracNeigh_ * (*(last - 1))))) {
     hitStripsTrim--;
     --last;
   }
   if (hitStripsTrim < std::floor(std::abs(hitPredPos) - maxTrimmedSizeDiffNeg_)) {
     filter_ = false;
   } else if (hitStripsTrim <= std::ceil(std::abs(hitPredPos) + maxTrimmedSizeDiffPos_)) {
     filter_ = true;
   } else {
     filter_ = peakFilter_;
   }
 }