#ifndef __CINT__
#include "TROOT.h"
#include "TH1.h"
#include "TObjArray.h"
#include "TSystem.h"
#include "TRegexp.h"
#include "TFile.h"
#include "TTree.h"
#include "TInterpreter.h"
#include <Stntuple/THBookStnEvent.hh>
#include <Stntuple/THBookStntuple.hh>
#include <Stntuple/TStntupleEvent.hh>
#include <Stntuple/TStntuple.hh>
#endif

//_____________________________________________________________________________
int book_histograms() {
  return 0;
}


//_____________________________________________________________________________
int end_job() {
  return 0;
}

//_____________________________________________________________________________
int process_event(TStntupleEvent* event) {
  event->InitJets();
  event->CorrectMissingEt();
  event->PrintAll();
  printf(" ********************* missing Et:\n");
  printf(" met[0:4], metjt, metc, corrmet: %11.4f %11.4f %11.4f %11.4f %11.4f %11.4f %11.4f\n",
	 event->fMet->fMet[0], 
	 event->fMet->fMet[1], 
	 event->fMet->fMet[2], 
	 event->fMet->fMet[3], 
	 event->fMet->fMet[4], 
	 event->fMet->fMetjt,
	 event->fMet->fMetc,
	 event->fCorrMet->Mod());
  printf(" phi[0:4], phijt, phic, corrmet: %11.4f %11.4f %11.4f %11.4f %11.4f %11.4f %11.4f\n",
	 event->fMet->fMetphi[0], 
	 event->fMet->fMetphi[1], 
	 event->fMet->fMetphi[2], 
	 event->fMet->fMetphi[3], 
	 event->fMet->fMetphi[4], 
	 event->fMet->fMetjtphi,
	 event->fMet->fMetcphi,
	 TVector2::Phi_0_2pi(event->fCorrMet->Phi()));
  return 0;
}