//----------------------------------
//  MC analysis driver script
//-----------------------------------
TStnAna* x   = 0;
TMcAnaModule* m_mc;

void mc_ana(const char* Filename = NULL, Int_t NEvents = 10) {
  char filename[200];

  if (Filename == 0) {
    sprintf(filename,"%s/%s",getenv("VAL_DATA_DIR"),
	    "isajet_qq_today_val_stntuple.root");
  }
  else {
    sprintf(filename,"%s",Filename);
  }

  x = new TStnAna(filename);
  m_mc = new TMcAnaModule();
  x->AddModule(m_mc);

  x->Run(NEvents);
  //  m_mc->PlotHistograms(-1,1);
}