// *******
// Macro for running analysis in a Task
// Take a look to tutorials/anatask/PndAnalysis.cxx
// Add the line add_subdirectory(tutorials/anatask)
// to pandaroot/CMakeLists.txt
// and do a 'make install' in build directory first
// *******


void makeTCands_stt(int nevts=10)
{
	;
  TStopwatch timer;
  timer.Start();
  gDebug=0;

  // Load basic libraries
  gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
  gSystem->Load("libAnalysisTools");
  gSystem->Load("libRho");

	TString parFile = "data/params_sttcombi.root";
	TString inSimuFile = "data/points_sttcombi.root";
	TString inDigiFile = "data/digi_sttcombi.root";
	TString inRecoFile = "data/reco_sttcombi.root";
	TString inKalmFile = "data/fit_sttcombi.root";
	TString inLheFile = "data/lhepidmaker.root";
	TString outFile = "data/tcands.root";
	
	FairRunAna *fRun= new FairRunAna();
	fRun->SetInputFile(inSimuFile);
	fRun->AddFriend(inDigiFile);
	fRun->AddFriend(inRecoFile);
	fRun->AddFriend(inKalmFile);
	//fRun->AddFriend(inLheFile);

	fRun->SetOutputFile(outFile);
  FairGeane *Geane = new FairGeane(inSimuFile);
	PndEmcMapper *emcMap = PndEmcMapper::Instance(2,inSimuFile);

	// -----  Parameter database   --------------------------------------------
	FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
	FairParRootFileIo* parInput1 = new FairParRootFileIo();
	parInput1->open(parFile.Data());
	rtdb->setFirstInput(parInput1);
	fRun->LoadGeometry();
	


  //append the analysis task
  PndMicroWriter* miwri = new PndMicroWriter("LheGenTrack");
  fRun->AddTask(miwri);
//  PndStoreTCands *tcandTask=new PndStoreTCands();
//  fRunA->AddTask(tcandTask);

  //PndSimpleAnalysis *anaTask=new PndSimpleAnalysis(cfgfile.Data());
  //fRunA->AddTask(anaTask);

  fRun->Init();
//  if (nevts==0) nevts=10;
  fRun->Run(0,nevts);


  timer.Stop();
  Double_t rtime = timer.RealTime();
  Double_t ctime = timer.CpuTime();
  printf("RealTime=%f seconds, CpuTime=%f seconds\n",rtime,ctime);
}