void Ratio()
{	
	TFile* fa = new TFile("n2s01.root", "READ");
	THStack* E = (THStack*) fa->Get("hs1");	
	
	THStack* F = (THStack*)E->Clone();
	
	TFile* fb = new TFile("n3s01.root", "READ"); 
	THStack* G = (THStack*) fb->Get("hs2");	
	
	
	gStyle->SetOptStat(0);
	auto c1 = new TCanvas("c1", "Ratio b/w n = 2 and n = 3");
	F->Sumw2();
	//B->Scale(1/1);
	F->SetTitle("Ratio b/w n = 2 and n = 3");
	F->GetXaxis()->SetTitle("Impact parameter b (in fm)");
	F->GetYaxis()->SetTitle("Eccentricity");
	auto rp = new TRatioPlot(F, G);
	c1->SetTicks(0, 1);
	rp->Draw();
	rp->GetLowYaxis()->SetNdivisions(505);
	c1->Update();
	
	
}