/*
 * Finds pixel hit pairs and stores them to pair File
 *
 *	Author: Roman Klasen, roklasen@uni-mainz.de or klasen@kph.uni-mainz.de
 *
 */
 
using namespace std;

/*
#include <TROOT.h>
#include <TString.h>
#include <TStopwatch.h>
#include <FairRunAna.h>
#include <FairRuntimeDb.h>
#include <FairParRootFileIo.h>
#include <FairParAsciiFileIo.h>
#include <LmdPairFinderTask.h>
#include <PndSdsGeoPar.h>
*/

void runLumiPixel2ePairFinder(const int nEvents=0, const int startEvent=00000, TString storePath="test/fullrun-1.5/digi", const int verboseLevel=0)
{
	// -----   Timer   --------------------------------------------------------
	TStopwatch timer;
	timer.Start();

	gSystem->Load("libLmd");
	gSystem->Load("libLmdPairFinder");
	gSystem->Load("libLmdSensorAligner");

	cout << "***********************\n";
	cout << "Running PairFinderTask.\n";
	cout << "***********************\n";

	// -----   Input File   ----------------------------------------------------
	TString inFile=storePath+"/Lumi_digi_";
	inFile += startEvent;
	inFile += ".root";

	// -----   Parameter Files   ------------------------------------------------
	TString parFile=storePath+"/Lumi_Params_";
	parFile += startEvent;
	parFile += ".root";
	TString digiparFile = "lumi.digi.par";
	std::string cutParameterfile = storePath + "/cutParameters.json";

	// -----   Output File   --------------------------------------------------
	TString outFile = storePath+"/Lumi_Pairs_";
	outFile += startEvent;
	outFile += ".root";
	std::cout << "DigiFileName: " << outFile.Data() << std::endl;

	// -----   Pair Finder / Cut Finder Runs   -------------------------------------------
	FairRunAna *fRun= new FairRunAna();
	fRun->SetInputFile(inFile);
	fRun->SetOutputFile(outFile);
	fRun->SetEventMeanTime(50);// TODO: 50 ???

	// -----  Parameter database   --------------------------------------------
	FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
	FairParRootFileIo* parInput1 = new FairParRootFileIo(kTRUE);
	parInput1->open(parFile.Data(),"UPDATE");
	rtdb->setFirstInput(parInput1);
	FairParAsciiFileIo* parInput2 = new FairParAsciiFileIo();
	parInput2->open(digiparFile.Data(),"in");
	rtdb->setSecondInput(parInput2);

	// =========================================================================
	// =====                 Start of PairFinder                           =====
	// -----   Actual Task   --------------------------------------------------

	//use dynamic cut parameters
	LmdPairFinderTask* lmdPairFinder = new LmdPairFinderTask();
	lmdPairFinder->useDynamicCut(true, cutParameterfile);
	fRun->AddTask(lmdPairFinder);

	rtdb->setOutput(parInput1);
	rtdb->print();
	// =====                 End of PairFinder                             =====
	// =========================================================================

	// -----   No Idea what this does   ----------------------------------------
	PndSdsGeoPar* geoPar  = (PndSdsGeoPar*)(rtdb->getContainer("PndSdsGeoPar"));

	// -----   Initialize and run PairFinder ----------------------------------
	fRun->Init();
	fRun->Run(0,nEvents);
	rtdb->saveOutput();
	rtdb->print();

	// -----   Finish   -------------------------------------------------------
	timer.Stop();
	Double_t rtime = timer.RealTime();
	Double_t ctime = timer.CpuTime();
	cout << endl << endl;
	cout << "Macro finished successfully." << endl;
	cout << "Real time " << rtime << " s, CPU time " << ctime << " s\a\a" << endl;
	cout << endl;

  return;
}