// Macro for running Panda simulation  with Geant3  or Geant4 (M. Al-Turany)
// This macro is supposed to run the full simulation of the panda detector
// to run the macro:
// root  sim_complete.C  or in root session root>.x  sim_complete.C
// to run with different options:(e.g more events, different momentum, Geant4)
// root  sim_complete.C"(100, "TGeant4",2)"

void sim_complete(TString prefix = "", Int_t nEvents = 100, TString  inputGenerator ="", Double_t BeamMomentum = 0., TString  SimEngine ="TGeant3")
{
  //-----User Settings:------------------------------------------------------
  TString  parAsciiFile   = "all_day1.par";
  // TString inputGenerator = 
  // EvtGen -> "xxxxxxxx.dec"
  // DPM    -> "dpm_xxxxx"
  // FTF    -> "ftf_xxxxx"
  // BoxGen -> "box:type(211,1):p(1):tht(10,120)"

  // TString  inputGenerator = "psi2s_Jpsi2pi_Jpsi_mumu.dec";
//  TString  inputDir = gSystem->Getenv("VMCWORKDIR");
//  inputDir += "/macro/qa/day1/";

  gRandom->SetSeed();
  //-------------------------------------------------------------------------
  // -----   Create the Simulation run manager ------------------------------
  PndMasterRunSim *fRun = new PndMasterRunSim();
  fRun->SetOptions("day1+gem");
  fRun->SetInput(inputGenerator);
//  fRun->SetInputDir(inputDir);
  fRun->SetName(SimEngine);
  fRun->SetParamAsciiFile(parAsciiFile);
  fRun->SetNumberOfEvents(nEvents);
  fRun->SetBeamMom(BeamMomentum);
  // -----  Initialization   ------------------------------------------------
  fRun->Setup(prefix);
  // -----   Geometry   -----------------------------------------------------
  fRun->CreateGeometry();
  // -----   Event generator   ----------------------------------------------
  fRun->SetGenerator();
  //------   Store trajectory -----------------------------------------------
  fRun->SetStoreTraj(kTRUE);
  // -----   Add tasks   ----------------------------------------------------
  fRun->AddSimTasks();
  // -----   Intialise and run   --------------------------------------------
  fRun->Init();
  fRun->Run(nEvents); 
  fRun->Finish();
  
  //exit(0);  
};