/* Eccentricity vs Impact Parameter b for n = 2 for both Shadowing = 0 and Shadowing = 1 */



void S()
{	
	TFile* fa = new TFile("2_S0.root", "READ");
	TProfile* A = (TProfile*) fa->Get("hprof1");	// for n = 2 for Shadowing = 0
	
	TFile* fb = new TFile("2_S1.root", "READ"); 
	TProfile* B = (TProfile*) fb->Get("hprof2");	// for n = 2 for Shadowing = 1
	
	TFile* fc = new TFile("3_S0.root", "READ");	
	TProfile* C = (TProfile*) fc->Get("hprof3");	// for n = 3 for Shadowing = 0
	
	TFile* fd = new TFile("3_S1.root", "READ");
	TProfile* D = (TProfile*) fd->Get("hprof4");	// for n = 3 for Shadowing = 1
	
	TCanvas *c1 = new TCanvas();
	c1->cd();
	
  	A->SetMarkerColor(kGreen);
  	B->SetMarkerColor(kRed);
  	 	
 
	THStack *hs1 = new THStack("hs","Eccentricity vs b for n = 2 for Shadowing = 0 and 1");
	hs1->Add(A);
	hs1->Add(B);
	hs1->Draw();
	
	TCanvas *c2 = new TCanvas();
	c2->cd();
	
	C->SetMarkerColor(kBlue);
   	D->SetMarkerColor(kBlack);
   	
	THStack *hs2 = new THStack("hs2","Eccentricity vs b for n = 3 for Shadowing = 0 and 1");
	hs2->Add(C);
	hs2->Add(D);
	hs2->Draw();
	
	
	c1->BuildLegend();
	c2->BuildLegend();
	
	/*
	TFile* f2 = new TFile("n2s01.root", "CREATE");
	hprof5->Write();
	TFile* f3 = new TFile("n3s01.root", "CREATE");
	hprof6->Write();*/

}