#include <GenericTruncatedAverageDeDxEstimator.h>

Inheritance diagram for GenericTruncatedAverageDeDxEstimator:

Public Member Functions
std::pair< float, float >	dedx (const reco::DeDxHitCollection &Hits) override

	GenericTruncatedAverageDeDxEstimator (const edm::ParameterSet &iConfig)

Public Member Functions inherited from BaseDeDxEstimator
virtual void	beginRun (edm::Run const &run, const edm::EventSetup &iSetup)

virtual	~BaseDeDxEstimator ()

Private Attributes
float	m_expo

float	m_fraction

Detailed Description

Definition at line 9 of file GenericTruncatedAverageDeDxEstimator.h.

Constructor & Destructor Documentation

◆ GenericTruncatedAverageDeDxEstimator()

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

inline

Definition at line 11 of file GenericTruncatedAverageDeDxEstimator.h.

                                                                        {
     m_fraction = iConfig.getParameter<double>("fraction");
     m_expo = iConfig.getParameter<double>("exponent");
   }

References edm::ParameterSet::getParameter(), m_expo, and m_fraction.

Member Function Documentation

◆ dedx()

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

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 16 of file GenericTruncatedAverageDeDxEstimator.h.

                                                                          {
     int first = 0, last = Hits.size();
     if (m_fraction > 0) {  // truncate high charge ones
       last -= int(Hits.size() * m_fraction);
     } else {
       first += int(Hits.size() * (-m_fraction));
     }
     double sumdedx = 0;
     for (int i = first; i < last; i++) {
       sumdedx += pow(Hits[i].charge(), m_expo);
     }
     double avrdedx = (last - first) ? pow(sumdedx / (last - first), 1.0 / m_expo) : 0.0;
     return std::make_pair(avrdedx, -1);
   }