#include "StandardLabels.C"
#include "GlueX_boundaries.C"


void cos_theta_inner(void)
{
	gROOT->Reset();
	gStyle->SetPadRightMargin(0.15);
	
	TCanvas *c1 = new TCanvas("c1");
	c1->SetTicks();
	c1->SetGrid();
	
	TFile *f = new TFile("bfield.root");
	TH2D *cos_theta_vs_r_vs_z = (TH2D*)gROOT->FindObject("cos_theta_vs_r_vs_z");

	cos_theta_vs_r_vs_z->GetXaxis()->SetRangeUser(0.0, 410.0);
	cos_theta_vs_r_vs_z->GetYaxis()->SetRangeUser(0.0, 90.0);
	cos_theta_vs_r_vs_z->GetZaxis()->SetRangeUser(-1.0, -0.90);
	cos_theta_vs_r_vs_z->GetZaxis()->SetTitleOffset(1.2);

	cos_theta_vs_r_vs_z->SetTitle("Magnetic Field Direction");
	cos_theta_vs_r_vs_z->SetZTitle("cos(#theta) w.r.t z-direction");
	cos_theta_vs_r_vs_z->Draw("colz");
	DrawGlueXBoundaries(kWhite, cos_theta_vs_r_vs_z->GetYaxis()->GetXmax());
	
	StandardLabels2D(cos_theta_vs_r_vs_z,"solenoid_1500_poisson_20090814_01");
	
	// Save
	c1->SaveAs("cos_theta_inner.pdf");
	c1->SaveAs("cos_theta_inner.gif");
}