//--------------------------
// DrawGlueXBoundaries
//--------------------------
void DrawGlueXBoundaries(int color=kBlack, double clip_tof=1000.0)
{
	DrawTargetBoundaries(color);
	DrawStartCounterBoundaries(color);
	DrawFDCBoundaries(color);
	DrawCDCBoundaries(color);
	DrawBCALBoundaries(color);
	DrawTOFBoundaries(color, clip_tof);

}

//--------------------------
// DrawTargetBoundaries
//--------------------------
void DrawTargetBoundaries(int color)
{
	double Zlo = 50.0;
	double Zhi = Zlo + 30.0;
	double Rlo = 0.0;
	double Rhi = Rlo + 1.5;
	
	TBox *box = new TBox(Zlo, Rlo, Zhi, Rhi);
	box->SetLineWidth(2.0);
	box->SetFillColor(0);
	box->SetFillStyle(0);
	box->SetLineColor(color);
	box->Draw();

	TLatex *lab = new TLatex(Zlo-5, Rhi-1, "target");
	lab->SetTextAlign(31);
	lab->SetTextSize(0.02);
	lab->SetTextColor(color);
	lab->Draw();
}

//--------------------------
// DrawStartCounterBoundaries
//--------------------------
void DrawStartCounterBoundaries(int color)
{
	double Zlo = 50.0;
	double Zhi = Zlo + 30.0;
	double Rlo = 0.0;
	double Rhi = Rlo + 1.5;
	
	// Values from Richard's old spreadsheet "start_geom.xls"
	int Npoints = 15;
	double r[] = {8.675, 8.675, 8.675, 7.759, 7.354, 4.353, 3.134,
	              2.290, 1.915, 5.593, 5.893, 6.907, 6.951, 6.951, 8.675};
	double z[] = {   0.0, 51.383, 51.726, 53.938, 54.914, 57.915, 59.134,
	              58.442, 58.290, 54.360, 54.238, 51.531, 51.488,    0.0, 0.0};
	
	// shift z to proper location in lab system
	for(int i=0; i<Npoints; i++)z[i] += 38.75;
	
	TPolyLine *pol = new TPolyLine(Npoints, z, r);
	pol->SetLineWidth(2.0);
	pol->SetFillColor(0);
	pol->SetFillStyle(0);
	pol->SetLineColor(color);
	pol->Draw();

	TLatex *lab = new TLatex(z[1]-1, r[1]-2, "start counter");
	lab->SetTextAlign(33);
	lab->SetTextSize(0.02);
	lab->SetTextColor(color);
	lab->Draw();
}

//--------------------------
// DrawFDCBoundaries
//--------------------------
void DrawFDCBoundaries(int color)
{
	double Rlo = 2.0; Rhi=48.5;

	double Zlo[7], Zhi[7];

	Zlo[0] = 176.1586;
	Zhi[0] = 187.1614;
	
	Zlo[1] = 233.7186;
	Zhi[1] = 244.7214;

	Zlo[2] = 291.2786;
	Zhi[2] = 302.2814;

	Zlo[3] = 348.8386;
	Zhi[3] = 359.8414;

	for(int sl=0; sl<4; sl++){
	
		TBox *box = new TBox(Zlo[sl], Rlo, Zhi[sl], Rhi);
		box->SetLineWidth(2.0);
		box->SetFillColor(0);
		box->SetFillStyle(0);
		box->SetLineColor(color);
		box->Draw();

		TLatex *lab = new TLatex((Zlo[sl]+Zhi[sl])/2.0+2.0, Rhi+0.5, "FDC");
		lab->SetTextAlign(21);
		lab->SetTextSize(0.025);
		lab->SetTextColor(color);
		lab->Draw();
	}	
}

//--------------------------
// DrawCDCBoundaries
//--------------------------
void DrawCDCBoundaries(int color)
{
	double Zlo=17.0; Zhi=Zlo+150.0;
	double Rlo[7], Rhi[7];

	Rlo[0] = 10.7219;
	Rhi[0] = 15.1621;
	
	Rlo[1] = 16.9321;
	Rhi[1] = sqrt(pow(21.8912,2.0) + pow(0.860106,2.0));

	Rlo[2] = 23.8544;
	Rhi[2] = sqrt(pow(28.5658,2.0) + pow(0.846871,2.0));

	Rlo[3] = 31.3799;
	Rhi[3] = 35.8301;

	Rlo[4] = 37.4446;
	Rhi[4] = sqrt(pow(41.9225,2.0) + pow(0.833676,2.0));

	Rlo[5] = 43.6152;
	Rhi[5] = sqrt(pow(48.0733,2.0) + pow(0.829899,2.0));

	Rlo[6] = 50.3747;
	Rhi[6] = 54.7617;
	
	TBox *box = new TBox(Zlo, Rlo[0], Zhi, Rhi[6]);
	box->SetLineWidth(2.0);
	box->SetFillColor(0);
	box->SetFillStyle(0);
	box->SetLineColor(color);
	box->Draw();
	
	TLatex *lab = new TLatex(Zlo+5, Rhi[6]-1, "CDC");
	lab->SetTextAlign(13);
	lab->SetTextSize(0.03);
	lab->SetTextColor(color);
	lab->Draw();

}

//--------------------------
// DrawTOFBoundaries
//--------------------------
void DrawTOFBoundaries(int color, double clip_tof=1000.0)
{
	double Zlo = 616.25;
	double Zhi = Zlo + 2.0*2.54;
	double Rlo = 3.0;
	double Rhi = 252.0/2.0;
	
	// Clip TOF at top of histogram area
	if(Rhi>clip_tof)Rhi = clip_tof;
	
	TBox *box = new TBox(Zlo, Rlo, Zhi, Rhi);
	box->SetLineWidth(2.0);
	box->SetFillColor(0);
	box->SetFillStyle(0);
	box->SetLineColor(color);
	box->Draw();

	TLatex *lab = new TLatex(Zlo-5, Rhi-1, "TOF");
	lab->SetTextAlign(31);
	lab->SetTextSize(0.03);
	lab->SetTextAngle(90.0);
	lab->SetTextColor(color);
	lab->Draw();
}

//--------------------------
// DrawBCALBoundaries
//--------------------------
void DrawBCALBoundaries(int color)
{
	double Zlo = 17.0;
	double Zhi = Zlo + 390.0;
	double Rlo = 64.2; // includes aluminum plate
	double Rhi = 90.0;
	
	TBox *box = new TBox(Zlo, Rlo, Zhi, Rhi);
	box->SetLineWidth(2.0);
	box->SetFillColor(0);
	box->SetFillStyle(0);
	box->SetLineColor(color);
	box->Draw();
	
	TLatex *lab = new TLatex(Zlo+5, Rhi-1, "BCAL");
	lab->SetTextAlign(13);
	lab->SetTextSize(0.03);
	lab->SetTextColor(color);
	lab->Draw();
}