#include "Stntuple/jet/TMetUtil.hh"
#include "Stntuple/jet/TJetUtil.hh"
#include "Stntuple/obj/TStnEvent.hh"
#include "Stntuple/obj/TStnHeaderBlock.hh"
#include "Stntuple/obj/TStnJetBlock.hh"
#include "Stntuple/obj/TStnMetBlock.hh"
#include "Stntuple/obj/TStnVertexBlock.hh"
#include "Stntuple/alg/TStntuple.hh"

//_____________________________________________________________________________
Int_t TMetUtil::MetCorr(TStnEvent* event) {

  TStnHeaderBlock* fHeaderBlock =  
    (TStnHeaderBlock*)event->UnpackDataBlock("HeaderBlock"); 

  TStnJetBlock* fJetBlock = 
    (TStnJetBlock*)event->UnpackDataBlock("JetBlock"); 
  if(!fJetBlock) return 1;

  TStnMetBlock* fMetBlock = 
    (TStnMetBlock*)event->UnpackDataBlock("MetBlock"); 
  if(!fMetBlock) return 1;

  int run = fHeaderBlock->RunNumber(); 
  int evt = fHeaderBlock->EventNumber(); 
  bool qMc = fHeaderBlock->McFlag();  

  TJetUtil::FillJetCorr(event,"JetBlock");

  for(int ij=0; ij<fJetBlock->NJets(); ij++) {
    TStnJet* jet = fJetBlock->Jet(ij);
    float l5corr = jet->Corfd();
    // ...
  }

  // save the answer in the met block
  //fMetBlock->fMetc = ;
  //fMetBlock->fMetcPhi = ;

  return 0;
}