#ifdef CDF
#include <BaBar/Experiment.hh>
#endif

#include <iostream>
#include <iomanip>

#include "Stntuple/obj/TStnPhotonBlock.hh"
#include "Stntuple/obj/TStnLinkBlock.hh"

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

   if (R__b.IsReading()) {
      Version_t R__v = R__b.ReadVersion(); if (R__v) { }
      R__b >> fNPhotons;
      R__b >> fNcencl;
      R__b >> fNphoidl;
      R__b >> fNphoisol;
      R__b >> fNphoid;
      R__b >> fNphoiso;
      fPhotonList->Streamer(R__b);
      if(R__v >=2 ) {
	fTowerLinkList->Streamer(R__b);
      }
   } else {
      R__b.WriteVersion(TStnPhotonBlock::IsA());
      R__b << fNPhotons;
      R__b << fNcencl;
      R__b << fNphoidl;
      R__b << fNphoisol;
      R__b << fNphoid;
      R__b << fNphoiso;
      fPhotonList->Streamer(R__b);
      fTowerLinkList->Streamer(R__b);
   }
}

//_____________________________________________________________________________
TStnPhotonBlock::TStnPhotonBlock() {
  fNPhotons   = 0;
  fPhotonList = new TClonesArray("TStnPhoton",10);
  fPhotonList->BypassStreamer(kFALSE);
  fTowerLinkList = new TStnLinkBlock();
}


//_____________________________________________________________________________
TStnPhotonBlock::~TStnPhotonBlock() {
  fPhotonList->Delete();
  delete fPhotonList;
  delete fTowerLinkList;
}


//_____________________________________________________________________________
void TStnPhotonBlock::Clear(const char* opt) {
  fNPhotons = 0;
  fNcencl   = 0;
  fNphoidl  = 0;
  fNphoisol = 0;
  fNphoid   = 0;
  fNphoiso  = 0;
  fPhotonList->Clear();
  fTowerLinkList->Clear();
}

//_____________________________________________________________________________
void TStnPhotonBlock::Print(const char* opt) const {
  std::cout << " **************** print_photons *********** Npho: " << fNPhotons
       << "--------------------- " 
       << std::endl
       << "  i     Etc     Eta      Phi     Px      Py    Status    "
       << " Co4  Had/Em    Lshr2        Z       Det   Time "
       << std::endl;
  std::cout.precision(5);
  std::cout.setf(std::ios::right,std::ios::adjustfield);
  for (int i=0; i<fNPhotons; i++) {
    TStnPhoton* photon = Photon(i);
    std::cout << std::setw(3) << i << " " 
	 << std::setw(7)  << photon->fEtc << " " 
	 << std::setw(7)  << photon->Momentum()->Eta() << " " 
	 << std::setw(7)  << photon->Momentum()->Phi() << " " 
	 << std::setw(7)  << photon->Momentum()->Px() << " " 
	 << std::setw(7)  << photon->Momentum()->Py() << " " 
	 << std::setw(7)  << photon->fStat << " " 
	 << std::setw(7)  << photon->fCo4 << " " 
	 << std::setw(7)  << photon->fHadem << " " 
	 << std::setw(7)  << photon->fLshr << " " 
	 << std::setw(7)  << photon->fZv << " " 
	 << std::setw(7)  << photon->fDetector << " " 
	 << std::setw(7)  << photon->fTime << " " 
	 << std::endl;
  }
}