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