LikelihoodFitDeDxEstimator Class Reference

#include <LikelihoodFitDeDxEstimator.h>

Inheritance diagram for LikelihoodFitDeDxEstimator:

Public Member Functions
std::pair< float, float >	dedx (const reco::DeDxHitCollection &Hits) override
	LikelihoodFitDeDxEstimator (const edm::ParameterSet &iConfig)
Public Member Functions inherited from BaseDeDxEstimator
virtual void	beginRun (edm::Run const &run, const edm::EventSetup &iSetup)
virtual	~BaseDeDxEstimator ()

Private Member Functions
void	calculate_wrt_epsilon (const reco::DeDxHit &, const double &, std::array< double, 3 > &)
double	estimate (const reco::DeDxHitCollection &, std::array< double, 2 > &)
void	functionEpsilon (const reco::DeDxHitCollection &, const double &, std::array< double, 3 > &)
double	minimizeAllSaturated (const reco::DeDxHitCollection &, std::array< double, 2 > &)
double	newtonMethodEpsilon (const reco::DeDxHitCollection &, std::array< double, 2 > &)

Detailed Description

Definition at line 7 of file LikelihoodFitDeDxEstimator.h.

Constructor & Destructor Documentation

LikelihoodFitDeDxEstimator()

LikelihoodFitDeDxEstimator::LikelihoodFitDeDxEstimator ( const edm::ParameterSet &  iConfig )

inline

Definition at line 9 of file LikelihoodFitDeDxEstimator.h.

{}

Member Function Documentation

calculate_wrt_epsilon()

void LikelihoodFitDeDxEstimator::calculate_wrt_epsilon ( const reco::DeDxHit &  h, const double &  epsilon, std::array< double, 3 > &  value  )

private

Definition at line 39 of file LikelihoodFitDeDxEstimator.h.

References ALPAKA_ACCELERATOR_NAMESPACE::brokenline::constexpr(), dumpMFGeometry_cfg::delta, PVValHelper::dy, geometryDiff::epsilon, h, isotrackApplyRegressor::k, relval_nano::ls, alignCSCRings::s, mathSSE::sqrt(), contentValuesCheck::ss, heppy_batch::val, and y.

Referenced by functionEpsilon().

                                                                                   {
  const auto& ls = h.pathLength();
  const auto& sn = h.error();      // energy sigma
  const auto y = h.charge() * ls;  // = g * y
  const auto sD = 2.E-3 + 0.095 * y;
  const auto ss = sD * sD + sn * sn;
  const auto s = std::sqrt(ss);
  const auto delta = epsilon * ls;
  const auto dy = delta - y;
  constexpr double nu(0.65);

  // calculate derivatives with respect to delta
  std::array<double, 3> val{{0.}};
  if (h.rawDetId() == 0) {  // normal
    if (dy < -nu * s) {
      val[0] = -2. * nu * dy / s - nu * nu;
      val[1] = -2. * nu / s;
      val[2] = 0.;
    } else {
      val[0] = dy * dy / ss;
      val[1] = 2. * dy / ss;
      val[2] = 2. / ss;
    }
  } else {  // saturated
    if (dy < s) {
      val[0] = -dy / s + 1.;
      val[1] = -1. / s;
      val[2] = 0.;
    } else {
      val[0] = 0.;
      val[1] = 0.;
      val[2] = 0.;
    }
  }

  // d/d delta -> d/d epsilon
  val[1] *= ls;
  val[2] *= ls * ls;

  // sum
  for (size_t k = 0; k < value.size(); k++)
    value[k] += val[k];
}

dedx()

std::pair<float, float> LikelihoodFitDeDxEstimator::dedx ( const reco::DeDxHitCollection &  Hits )

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 11 of file LikelihoodFitDeDxEstimator.h.

References isoTrack_cff::chi2, estimate(), create_idmaps::n, mathSSE::sqrt(), findQualityFiles::v, and relativeConstraints::value.

                                                                         {
    if (Hits.empty())
      return {0., 0.};
    // compute original
    std::array<double, 2> value;
    const auto& chi2 = estimate(Hits, value);
    // try to remove lowest dE/dx measurement
    const auto& n = Hits.size();
    if (n >= 3 && (chi2 > 1.3 * n + 4 * std::sqrt(1.3 * n))) {
      auto hs = Hits;
      hs.erase(std::min_element(hs.begin(), hs.end()));
      // if got better, accept
      std::array<double, 2> v;
      if (estimate(hs, v) < chi2 - 12)
        value = v;
    }
    return {value[0], std::sqrt(value[1])};
  }

estimate()

double LikelihoodFitDeDxEstimator::estimate ( const reco::DeDxHitCollection &  Hits, std::array< double, 2 > &  value  )

private

Definition at line 142 of file LikelihoodFitDeDxEstimator.h.

References h, minimizeAllSaturated(), and newtonMethodEpsilon().

Referenced by dedx().

                                                                                                         {
  // use newton method if at least one hit is not saturated
  for (const auto& h : Hits)
    if (h.rawDetId() == 0)
      return newtonMethodEpsilon(Hits, value);
  // else use minimize all saturated
  return minimizeAllSaturated(Hits, value);
}

functionEpsilon()

void LikelihoodFitDeDxEstimator::functionEpsilon ( const reco::DeDxHitCollection &  Hits, const double &  epsilon, std::array< double, 3 > &  val  )

private

Definition at line 86 of file LikelihoodFitDeDxEstimator.h.

References calculate_wrt_epsilon(), geometryDiff::epsilon, h, and heppy_batch::val.

Referenced by minimizeAllSaturated(), and newtonMethodEpsilon().

                                                                           {
  val = {{0, 0, 0}};
  for (const auto& h : Hits)
    calculate_wrt_epsilon(h, epsilon, val);
}

minimizeAllSaturated()

double LikelihoodFitDeDxEstimator::minimizeAllSaturated ( const reco::DeDxHitCollection &  Hits, std::array< double, 2 > &  value  )

private

Definition at line 95 of file LikelihoodFitDeDxEstimator.h.

References MillePedeFileConverter_cfg::e, functionEpsilon(), and heppy_batch::val.

Referenced by estimate().

                                                                                    {
  int nStep(0);
  double par(3.0);  // input MeV/cm

  std::array<double, 3> val{{0}};
  do {
    functionEpsilon(Hits, par, val);
    if (val[1] != 0)
      par += -val[0] / val[1];
    nStep++;
  } while (val[0] > 1e-3 && val[1] != 0 && nStep < 10);

  value[0] = par * 1.1;
  value[1] = par * par * 0.01;

  return val[0];
}

newtonMethodEpsilon()

double LikelihoodFitDeDxEstimator::newtonMethodEpsilon ( const reco::DeDxHitCollection &  Hits, std::array< double, 2 > &  value  )

private

Definition at line 115 of file LikelihoodFitDeDxEstimator.h.

References funct::abs(), MillePedeFileConverter_cfg::e, functionEpsilon(), and heppy_batch::val.

Referenced by estimate().

                                                                                   {
  int nStep(0);
  double par(3.0);  // input MeV/cm
  double dpar(0);

  std::array<double, 3> val{{0}};
  do {
    functionEpsilon(Hits, par, val);
    if (val[2] != 0.)
      dpar = -val[1] / std::abs(val[2]);
    else
      dpar = 1.;  // step up, for epsilon
    if (par + dpar > 0)
      par += dpar;  // ok
    else
      par /= 2.;  // half
    nStep++;
  } while (std::abs(dpar) > 1e-3 && nStep < 50);

  value[0] = par;
  value[1] = 2. / val[2];

  return val[0];
}