#include<iostream>
#include<cmath>
#include <iomanip>
#include <fstream>
#include <string>
#include <sstream>

using namespace std;

void main()
{
	Int_t nb = 40;
	Emin = 0; 
	Emax = 10;
	TH1D *I_exp = new TH1D("I_exp","Igg experiment",nb,Emin,Emax);
	ifstream if1; 
	double t[9];
	// if1.open("resolve.in");
	// while (if1 >> t[0] >> t[1] >> t[2] >> t[3] >> t[4] >> t[5] >> t[6] >> t[7] >> t[8])
	// {
		// I_exp->SetBinContent(I_exp->GetXaxis()->FindBin(t[0]/1000),I_exp->GetBinContent(I_exp->GetXaxis()->FindBin(t[0]/1000))+t[7]);
		// I_exp->SetBinError(I_exp->GetXaxis()->FindBin(t[0]/1000),I_exp->GetBinError(I_exp->GetXaxis()->FindBin(t[0]/1000))+t[8]);
	// }
	// if1.close();
	
	if1.open("unresolve.in");
	while (if1 >> t[0] >> t[1])
	{
		I_exp->SetBinContent(I_exp->GetXaxis()->FindBin(t[0]/1000),I_exp->GetBinContent(I_exp->GetXaxis()->FindBin(t[0]/1000))+t[1]);
		I_exp->SetBinError(I_exp->GetXaxis()->FindBin(t[0]/1000),I_exp->GetBinError(I_exp->GetXaxis()->FindBin(t[0]/1000))+ sqrt(abs(t[1])));
	}
	if1.close();
	I_exp->Draw("");
	I_exp->Draw("HIST SAME");
	
}