#ifndef FIT_H
#define FIT_H
#include <TMath.h>
#include <Math/PdfFuncMathCore.h>
#include <Math/ProbFuncMathCore.h>



double model(double *x,double *par){
        //Cruijff
        double arg = 0;
        double arg2 = (x[0] - par[2]);
        if (par[0] != 0 && par[1] != 0){
                if( arg2 <= 0){
                        arg = par[5]*TMath::Exp(-TMath::Power(arg2,2) / (2*par[0]*par[0] + par[3] *TMath::Power(arg2,2)));
                }
                else{
                        arg = par[5]*TMath::Exp(-TMath::Power(arg2,2) / (2*par[1]*par[1] + par[4] *TMath::Power(arg2,2)));
                }
        }
        //4body
        double xx = x[0] - par[9];
        double arg1 = 0;
        if(xx > 0 ){
                arg1 = xx*TMath::Power(1 - TMath::Power(xx/par[6],2), par[8]) * TMath::Exp(par[7] * (1 - TMath::Power(xx/par[6],2)));
        }
        //exponential
        double arg3 = 0;
        arg3 = par[12]*par[11]*TMath::Exp(-par[11]*xx);

        return (arg + par[10]*arg1 +arg3);
}

#endif