#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	expo_

float	fraction_

bool	truncate_

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.

References expo_, fraction_, edm::ParameterSet::getParameter(), and truncate_.

                                                                        {
     fraction_ = iConfig.getParameter<double>("fraction");
     expo_ = iConfig.getParameter<double>("exponent");
     truncate_ = iConfig.getParameter<bool>("truncate");
   }

Member Function Documentation

◆ dedx()

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

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 17 of file GenericTruncatedAverageDeDxEstimator.h.

References ALCARECOTkAlJpsiMuMu_cff::charge, expo_, dqmdumpme::first, fraction_, mps_fire::i, createfilelist::int, dqmdumpme::last, conifer::pow(), and truncate_.

Referenced by DeDxHitInfoProducer::produce().

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