#include <iostream>
#include <cmath>
#include <iomanip>
#include <fstream>
#include <string>
#include <sstream>
#include <TH1.h>
#include <TFile.h>
#include <TF1.h>
#include <TCanvas.h>
#include <TRoot.h>
#include <TGraphErrors.h>
#include <TRandom.h>

using namespace std;

	Double_t Emin = 0; 
	Double_t Emax = 10;
	Double_t bin_size = 0.25;
	Int_t nb = (Emax-Emin)/bin_size;
	
	Double_t Bn = 8.020;
	Double_t E_thres = 0.52; 
	
	Double_t B = 1;
	Double_t Bcal = 0.0;
	
	
	TH1D *LDT = new TH1D("LDT","Total Level density",nb,Emin,Emax);
	TH1D *LDJ = new TH1D("LDJ","0<J<2 Level density",nb,Emin,Emax);
	TH1D *RSF = new TH1D("RSF","Partial gamma width via E gamma",nb,Emin,Emax);
	TH1D *rsf = new TH1D("rsf","Gamma strength function",nb,Emin,Emax);
	TH1D *gRSF = new TH1D("gRSF","Total Gamma",nb,Emin,Emax);
	TH1D *mRSF = new TH1D("mRSF","Partial gamma width for intermediate",nb,Emin,Emax);
	
	TH1D *SD = new TH1D("SD","spin distribution",nb,Emin,Emax);
	TH1D *Igg = new TH1D("Igg","Cascade intensity",nb,Emin,Emax);
	TH1D *IggC = new TH1D("IggC","Calculated Cascade intensity",nb,Emin,Emax);
	
	TH1D *LD_RIPL3 = new TH1D("LD_RIPL3","Counting LD", nb, Emin, Emax);
	TH1D *LD_RIPL3_LS = new TH1D("LD_RIPL3_LS","Counting LD J<2", nb, Emin, Emax);
	
	TH1D *tmpld = new TH1D("tmpld","tmp ld",nb,Emin,Emax);
	TH1D *tmprsf = new TH1D("tmprsf","tmp rsf",nb,Emin,Emax);
	
	TH1D *addLD = new TH1D("addLD","adding ld logarithm",nb,Emin,Emax);
	TH1D *addrsf = new TH1D("addrsf","adding rsf logarithm",nb,Emin,Emax);
	
	TH1D *simLD = new TH1D("simLD","LD simulation", nb, Emin, Emax);
	TH1D *simrsf = new TH1D("simrsf","RSF simulation", nb, Emin, Emax);

// void sim_igg()
// {
	// for (int i = 1; i<=nb; i++)
	// {
		// Double_t E = simLD->GetXaxis()->GetBinCenter(i);
		// simrsf->SetBinContent(i, pow(10.,-9)*2.5*E);
		// simLD->SetBinContent(i,exp(E*1.2));
		// if (E<=2) simLD->SetBinContent(i,10.);
	// }
	// for (int i = 1; i<=nb; i++)
	// {
		// LDJ->SetBinContent(i,simLD->GetBinContent(i));
		// rsf->SetBinContent(i,simrsf->GetBinContent(i));
	// }
	// TFile f("kq.root","new");
	// f.cd();
	// LDJ->Write();
	// f.Close();
	// //f.Delete();
// }






void main()
{
	//sim_igg();
	
	TFile f("kq.root");
	f.ls();
	LDJ = (TH1D*)f.Get("LDJ");
	LDJ->Draw();
	f.Close();
	
	
}