median.h File Reference

#include <vector>
#include <algorithm>
#include <cmath>

Go to the source code of this file.

Functions
template<class T >
T	median (std::vector< T > values)

Function Documentation

median()

template<class T >

T median

(

std::vector< T >

values

)

Median in one dimension. T must be sortable, 'addable', and dividable by 2.

Definition at line 16 of file median.h.

References reco::ceil(), createfilelist::int, runTheMatrix::ret, findQualityFiles::size, and contentValuesCheck::values.

Referenced by SiStripAPVRestorer::baselineFollower(), SiStripAPVRestorer::baselineFollowerRestore(), BoostedJetONNXJetTagsProducer::BoostedJetONNXJetTagsProducer(), SiStripMonitorDigi::FillApvShotsMap(), PVValHelper::getMAD(), TrackerOfflineValidationSummary::getMedian(), PVValHelper::getMedian(), TrackerOfflineValidation::getMedian(), SiStripHotStripAlgorithmFromClusterOccupancy::iterativeSearch(), SiStripBadAPVandHotStripAlgorithmFromClusterOccupancy::iterativeSearch(), KDTreeLinkerAlgo< reco::PFRecHit const *>::medianSearch(), PuppiContainer::puppiWeights(), KDTreeLinkerAlgo< reco::PFRecHit const *>::recSearch(), CompareToMedian::runTest(), and SiStripZeroSuppression::storeBaseline().

                              {
  switch (values.size()) {
    case 0:
      //      throw LogicError ("Median of empty vector");
    case 1:
      return values[0];
    case 2:
      return (values[0] + values[1]) / 2.;
  };
  T ret;
  const int size = values.size();
  const int half = (int)ceil(size / 2.);
  partial_sort(values.begin(), values.begin() + half + 1, values.end());
  if (size & 1) {
    ret = values[half - 1];
  } else {
    ret = (values[half - 1] + values[half]) / 2.;
  };
  return ret;
}