#include <TSystem.h>
#include <TRandom3.h>
#include <TCanvas.h>
#include <TF1.h>

//double Binomial(double k, double mu, double n);

void EnergyProductionFit() {
	gSystem->Load("libMathCore");
	


	
	//TF1 * NBD = new TF1("NBD", "(ROOT::Math::tgamma(x+[0])/(ROOT::Math::tgamma(x+1)*ROOT::Math::tgamma([0])))*(pow([1]/[0], x)/pow([1]/([0]+1), x+[0]))", 0,5);

	TF1 * NBD = new TF1("NBD", "[0]*ROOT::Math::negative_binomial_pdf([1],[2], x)", 0,5);

	NBD->SetParameter(0, 1);
	NBD->SetParameter(1, 2);
	NBD->SetParameter(2, 0.57);

	TCanvas * c1 = new TCanvas("c1", "Match Fit", 800, 800);
	NBD->Draw("");
}

/*double Binomial(double k, double mu, double n) {
namespace ROOT {
namespace Math {
	double a = (tgamma(n+k)/(tgamma(n+1)*tgamma(k)))*(pow(mu/k, n)/pow(mu/(k+1), n+k));
}
}
return a;
}   */