#include "StandardLabels.C"


void dphi_vs_p_vs_theta(const char *suffix="")
{
	gROOT->Reset();
	gStyle->SetPalette(1);
	gStyle->SetPadRightMargin(0.15);
	
	TCanvas *c1 = new TCanvas("c1");
	c1->SetTicks();
	c1->SetLogz();
	
	char fname[256];
	sprintf(fname, "hd_res_charged%s.root", suffix);
	TFile *f = new TFile(fname);
	TH2D *h = (TH2D*)gROOT->FindObject("dphi_sigma");
	h->SetStats(0);
	h->SetTitle("#phi angle resolution vs. p_{tot} and #theta");
	h->GetXaxis()->SetRangeUser(0.0, 170.0);
	h->GetZaxis()->SetRangeUser(1.0, 100.0);
	h->SetZTitle("Azimuthal angle #phi resolution (mrad)");
	h->Draw("colz");
	
	StandardLabels2D(h, suffix);
	
	char fname[256];
	sprintf(fname, "dphi_vs_p_vs_theta%s.gif", suffix);
	c1->SaveAs(fname);
	sprintf(fname, "dphi_vs_p_vs_theta%s.pdf", suffix);
	c1->SaveAs(fname);
}