//==================================================================
//
// Run TEleCand.cc
//
//	Version 0.0:      P. Koehn    05/05/01
//	Version 1.0:      E. Thomson  02/06/02
//	Version 2.0:      E. Thomson  02/22/02
//
// Run TopEventModule to classify electrons as tight and loose
// Run Stntuple
// Run this to make histograms (done-02/22/02) or perform analysis (not yet!)
// Run ele_display.C to make a postscript file of 1D histograms
//==================================================================
{

  //timer start
  gROOT->Time();


  // Instantiate a StnAna object
  //----------------------------
  TStnAna x("/nosu6a/thomsone/stntuple_topevent_ttbar.root"); // your input Stntuple file
  //TChain* chain = new TChain("STNTUPLE");
  //chain->Add("/nosu6a/thomsone/topevent/topevent_4166.root");
  //chain->Add("/nosu6a/thomsone/topevent/topevent_4167.root");
  //TStnAna x(chain); 

  // Make the output histogram file
  //--------------------------------
  TFile* HistFile = new TFile("ele_cand.root","RECREATE","test file");
  
  // Compile modules needed below
  //------------------------------
  //gSystem->CompileMacro("Stntuple/top/TConvEleCand.cc","k");
  gSystem->CompileMacro("Stntuple/top/TEleCand.cc","k");

  // Construct modules, store pointers
  //----------------------------------------
  //TConvEleCand* conv  = new TConvEleCand(HistFile,"ConvEleC","ConvEleC");
  TEleCand* ec  = new TEleCand(HistFile,"ElectronBlock","EleCandModule ");
  
  // User-defined options for TEleCand
  //----------------------------------
  ec->enableEleHist(true);
  ec->enableEleTopEvent(true);
  ec->setDebugLevel(1);
  ec->enableEleDelXCut(false);
  ec->enableEleDelZCut(false);
  ec->enableEleChi2StripCut(false);
  ec->enableEleLshrCut(false);

  
  // Add module to TStnAna object
  //-----------------------------
  //x.AddModule(conv);
  x.AddModule(ec);
  
  // Run the StnAna object over the root file specified at instantiation
  //--------------------------------------------------------------------
  x.Run(); 

  printf("Job Finished.\n");
  printf("Output histograms are stored in: %s\n",gFile->GetName());
  
  // Close the histogram file
  //--------------------------
  HistFile->Write();
  HistFile->Close();

  //timer stop
  gROOT->Time();
}