GausRmsFitter Class Reference

#include <GausRmsFitter.h>

Public Types
using	Name = std::string
using	Index = unsigned int

Public Member Functions
	GausRmsFitter (Name fnam)
	GausRmsFitter (double sigma0, double nsigma, Name fnam)
int	fit (TH1 *ph, double mean0) const
void	setLogLevel (Index lev)

Private Attributes
double	m_sigma0
double	m_nsigma
Name	m_fnam
Index	m_LogLevel

Detailed Description

Definition at line 35 of file GausRmsFitter.h.

Member Typedef Documentation

using GausRmsFitter::Index = unsigned int

Definition at line 40 of file GausRmsFitter.h.

using GausRmsFitter::Name = std::string

Definition at line 39 of file GausRmsFitter.h.

Constructor & Destructor Documentation

GausRmsFitter::GausRmsFitter

(

Name

fnam

)

Definition at line 18 of file GausRmsFitter.cxx.

: m_sigma0(0.0), m_nsigma(0.0), m_fnam(fnam), m_LogLevel(0) { }

GausRmsFitter::GausRmsFitter	(	double	sigma0,
		double	nsigma,
		Name	fnam
	)

Definition at line 23 of file GausRmsFitter.cxx.

: m_sigma0(sigma0), m_nsigma(nsigma), m_fnam(fnam), m_LogLevel(0) { }

Member Function Documentation

int GausRmsFitter::fit	(	TH1 *	ph,
		double	mean0
	)		const

Definition at line 28 of file GausRmsFitter.cxx.

                                                  {
  Name myname = "GausRmsFitter::fit: ";
  if ( ph->GetEntries() == 0 ) {
    if ( m_LogLevel >= 2 ) cout << "Histogram has no entries." << endl;
    return 101;  // E.g. empty histo.
  }
  if ( m_LogLevel >= 2 ) cout << myname << "Fitting with sigma " << m_sigma0
                              << " and nsigma " << m_nsigma << endl;
  double mean = mean0;
  double sigma = ph->GetRMS();
  double height = ph->GetMaximum();
  int nbin = ph->GetXaxis()->GetNbins();
  int ibin1Start = ph->GetXaxis()->GetFirst();
  int ibin2Start = ph->GetXaxis()->GetLast();
  double x1 = ph->GetXaxis()->GetBinLowEdge(ibin1Start);
  double x2 = ph->GetXaxis()->GetBinUpEdge(ibin2Start);
  // Evaluate mean and sigma.
  if ( m_sigma0 > 0.0 && m_nsigma > 0.0 ) {
    sigma = 0.0;
    double sigmaOld = m_sigma0;
    if ( m_LogLevel >= 2 ) {
      cout << myname << "Start loop with mean, sigma: " << mean << ", " << sigmaOld << endl;
    }
    const Index maxtry = 100;
    Index itry = 0;
    int ibin1 = ibin1Start;
    int ibin2 = ibin2Start;
    for ( ; itry<maxtry; ++itry ) {
      double dx = m_nsigma*sigmaOld;
      x1 = mean - dx;
      ibin1 = ph->GetXaxis()->FindFixBin(x1);
      if ( ibin1 == 0 ) ibin1 = 1;
      x1 = ph->GetXaxis()->GetBinLowEdge(ibin1);
      x2 = mean + dx;
      ibin2 = ph->GetXaxis()->FindFixBin(x2);
      if ( ibin2 == nbin+1 ) ibin2 = nbin;
      x2 = ph->GetXaxis()->GetBinUpEdge(ibin2);
      ph->GetXaxis()->SetRange(ibin1, ibin2);
      sigma = ph->GetRMS();
      mean = ph->GetMean();
      if ( m_LogLevel >= 2 ) {
        cout << myname << setw(5) << itry << ": (" << x1 << ", " << x2 << "): "
             << mean << ", " << sigma << endl;
      }
      if ( sigma < 1.001*sigmaOld ) {
        break;
      }
      sigmaOld = sigma;
    }
    ph->GetXaxis()->SetRange(ibin1Start, ibin2Start);
    if ( itry >= 100 ) cout << myname << "WARNING: too many iterations." << endl;
  // Use input mean and configured sigma.
  } else if ( m_sigma0 > 0.0 ) {
    sigma = m_sigma0;
  // Use histogram mean and RMS.
  } else {
    mean = ph->GetMean();
  }
  // Protect against sigma == 0.
  if ( sigma <= 0.0 ) {
    sigma = ph->GetXaxis()->GetBinWidth(1)/sqrt(12.0);
    if ( m_LogLevel >= 1 ) cout << myname << "Reset sigma = 0 --> " << sigma << endl;
  }
  TF1* pffix = gausTF1(height, mean, sigma, m_fnam);
  pffix->FixParameter(1, mean);
  pffix->FixParameter(2, sigma);
  pffix->SetRange(x1, x2);
  string fopt = "SR";
  if ( m_LogLevel <= 1 ) fopt += "Q";
  int fstat = quietHistFit(ph, pffix, fopt);
  if ( m_LogLevel >= 1 ) cout << myname << "  status " << fstat << endl;
  if ( fstat == 0 ) {
    if ( m_LogLevel >= 1 ) cout << myname << "Fit suceeded." << endl;
    return 0;
  }
  // Otherwise we may get a crash when we try to view saved copy of histo:
  ph->GetListOfFunctions()->Clear();
  delete pffix;
  if ( m_LogLevel >= 1 ) cout << myname << "Fit failed." << endl;
  return 1;
}

void GausRmsFitter::setLogLevel ( Index  lev )

inline

Definition at line 55 of file GausRmsFitter.h.

{ m_LogLevel = lev; }

Member Data Documentation

Name GausRmsFitter::m_fnam

private

Definition at line 61 of file GausRmsFitter.h.

Index GausRmsFitter::m_LogLevel

private

Definition at line 62 of file GausRmsFitter.h.

double GausRmsFitter::m_nsigma

private

Definition at line 60 of file GausRmsFitter.h.

double GausRmsFitter::m_sigma0

private

Definition at line 59 of file GausRmsFitter.h.

---

The documentation for this class was generated from the following files:

• dunecore/dunecore/DuneCommon/Utility/GausRmsFitter.h
• dunecore/dunecore/DuneCommon/Utility/GausRmsFitter.cxx