#include "StandardLabels.C"

void PiPi_background(void)
{
	TFile *f = new TFile("hd_root.root");
	TH1D *Wplus_minusPiPi_exclusive = (TH1D*)gROOT->FindObject("Wplus_minusPiPi_exclusive");
	TH1D *Wplus_minusPiPi_inclusive = (TH1D*)gROOT->FindObject("Wplus_minusPiPi_inclusive");

	//gStyle->SetPadBottomMargin(0.15);
	//gStyle->SetPadRightMargin(0.2);
	
	Wplus_minusPiPi_exclusive->Rebin(4);
	Wplus_minusPiPi_inclusive->Rebin(4);

	TCanvas *c1 = new TCanvas("c1");
	c1->SetTicks();
	c1->SetGridy();
	
	Wplus_minusPiPi_inclusive->SetLineColor(kRed);

	Wplus_minusPiPi_inclusive->SetLineWidth(3);
	Wplus_minusPiPi_exclusive->SetLineWidth(3);
	
	Wplus_minusPiPi_inclusive->SetTitle("PYTHIA generated #pi^{+}#pi^{-}");
	Wplus_minusPiPi_inclusive->SetXTitle("Invariant mass W of #pi^{+}#pi^{-} (GeV)");
	
	Wplus_minusPiPi_inclusive->Draw();
	Wplus_minusPiPi_exclusive->Draw("same");
	
	TLegend *leg = new TLegend(0.5, 0.7, 0.89, 0.85);
	leg->SetFillColor(kWhite);
	leg->AddEntry(Wplus_minusPiPi_exclusive, "exclusive #pi^{+}#pi^{-}");
	leg->AddEntry(Wplus_minusPiPi_inclusive, "inclusive #pi^{+}#pi^{-}");
	leg->Draw();
	
	StandardLabels1D(Wplus_minusPiPi_inclusive, "0.5^{o} #leq #theta_{#pm} #leq 10^{o}", "#gamma p #rightarrow p #pi^{+} #pi^{-} (?)", "generated events only", "");
	
	c1->SaveAs("PiPi_background.pdf");
	c1->SaveAs("PiPi_background.png");
}