#include <UnbinnedFitDeDxEstimator.h>

Inheritance diagram for UnbinnedFitDeDxEstimator:

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

	UnbinnedFitDeDxEstimator (const edm::ParameterSet &iConfig)

	~UnbinnedFitDeDxEstimator () override

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

virtual	~BaseDeDxEstimator ()

Public Attributes
double	data [50]

TF1 *	f1

UnbinnedLikelihoodFit	fitter

Detailed Description

Definition at line 14 of file UnbinnedFitDeDxEstimator.h.

Constructor & Destructor Documentation

◆ UnbinnedFitDeDxEstimator()

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

inline

Definition at line 16 of file UnbinnedFitDeDxEstimator.h.

References f1, fitter, and UnbinnedLikelihoodFit::setFunction().

                                                            {
     fitter.setFunction((f1 = new TF1("myLandau", "TMath::Landau(x,[0],[1],1)", 0, 255)));
   }

◆ ~UnbinnedFitDeDxEstimator()

UnbinnedFitDeDxEstimator::~UnbinnedFitDeDxEstimator ( )

inlineoverride

Definition at line 20 of file UnbinnedFitDeDxEstimator.h.

References f1.

                                        {
     // clean up everything
     delete f1;
   }

Member Function Documentation

◆ dedx()

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

inlineoverridevirtual

Implements BaseDeDxEstimator.

Definition at line 25 of file UnbinnedFitDeDxEstimator.h.

References data, f1, UnbinnedLikelihoodFit::fit(), fitter, and mps_fire::i.

                                                                          {
     // if there is no hit, returns invalid.
     if (Hits.empty())
       return std::make_pair(-1, -1);
     // sets the initial parameters
     f1->SetParameters(3.0, 0.3);
     // fills a temporary array and performs the fit
     uint32_t i = 0;
     for (reco::DeDxHitCollection::const_iterator hit = Hits.begin(); hit != Hits.end(); ++hit, ++i) {
       data[i] = hit->charge();
     }
     // fit !
     fitter.fit(Hits.size(), data);
     // returns the mpv and its error
     return std::make_pair(f1->GetParameter(0), f1->GetParError(0));
   }