#include <iostream.h>

//Type .L xft_hotcold.C then .x xft_problemhunt.C("122621.hbk") > ouput.txt
int xft_problemhunt(Char_t* filename="119632.hbk")
{
	//Clear out variables on the stack (avoids much grief)
	gROOT->Reset();
	gROOT->LoadMacro("xft_hotcold.C");
	//Make a canvas, on which plots will flash by
   myc = new TCanvas("myc","Test Canvas",700,700);
	myc->Divide(1,1);
	myc->cd(1);

	//Open file	
	cout << "filename is "<< filename << endl;	
 	TFile* f= new TFile(filename);

	//Look for hot/cold channels 
	//Linker tracks	
	xft_hotcold(xftAbsPhi,12,0,3.); myc->Update();

	//Finder track segments per chip
	xft_hotcold(SL1PixelChip,2); myc->Update();
	xft_hotcold(SL2PixelChip,3); myc->Update();
	xft_hotcold(SL3PixelChip,4); myc->Update();
	xft_hotcold(SL4PixelChip,5); myc->Update();

	//Finder track segments per cell
	xft_hotcold(SL1PixelCell,2,1); myc->Update();
	xft_hotcold(SL2PixelCell,3,1); myc->Update();
	xft_hotcold(SL3PixelCell,4,1); myc->Update();
	xft_hotcold(SL4PixelCell,5,1); myc->Update();

	//XTC prompt/delayed/PD hits per cell
	xft_hotcoldhits(cellsPrompt0,cellsDelayed0, cellsPD0,2,1); myc->Update();
	xft_hotcoldhits(cellsPrompt1,cellsDelayed1, cellsPD1,3,1); myc->Update();
	xft_hotcoldhits(cellsPrompt2,cellsDelayed2, cellsPD2,4,1); myc->Update();
	xft_hotcoldhits(cellsPrompt3,cellsDelayed3, cellsPD3,5,1); myc->Update();
	
	//XTC prompt/delayed/PD hits per wire
	xft_hotcoldhits(wiresPrompt0,wiresDelayed0, wiresPD0,2,2,3.); myc->Update();
	xft_hotcoldhits(wiresPrompt1,wiresDelayed1, wiresPD1,3,2,3.); myc->Update();
	xft_hotcoldhits(wiresPrompt2,wiresDelayed2, wiresPD2,4,2,3.); myc->Update();
	xft_hotcoldhits(wiresPrompt3,wiresDelayed3, wiresPD3,5,2,3.); myc->Update();
	
//	delete f;
	return 0;
}