// Lmd FTF Sim macro
void runLumiPixel0SimFTF(const int nEvents = 10, const int startEvent = 0,
		const double mom = 15, TString input, TString storePath = "tmpOutputDPM",
		const double beam_X0 = 0.0, const double beam_Y0 = 0.0,
		const double target_Z0 = 0.0, const double beam_width_sigma_X = 0.0,
		const double beam_width_sigma_Y = 0.0,
		const double target_width_Z = 0.0, // beam offset and smearing parameters
		const double beam_grad_X = 0.0, const double beam_grad_Y = 0.0,
		const double beam_grad_sigma_X = 0.0, const double beam_grad_sigma_Y = 0.0, // beam gradiant parameters
		const int verboseLevel = 0) {
	// gRandom->SetSeed(seed);
	Int_t mode = 1;
	TStopwatch timer;
	timer.Start();
	gDebug = 0;

	//output1
	TString simOutput = storePath + "/Lumi_MC_";
	simOutput += startEvent;
	simOutput += ".root";
	TString parOutput = storePath + "/Lumi_Params_";
	parOutput += startEvent;
	parOutput += ".root";
	//Load basic libraries
	gROOT->Macro("$VMCWORKDIR/gconfig/rootlogon.C");
	gSystem->Load("libSds");
	gSystem->Load("libLmd");
	FairRunSim *fRun = new FairRunSim();
	cout << "All libraries succsesfully loaded!" << endl;

	//set the MC version used
	fRun->SetName("TGeant4");

	fRun->SetOutputFile(simOutput);

	//set material
	fRun->SetMaterials("media_pnd.geo");

	//  //create and add detectors
// //-------------------------  CAVE      -----------------

	FairModule *Cave = new PndCave("CAVE");
	//	Cave->SetGeometryFileName("pndcave.geo");
	Cave->SetGeometryFileName("pndcaveVAC.geo"); //vac cave is not really necessary 
	fRun->AddModule(Cave);
	//-------------------------  Magnet   -----------------
	FairModule *Magnet = new PndMagnet("MAGNET");
	Magnet->SetGeometryFileName("FullSuperconductingSolenoid_v831.root");
	fRun->AddModule(Magnet);
	FairModule *Dipole = new PndMagnet("MAGNET");
	Dipole->SetGeometryFileName("dipole.geo");
	fRun->AddModule(Dipole);
	//-------------------------  Pipe     -----------------
	FairModule *Pipe = new PndPipe("PIPE");
	Pipe->SetGeometryFileName("beampipe_201309.root");
	fRun->AddModule(Pipe);


	PndLmdDetector *Lum = new PndLmdDetector("LUM", kTRUE);
	Lum->SetExclusiveSensorType("LumActive"); //ignore MVD
	Lum->SetGeometryFileName("Luminosity-Detector.root"); // new sensors
	Lum->SetVerboseLevel(verboseLevel);
	fRun->AddModule(Lum);

	//particle generator
	FairPrimaryGenerator* primGen = new FairPrimaryGenerator();
	if (beam_X0 != 0.0 || beam_Y0 != 0.0 || beam_width_sigma_X > 0.0 || beam_width_sigma_Y > 0.0) {
		primGen->SmearGausVertexXY(true);
		primGen->SetBeam(beam_X0, beam_Y0, beam_width_sigma_X, beam_width_sigma_Y);
	}
	if(target_Z0 != 0.0 || target_width_Z > 0.0) {
		primGen->SmearGausVertexZ(true);
		primGen->SetTarget(target_Z0, target_width_Z);
	}
	if (beam_grad_X != 0.0 || beam_grad_Y != 0.0 || beam_grad_sigma_X > 0.0 || beam_grad_sigma_Y > 0.0) {
		primGen->SetBeamAngle(beam_grad_X, beam_grad_Y, beam_grad_sigma_X,
				beam_grad_sigma_Y);
	}

	fRun->SetGenerator(primGen);

	// FTF Generator
	PndFtfGenerator* ftfGen = new PndFtfGenerator(input);
	primGen->AddGenerator(ftfGen);


//reading the new field map in the old format
	fRun->SetBeamMom(mom);
	PndMultiField *fField= new PndMultiField("AUTO");
	fRun->SetField(fField);

	fRun->SetField(fField);

	if (nEvents < 450)
		fRun->SetStoreTraj(kTRUE); // toggle this for use with EVE
	else
		fRun->SetStoreTraj(kFALSE);

// //FairLogger
// // get handle
// 	FairLogger *logger = FairLogger::GetLogger();
// //
// // log to screen and to file
// 	logger->SetLogToScreen(kTRUE);
// 	logger->SetLogToFile(kFALSE);

// 	logger->SetLogVerbosityLevel("LOW");

// // Set different levels of verbosity. In the example everything >=INFO goes to the
// // file and everything >= ERROR is printed on the screen
// // LogLevels are (FATAL, ERROR, WARNING, INFO, DEBUG, DEBUG1, DEBUG2, DEBUG3, DEBUG4)
// 	logger->SetLogScreenLevel("ERROR"); //Only FATAL and ERROR to screen

	fRun->Init();
	//	((TGeant4*)gMC)->ProcessGeantCommand("/mcVerbose/eventAction 0");

// // Fill the Parameter containers for this run
// //-------------------------------------------
	FairRuntimeDb *rtdb = fRun->GetRuntimeDb();
	Bool_t kParameterMerged = kTRUE;
	FairParRootFileIo* output = new FairParRootFileIo(kParameterMerged);
	output->open(parOutput.Data(), "RECREATE");
	rtdb->setOutput(output);

// PndMultiFieldPar* Par = (PndMultiFieldPar*) rtdb->getContainer("PndMultiFieldPar");
// if (fField) {  Par->SetParameters(fField); }
// Par->setInputVersion(fRun->GetRunId(),1);
// Par->setChanged();

// Transport nEvents
// -----------------

	fRun->Run(nEvents);

	rtdb->saveOutput();
	rtdb->print();

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