void Test()
{
TH1F * h[15];

for (int i=0;i<sizeof(h)/sizeof(*h);i++)
  h[i]=new TH1F();

auto canvas = new TCanvas("canvas","Fitting",4000,2000);
canvas->Divide(5,3);
TFile f("resolution_plots_HCALOUT.root");
TTree * t;
f.GetObject("t",t);
Float_t  e,del_e;
t->SetBranchAddress("truth_e",&e);
t->SetBranchAddress("delta_e",&del_e);

vector<vector<float> > M;
vector<float> m;
for (int i=0;i<15;i++)
{
for (int j=0;j<3130;j++)
{
t->GetEntry(i*j);
m.push_back(e);
h[j]->Fill(del_e);
}
M.push_back(m);
}
Double_t energy[15],mean[15],sigma[15];
for (int i=0;i<15;i++)
{
canvas->cd(i);
h[i]->Fit("gaus");
TF1 *gaus= h[i]->GetFunction("gaus");
mean[i]=gaus->GetParameter("Mean");
sigma[i]=gaus->GetParameter("Sigma");
delete gaus;
//energy[i]=(*std::max_element(M[i],M[i]+3130)+*std::min_element(M[i],M[i]+3130))/2;
}
std::ofstream mean_data("e vs mean.txt");
std::ofstream sigma_data("e vs sigma.txt");
if(mean_data.is_open())
{
	for(int r=0;r<15;r++)
{
	mean_data<< energy[r]<< "         ";
	mean_data<<mean[r]<< endl;
}
mean_data.close();
}
if(sigma_data.is_open())
{
	for(int r=0;r<15;r++)
{
	sigma_data<< energy[r]<< "         ";
	sigma_data<<sigma[r]<< endl;
}
sigma_data.close();
}
}