//  -------------------------------------------------------------------------
//
//   ----- General Macro for R3B Tracker Analysis
//         Author: Nick Ashwood <n.i.ashwood@bham.ac.uk>
//         Last Update: 01/06/11
//         Comments: Adapted from califaAna.C
//			
//
//  -------------------------------------------------------------------------
//
//   Usage: 
//        > root -l trackAna.C
//        use in /macros/r3b directory                
//
//  -------------------------------------------------------------------------

void trackerAna(Int_t nEvents = 1) {
	
	/*TString dir = getenv("VMCWORKDIR");
	
	TString r3b_geomdir = dir + "/geometry";
	gSystem->Setenv("GEOMPATH",r3b_geomdir.Data());
	
	TString r3b_confdir = dir + "gconfig";
	gSystem->Setenv("CONFIG_DIR",r3b_confdir.Data());
	*/
	
	
	// In general, the following parts need not be touched
	// ========================================================================
	
	// ----    Debug option   -------------------------------------------------
	gDebug = 0;
	// ------------------------------------------------------------------------
	
	// -----   Timer   --------------------------------------------------------
	TStopwatch timer;
	timer.Start();
	// ------------------------------------------------------------------------
	
	
	// -----   Create analysis run   ----------------------------------------
	FairRunAna* fRun = new FairRunAna();
	
	
	TFile* file = new TFile("r3bpar.root");
	file->Get("FairBaseParSet"); 
	
	// -----   Runtime database   ---------------------------------------------
	/*
	FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
	FairParRootFileIo* parIn = new FairParRootFileIo();
	parIn->open("r3bpar.root");
	rtdb->print();
	*/
	// -----  Analysis routines for Si Tracker	
	
	fRun->SetInputFile("r3bsim.root");
	fRun->SetOutputFile("tra_digi.root");
	
	//fRun->LoadGeometry(); //Needed for what??
	
	R3BTraHitFinder* traHF = new R3BTraHitFinder();

	traHF->SetDetectionThreshold(0.000050);   // 50 keV
	traHF->SetExperimentalResolution(0.000050);    // 50 keV 
	fRun->AddTask(traHF);
	
	// Number of events to process
	Int_t nEvents = 50000;

	fRun->Init(); 
	fRun->Run(0, nEvents);
	
	// -----   Finish   -------------------------------------------------------
	timer.Stop();
	Double_t rtime = timer.RealTime();
	Double_t ctime = timer.CpuTime();
	cout << endl << endl;
	cout << "Macro finished succesfully." << endl;
	cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
	cout << endl;
	// ------------------------------------------------------------------------
	
	
	
}