void phys_UrQMD() {

  TStopwatch timer;
  timer.Start();

 
    // ----  Load libraries   -------------------------------------------------
  gROOT->LoadMacro("$VMCWORKDIR/gconfig/basiclibs.C");
  basiclibs();
  gSystem->Load("libGeoBase");
  gSystem->Load("libParBase");
  gSystem->Load("libBase");
  gSystem->Load("libCbmBase");
  gSystem->Load("libCbmData");
  gSystem->Load("libField");
  gSystem->Load("libGen");
  gSystem->Load("libPassive");
  gSystem->Load("libMvd");
  gSystem->Load("libSts");
  gSystem->Load("libRich");
  gSystem->Load("libTrd");
  gSystem->Load("libTof");
  gSystem->Load("libEcal");
  gSystem->Load("libGlobal");
  gSystem->Load("libKF");
  gSystem->Load("libL1");
  gSystem->Load("libMinuit2"); // Nedded for rich ellipse fitter
  gSystem->Load("libAnalysis"); // D0
  // ------------------------------------------------------------------------
   TString name ="auau.25gev.centr.v11b4"; //->
  TString inFile  = name + ".mc.root";
  TString parFile = name + ".params.root";
  TString bgFile  =  name + ".mc.root";
//  TString recFile =  name + ".rec.root";
  TString recFile =  name + ".rec.root";
  char *outFile = "out.root";
  char *version = "v05a";

  // -----   Reconstruction run   -------------------------------------------
  FairRunAna *fRun= new FairRunAna();
  fRun->SetInputFile(inFile);
  fRun->AddFriend(recFile);
  fRun->SetOutputFile(outFile);
  // ------------------------------------------------------------------------
    TString stsDigiFile = "sts_v11b.digi.par";
  
  // -----  Parameter database   --------------------------------------------
  TString stsDigi = gSystem->Getenv("VMCWORKDIR");
  stsDigi += "/parameters/sts/";
  stsDigi += stsDigiFile;
  FairRuntimeDb* rtdb = fRun->GetRuntimeDb();
  FairParRootFileIo* parIo1 = new FairParRootFileIo();
  FairParAsciiFileIo* parIo2 = new FairParAsciiFileIo();
  parIo1->open(parFile.Data());
  parIo2->open(stsDigi.Data(),"in");
  rtdb->setFirstInput(parIo1);
  rtdb->setSecondInput(parIo2);
  rtdb->setOutput(parIo1);
  rtdb->saveOutput();

//  fRun->LoadGeometry();
  // ------------------------------------------------------------------------


  int iVerbose = 0;
  // ------------------------------------------------------------------------

 
  // Kalman Filter

  CbmKF *KF = new CbmKF();
  fRun->AddTask(KF);

  CbmV0Analysis *V0 = new CbmV0Analysis();
  fRun->AddTask(V0);



  cout << " INITIALISATION " << endl;

  fRun->Init();
//  fRun->Run(0,100);
  fRun->Run(0,2);
//  fRun->Run(0,1000);

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