Public Types | Public Member Functions | Private Attributes | List of all members
GausStepFitter Class Reference

#include <GausStepFitter.h>

Public Types

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

Public Member Functions

 GausStepFitter (double pos, double sigma, double height, Name fnam, Name fopt="WWS")
 
int fit (TH1 *ph) const
 

Private Attributes

double m_pos
 
double m_sigma
 
double m_height
 
Name m_fnam
 
Name m_fopt
 
Index m_LogLevel
 

Detailed Description

Definition at line 15 of file GausStepFitter.h.

Member Typedef Documentation

using GausStepFitter::Index = unsigned int

Definition at line 20 of file GausStepFitter.h.

using GausStepFitter::Name = std::string

Definition at line 19 of file GausStepFitter.h.

Constructor & Destructor Documentation

GausStepFitter::GausStepFitter ( double  pos,
double  sigma,
double  height,
Name  fnam,
Name  fopt = "WWS" 
)

Definition at line 16 of file GausStepFitter.cxx.

17 : m_pos(pos), m_sigma(sigma), m_height(height), m_fnam(fnam), m_fopt(fopt),
18  m_LogLevel(0) { }

Member Function Documentation

int GausStepFitter::fit ( TH1 *  ph) const

Definition at line 22 of file GausStepFitter.cxx.

22  {
23  Name myname = "GausStepFitter::fit: ";
24  double sigma = m_sigma;
25  double height = m_height;
26  double sigfac = 2.0;
27  for ( int ifit=0; ifit<5; ++ifit ) {
28  TF1* pffix = gausTF1(height, m_pos, sigma, m_fnam);
29  double sigmax = 1.1*sigfac*sigma;
30  double sigmin = 0.9*sigma/sigfac;
31  if ( m_LogLevel >= 1 ) cout << myname << " Doing constrained fit " << ifit
32  << " with pos=" << m_pos << ", sigma=" << sigma
33  << " (" << sigmin << ", " << sigmax << ")" << endl;
34  pffix->SetParameter(2, sigma);
35  pffix->SetParLimits(2, sigmin, sigmax);
36  if ( m_height == 0.0 ) {
37  pffix->SetParLimits(0, 0.0, 1.e10); // Don't let height go negative.
38  } else if ( height > 0.0 ) {
39  pffix->SetParLimits(0, 0.1*height, 2.0*height); // Don't let height go negative.
40  } else if ( height < 0.0 ) {
41  pffix->SetParLimits(0, 2.0*height, 0.1*height); // Don't let height go positive.
42  }
43  Name fopt = m_fopt;
44  if ( m_LogLevel >= 2 ) fopt += "Q";
45  int fstat = quietHistFit(ph, pffix, fopt.c_str());
46  double signew = pffix->GetParameter(2);
47  bool atHiLimit = signew > 0.999*sigmax;
48  bool atLoLimit = signew < 1.001*sigmin;
49  if ( m_LogLevel >= 1 ) cout << myname << " status " << fstat << ", fit sigma=" << signew << endl;
50  if ( fstat == 0 && !atHiLimit && !atLoLimit ) {
51  if ( m_LogLevel >= 1 ) cout << myname << "Fit suceeded." << endl;
52  return 0;
53  }
54  // Otherwise we may get a crash when we try to view saved copy of histo:
55  ph->GetListOfFunctions()->Clear();
56  delete pffix;
57  if ( atLoLimit ) sigma /= sigfac;
58  else sigma *= sigfac;
59  }
60  if ( m_LogLevel >= 1 ) cout << myname << "Fit failed." << endl;
61  return 1;
62 }
Name
ChannelGroupService::Name Name
Definition: FixedChannelGroupService_service.cc:19
quietHistFit
int quietHistFit(TH1 *ph, std::string fname, std::string fopt)
Definition: quietHistFit.cxx:25
gausTF1
TF1 * gausTF1(double heightIn, double meanIn, double sigmaIn, std::string fname)
Definition: gausTF1.cxx:28
endl
QTextStream & endl(QTextStream &s)
Definition: qtextstream.cpp:2030

Member Data Documentation

Name GausStepFitter::m_fnam
private

Definition at line 39 of file GausStepFitter.h.

Name GausStepFitter::m_fopt
private

Definition at line 40 of file GausStepFitter.h.

double GausStepFitter::m_height
private

Definition at line 38 of file GausStepFitter.h.

Index GausStepFitter::m_LogLevel
private

Definition at line 41 of file GausStepFitter.h.

double GausStepFitter::m_pos
private

Definition at line 36 of file GausStepFitter.h.

double GausStepFitter::m_sigma
private

Definition at line 37 of file GausStepFitter.h.

The documentation for this class was generated from the following files:
Generated by   doxygen 1.8.11