#include <Stntuple/obj/TDcasDataBlock.hh>
#include <iostream>

ClassImp(TDcasDataBlock)

//______________________________________________________________________________
void TDcasDataBlock::Streamer(TBuffer &R__b)
{
   // Stream an object of class TDcasDataBlock - this is the only 
  // non-trivial method of this class

  if (R__b.IsReading()) {
    Version_t R__v = R__b.ReadVersion(); if (R__v) { }
    R__b.ReadFastArray((Short_t*) fEmEnergy,24*24);
    R__b.ReadFastArray((Short_t*) fHdEnergy,24*24);
    R__b.ReadFastArray((Short_t*) fTowUsage,24*24);
  } 
  else {
    R__b.WriteVersion(TDcasDataBlock::IsA());
    R__b.WriteFastArray((Short_t*) fEmEnergy,24*24);
    R__b.WriteFastArray((Short_t*) fHdEnergy,24*24);
    R__b.WriteFastArray((Short_t*) fTowUsage,24*24);
  }
}


//_____________________________________________________________________________
TDcasDataBlock::TDcasDataBlock() {
}

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


//_____________________________________________________________________________
void  TDcasDataBlock::Clear(Option_t* opt)  {
  memset(fEmEnergy,0,24*24*sizeof(Short_t));
  memset(fHdEnergy,0,24*24*sizeof(Short_t));
  memset(fTowUsage,0,24*24*sizeof(Short_t));
}


//_____________________________________________________________________________
void  TDcasDataBlock::Print(Option_t* opt)  const {

  if (strcmp(opt,"") == 0) {
    for(int i=0; i<3; i++) std::cout << " eta phi  Em Had Use   ";
    std::cout << std::endl;

    for (int ieta=0; ieta<24; ieta++) {
      for (int iphi=0; iphi<24; iphi++) {
	printf("%3i %3i %4i %3i %3x   ",
	       ieta, iphi, fEmEnergy[iphi][ieta],fHdEnergy[iphi][ieta],
	       fTowUsage[iphi][ieta]);
	if(iphi%3==2) std::cout << std::endl;
      }
    }
  }
  else if (strcmp(opt,"tot") == 0) {
//-----------------------------------------------------------------------------
//  print 24x24 matrix with the total energies to cross check MET
//-----------------------------------------------------------------------------
    printf("iphi:ieta ");
    for (int i=0; i<24; i++) printf("%4i ",i);
    printf("\n");

    for (int iphi=0; iphi<24; iphi++) {
      printf("%5i     ",iphi);
      for (int ieta=0; ieta<24; ieta++) {
	printf("%5i",GetEnergy(iphi,ieta));
      }
      printf("\n");
    }
  }
}