FitUtils.h

#ifndef FitUtils_hh
#define FitUtils_hh 1
 
#include <RtypesCore.h>
#include <TF1.h>
#include <TVector2.h>
 
namespace FitUtils {
 
  class CrystalBallFunctionObject {
    public:
      CrystalBallFunctionObject(Bool_t _isInverted = kFALSE);
      Double_t operator()(double *_x, double *par);
 
 
    private:
      Bool_t isReversed;
  };
 
  TF1* getCrystalBallFunction(TH1 *hist, Bool_t reversed);
 
  TF1* getCrystalBallFunction(TH1 *hist);
 
  TVector2 getCrystalBallMean(TF1* cball);
 
  TVector2 getCrystalBallDispersion(TF1* cball);
 
  TF1* getCrystalBallGaussComponent(TF1* cball);
 
  enum ResolutionType {
    CBFunctionMomenta,    
    CBGaussianPart        
  };
 
  TVector2 getCrystalBallResolution(TF1* cball, ResolutionType rt = CBFunctionMomenta);
 
  TVector2 evalResolution(Double_t mean, Double_t meanErr, Double_t std, Double_t stdErr);
}
 
#endif