void ana_etac_task(TString InFile="test_fast.root", int nevts=0, double pbarmom = 6.232069, int decayMode = 0, int recoilMode = 0, int dataMode = 0, TString pidQuality = "All")
{
	// *** the output file for FairRunAna
	TString OutFile = InFile;
	OutFile.ReplaceAll(".root","");
	OutFile += Form("_%s_ntp.root",pidQuality.Data());
	// *** initialization
	FairLogger::GetLogger()->SetLogToFile(kFALSE);

	FairRunAna* fRun = new FairRunAna();
	fRun->SetWriteRunInfoFile(kFALSE);
	fRun->SetInputFile(InFile);
	fRun->SetOutputFile("data/fairout_dummy.root");

	// *** take constant field; needed for PocaVtx
	RhoCalculationTools::ForceConstantBz(20.0);

	// *** HERE OUR ANALYSIS TASK GOES!
	PndScrutJpsiEtacTask *scrutTask = new PndScrutJpsiEtacTask(pbarmom, OutFile, decayMode, recoilMode, dataMode, pidQuality);
	fRun->AddTask(scrutTask);

	// *** and run analysis
	fRun->Init();
	fRun->Run(0,nevts);
}