//Figure 1 IZV RAS

#if !defined(__CLING__)

#include "TFile.h"
#include "TCanvas.h"
#include "TH1F.h"
#include "TText.h"

#endif

void figSpectraIZV(const Bool_t saveFigs = 0)
{

//	TFile *fr = new TFile("figures/spectra.root", "READ");
//	TFile *fr = new TFile("figures/spectra100new.root", "READ");
	TFile *fr = new TFile("figures/spectra50new.root", "READ");

//	fr->ls();

	TH1F *hExp[4];
	TH1F *hSim[4];

	TH1F *hCurr = 0;

	TString heName;
	TString hsName;

	hCurr = (TH1F*)fr->Get("heSpectrum4560");
	hExp[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("heSpectrum6075");
	hExp[1] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("heSpectrum7590");
	hExp[2] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("heSpectrum90120");
	hExp[3] = new TH1F(*hCurr);


	hCurr = (TH1F*)fr->Get("hsSpectrum4560");
	hSim[0] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("hsSpectrum6075");
	hSim[1] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("hsSpectrum7590");
	hSim[2] = new TH1F(*hCurr);
	hCurr = (TH1F*)fr->Get("hsSpectrum90120");
	hSim[3] = new TH1F(*hCurr);



	TCanvas *c1 = new TCanvas("c1", "Fig. 1, IZV RAS");
	c1->ToggleEditor();
//	c1->SetWindowSize(800, 1200);
//	c1->SetWindowSize(800, 800);
	c1->SetWindowSize(1333, 787);
	c1->Divide(2, 2, .001, .001);

	TText *t1[4];
	TString panel[4] = {"a)", "b)", "c)", "d)"};


	for (Int_t i = 0; i < 4; i++) {
		c1->cd(i+1);

		hExp[i]->SetTitle("");
		hSim[i]->SetTitle("");

		hSim[i]->SetXTitle("#it{E}, MeV");
		hSim[i]->Smooth();
		hExp[i]->SetLineWidth(2);
		hSim[i]->SetLineColor(kGray+1);



		hSim[i]->Draw("");
		hExp[i]->Draw("E same");

		if (i==0) hSim[i]->GetYaxis()->SetRangeUser(0,13000);

		t1[i] = new TText(0.18,0.87,panel[i].Data());
		t1[i]->SetNDC();
		t1[i]->SetTextFont(132);
		t1[i]->SetTextSize(0.07);
		t1[i]->Draw();
	}

	if (saveFigs) {
//		c2->SaveAs("figures/PresGSepsilonTright7590.png");

		c1->SaveAs("figures/Fig1izvSpectra.png");
		c1->SaveAs("figures/Fig1izvSpectra.eps");
	}

	return;
}