///////////////////////////////////////////////////////////////////////////////
//  Missing Et block (P.Murat)
//  Description of the variables:
//  -----------------------------
// - fMet, fMetphi - magnitude and angl of missing Et, corrected in different 
//   ways:
//   [0] - raw values
//   [1] - calculated at best Z
//   ....
// - fSumet - sum(Et) for all the calorimeter towers
//
///////////////////////////////////////////////////////////////////////////////
#ifdef __GNUG__
#pragma implementation
#endif

#include "Stntuple/obj/TStnMetBlock.hh"

ClassImp(TStnMetBlock)


//_____________________________________________________________________________
void TStnMetBlock::ReadV1(TBuffer& R__b) {
  //  read in version 1 of MET block

  Float_t* emoe = fMet;
  R__b.ReadArray(emoe);
  emoe = fMetphi;
  R__b.ReadArray(emoe);
  R__b >> fMetjt;
  R__b >> fMetjtphi;
  R__b >> fMetj1;
  R__b >> fMetj1phi;
  R__b >> fMetc;
  R__b >> fMetcphi;
  R__b >> fHtc;
  R__b >> fMetcf;
  R__b >> fMetcfphi;
  R__b >> fHtcf;
  R__b >> fEtatrks;
  R__b >> fPhitrks;
  R__b >> fPttrks;
  R__b >> fSeedpt;
  R__b >> fSumet;
  R__b >> fSumetjet;
  R__b >> fMetsig;
}


//______________________________________________________________________________
void TStnMetBlock::Streamer(TBuffer &R__b)
{
   // Stream an object of class TStnMetBlock.

  int nwf = &fEOR-fMet;

   if (R__b.IsReading()) {
     Version_t R__v = R__b.ReadVersion(); 
     if (R__v == 1) { 
       ReadV1(R__b);
     }
     else {
					// v2=v3 (v3 is current for the moment)
       R__b.ReadFastArray(fMet,nwf);
     }
   } 
   else {
      R__b.WriteVersion(TStnMetBlock::IsA());
      R__b.WriteFastArray(fMet,nwf);
   }
}

//_____________________________________________________________________________
TStnMetBlock::TStnMetBlock() {
  Clear();
}

//_____________________________________________________________________________
TStnMetBlock::~TStnMetBlock() {
}


//_____________________________________________________________________________
void TStnMetBlock::Clear(const char* opt) {
  // fill MET block with non-physical values

  fLinksInitialized = 0;

  for (int i=0; i<5; i++) {
    fMet   [i] = -1.;
    fMetphi[i] = 100.;
  }
  fMetjt       = -1.;
  fMetjtphi    = 100.;
  fMetj1       = -1.;
  fMetj1phi    = 100.;
  fMetc        = -1.;
  fMetcphi     = 100.;
  fHtc         = -1.;
  fMetcf       = -1.;
  fMetcfphi    = 100.;
  fHtcf        = -1.;
  fEtatrks     = 100.;
  fPhitrks     = 100.;
  fPttrks      = -1.;
  fSeedpt      = -1.;
  fSumet       = -1;
  fSumetjet    = -1.;
  fMetsig      = -1;

  fZ0          = 0.;
}

//_____________________________________________________________________________
void TStnMetBlock::Print(const char* opt) const {
  printf(" met[0], phi[0], met[1], phi[1] = %10.3f %10.3f %10.3f %10.3f\n",
	 fMet[0], fMetphi[0],
	 fMet[1], fMetphi[1]);
}