figSpectraIZV.cxx

#if !defined(__CLING__)

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

#endif

void figSpectraIZV()
{

//	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->ToggleEditor();
//	c1->SetWindowSize(800, 1200);
//	c1->SetWindowSize(800, 800);
	c1->SetWindowSize(1333, 787);
	c1->Divide(2, 2, .001, .001);


	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);
	}

//	c1->SaveAs("figures/spectra.eps");

}