//////////////////////////////////////////////////////////////////////////
//                                                                      //
// TMBTree Jets class                                                   //
//                                                                      //
//////////////////////////////////////////////////////////////////////////

#include "tmb_tree/TMBJets.hpp"
#include "tmb_tree/TMBTrks.hpp"
#include "tmb_tree/TMBVrts.hpp"


ClassImp(TMBJets)


//_____________________________________________________________________________
 TMBJets::TMBJets      (Float_t px,   Float_t py,   Float_t pz,   Float_t E,
                   Float_t pT,   Float_t phi,  Float_t eta,   Float_t q)
{
Set0   (px, py, pz, E, pT, phi, eta, q);

}

 void TMBJets::Set0    (Float_t px,   Float_t py,   Float_t pz,   Float_t E,
                   Float_t pT,   Float_t phi,  Float_t eta,   Float_t q)
       {

               _px    = px; 
               _py    = py; 
               _pz    = pz; 
               _E     = E; 
               _pT    = pT; 
               _phi   = phi; 
               _eta   = eta; 
               _q     = q;
	       //  alwyas initialize corrections to 0 
	       _jes_C=0;
	       _jes_dC_stat=0;
	       _jes_dC_sys=0;
	       _jes_metC=0;
	       _jes_metdC_stat=0;
	       _jes_metdC_sys=0;

 
       }


 void TMBJets::Set1    (Float_t dphi, Float_t deta, 
                   Float_t emf,  Float_t em1f, Float_t em2f, Float_t em3f,
                   Float_t ccmg, Float_t icdf, Float_t ecmg, Float_t icrf,
                   Float_t fh1f, Float_t fh2f, Float_t fh3f, Float_t chf,
                   Float_t emcc, Float_t hadcc,Float_t emec, Float_t hadec,
                   Float_t hot,  Float_t mxET,
                   Float_t cpsE, Float_t etaW, Float_t phiW, 
                   Float_t sET,  Float_t vPT,  Float_t iPT,  Float_t seedET)
       {

               _dphi  = dphi; 
               _deta  = deta; 
               _emf   = emf; 
               _em1f  = em1f; 
               _em2f  = em2f; 
               _em3f  = em3f;
               _ccmg  = ccmg;
               _icdf  = icdf; 
               _ecmg  = ecmg;
               _icrf  = icrf; 
               _fh1f  = fh1f; 
               _fh2f  = fh2f; 
               _fh3f  = fh3f; 
               _emcc  = emcc;
               _hadcc = hadcc;
               _emec  = emec;
               _hadec = hadec;
               _chf   = chf; 
               _hot   = hot; 
               _mxET  = mxET; 
               _cpsE  = cpsE; 
               _etaW  = etaW; 
               _phiW  = phiW; 
               _sET   = sET; 
               _vPT   = vPT; 
               _iPT   = iPT; 
               _seedET= seedET; 
 
       }


 void TMBJets::Set2    (Int_t ntrk, Int_t nps,  Int_t Nitems, 
                   Int_t n90,  Int_t split_merge_word, 
                   const char *algoname, const TRefArray *tracks, TRef vtxref)
       {

               _ntrk  = ntrk; 
               _nps   = nps; 
               _Nitems= Nitems; 
               _n90   = n90;
               _split_merge_word =  split_merge_word;
               _algoname = algoname;              
               _vtxref   = vtxref;

	      _tracks.Clear();
	      if (tracks) {
		 for (Int_t i=0; i<tracks->GetLast()+1; i++) {
		    _tracks.Add( tracks->At(i) );
		 }
	      }
	       
       }

 void TMBJets::SetL1(Float_t l1set)
{
  _l1set = l1set;
}

 void TMBJets::SetCorr (Float_t jes_C, Float_t jes_dC_stat, Float_t jes_dC_sys,
		   Float_t jes_metC, Float_t jes_metdC_stat, Float_t jes_metdC_sys){

	       _jes_C=jes_C;
	       _jes_dC_stat=jes_dC_stat;
	       _jes_dC_sys=jes_dC_sys;
	       _jes_metC=jes_metC;
	       _jes_metdC_stat=jes_metdC_stat;
	       _jes_metdC_sys=jes_metdC_sys;
}

//______________________________________________________________________________