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


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

	Btot_vs_r_vs_z->GetZaxis()->SetRangeUser(0.0, 2.8);

	Btot_vs_r_vs_z->SetTitle("Magnetic Field Magnitude");
	Btot_vs_r_vs_z->SetZTitle("B-field magnitude (Tesla)");
	Btot_vs_r_vs_z->Draw("colz");
	DrawGlueXBoundaries();
	
	StandardLabels2D(Btot_vs_r_vs_z,"solenoid_1500_poisson_20090814_01");
	
	// Save
	c1->SaveAs("Bmag.pdf");
	c1->SaveAs("Bmag.gif");
}