void runMain//(int iseq, int iseed)
  (char* seq = "", char* seed = "", char* fileGEN="", char* fileFIT="")
{
  //load the shared library
  gSystem->Load("library.so");
  //set up main, eg command line opts
  TList* args = new TList();
  args->Add(new TObjString("Gen_byEv"));//program name
  if(strlen(seq)>0) {
    args->Add(new TObjString("-seq"));    args->Add(new TObjString(seq));
  }
  
  if(strlen(seed)>0) {
    args->Add(new TObjString("-seed"));    args->Add(new TObjString(seed));
  }  
  args->Add(new TObjString("-ROOTfile"));  args->Add(new TObjString("PTerror_uguali.root"));

  if(strlen(fileGEN)>0) {
    args->Add(new TObjString("-genfile"));    args->Add(new TObjString(fileGEN));
  }
  
  //run the code
  args->Print();
  Main m(args);
  int returnCode = m.run_Gen_byEv();
///   ???? 
  gROOT->Reset();
////
  TList* args_fit = new TList();
  args_fit->Add(new TObjString("Fit_fix"));//program name
  if(strlen(seq)>0) {
    args_fit->Add(new TObjString("-seq"));    args_fit->Add(new TObjString(seq));
  }
  
  if(strlen(seed)>0) {
    args_fit->Add(new TObjString("-seed"));    args_fit->Add(new TObjString(seed));
  }
  if(strlen(fileFIT)>0) {
    args_fit->Add(new TObjString("-fitfile"));    args_fit->Add(new TObjString(fileFIT));
  }
  args_fit->Print();
  Main mfit(args_fit);
  int returnCodefit = mfit.run_Fit_fix();
  
}