{
	// Load basic libraries
	gROOT->LoadMacro("$VMCWORKDIR/gconfig/rootlogon.C");
	rootlogon();

	Int_t iVerbose = 0;
	Int_t nEvents = 0;

	// -----   Timer   --------------------------------------------------------
	TStopwatch timer;
	timer.Start();
	// --------------
	TString suffix="_sm_1250";//dpm_0370";

	TString simFile = "sim_" + suffix +".root";
	TString digiFile = "digi_"+ suffix +".root";

	TString parFile = "params_"+ suffix +".root";

	TString digiparFile = gSystem->Getenv("VMCWORKDIR");
 	digiparFile += "/macro/params/lumi.digi.par";
    FairRunAna *fRun = new FairRunAna();
    fRun->SetInputFile(simFile);
    fRun->SetOutputFile(digiFile);

    // -----  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);

    fRun->LoadGeometry();

    /** Run Digitization */

    PndLumiDigiProducer* lumiDigiProducer = new PndLumiDigiProducer(iVerbose);
    fRun->AddTask(lumiDigiProducer);


    rtdb->setOutput(parInput1);
    rtdb->print();

      // -----   Intialise and run   --------------------------------------------
    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 succesfully." << endl;
    cout << "Output file is "    << digiFile << endl;
    cout << "Parameter file is " << parFile << endl;
    cout << "Real time " << rtime << " s, CPU time " << ctime << " s\a\a" << endl;
    cout << endl;
    exit();
}