gausTF1.h File Reference

#include "TF1.h"
#include <string>

Go to the source code of this file.

Functions
double	gausFunction (double *x, double *pars)
TF1 *	gausTF1 (double heightIn=0.0, double meanIn=0.0, double sigmaIn=0.0, std::string fname="dunegaus")

Function Documentation

double gausFunction	(	double *	x,
		double *	pars
	)

Definition at line 18 of file gausTF1.cxx.

                                             {
  double pos = x[0];
  double height = pars[0];
  double mean = pars[1];
  double sigma = pars[2];
  if ( sigma <= 0.0 ) return 0;
  double rat = (pos-mean)/sigma;
  return height*exp(-0.5*rat*rat);
}

TF1* gausTF1	(	double	heightIn = 0.0,
		double	meanIn = 0.0,
		double	sigmaIn = 0.0,
		std::string	fname = "dunegaus"
	)

Definition at line 28 of file gausTF1.cxx.

                                                                              {
  double height = heightIn != 0.0 ? heightIn : 1.0;
  double sigma = sigmaIn > 0.0 ? sigmaIn : 1.0;
  double mean = meanIn;
  string sform = "[0]*exp(-0.5*((x-[1])/[2])*((x-[1])/[2]))";
  TF1* pf = new TF1(fname.c_str(), sform.c_str(), -10, 10);
  //TF1* pf = new TF1(fname.c_str(), gausFunction, -10, 10, 3);
  pf->SetParName(0, "Height");
  pf->SetParName(1, "Mean");
  pf->SetParName(2, "Sigma");
  pf->SetParameter(0, height);
  pf->SetParameter(1, mean);
  pf->SetParameter(2, sigma);
  double sigmaMin = 0.0;
  double sigmaMax = 1.e10;
  if ( sigmaIn > 0.0 ) {
    sigmaMin =  0.1*sigmaIn;
    sigmaMax = 10.0*sigmaIn;
  }
  pf->SetParLimits(2, sigmaMin, sigmaMax);
  if ( sigmaIn > 0.0 ) pf->SetParLimits(2, 0.1*sigma, 10*sigma);
  return pf;
}