//==================================================================
//
//	jet_cand.C Macro
//
//	Version 0.0:      P. Koehn    5/5/01
//
// root[] .x jet_cand.C
//
//==================================================================
{

  //timer start
  gROOT->Time();

  // Instantiate a StnAna object
  //----------------------------
  TStnAna x("/nosu3a/pkoehn/mc/top/new_pyt_ttbar_dilep_ee_stntuple_440.root"); // data

  // Compile required code
  //-----------------------
  gSystem->CompileMacro("TJetCand.cc","+");
  
  //
  // Make the output histogram file
  TFile* HistFile = new TFile("jet_cand.root","RECREATE","test file");

  // Create pointer to a JetCand object
  //----------------------------------------
  cout << ">>> Making JetCand object <<<" << endl;
  TJetCand* myjet  = new TJetCand(HistFile,"JetC","JetC");
  
  //
  // Turn the histograms on
  myjet->enableJetHist(true);  // default jets
  myjet->enableTJetHist(true); // tight top jets
  myjet->enableLJetHist(true); // loose top jets
	
  //
  // enable and set jet cuts		  
  myjet->enableJetEtCut(true); myjet->setJetEtCut(10.0);
  myjet->enableJetEmfCut(true); myjet->setJetEmfCut(0.8);
  myjet->enableJetDtetaCut(true); myjet->setJetDtetaCut(2.0);
  
  //  The module must be added to TStnAna object the for 
  //  events to be processed by the TopCand module properly.
  x.AddModule(myjet);

  //
  // Run the StnAna object over the input file specified above
  x.Run(500); 

  myjet->PlotHistograms();

  printf("\n Job Finished.  Output histograms are stored in: %s\n",gFile->GetName());

  HistFile->Write();
  //HistFile->Close();

}