#include <iostream.h>

int xft_spacer(Char_t* filename="119632.hbk")
{
 
	//Clear out variables on the stack (avoids much grief)
	gROOT->Reset();
	gROOT->LoadMacro("xft_effCalc.C");
	gStyle->SetHistLineWidth(1);	
	gStyle->SetPadGridX(true);
	gStyle->SetPadGridY(true);
	gStyle->SetLabelSize(0.07,"XY");
	gStyle->SetTitleSize(0.08, "XY");
	gStyle->SetTitleOffset(0.3, "XY");
	gStyle->SetFuncColor(6);
	gStyle->SetFuncWidth(6.);
	
	gStyle->SetOptStat(0); //Integral, Overflow, Underflow, RMS, Mean, Nent, Name
	gStyle->SetOptFit(11);     //probability, Chi2, errors, name/values of parameters
	//Open histogram file	
	cout << "filename is "<< filename << endl;	
 	TFile* f= new TFile(filename);

	//Open file in which to store pictures
	TPostScript *ps = new TPostScript("XFT_spacer.ps",111);

	//Open Canvas
	TCanvas *myc = new TCanvas("MyC","Test Canvas",1000,700);

	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,3);

	CotHitsTotZSl3Found->ProfileX("CotHitsTotZSl3Found_px");
	CotHitsAxZSl3Found->ProfileX("CotHitsAxZSl3Found_px");
	CotHitsStZSl3Found->ProfileX("CotHitsStZSl3Found_px");
	CotHitsTotZSl3Missed->ProfileX("CotHitsTotZSl3Missed_px");
	CotHitsAxZSl3Missed->ProfileX("CotHitsAxZSl3Missed_px");
	CotHitsStZSl3Missed->ProfileX("CotHitsStZSl3Missed_px");
	myc->cd(1);
	CotHitsTotZSl3Found_px->Draw();
	CotHitsTotZSl3Found_px->SetXTitle("Z in SL8 (cm)");
	CotHitsTotZSl3Missed_px->SetMarkerColor(2);
	CotHitsTotZSl3Missed_px->SetLineColor(2);
	CotHitsTotZSl3Missed_px->Draw("same");
	leg= new TLegend(.6,0.1,.9,0.5,""); 
	leg->AddEntry(CotHitsTotZSl3Found_px,"Mean hits attached with XFT");
	leg->AddEntry(CotHitsTotZSl3Missed_px,"Mean hits attached without XFT");
	leg->SetFillColor(0);    
   leg->Draw();

	myc->cd(2);
	CotHitsAxZSl3Found_px->Draw();
	CotHitsAxZSl3Found_px->SetXTitle("Z in SL8 (cm)");
	CotHitsAxZSl3Missed_px->SetMarkerColor(2);
	CotHitsAxZSl3Missed_px->SetLineColor(2);
	CotHitsAxZSl3Missed_px->Draw("same");

	myc->cd(3);
	CotHitsStZSl3Found_px->Draw();
	CotHitsStZSl3Found_px->SetXTitle("Z in SL8 (cm)");
	CotHitsStZSl3Missed_px->SetMarkerColor(2);
	CotHitsStZSl3Missed_px->SetLineColor(2);
	CotHitsStZSl3Missed_px->Draw("same");


	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	myc->cd(1);
	foundHitsZSl0->Draw();
	foundHitsZSl0->SetXTitle("Z in SL2 (cm)");

	myc->cd(2);
	foundHitsZSl1->Draw();
	foundHitsZSl1->SetXTitle("Z in SL4 (cm)");

	myc->cd(3);
	foundHitsZSl2->Draw();
	foundHitsZSl2->SetXTitle("Z in SL6 (cm)");

	myc->cd(4);
	foundHitsZSl3->Draw();
	foundHitsZSl3->SetXTitle("Z in SL8 (cm)");

	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	myc->cd(1);
	foundHitsIsoZSl0->Draw();
	foundHitsIsoZSl0->SetXTitle("Z in SL2 (cm)");

	myc->cd(2);
	foundHitsIsoZSl1->Draw();
	foundHitsIsoZSl1->SetXTitle("Z in SL4 (cm)");

	myc->cd(3);
	foundHitsIsoZSl2->SetXTitle("Z in SL6 (cm)");
	foundHitsIsoZSl2->Draw();

	myc->cd(4);
	foundHitsIsoZSl3->Draw();
	foundHitsIsoZSl3->SetXTitle("Z in SL8 (cm)");

	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	myc->cd(1);
	foundHitsIsoUMZSl0->SetXTitle("Z in SL2 (cm)");
	foundHitsIsoUMZSl0->Draw();

	myc->cd(2);
	foundHitsIsoUMZSl1->SetXTitle("Z in SL4 (cm)");
	foundHitsIsoUMZSl1->Draw();

	myc->cd(3);
	foundHitsIsoUMZSl2->SetXTitle("Z in SL6 (cm)");
	foundHitsIsoUMZSl2->Draw();

	myc->cd(4);
	foundHitsIsoUMZSl3->SetXTitle("Z in SL8 (cm)");
	foundHitsIsoUMZSl3->Draw();

	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	foundHitsZSl0->ProfileX("foundHitsZSl0_px");
	foundHitsZSl1->ProfileX("foundHitsZSl1_px");
	foundHitsZSl2->ProfileX("foundHitsZSl2_px");
	foundHitsZSl3->ProfileX("foundHitsZSl3_px");
	myc->cd(1);
	foundHitsZSl0_px->SetXTitle("Z in SL2 (cm)");
	foundHitsZSl0_px->Draw();

	myc->cd(2);
	foundHitsZSl1_px->SetXTitle("Z in SL4 (cm)");
	foundHitsZSl1_px->Draw();

	myc->cd(3);
	foundHitsZSl2_px->SetXTitle("Z in SL6 (cm)");
	foundHitsZSl2_px->Draw();

	myc->cd(4);
	foundHitsZSl3_px->SetXTitle("Z in SL8 (cm)");
	foundHitsZSl3_px->Draw();

	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	foundHitsIsoZSl0->ProfileX("foundHitsIsoZSl0_px");
	foundHitsIsoZSl1->ProfileX("foundHitsIsoZSl1_px");
	foundHitsIsoZSl2->ProfileX("foundHitsIsoZSl2_px");
	foundHitsIsoZSl3->ProfileX("foundHitsIsoZSl3_px");
	myc->cd(1);
	foundHitsIsoZSl0_px->SetXTitle("Z in SL2 (cm)");
	foundHitsIsoZSl0_px->Draw();

	myc->cd(2);
	foundHitsIsoZSl1_px->SetXTitle("Z in SL4 (cm)");
	foundHitsIsoZSl1_px->Draw();

	myc->cd(3);
	foundHitsIsoZSl2_px->SetXTitle("Z in SL6 (cm)");
	foundHitsIsoZSl2_px->Draw();

	myc->cd(4);
	foundHitsIsoZSl3_px->SetXTitle("Z in SL8 (cm)");
	foundHitsIsoZSl3_px->Draw();

	myc->Update();
	char *ch = new char[1]; printf("Hit Enter to continue \n");	gets(ch);
	ps->NewPage();	
	myc->Clear();
 	myc->Divide(1,4);

	foundHitsIsoUMZSl0->ProfileX("foundHitsIsoUMZSl0_px");
	foundHitsIsoUMZSl1->ProfileX("foundHitsIsoUMZSl1_px");
	foundHitsIsoUMZSl2->ProfileX("foundHitsIsoUMZSl2_px");
	foundHitsIsoUMZSl3->ProfileX("foundHitsIsoUMZSl3_px");
	myc->cd(1);
	foundHitsIsoUMZSl0_px->SetXTitle("Z in SL2 (cm)");
	foundHitsIsoUMZSl0_px->Draw();

	myc->cd(2);
	foundHitsIsoUMZSl1_px->SetXTitle("Z in SL4 (cm)");
	foundHitsIsoUMZSl1_px->Draw();

	myc->cd(3);
	foundHitsIsoUMZSl2_px->SetXTitle("Z in SL6 (cm)");
	foundHitsIsoUMZSl2_px->Draw();

	myc->cd(4);
	foundHitsIsoUMZSl3_px->SetXTitle("Z in SL8 (cm)");
	foundHitsIsoUMZSl3_px->Draw();

	myc->Update();
	ps->Close();
	char *ch = new char[1]; printf(" Hit Enter to continue \n");	gets(ch);	



	return 0;	
}