#include "StandardLabels.C"
#include "GlueX_boundaries.C"
#include <TH2D.h>
#include <TCanvas.h>

void MagnetDetector()
{
	TColor::CreateColorWheel();

	TCanvas *c1 = new TCanvas("boundaries canvas");
	c1->SetTicks();

	Int_t zmin = -100;
	Int_t zmax=700, rmax=200;
	TH2D *boundaries = new TH2D("boundaries", "GlueX Detectors and Magnet",
                              zmax-zmin, zmin, zmax, rmax*2, -rmax, rmax);
	boundaries->SetStats(0);
	boundaries->SetXTitle("Z (cm)");
	boundaries->SetYTitle("X (cm)");
	
	boundaries->Draw();

	TFile *f = new TFile("bfield.root");
	TH2D *Btot_vs_x_vs_z = (TH2D*)gROOT->FindObject("Btot_vs_x_vs_z");

	Btot_vs_x_vs_z->GetZaxis()->SetRangeUser(0.0, 2.5);
	Btot_vs_x_vs_z->Draw("colz same");

	DrawGlueXBoundaries(true, true);
	DrawGlueXBoundaries(true, false, 100000, kBlack);
	StandardLabels(boundaries,"","","solenoid_1500_poisson_20090814_01");
	
	c1->SaveAs("MagnetDetector.png");
	c1->SaveAs("MagnetDetector.pdf");
}