#include "StandardLabels.C"

void eff_vs_p_vs_theta_wb(const char *suffix="")
{
	gROOT->Reset();
	gStyle->SetPalette(1);
	gStyle->SetPadRightMargin(0.16);
	
	TCanvas *c1 = new TCanvas("c1");
	c1->SetTicks();
	
	char fname[256];
	sprintf(fname, "hd_res_charged%s.root", suffix);
	TFile *f = new TFile(fname);
	TH2D *h = (TH2D*)gROOT->FindObject("eff_vs_p_vs_theta_wb");
	string str = string(ParticleName())+" Wire-based efficiency vs. p_{tot} and #theta";
	h->SetTitle(str.c_str());
	h->GetXaxis()->SetRangeUser(0.0, 165.0);
	h->GetZaxis()->SetRangeUser(0.80, 1.0);
	h->GetZaxis()->SetTitleOffset(1.35);
	h->SetXTitle("Polar angle #theta (degrees)");
	h->SetYTitle("Total Momentum (GeV/c)");
	h->SetZTitle("Wire-based tracking Efficiency (#chi^{2}/N_{dof} < 5.0)");
	h->Draw("colz");
	
	StandardLabels2D(h,suffix);
	
	char fname[256];
	sprintf(fname, "eff_vs_p_vs_theta%s_wb.gif", suffix);
	c1->SaveAs(fname);
	sprintf(fname, "eff_vs_p_vs_theta%s_wb.pdf", suffix);
	c1->SaveAs(fname);
}