16 arg = (x[0]-par[1])/par[2];
18 return par[0]*TMath::Exp(-0.5*arg*arg);
24 funcType_ (funcType ),
25 funcName_ (funcName ),
26 lowerBound_(lowerBound),
27 upperBound_(upperBound)
31 func_->SetParNames(
"Const",
"Mean",
"Sigma" );
34 std::cout <<
"Sorry: not yet implemented" << std::endl;
43 double maxValue=hist.GetBinCenter(hist.GetMaximumBin());
44 func_->SetParameter(1, maxValue);
45 func_->SetParameter(2, hist.GetRMS());
51 func_->SetParLimits(2, 0., 5.*hist.GetRMS());
59 if(hist.GetFunction(
"func")){
62 mean_ = hist.GetFunction(
"func")->GetParameter(1);
63 sigma_= hist.GetFunction(
"func")->GetParameter(2);
72 std::cout <<
"sorry... no fit function found..." << std::endl;
81 int idx=hist.GetMaximumBin(), jdx=hist.GetMaximumBin();
82 if(0<=frac && frac<=1){
83 while( hist.GetBinContent(idx)/hist.GetMaximum()>
frac) --idx;
84 while( hist.GetBinContent(jdx)/hist.GetMaximum()>
frac) ++jdx;
86 return std::pair<int, int>(idx, jdx);
96 return std::fabs(outer-inner)/2;
std::pair< int, int > contour(TH1F &, double &)
double distance(TH1F &hist)
double Gauss(double *x, double *par)
double spreadError(TH1F &hist)
const double par[8 *NPar][4]
void quantiles(TH1F &hist, double err)