void runLumiPixel5BackProp(const int nEvents = 10, const int startEvent = 0,
		TString storePath = "tmpOutput", const int verboseLevel = 0,
		const TString Method = "Geane", const bool mergedHits = true,
		const double mom = 15, double ipc_x = 0.0, double ipc_y = 0.0,
		double ipc_z = 0.0) {
	// ========================================================================
	// Input file (MC events)
	TString MCFile = storePath + "/Lumi_MC_";
	MCFile += startEvent;
	MCFile += ".root";
	TString DigiFile = storePath + "/Lumi_digi_";
	DigiFile += startEvent;
	DigiFile += ".root";
	// Digi file
	TString RecoFile = storePath;
	if (mergedHits)
		RecoFile += "/Lumi_recoMerged_";
	else
		RecoFile += "/Lumi_reco_";
	RecoFile += startEvent;
	RecoFile += ".root";
	// TCand file
	TString CandFile = storePath + "/Lumi_TCand_";
	CandFile += startEvent;
	CandFile += ".root";
	// Parameter file
	TString parFile = storePath + "/Lumi_Params_";
	parFile += startEvent;
	parFile += ".root";
	// Track file
	TString TrkFile = storePath + "/Lumi_Track_";
	TrkFile += startEvent;
	TrkFile += ".root";

	// // ----  Load libraries   -------------------------------------------------
	// gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
	// gSystem->Load("libSds");
	// gSystem->Load("libSdsReco");
	// gSystem->Load("libLmd");
	// gSystem->Load("libLmdReco");
	// gSystem->Load("libLmdTrk");
	// // ------------------------------------------------------------------------

	// Output file
	TString outFile = storePath + "/Lumi_Geane_";
	outFile += startEvent;
	outFile += ".root";

	std::cout << "MCFile  : " << MCFile.Data() << std::endl;
	std::cout << "DigiFile: " << DigiFile.Data() << std::endl;
	std::cout << "RecoFile: " << RecoFile.Data() << std::endl;
	std::cout << "TCandFile: " << CandFile.Data() << std::endl;
	std::cout << "TrackFile: " << TrkFile.Data() << std::endl;
	std::cout << "GeaneFile: " << outFile.Data() << std::endl;
	// ---  Now choose concrete engines for the different tasks   -------------
	// ------------------------------------------------------------------------

	// In general, the following parts need not be touched
	// ========================================================================

	// -----   Timer   ------ --------------------------------------------------
	TStopwatch timer;
	timer.Start();
	// ------------------------------------------------------------------------

	// -----   Reconstruction run   -------------------------------------------
	FairRunAna *fRun = new FairRunAna();
	fRun->SetInputFile(MCFile);
	fRun->AddFriend(DigiFile);
	fRun->AddFriend(RecoFile);
	fRun->AddFriend(CandFile);
	fRun->AddFriend(TrkFile);
	fRun->SetOutputFile(outFile);
	// ------------------------------------------------------------------------

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

	// =========================================================================
	// ======                 Back-Propagating                     ======
	// =========================================================================

	Double_t fpBeam = mom;
	TVector3 IP(ipc_x, ipc_y, ipc_z);
	if (Method == "Geane") {
		FairGeane *Geane = new FairGeane();
		fRun->AddTask(Geane);
		PndLmdGeaneTask* lmdgeane = new PndLmdGeaneTask(fpBeam, IP);
		lmdgeane->SetVerbose(verboseLevel);
		fRun->AddTask(lmdgeane);
	} else {
	  if (Method == "RK") {
		  //PndLmdBPRungeKuttaTask have dependences from genfit -> switched off until movement to gentfit2
		  cout<<"I'm very sorry, but BPRungeKuttaTask  option for LMD is currently switched off :("<<endl;
			// PndLmdBPRungeKuttaTask* lmdbp = new PndLmdBPRungeKuttaTask(fpBeam, IP);
			// lmdbp->SetVerbose(verboseLevel);
			// fRun->AddTask(lmdbp);
		} else {
			cout << "Method " << Method.Data() << " doesn't exist!" << endl;
			break;
		}
	}
	rtdb->setOutput(parInput1);
	rtdb->print();
	// =====                        End of Geane                           =====
	// =========================================================================

	// -----   Intialise and run   --------------------------------------------
	fRun->Init();
	//  // PndEmcMapper *emcMap = PndEmcMapper::Instance(6);
	//  PndEmcMapper *emcMap = PndEmcMapper::Instance();
	// //Geane->SetField(fRun->GetField());
	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 succesfully." << endl;
	cout << "Output file is " << outFile << endl;
	cout << "Parameter file is " << parFile << endl;
	cout << "Real time " << rtime << " s, CPU time " << ctime << " s" << endl;
	cout << endl;
	// ------------------------------------------------------------------------

}