#include "TCanvas.h"
#include "TF1.h"




using namespace std;


void test(){


c1 = new TCanvas("c1","Double Beta Decay Distributions",1000,1000);


TF1 *f = new TF1("f","1/10423.1*x*(4.8-x)**5*(1.+2.*x+4./3.*x**2+1./3.*x**3+1./30.*x**4)",0,10./511.707);
f->SetLineColor(kBlue);
f->GetHistogram()->GetXaxis()->SetTitle("E [units of e mass]");
c1->Modified();

 	
TF1 *f1 = new TF1("f1", "63.9532*x*(2447.8-x)**5*((7.22776e-23)+(1.54270e-24)*x+(8.93146e-27)*x**2+(-2.14760e-28)*x**3+(-6.35217e-31)*x**4+(7.34418e-33)*x**5+(2.89105e-35)*x**6+(-5.25903e-38) *x**7+(-7.51820e-39)*x**8+(3.49659e-41)*x**9+(1.24854e-42)*x**10+(5.29055e-44)*x**11+(-1.98767e-46)*x**12+(-4.11338e-48)*x**13)", 0, 10.);
f1->Draw();
f->Draw("same");
f1->GetHistogram()->GetXaxis()->SetTitle("E [keV]");
c1->Modified();
	


}