#ifdef USE_CDFEDM2
//////////////////////////////////////////////////////////////////////////
//
// Component: TowerType.cc
// Purpose: Implementation of TowerType classes
//          
// Created: 05/05/99  Pierre Savard
// History: 
//
//////////////////////////////////////////////////////////////////////////
#include <iostream>
#include "RawDataBanks/RawQieStorableBank.hh"
#include "CalorObjects/TowerType.hh"
#include "CalorObjects/getScintillatorEnergy.hh"

using namespace std;

// namespace calor {

  //
  // TowerType 0
  //
    
  // Equality test

bool TowerType0::operator==(const TowerType0& rhs) const
{
  if (this == &rhs) return true;
  return (_cemEnergy[0] == rhs._cemEnergy[0] &&
	  _cemEnergy[1] == rhs._cemEnergy[1] &&
	  _chaEnergy[0] == rhs._chaEnergy[0] &&
	  _chaEnergy[1] == rhs._chaEnergy[1] &&
	  _stripEnergy == rhs._stripEnergy &&
	  _crackEnergy == rhs._crackEnergy &&
	  _iChaTime[1] == rhs._iChaTime[1] &&
	  _iChaTime[0] == rhs._iChaTime[0]);
}

// EDM stuff

void TowerType0::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >> _iCemEnergy[0];
  iobuffer >> _iCemEnergy[1];
  iobuffer >> _iChaEnergy[0]; 
  iobuffer >> _iChaEnergy[1]; 

  if (version > 6){
    iobuffer >> _nChaHits;
    for (int i=0; i< _nChaHits; i++) {
      if(_nChaHits > i) iobuffer >> _iChaTime[i]; 
    }
  } else if(version > 4){
    iobuffer >> _nChaHits;
    if(_nChaHits > 0) iobuffer >> _iChaTime[0]; 
    if(_nChaHits > 1) iobuffer >> _iChaTime[1]; 
  }
  else{
    iobuffer >> _iChaTime[0]; 
  }
  // we store now the timing with a factor 10 better precision and have to divide by 10
  // such that it is still in ns
  for (int i=0; i< _nChaHits; i++) {
    if (version>8){
      _chaTime[i]=float(_iChaTime[i])/10-float(HTDC_Offset);
    } else {
      _chaTime[i]=float(_iChaTime[i])-float(HTDC_Offset);
    }
  }

  if (version>10){
    iobuffer >> _nCemHits;
    for (int i=0; i< _nCemHits; i++) {
      if(_nCemHits > i) iobuffer >> _iCemTime[i]; 
      _cemTime[i]=float(_iCemTime[i])/10-float(EMTDC_Offset);
    }   
  }
  
}

void TowerType0::writeToBuffer(TBuffer& iobuffer,Version_t version){
  iobuffer << _iCemEnergy[0];
  iobuffer <<  _iCemEnergy[1];
  iobuffer << _iChaEnergy[0]; 
  iobuffer << _iChaEnergy[1]; 

  iobuffer << _nChaHits;
  for (int i=0; i< _nChaHits; i++) {
    if(_nChaHits > i) iobuffer << _iChaTime[i]; 
  }

  iobuffer << _nCemHits;
  for (int i=0; i< _nCemHits; i++) {
    if(_nCemHits > i) iobuffer << _iCemTime[i]; 
  }

}


//
// Printing    (incomplete)
//
 
void TowerType0::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType0& t)
{ 
  os << "CEM energy (0): " << setw(10) <<setprecision(5) << t._cemEnergy[0] << "   ";
  os << "CEM energy (1): " << setw(10) <<setprecision(5) << t._cemEnergy[1] << "   ";
  os << "CHA energy (0): " << setw(10) <<setprecision(5) << t._chaEnergy[0] << "   ";
  os << "CHA energy (1): " << setw(10) <<setprecision(5) << t._chaEnergy[1] << '\n';
  os << "CHA time   (0): " << setw(10) <<setprecision(5) << t._chaTime[0]   << "   ";
  os << "CHA time   (1): " << setw(10) <<setprecision(5) << t._chaTime[1]   << "   ";
  os << "CHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nChaHits)<< '\n';
  os << "CEM time   (0): " << setw(10) <<setprecision(5) << t._cemTime[0]   << "   ";
  os << "CEM time   (1): " << setw(10) <<setprecision(5) << t._cemTime[1]   << "   ";
  os << "CEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nCemHits)<< '\n';
  return os;
}

//
// Manipulation and access
//
  

float TowerType0::pmtEnergy(CellKey& cKey, size_t ipmt) const{

  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _cemEnergy[ipmt];
  case CHA:
    return _chaEnergy[ipmt];
  case STRIP:
    return _stripEnergy;
  case CRACK:
    return _crackEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType0::pmtData(CellKey& cKey, size_t ipmt) const{

  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _iCemEnergy[ipmt];
  case CHA:
    return _iChaEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }
}
  
 
int TowerType0::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::CEM:
    _cemEnergy[ipmt]+=energy;
    _iCemEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::CHA:
    _chaEnergy[ipmt]+=energy;
    _iChaEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType0::addPmtEnergy(Detector det , size_t ipmt,float& energy){

  switch ( det) {
  case CEM:
    _cemEnergy[ipmt]+=energy;
    return 1;
  case CHA:
    _chaEnergy[ipmt]+=energy;
    return 1;
  case STRIP:
    _stripEnergy += energy;
    return 1;
  case CRACK:
    _crackEnergy += energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType0::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _nChaHits;
  case CEM:
    return _nCemHits;
  default:
    return 0;
  }
}
float TowerType0::time(CellKey& cKey) const{ 
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _chaTime[0];
  case CEM:
    return _cemTime[0];
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }  
}

float TowerType0::time(CellKey& cKey, int ihit) const{ 
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _chaTime[ihit];
  case CEM:
    return _cemTime[ihit];
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }  
}



int TowerType0::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case CHA:
    _iChaTime[ihit] = time;
    _chaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  } 
}
int TowerType0::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    _iChaTime[ihit] = time;
    _chaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1; 
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType0::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case CHA:
    _nChaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1;  
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType0::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    _nChaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1; 
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType0 for detector " << d <<std::endl;
    return 0;
  }
}

TowerType0* TowerType0::clone() const{
  return new TowerType0(*this);
}

void TowerType0::puffTower(const CalDataCalib* calib,float CalOffLer[]){

  _cemEnergy[0] = (RawQieStorableBank::unpack_energy(_iCemEnergy[0]) -calib->pedCem()) * calib->calibCem()*CalOffLer[0];
  _cemEnergy[1] = (RawQieStorableBank::unpack_energy(_iCemEnergy[1]) -calib->pedCem()) * calib->calibCem()*CalOffLer[1];
  _chaEnergy[0] = (RawQieStorableBank::unpack_energy(_iChaEnergy[0]) -calib->pedCha()) * calib->calibCha()*CalOffLer[4];
  _chaEnergy[1] = (RawQieStorableBank::unpack_energy(_iChaEnergy[1]) -calib->pedCha()) * calib->calibCha()*CalOffLer[5];
  // just convert to float for now
  //  _chaTime[0] = float(_iChaTime[0]/10)=HTDC_Offset;
}
int TowerType0::iTower() const {
  return -1;
}

//
// TowerType 1
//

  // Equality test

bool TowerType1::operator==(const TowerType1& rhs) const
{
  if (this == &rhs) return true;
  return (_cemEnergy[0] == rhs._cemEnergy[0] &&
	  _cemEnergy[1] == rhs._cemEnergy[1] &&
	  _chaEnergy[0] == rhs._chaEnergy[0] &&
	  _chaEnergy[1] == rhs._chaEnergy[1] &&
	  _whaEnergy[0] == rhs._whaEnergy[0] &&
	  _whaEnergy[1] == rhs._whaEnergy[1] &&
	  _stripEnergy == rhs._stripEnergy &&
	  _crackEnergy == rhs._crackEnergy &&
	  _iChaTime[0] == rhs._iChaTime[0] &&
	  _iWhaTime[0] == rhs._iWhaTime[0]);
}

// EDM stuff

void TowerType1::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iCemEnergy[0];
  iobuffer >>  _iCemEnergy[1];
  iobuffer >>  _iChaEnergy[0]; 
  iobuffer >>  _iChaEnergy[1]; 
  iobuffer >>  _iWhaEnergy[0]; 
  iobuffer >>  _iWhaEnergy[1]; 
  //  iobuffer >>  _stripEnergy;
  //  iobuffer >>  _crackEnergy;

  if (version > 6){
    iobuffer >> _nChaHits;
    for (int i=0; i< _nChaHits; i++) {
      uint2 dummy=0;
      if(_nChaHits > i) {
	iobuffer >> dummy;
	if (i<8) _iChaTime[i]=dummy; 
      }
    }
    iobuffer >> _nWhaHits;
    for (int i=0; i<_nWhaHits; i++) {
      uint2 dummy=0;
      if(_nWhaHits > i) {
	iobuffer >> dummy;
	if (i<8) _iWhaTime[i]=dummy; 
      }
    }
  } else if(version > 4){
    iobuffer >> _nChaHits;
    if(_nChaHits > 0) iobuffer >> _iChaTime[0]; 
    if(_nChaHits > 1) iobuffer >> _iChaTime[1];
    iobuffer >> _nWhaHits;
    if(_nWhaHits > 0) iobuffer >> _iWhaTime[0]; 
    if(_nWhaHits > 1) iobuffer >> _iWhaTime[1]; 
  }
  else{
    iobuffer >>  _iChaTime[0]; 
    iobuffer >>  _iWhaTime[0]; 
  }
  for (int i=0; i< _nChaHits; i++) {
    if (version>8){
      _chaTime[i]=float(_iChaTime[i])/10-float(HTDC_Offset);
    } else {
      _chaTime[i]=float(_iChaTime[i])-float(HTDC_Offset);
    }
  }
  for (int i=0; i< _nWhaHits; i++) {
    if (version>8){
      _whaTime[i]=float(_iWhaTime[i])/10-float(HTDC_Offset);
    } else {
      _whaTime[i]=float(_iWhaTime[i])-float(HTDC_Offset);
    }
  }
  if (version>10){
    iobuffer >> _nCemHits;
    for (int i=0; i< _nCemHits; i++) {
      if(_nCemHits > i) iobuffer >> _iCemTime[i]; 
      _cemTime[i]=float(_iCemTime[i])/10-float(EMTDC_Offset);
    }   
  }
}


void TowerType1::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iCemEnergy[0];
  iobuffer <<  _iCemEnergy[1];
  iobuffer <<  _iChaEnergy[0]; 
  iobuffer <<  _iChaEnergy[1]; 
  iobuffer <<  _iWhaEnergy[0]; 
  iobuffer <<  _iWhaEnergy[1]; 
  //  iobuffer <<  _stripEnergy;
  //  iobuffer <<  _crackEnergy;
  iobuffer << _nChaHits;
  for (int i=0; i< _nChaHits; i++) {
    if(_nChaHits > i) iobuffer << _iChaTime[i]; 
  }
  iobuffer << _nWhaHits;
  for (int i=0; i< _nWhaHits; i++) {
    if(_nWhaHits > i) iobuffer << _iWhaTime[i]; 
  }


  iobuffer << _nCemHits;
  for (int i=0; i< _nCemHits; i++) {
    if(_nCemHits > i) iobuffer << _iCemTime[i]; 
  }
}


//
// Printing  (incomplete)
//

void TowerType1::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType1& t)
{
  os << "CEM energy (0): " << setw(10) <<setprecision(5) << t._cemEnergy[0] << "   ";
  os << "CEM energy (1): " << setw(10) <<setprecision(5) << t._cemEnergy[1] << "   ";
  os << "CHA energy (0): " << setw(10) <<setprecision(5) << t._chaEnergy[0] << "   ";
  os << "CHA energy (1): " << setw(10) <<setprecision(5) << t._chaEnergy[1] << "   ";
  os << "WHA energy (0): " << setw(10) <<setprecision(5) << t._whaEnergy[0] << "   ";
  os << "WHA energy (1): " << setw(10) <<setprecision(5) << t._whaEnergy[1] <<  '\n';
  os << "CHA time   (0): " << setw(10) <<setprecision(5) << t._chaTime[0]   << "   ";
  os << "CHA time   (1): " << setw(10) <<setprecision(5) << t._chaTime[1]   << "   ";
  os << "CEM time   (0): " << setw(10) <<setprecision(5) << t._cemTime[0]   << "   ";
  os << "CEM time   (1): " << setw(10) <<setprecision(5) << t._cemTime[1]   << "   ";
  os << "WHA time   (0): " << setw(10) <<setprecision(5) << t._whaTime[0]   << "   ";
  os << "WHA time   (1): " << setw(10) <<setprecision(5) << t._whaTime[1]   << "   ";
  os << "CHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nChaHits)<< "   ";
  os << "CEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nCemHits)<< "   ";
  os << "WHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nWhaHits)<< '\n';
  return os;
}
   
//
// Manipulation and access
//

float TowerType1::time(CellKey& cKey) const{ 

  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _chaTime[0];
  case WHA:
    return _whaTime[0];
  case CEM:
    return _cemTime[0];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0.0;
  }
}

float TowerType1::time(CellKey& cKey,int ihit) const{ 

  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _chaTime[ihit];
  case WHA:
    return _whaTime[ihit];
  case CEM:
    return _cemTime[ihit];
  default:
    // should throw an error here...
    std::cerr << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0.0;
  }
}
	
float TowerType1::pmtEnergy(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _cemEnergy[ipmt];
  case CHA:
    return _chaEnergy[ipmt];
  case WHA:
    return _whaEnergy[ipmt];
  case STRIP:
    return _stripEnergy;
  case CRACK:
    return _crackEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType1::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _iCemEnergy[ipmt];
  case CHA:
    return _iChaEnergy[ipmt];
  case WHA:
    return _iWhaEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0;
  }
}
    

int TowerType1::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){
  
  switch ( d) {
  case CalorDetInfo::CEM:
    _cemEnergy[ipmt]+=energy;
    _iCemEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::CHA:
    _chaEnergy[ipmt]+=energy;
    _iChaEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::WHA:
    _whaEnergy[ipmt]+=energy;
    _iWhaEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}



int TowerType1::addPmtEnergy(Detector det , size_t ipmt,float& energy){
  
  switch ( det ) {
  case CEM:
    _cemEnergy[ipmt]+=energy;
    return 1;
  case CHA:
    _chaEnergy[ipmt]+=energy;
    return 1;
  case WHA:
    _whaEnergy[ipmt]+=energy;
    return 1;
  case STRIP:
    _stripEnergy += energy;
    return 1;
  case CRACK:
    _crackEnergy += energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType1::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    return _nChaHits;
  case WHA:
    return _nWhaHits;
  case CEM:
    return _nCemHits;
  default:
    return 0;
  }
}

int TowerType1::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case CHA:
    _iChaTime[ihit] = time;
    _chaTime[ihit]=float(time)/10-float(HTDC_Offset);
    return 1;
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit]=float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit]=float(time)/10-float(EMTDC_Offset);
    return 1;  
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType1::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    _iChaTime[ihit] = time;
    _chaTime[ihit]  = float(time)/10-float(HTDC_Offset);
    return 1;
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit]  = float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit]  = float(time)/10-float(EMTDC_Offset);
    return 1; 
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType1::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case CHA:
    _nChaHits = nhits;
    return 1;
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1;  
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType1::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d) {
  case CHA:
    _nChaHits = nhits;
    return 1;
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1;     
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType1 for detector " << d <<std::endl;
    return 0;
  }
}

TowerType1* TowerType1::clone() const{
  return new TowerType1(*this);
}
void TowerType1::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _cemEnergy[0] = (RawQieStorableBank::unpack_energy(_iCemEnergy[0]) -calib->pedCem())* calib->calibCem()*CalOffLer[0];
  _cemEnergy[1] = (RawQieStorableBank::unpack_energy(_iCemEnergy[1]) -calib->pedCem())* calib->calibCem()*CalOffLer[1];
  _chaEnergy[0] = (RawQieStorableBank::unpack_energy(_iChaEnergy[0]) -calib->pedCha())* calib->calibCha()*CalOffLer[4];
  _chaEnergy[1] = (RawQieStorableBank::unpack_energy(_iChaEnergy[1]) -calib->pedCha())* calib->calibCha()*CalOffLer[5];
  _whaEnergy[0] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[0]) -calib->pedWha())* calib->calibWha()*CalOffLer[6];
  _whaEnergy[1] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[1]) -calib->pedWha())* calib->calibWha()*CalOffLer[7];
  // just convert to float for now
  //_chaTime[0] = float(_iChaTime[0]);
  //_whaTime[0] = float(_iWhaTime[0]);
}
int TowerType1::iTower() const {
  return -1;
}
//
// TowerType 2
//

  // Equality test

bool TowerType2::operator==(const TowerType2& rhs) const
{
  if (this == &rhs) return true;
  return (_cemEnergy[0] == rhs._cemEnergy[0] &&
	  _cemEnergy[1] == rhs._cemEnergy[1] &&
	  _whaEnergy[0] == rhs._whaEnergy[0] &&
	  _whaEnergy[1] == rhs._whaEnergy[1] &&
	  _stripEnergy == rhs._stripEnergy &&
	  _crackEnergy == rhs._crackEnergy &&
	  _iWhaTime[0] == rhs._iWhaTime[0]);
}

// EDM stuff

void TowerType2::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iCemEnergy[0];
  iobuffer >>  _iCemEnergy[1];
  iobuffer >>  _iWhaEnergy[0]; 
  iobuffer >>  _iWhaEnergy[1]; 
  //  iobuffer >>  _stripEnergy;
  //  iobuffer >>  _crackEnergy;


  if (version > 6){
    iobuffer >> _nWhaHits;
    for (int i=0; i< _nWhaHits; i++) {
      if(_nWhaHits > i) iobuffer >> _iWhaTime[i]; 
    }
  } else if(version > 4){
    iobuffer >> _nWhaHits;
    if(_nWhaHits > 0) iobuffer >> _iWhaTime[0]; 
    if(_nWhaHits > 1) iobuffer >> _iWhaTime[1]; 
  }
  else{
    iobuffer >>  _iWhaTime[0]; 
  }
  for (int i=0; i< _nWhaHits; i++) {
    if (version>8){
      _whaTime[i]=float(_iWhaTime[i])/10-float(HTDC_Offset);
    } else {
      _whaTime[i]=float(_iWhaTime[i])-float(HTDC_Offset);
    }
  }
  if (version>10){
    iobuffer >> _nCemHits;
    for (int i=0; i< _nCemHits; i++) {
      if(_nCemHits > i) iobuffer >> _iCemTime[i]; 
      _cemTime[i]=float(_iCemTime[i])/10-float(EMTDC_Offset);
    }   
  }
}

void TowerType2::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iCemEnergy[0];
  iobuffer <<  _iCemEnergy[1];
  iobuffer <<  _iWhaEnergy[0]; 
  iobuffer <<  _iWhaEnergy[1]; 
  //  iobuffer <<  _stripEnergy;
  //  iobuffer <<  _crackEnergy;
 
  iobuffer << _nWhaHits;
  for (int i=0; i< _nWhaHits; i++) {
    if(_nWhaHits > i) iobuffer << _iWhaTime[i]; 
  }


  iobuffer << _nCemHits;
  for (int i=0; i< _nCemHits; i++) {
    if(_nCemHits > i) iobuffer << _iCemTime[i]; 
  }

}
  
//
// Printing  (incomplete)
//

void TowerType2::printTower(std::ostream& os) const{os << *this;}  
ostream& operator<<(std::ostream& os, const TowerType2& t)
{
  
  os << "CEM energy (0): " << setw(10) <<setprecision(5) << t._cemEnergy[0] << "   ";
  os << "CEM energy (1): " << setw(10) <<setprecision(5) << t._cemEnergy[1] << "   ";
  os << "WHA energy (0): " << setw(10) <<setprecision(5) << t._whaEnergy[0] << "   ";
  os << "WHA energy (1): " << setw(10) <<setprecision(5) << t._whaEnergy[1] <<  '\n';
  os << "CEM time   (0): " << setw(10) <<setprecision(5) << t._cemTime[0]   << "   ";
  os << "CEM time   (1): " << setw(10) <<setprecision(5) << t._cemTime[1]   << "   ";
  os << "WHA time   (0): " << setw(10) <<setprecision(5) << t._whaTime[0]   << "   ";
  os << "WHA time   (1): " << setw(10) <<setprecision(5) << t._whaTime[1]   << "   ";
  os << "CEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nCemHits)<< '\n';
  os << "WHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nWhaHits)<< '\n';
  return os;
}  

//
// Manipulation and access
//
  

float TowerType2::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _cemEnergy[ipmt];
  case WHA:
    return _whaEnergy[ipmt];
  case STRIP:
    return _stripEnergy;
  case CRACK:
    return _crackEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType2::pmtData(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case CEM:
    return _iCemEnergy[ipmt];
  case WHA:
    return _iWhaEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType2::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){
  
  switch ( d) {
  case CalorDetInfo::CEM:
    _cemEnergy[ipmt]+=energy;
    _iCemEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::WHA:
    _whaEnergy[ipmt]+=energy;
    _iWhaEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType2::addPmtEnergy(Detector det , size_t ipmt,float& energy){
  
  switch ( det ) {
  case CEM:
    _cemEnergy[ipmt]+=energy;
    return 1;
  case WHA:
    _whaEnergy[ipmt]+=energy;
    return 1;
  case STRIP:
    _stripEnergy += energy;
    return 1;
  case CRACK:
    _crackEnergy += energy;
    return 1;
  default:
    return 0;
  }
}

float TowerType2::time(CellKey& cKey) const{

  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _whaTime[0];
  case CEM:
    return _cemTime[0];
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }  
}
float TowerType2::time(CellKey& cKey,int ihit) const{

  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _whaTime[ihit];
  case CEM:
    return _cemTime[ihit];  
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }  
}
  
int TowerType2::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _nWhaHits;
  case CEM:
    return _nCemHits; 
  default:
    return 0;
  }  
}
int TowerType2::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;  
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  } 
}
int TowerType2::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case CEM:
    _iCemTime[ihit] = time;
    _cemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;  
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType2::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1; 
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType2::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case CEM:
    _nCemHits = nhits;
    return 1;  
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }
}


TowerType2* TowerType2::clone() const{
  return new TowerType2(*this);
}
void TowerType2::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _cemEnergy[0] = (RawQieStorableBank::unpack_energy(_iCemEnergy[0]) -calib->pedCem()) * calib->calibCem()*CalOffLer[0];
  _cemEnergy[1] = (RawQieStorableBank::unpack_energy(_iCemEnergy[1]) -calib->pedCem()) * calib->calibCem()*CalOffLer[1];
  _whaEnergy[0] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[0]) -calib->pedWha()) * calib->calibWha()*CalOffLer[6];
  _whaEnergy[1] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[1]) -calib->pedWha()) * calib->calibWha()*CalOffLer[7];
  // just convert to float for now
  // _whaTime[0] = float(_iWhaTime[0]);
}
int TowerType2::iTower() const {
  return -1;
}


//
// TowerType 3
//


  // Equality test

bool TowerType3::operator==(const TowerType3& rhs) const
{
  if (this == &rhs) return true;
  return (_pemEnergy[0] == rhs._pemEnergy[0] &&
	  _pemEnergy[1] == rhs._pemEnergy[1] &&
	  _whaEnergy[0] == rhs._whaEnergy[0] &&
	  _whaEnergy[1] == rhs._whaEnergy[1] &&
	  _pprEnergy[0] == rhs._pprEnergy[0] &&
	  _pprEnergy[1] == rhs._pprEnergy[1] &&
	  _iWhaTime[0] == rhs._iWhaTime[0]);
}

// EDM stuff

void TowerType3::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iPemEnergy[0];
  iobuffer >>  _iPemEnergy[1];
  iobuffer >>  _iWhaEnergy[0]; 
  iobuffer >>  _iWhaEnergy[1]; 
  iobuffer >>  _iPprEnergy[0]; 
  iobuffer >>  _iPprEnergy[1];


  if (version > 6){
    iobuffer >> _nWhaHits;
    for (int i=0; i< _nWhaHits; i++) {
      if(_nWhaHits > i) iobuffer >> _iWhaTime[i]; 
    }
  } else if(version > 4){
    
    iobuffer >> _nWhaHits;
    if(_nWhaHits > 0) iobuffer >> _iWhaTime[0]; 
    if(_nWhaHits > 1) iobuffer >> _iWhaTime[1]; 
  }
  else{
    iobuffer >>  _iWhaTime[0]; 
  } 
  for (int i=0; i< _nWhaHits; i++) {
    if (version>8){
      _whaTime[i]=float(_iWhaTime[i])/10-float(HTDC_Offset);
    } else {
      _whaTime[i]=float(_iWhaTime[i])-float(HTDC_Offset);
    }
  }
  if (version>10){
    iobuffer >> _nPemHits;
    for (int i=0; i< _nPemHits; i++) {
      if(_nPemHits > i) iobuffer >> _iPemTime[i]; 
      _pemTime[i]=float(_iPemTime[i])/10-float(EMTDC_Offset);
    }   
  }
  
}

void TowerType3::writeToBuffer(TBuffer& iobuffer,Version_t version){
  
  iobuffer <<  _iPemEnergy[0];
  iobuffer <<  _iPemEnergy[1];
  iobuffer <<  _iWhaEnergy[0]; 
  iobuffer <<  _iWhaEnergy[1]; 
  iobuffer <<  _iPprEnergy[0]; 
  iobuffer <<  _iPprEnergy[1];

  iobuffer << _nWhaHits;
  for (int i=0; i< _nWhaHits; i++) {
    if(_nWhaHits > i) iobuffer << _iWhaTime[i]; 
  }

  iobuffer << _nPemHits;
  for (int i=0; i< _nPemHits; i++) {
    if(_nPemHits > i) iobuffer << _iPemTime[i]; 
  }
}

//
// Printing  (incomplete)
//

void TowerType3::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType3& t)
{

  os << "PPR energy (0): " << setw(10) <<setprecision(5) << t._pprEnergy[0] << "   ";
  os << "PPR energy (1): " << setw(10) <<setprecision(5) << t._pprEnergy[1] << "   ";
  os << "PEM energy (0): " << setw(10) <<setprecision(5) << t._pemEnergy[0] << "   ";
  os << "PEM energy (1): " << setw(10) <<setprecision(5) << t._pemEnergy[1] << "   ";
  os << "WHA energy (0): " << setw(10) <<setprecision(5) << t._whaEnergy[0] << "   ";
  os << "WHA energy (1): " << setw(10) <<setprecision(5) << t._whaEnergy[1] << '\n';
  os << "WHA time   (0): " << setw(10) <<setprecision(5) << t._whaTime[0]   << "   ";
  os << "WHA time   (1): " << setw(10) <<setprecision(5) << t._whaTime[1]   << "   ";
  os << "PEM time   (0): " << setw(10) <<setprecision(5) << t._pemTime[0]   << "   ";
  os << "PEM time   (1): " << setw(10) <<setprecision(5) << t._pemTime[1]   << "   ";
  os << "WHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nWhaHits)<< '\n';
  os << "PEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPemHits)<< '\n';

  return os;
}
  
  
 
//
// Manipulation and access
//
  
	
float TowerType3::pmtEnergy(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _pemEnergy[ipmt];
  case WHA:
    return _whaEnergy[ipmt];
  case PPR:
    return _pprEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0.0;
  }
}
  
uint2 TowerType3::pmtData(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _iPemEnergy[ipmt];
  case WHA:
    return _iWhaEnergy[ipmt];
  case PPR:
    return _iPprEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType3::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

    
  switch ( d) {
  case CalorDetInfo::PEM:
    _pemEnergy[ipmt]+=energy;
    _iPemEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::WHA:
    _whaEnergy[ipmt]+=energy;
    _iWhaEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::PPR:
    _pprEnergy[ipmt] += energy;
    _iPprEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}


int TowerType3::addPmtEnergy(Detector det , size_t ipmt,float& energy){

  switch ( det ) {
  case PEM:
    _pemEnergy[ipmt]+=energy;
    return 1;
  case WHA:
    _whaEnergy[ipmt]+=energy;
    return 1;
  case PPR:
    _pprEnergy[ipmt] += energy;
    return 1;
  default:
    return 0;
  }
}

float TowerType3::time(CellKey& cKey) const{ 
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _whaTime[0];
  case PEM:
    return _pemTime[0]; 
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType2 for detector " << d <<std::endl;
    return 0;
  }  
}
float TowerType3::time(CellKey& cKey,int ihit) const{

  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _whaTime[ihit];
  case PEM:
    return _pemTime[ihit];  
  default:
    // should throw an error here...
    std::cout << "time: Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  }  
}

int TowerType3::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    return _nWhaHits;
  case PEM:
    return _nPemHits;     
  default:
    return 0;
  }  
}
int TowerType3::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case PEM:
    _iPemTime[ihit] = time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  } 
}
int TowerType3::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case PEM:
    _iPemTime[ihit] = time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;
  default:
    // should throw an error here...
    std::cout << "setTime: Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType3::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1;   
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType3::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d) {
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1;  
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType3 for detector " << d <<std::endl;
    return 0;
  }
}

TowerType3* TowerType3::clone() const{
  return new TowerType3(*this);
}
void TowerType3::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _pemEnergy[0] = (RawQieStorableBank::unpack_energy(_iPemEnergy[0]) -calib->pedPem()) * calib->calibPem()*CalOffLer[2];
  _pemEnergy[1] = (RawQieStorableBank::unpack_energy(_iPemEnergy[1]) -calib->pedPem()) * calib->calibPem()*CalOffLer[3];
  _pprEnergy[0] = (RawQieStorableBank::unpack_energy(_iPprEnergy[0]) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[10];
  _pprEnergy[1] = (RawQieStorableBank::unpack_energy(_iPprEnergy[1]) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[11];
  _whaEnergy[0] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[0]) -calib->pedWha()) * calib->calibWha()*CalOffLer[6];
  _whaEnergy[1] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[1]) -calib->pedWha()) * calib->calibWha()*CalOffLer[7];
  // just convert to float for now
  // _whaTime[0] = float(_iWhaTime[0]);
}
int TowerType3::iTower() const {
  return -1;
}

//
// TowerType 4
//


  // Equality test

bool TowerType4::operator==(const TowerType4& rhs) const
{
  if (this == &rhs) return true;
  return (_pemEnergy[0] == rhs._pemEnergy[0] &&
	  _pemEnergy[1] == rhs._pemEnergy[1] &&
	  _whaEnergy[0] == rhs._whaEnergy[0] &&
	  _whaEnergy[1] == rhs._whaEnergy[1] &&
	  _phaEnergy[0] == rhs._phaEnergy[0] &&
	  _phaEnergy[1] == rhs._phaEnergy[1] &&
	  _pprEnergy[0] == rhs._pprEnergy[0] &&
	  _pprEnergy[1] == rhs._pprEnergy[1] &&
	  _iPhaTime[0] == rhs._iPhaTime[0] &&
	  _iWhaTime[0] == rhs._iWhaTime[0] );
}

// EDM stuff

void TowerType4::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iPemEnergy[0];
  iobuffer >>  _iPemEnergy[1];
  iobuffer >>  _iPhaEnergy[0];
  iobuffer >>  _iPhaEnergy[1];
  iobuffer >>  _iWhaEnergy[0]; 
  iobuffer >>  _iWhaEnergy[1]; 
  iobuffer >>  _iPprEnergy[0]; 
  iobuffer >>  _iPprEnergy[1];

  if (version > 6){
    iobuffer >> _nWhaHits;
    for (int i=0; i< _nWhaHits; i++) {
      uint2 dummy=0;
      if(_nWhaHits > i) {
	iobuffer >> dummy;
	if (i<8) _iWhaTime[i]=dummy; 
      }
    }
    iobuffer >> _nPhaHits;
    for (int i=0; i< _nPhaHits; i++) {
      uint2 dummy=0;
      if(_nPhaHits > i) {
	iobuffer >> dummy;
	if (i<8) _iPhaTime[i]=dummy; 
      }
    }
  } else if(version > 4){
   
    iobuffer >> _nWhaHits;
    if(_nWhaHits > 0) iobuffer >> _iWhaTime[0]; 
    if(_nWhaHits > 1) iobuffer >> _iWhaTime[1]; 
    iobuffer >> _nPhaHits;
    if(_nPhaHits > 0) iobuffer >> _iPhaTime[0]; 
    if(_nPhaHits > 1) iobuffer >> _iPhaTime[1];

  }
  else{
    iobuffer >>  _iWhaTime[0]; 
    if(version > 3) iobuffer >> _iPhaTime[0];
  }
  for (int i=0; i< _nWhaHits; i++) {
    if (version>8){
      _whaTime[i]=float(_iWhaTime[i])/10-float(HTDC_Offset);
    } else {
      _whaTime[i]=float(_iWhaTime[i])-float(HTDC_Offset);
    }
  }
  for (int i=0; i< _nPhaHits; i++) {
    if (version>8){
      _phaTime[i]=float(_iPhaTime[i])/10-float(HTDC_Offset);
    } else {
      _phaTime[i]=float(_iPhaTime[i])-float(HTDC_Offset);
    }
  }
  if (version>10){
    iobuffer >> _nPemHits;
    for (int i=0; i< _nPemHits; i++) {
      if(_nPemHits > i) iobuffer >> _iPemTime[i]; 
      _pemTime[i]=float(_iPemTime[i])/10-float(EMTDC_Offset);
    }   
  }
}

void TowerType4::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iPemEnergy[0];
  iobuffer <<  _iPemEnergy[1];
  iobuffer <<  _iPhaEnergy[0];
  iobuffer <<  _iPhaEnergy[1];
  iobuffer <<  _iWhaEnergy[0]; 
  iobuffer <<  _iWhaEnergy[1]; 
  iobuffer <<  _iPprEnergy[0]; 
  iobuffer <<  _iPprEnergy[1];
  iobuffer << _nWhaHits;
  for (int i=0; i< _nWhaHits; i++) {
    if(_nWhaHits > i) iobuffer << _iWhaTime[i]; 
  }
  iobuffer << _nPhaHits;
  for (int i=0; i< _nPhaHits; i++) {
    if(_nPhaHits > i) iobuffer << _iPhaTime[i]; 
  }

  iobuffer << _nPemHits;
  for (int i=0; i< _nPemHits; i++) {
    if(_nPemHits > i) iobuffer << _iPemTime[i]; 
  }
}

//
// Printing  (incomplete)
//

void TowerType4::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType4& t)
{
  os << "PPR energy (0): " << setw(10) <<setprecision(5) << t._pprEnergy[0] << "   ";
  os << "PPR energy (1): " << setw(10) <<setprecision(5) << t._pprEnergy[1] << "   ";
  os << "WHA energy (0): " << setw(10) <<setprecision(5) << t._whaEnergy[0] << "   ";
  os << "WHA energy (1): " << setw(10) <<setprecision(5) << t._whaEnergy[1] << '\n';
  os << "PEM energy (0): " << setw(10) <<setprecision(5) << t._pemEnergy[0] << "   ";
  os << "PEM energy (1): " << setw(10) <<setprecision(5) << t._pemEnergy[1] << "   ";
  os << "PHA energy (0): " << setw(10) <<setprecision(5) << t._phaEnergy[0] << "   ";
  os << "PHA energy (1): " << setw(10) <<setprecision(5) << t._phaEnergy[1] <<  '\n';
  os << "WHA time   (0): " << setw(10) <<setprecision(5) << t._whaTime[0]   << "   ";
  os << "WHA time   (1): " << setw(10) <<setprecision(5) << t._whaTime[1]   << "   ";
  os << "PHA time   (0): " << setw(10) <<setprecision(5) << t._phaTime[0]   << "   ";
  os << "PHA time   (1): " << setw(10) <<setprecision(5) << t._phaTime[1]   << "   ";
  os << "PEM time   (0): " << setw(10) <<setprecision(5) << t._pemTime[0]   << "   ";
  os << "PEM time   (1): " << setw(10) <<setprecision(5) << t._pemTime[1]   << "   ";
  os << "WHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nWhaHits)<< "   ";
  os << "PHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPhaHits)<< '\n';
  os << "PEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPemHits)<< '\n';
  return os;
}  
  
//
// Manipulation and access
//
  

float TowerType4::pmtEnergy(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _pemEnergy[ipmt];
  case PHA:
    return _phaEnergy[ipmt];
  case WHA:
    return _whaEnergy[ipmt];
  case PPR:
    return _pprEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType4::pmtData(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _iPemEnergy[ipmt];
  case PHA:
    return _iPhaEnergy[ipmt];
  case WHA:
    return _iWhaEnergy[ipmt];
  case PPR:
    return _iPprEnergy[ipmt];
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0;
  }
}
    

int TowerType4::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::PEM:
    _pemEnergy[ipmt]+=energy;
    _iPemEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::PHA:
    _phaEnergy[ipmt]+=energy;
    _iPhaEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::WHA:
    _whaEnergy[ipmt]+=energy;
    _iWhaEnergy[ipmt] = iData;
    return 1;
  case CalorDetInfo::PPR:
    _pprEnergy[ipmt] += energy;
    _iPprEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType4::addPmtEnergy(Detector det , size_t ipmt,float& energy){

  switch ( det ) {
  case PEM:
    _pemEnergy[ipmt] += energy;
    return 1;
  case PHA:
    _phaEnergy[ipmt] += energy;
    return 1;
  case WHA:
    _whaEnergy[ipmt] += energy;
    return 1;
  case PPR:
    _pprEnergy[ipmt] += energy;
    return 1;
  default:
    return 0;
  }
}

float TowerType4::time(CellKey& cKey) const{
  Detector d = cKey.detector();
  switch ( d) {
  case PHA:
    return _phaTime[0];
  case WHA:
    return _whaTime[0];
  case PEM:
    return _pemTime[0];
  default:
    // should throw an error here...
    std::cerr << "time: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0.0;
  }
}
float TowerType4::time(CellKey& cKey,int ihit) const{ 
  Detector d = cKey.detector();
  switch ( d) {
  case PHA:
    return _phaTime[ihit];
  case WHA:
    return _whaTime[ihit];
  case PEM:
    return _pemTime[ihit]; 
  default:
    // should throw an error here...
    std::cerr << "time: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0.0;
  }
}	

int TowerType4::nTDCHits(CellKey& cKey) const{                                

  Detector d = cKey.detector();
  switch ( d) {
  case PHA:
    return _nPhaHits;
  case WHA:
    return _nWhaHits;
  case PEM:
    return _nPemHits;  
  default:
    return 0;
  }
}

int TowerType4::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case PHA:
    _iPhaTime[ihit] = time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case PEM:
    _iPemTime[ihit] = time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;
  default:
    std::cerr << "setTime: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0;
  }

}

int TowerType4::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d) {
  case PHA:
    _iPhaTime[ihit] = time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    break;
  case WHA:
    _iWhaTime[ihit] = time;
    _whaTime[ihit] = float(time)/10-float(HTDC_Offset);
    break;
  case PEM:
    _iPemTime[ihit] = time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    break; 
  default:
    std::cerr << "setTime: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0;
  }
  
  return 1;
}
int TowerType4::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1; 
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType4::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  case WHA:
    _nWhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1; 
  default:
    // should throw an error here...
    std::cout << "setTDCHits: Bad Key in TowerType4 for detector " << d <<std::endl;
    return 0;
  }
}
TowerType4* TowerType4::clone() const{
  return new TowerType4(*this);
}
void TowerType4::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _pemEnergy[0] = (RawQieStorableBank::unpack_energy(_iPemEnergy[0]) -calib->pedPem()) * calib->calibPem()*CalOffLer[2];
  _pemEnergy[1] = (RawQieStorableBank::unpack_energy(_iPemEnergy[1]) -calib->pedPem()) * calib->calibPem()*CalOffLer[3];
  _pprEnergy[0] = (RawQieStorableBank::unpack_energy(_iPprEnergy[0]) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[10];
  _pprEnergy[1] = (RawQieStorableBank::unpack_energy(_iPprEnergy[1]) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[11];
  _phaEnergy[0] = (RawQieStorableBank::unpack_energy(_iPhaEnergy[0]) -calib->pedPha()) * calib->calibPha()*CalOffLer[8];
  _phaEnergy[1] = (RawQieStorableBank::unpack_energy(_iPhaEnergy[1]) -calib->pedPha()) * calib->calibPha()*CalOffLer[9];
  _whaEnergy[0] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[0]) -calib->pedWha()) * calib->calibWha()*CalOffLer[6];
  _whaEnergy[1] = (RawQieStorableBank::unpack_energy(_iWhaEnergy[1]) -calib->pedWha()) * calib->calibWha()*CalOffLer[7];
  // just convert to float for now
  // _whaTime[0] = float(_iWhaTime[0]);
}
int TowerType4::iTower() const {
  return -1;
}
//
// TowerType 5
//

  // Equality test

bool TowerType5::operator==(const TowerType5& rhs) const
{
  if (this == &rhs) return true;
  return (_pemEnergy == rhs._pemEnergy &&
	  _phaEnergy == rhs._phaEnergy &&
	  _pprEnergy == rhs._pprEnergy);
}

// EDM stuff

void TowerType5::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iPemEnergy;
  iobuffer >>  _iPhaEnergy;
  iobuffer >>  _iPprEnergy;
  
  if (version > 6){
    iobuffer >> _nPhaHits;
    for (int i=0; i< _nPhaHits; i++) {
      if(_nPhaHits > i) iobuffer >> _iPhaTime[i]; 
    }
  } else if(version > 4){   
    iobuffer >> _nPhaHits;
    if(_nPhaHits > 0) iobuffer >> _iPhaTime[0]; 
    if(_nPhaHits > 1) iobuffer >> _iPhaTime[1];
    
  }
  else{
    if(version > 3) iobuffer >> _iPhaTime[0];
  }
  for (int i=0; i< _nPhaHits; i++) {
    if (version>8){
      _phaTime[i]=float(_iPhaTime[i])/10-float(HTDC_Offset);
    } else {
      _phaTime[i]=float(_iPhaTime[i])-float(HTDC_Offset);
    }
  }
  if (version>10){
    iobuffer >> _nPemHits;
    for (int i=0; i< _nPemHits; i++) {
      if(_nPemHits > i) iobuffer >> _iPemTime[i]; 
      _pemTime[i]=float(_iPemTime[i])/10-float(EMTDC_Offset);
    }   
  }

}

void TowerType5::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iPemEnergy;
  iobuffer <<  _iPhaEnergy;
  iobuffer <<  _iPprEnergy;
  iobuffer << _nPhaHits;
  for (int i=0; i< _nPhaHits; i++) {
    if(_nPhaHits > i) iobuffer << _iPhaTime[i]; 
  }

  iobuffer << _nPemHits;
  for (int i=0; i< _nPemHits; i++) {
    if(_nPemHits > i) iobuffer << _iPemTime[i]; 
  }
}

//
// Printing  (incomplete)
//

void TowerType5::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType5& t)
{
  //    float time0=float(t._iPhaTime[0])/10-float(HTDC_Offset);
  //    float time1=float(t._iPhaTime[0])/10-float(HTDC_Offset);
  //    os << "my time :" << t._iPhaTime[0] << "    " << float(t._iPhaTime[0])/10 << "  " << float(t._iPhaTime[0])/10-float(HTDC_Offset) <<'\n';
  os << "PPR energy    : " << setw(10) <<setprecision(5) << t._pprEnergy << "   ";
  os << "PEM energy    : " << setw(10) <<setprecision(5) << t._pemEnergy << "   ";
  os << "PHA energy    : " << setw(10) <<setprecision(5) << t._phaEnergy << '\n';
  os << "PHA time   (0): " << setw(10) <<setprecision(5) << t._phaTime[0]   << "   ";
  os << "PHA time   (1): " << setw(10) <<setprecision(5) << t._phaTime[1]   << "   ";
  os << "PHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPhaHits)<< '\n';
  os << "PEM time   (0): " << setw(10) <<setprecision(5) << t._pemTime[0]   << "   ";
  os << "PEM time   (1): " << setw(10) <<setprecision(5) << t._pemTime[1]   << "   ";
  os << "PEM TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPemHits)<< '\n';
  return os;
}

//
// Manipulation and access
//
  
	
float TowerType5::pmtEnergy(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _pemEnergy;
  case PHA:
    return _phaEnergy;
  case PPR:
    return _pprEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType5::pmtData(CellKey& cKey, size_t ipmt) const{
    
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _iPemEnergy;
  case PHA:
    return _iPhaEnergy;
  case PPR:
    return _iPprEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType5::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::PEM:
    _pemEnergy += energy;
    _iPemEnergy = iData;
    return 1;
  case CalorDetInfo::PHA:
    _phaEnergy += energy;
    _iPhaEnergy = iData;
    return 1;
  case CalorDetInfo::PPR:
    _pprEnergy += energy;
    _iPprEnergy = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType5::addPmtEnergy(Detector det , size_t ipmt,float& energy){
  switch ( det ) {
  case PEM:
    _pemEnergy += energy;
    return 1;
  case PHA:
    _phaEnergy += energy;
    return 1;
  case PPR:
    _pprEnergy += energy;
    return 1;
  default:
    return 0;
  }
}
  
float TowerType5::time(CellKey& cKey) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[0];
  case PEM:
    return _pemTime[0]; 
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}
float TowerType5::time(CellKey& cKey,int ihit) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[ihit];
  case PEM:
    return _pemTime[ihit];  
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType5::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _nPhaHits;
  case PEM:
    return _nPemHits;   
  default:
    return 0;
  }
}

int TowerType5::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] = time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case PEM:
    _iPemTime[ihit] = time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1; 
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType5::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] =time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  case PEM:
    _iPemTime[ihit] =time;
    _pemTime[ihit] = float(time)/10-float(EMTDC_Offset);
    return 1;  
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType5::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1;  
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType5::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  case PEM:
    _nPemHits = nhits;
    return 1; 
  default:
    std::cout << "Bad Key in TowerType5 for detector " << d <<std::endl;
    return 0;
  }
}


TowerType5* TowerType5::clone() const{
  return new TowerType5(*this);
}
void TowerType5::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _pemEnergy = (RawQieStorableBank::unpack_energy(_iPemEnergy) -calib->pedPem()) * calib->calibPem()*CalOffLer[2];
  _pprEnergy = (RawQieStorableBank::unpack_energy(_iPprEnergy) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[10];
  _phaEnergy = (RawQieStorableBank::unpack_energy(_iPhaEnergy) -calib->pedPha()) * calib->calibPha()*CalOffLer[8];
}
int TowerType5::iTower() const {
  return -1;
}
//
// TowerType 6
//


  // Equality test

bool TowerType6::operator==(const TowerType6& rhs) const
{
  if (this == &rhs) return true;
  return (_pemEnergy == rhs._pemEnergy &&
	  _phaEnergy == rhs._phaEnergy &&
	  _pprEnergy == rhs._pprEnergy);
}

// EDM stuff

void TowerType6::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iPemEnergy;
  iobuffer >>  _iPhaEnergy;
  iobuffer >>  _iPprEnergy;


  if (version > 6){
    iobuffer >> _nPhaHits;
    for (int i=0; i< _nPhaHits; i++) {
      if(_nPhaHits > i) iobuffer >> _iPhaTime[i]; 
    }
  } else if(version > 4){   
    iobuffer >> _nPhaHits;
    if(_nPhaHits > 0) iobuffer >> _iPhaTime[0]; 
    if(_nPhaHits > 1) iobuffer >> _iPhaTime[1];
    
  }
  else{
    if(version > 3) iobuffer >> _iPhaTime[0];
  }
  for (int i=0; i< _nPhaHits; i++) {
    if (version>8){
      _phaTime[i]=float(_iPhaTime[i])/10-float(HTDC_Offset);
    } else {
      _phaTime[i]=float(_iPhaTime[i])-float(HTDC_Offset);
    }
  }
}

void TowerType6::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iPemEnergy;
  iobuffer <<  _iPhaEnergy;
  iobuffer <<  _iPprEnergy;
  iobuffer << _nPhaHits;
  for (int i=0; i< _nPhaHits; i++) {
    if(_nPhaHits > i) iobuffer << _iPhaTime[i]; 
  }
}

//
// Printing  (incomplete)
//

void TowerType6::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType6& t){
 
  os << "PPR energy    : " << setw(10) <<setprecision(5) << t._pprEnergy << "   ";
  os << "PEM energy    : " << setw(10) <<setprecision(5) << t._pemEnergy << "   ";
  os << "PHA energy    : " << setw(10) <<setprecision(5) << t._phaEnergy << '\n';
  os << "PHA time   (0): " << setw(10) <<setprecision(5) << t._phaTime[0]   << "   ";
  os << "PHA time   (1): " << setw(10) <<setprecision(5) << t._phaTime[1]   << "   ";
  os << "PHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPhaHits)<< '\n';
  return os;
}
   
//
// Manipulation and access
//
  
float TowerType6::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _pemEnergy;
  case PHA:
    return _phaEnergy;
  case PPR:
    return _pprEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0.0;
  }
}
	
uint2 TowerType6::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _iPemEnergy;
  case PHA:
    return _iPhaEnergy;
  case PPR:
    return _iPprEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType6::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){
  switch ( d) {
  case CalorDetInfo::PEM:
    _pemEnergy+=energy;
    _iPemEnergy = iData;

    return 1;
  case CalorDetInfo::PHA:
    _phaEnergy+=energy;
    _iPhaEnergy = iData;

    return 1;
  case CalorDetInfo::PPR:
    _pprEnergy+= energy;
    _iPprEnergy = iData;

    return 1;
  default:
    return 0;
  }
}


int TowerType6::addPmtEnergy(Detector det , size_t ipmt,float& energy){
  switch ( det ) {
  case PEM:
    _pemEnergy+=energy;
    return 1;
  case PHA:
    _phaEnergy+=energy;
    return 1;
  case PPR:
    _pprEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

float TowerType6::time(CellKey& cKey) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[0];
  default:
    std::cout << "time: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}
float TowerType6::time(CellKey& cKey,int ihit) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[ihit];
  default:
    std::cout << "time: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType6::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _nPhaHits;
  default:
    return 0;
  }
}

int TowerType6::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] = time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  default:
    std::cout << "setTime: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType6::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] =time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  default:
    std::cout << "setTime: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType6::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  default:
    std::cout << "setTDCHits: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType6::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  default:
    std::cout << "setTDCHits: Bad Key in TowerType6 for detector " << d <<std::endl;
    return 0;
  }
}

TowerType6* TowerType6::clone() const{
  return new TowerType6(*this);
}
void TowerType6::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  _pemEnergy = (RawQieStorableBank::unpack_energy(_iPemEnergy) -calib->pedPem()) * calib->calibPem()*CalOffLer[2];
  _pprEnergy = (RawQieStorableBank::unpack_energy(_iPprEnergy) -calib->pedPpr()) * calib->calibPpr()*CalOffLer[10];
  _phaEnergy = (RawQieStorableBank::unpack_energy(_iPhaEnergy) -calib->pedPha()) * calib->calibPha()*CalOffLer[8];

}
int TowerType6::iTower() const {
  return -1;
}  
//
// TowerType 7
//

// Equality test

bool TowerType7::operator==(const TowerType7& rhs) const
{
  if (this == &rhs) return true;
  return (_pemEnergy == rhs._pemEnergy &&
	  _phaEnergy == rhs._phaEnergy &&
	  _pprEnergy == rhs._pprEnergy &&
	  _mpaEnergy == rhs._mpaEnergy );
}

// EDM stuff

void TowerType7::readFromBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer >>  _iPemEnergy;
  iobuffer >>  _iPhaEnergy;
  iobuffer >>  _iPprEnergy;
  if (version<8) {
    iobuffer >> _iMpaEnergy; 
    uint2 dummy;
    iobuffer >> dummy;
  } else {
    iobuffer >> _iMpaEnergy; 
  }
  if (version > 6){
    iobuffer >> _nPhaHits;
    for (int i=0; i< _nPhaHits; i++) {
      if(_nPhaHits > i) iobuffer >> _iPhaTime[i]; 
    }
  } else if(version > 4){   
    iobuffer >> _nPhaHits;
    if(_nPhaHits > 0) iobuffer >> _iPhaTime[0]; 
    if(_nPhaHits > 1) iobuffer >> _iPhaTime[1];
    
  }
  else{
    if(version > 3) iobuffer >> _iPhaTime[0];
  }
  for (int i=0; i< _nPhaHits; i++) {
    if (version>8){
      _phaTime[i]=float(_iPhaTime[i])/10-float(HTDC_Offset);
    } else {
      _phaTime[i]=float(_iPhaTime[i])-float(HTDC_Offset);
    }
  }
}

void TowerType7::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer <<  _iPemEnergy;
  iobuffer <<  _iPhaEnergy;
  iobuffer <<  _iPprEnergy;
  iobuffer <<  _iMpaEnergy; 
  iobuffer <<  _nPhaHits;
  for (int i=0; i< _nPhaHits; i++) {
    if(_nPhaHits > i) iobuffer << _iPhaTime[i]; 
  }
}

//
// Printing  (incomplete)
//

void TowerType7::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType7& t)
{

  os << "PPR energy    : " << setw(10) <<setprecision(5) << t._pprEnergy << "   ";
  os << "PEM energy    : " << setw(10) <<setprecision(5) << t._pemEnergy << "   ";
  os << "PHA energy    : " << setw(10) <<setprecision(5) << t._phaEnergy << '\n';
  os << "PHA time   (0): " << setw(10) <<setprecision(5) << t._phaTime[0]   << "   ";
  os << "PHA time   (1): " << setw(10) <<setprecision(5) << t._phaTime[1]   << "   ";
  os << "PHA TDC Hits  : " << setw(10) <<setprecision(0) << uint2(t._nPhaHits)<< '\n';
  return os;
}
   
 
//
// Manipulation and access
//
  
float TowerType7::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _pemEnergy;
  case PHA:
    return _phaEnergy;
  case PPR:
    return _pprEnergy;
  case MPA:
    return _mpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType7::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case PEM:
    return _iPemEnergy;
  case PHA:
    return _iPhaEnergy;
  case PPR:
    return _iPprEnergy;
  case MPA:
    return _iMpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType7::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::PEM:
    _pemEnergy+=energy;
    _iPemEnergy = iData;
    return 1;
  case CalorDetInfo::PHA:
    _phaEnergy+=energy;
    _iPhaEnergy = iData;
    return 1;
  case CalorDetInfo::PPR:
    _pprEnergy+= energy;
    _iPprEnergy = iData;
    return 1;
    // do not sum MiniPlug energy in for TowerType 7 yet (Beate, 05/06/02)
  case MPA:
    //	_mpaEnergy[ipmt]+= energy;
    //	_iMpaEnergy[ipmt] = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType7::addPmtEnergy(Detector det , size_t ipmt,float& energy){
  switch ( det) {
  case PEM:
    _pemEnergy+=energy;
    return 1;
  case PHA:
    _phaEnergy+=energy;
    return 1;
  case PPR:
    _pprEnergy+= energy;
    return 1;
  case MPA:
    //	_mpaEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

float TowerType7::time(CellKey& cKey) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[0];
  default:
    std::cout << "time: Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}
float TowerType7::time(CellKey& cKey,int ihit) const{ 
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _phaTime[ihit];
  default:
    std::cout << "time : Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType7::nTDCHits(CellKey& cKey) const{                                
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    return _nPhaHits;
  default:
    return 0;
  }
}

int TowerType7::setTime(CalorKey& bKey,int time, int ihit){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] = time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  default:
    std::cout << "setTime: Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}

int TowerType7::setTime(CellKey& cKey,int time, int ihit){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _iPhaTime[ihit] =time;
    _phaTime[ihit] = float(time)/10-float(HTDC_Offset);
    return 1;
  default:
    std::cout << "setTime: Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType7::setTDCHits(CalorKey& bKey,uint1 nhits){
  CalorDetInfo::detectorTag d = bKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  default:
    std::cout << "setTDCHits: Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}
int TowerType7::setTDCHits(CellKey& cKey,uint1 nhits){
  Detector d = cKey.detector();
  switch ( d ) {
  case PHA:
    _nPhaHits = nhits;
    return 1;
  default:
    std::cout << "setTDCHits: Bad Key in TowerType7 for detector " << d <<std::endl;
    return 0;
  }
}

TowerType7* TowerType7::clone() const{
  return new TowerType7(*this);
}
void TowerType7::puffTower(const CalDataCalib* calib,float CalOffLer[]){

 
  _pemEnergy = (RawQieStorableBank::unpack_energy(_iPemEnergy) -calib->pedPem()) * calib->calibPem()*CalOffLer[2];
  _pprEnergy = (RawQieStorableBank::unpack_energy(_iPprEnergy) -calib->pedPpr() )* calib->calibPpr()*CalOffLer[10];
  _phaEnergy = (RawQieStorableBank::unpack_energy(_iPhaEnergy) -calib->pedPha()) * calib->calibPha()*CalOffLer[8];
  _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * calib->calibMpa();  

}
int TowerType7::iTower() const {
  return -1;
}  
//
// TowerType 8
//


  // Equality test

bool TowerType8::operator==(const TowerType8& rhs) const
{
  if (this == &rhs) return true;
  return (_mpaEnergy == rhs._mpaEnergy);
}

// EDM stuff

void TowerType8::readFromBuffer(TBuffer& iobuffer,Version_t version){
  
  if (version<8) {
    iobuffer >> _iMpaEnergy; 
    uint2 dummy;
    iobuffer >> dummy;
  } else {
    iobuffer >> _iMpaEnergy; 
  }
}


void TowerType8::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer << _iMpaEnergy; 
}

//
// Printing  (incomplete)
//

void TowerType8::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType8& t)
{

  os << "MPA energy     : " <<setw(10) <<setprecision(5) << t._mpaEnergy << '\n';

  return os;
}
    
//
// Manipulation and access
//
  
float TowerType8::time(CellKey& cKey) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}

float TowerType8::time(CellKey& cKey, int ihit) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}
	
float TowerType8::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _mpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType8 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType8::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _iMpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType8 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType8::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::MPA:
    _mpaEnergy+= energy;
    _iMpaEnergy = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType8::addPmtEnergy(Detector det , size_t ipmt,float& energy){  

  switch ( det) {
  case MPA:
    _mpaEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType8::nTDCHits(CellKey& cKey) const{                                
  return 0;
}

int TowerType8::setTime(CalorKey& bKey,int time, int ihit){
  
  return 0;
}
int TowerType8::setTime(CellKey& cKey,int time, int ihit){
  
  return 0;
}
int TowerType8::setTDCHits(CalorKey& bKey,uint1 nhits){

  return 0;
}
int TowerType8::setTDCHits(CellKey& cKey,uint1 nhits){

  return 0;
}
TowerType8* TowerType8::clone() const{
  return new TowerType8(*this);
}

void TowerType8::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  if (CalOffLer[12] == -1.) 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * calib->calibMpa();  
  else 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * CalOffLer[12];  
  if ( _mpaEnergy < 0.)  _mpaEnergy = 0.;
}
int TowerType8::iTower() const {
  int phi = this->iPhi();
  int wedge = phi / 5;
  int phiInWedge = phi % 5;
  int tower = 9 + 14 * wedge + phiInWedge;
  return tower;
}  
//
// TowerType 9
//


  // Equality test

bool TowerType9::operator==(const TowerType9& rhs) const
{
  if (this == &rhs) return true;
  return (_mpaEnergy == rhs._mpaEnergy);
}

// EDM stuff

void TowerType9::readFromBuffer(TBuffer& iobuffer,Version_t version){
  
  if (version<8) {
    iobuffer >> _iMpaEnergy; 
    uint2 dummy;
    iobuffer >> dummy;
  } else {
    iobuffer >> _iMpaEnergy; 
  }
}


void TowerType9::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer << _iMpaEnergy; 
}

//
// Printing  (incomplete)
//

void TowerType9::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType9& t)
{
  os << "MPA energy     : " << setw(10) <<setprecision(5) << t._mpaEnergy << '\n';


  return os;
}
    
//
// Manipulation and access
//
  
float TowerType9::time(CellKey& cKey) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}
 
float TowerType9::time(CellKey& cKey, int ihit) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}
	
float TowerType9::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _mpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType9 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType9::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _iMpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType9 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType9::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::MPA:
    _mpaEnergy+= energy;
    _iMpaEnergy = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType9::addPmtEnergy(Detector det , size_t ipmt,float& energy){  

  switch ( det) {
  case MPA:
    _mpaEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType9::nTDCHits(CellKey& cKey) const{                                
  return 0;
}

int TowerType9::setTime(CalorKey& bKey,int time, int ihit){
  
  return 0;
}
int TowerType9::setTime(CellKey& cKey,int time, int ihit){
  
  return 0;
}
int TowerType9::setTDCHits(CalorKey& bKey,uint1 nhits){

  return 0;
}
int TowerType9::setTDCHits(CellKey& cKey,uint1 nhits){

  return 0;
}
TowerType9* TowerType9::clone() const{
  return new TowerType9(*this);
}

void TowerType9::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  if (CalOffLer[12] == -1.) 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * calib->calibMpa();  
  else 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * CalOffLer[12];  
  if ( _mpaEnergy < 0.)  _mpaEnergy = 0.;
}
int TowerType9::iTower() const {
  int phi = this->iPhi();
  int wedge = phi / 4;
  int phiInWedge = phi % 4;
  int tower = 5 + 14 * wedge + phiInWedge;
  return tower;
}  
  
//
// TowerType 10
//


  // Equality test

bool TowerType10::operator==(const TowerType10& rhs) const
{
  if (this == &rhs) return true;
  return (_mpaEnergy == rhs._mpaEnergy);
}

// EDM stuff

void TowerType10::readFromBuffer(TBuffer& iobuffer,Version_t version){
  
  if (version<8) {
    iobuffer >> _iMpaEnergy; 
    uint2 dummy;
    iobuffer >> dummy;
  } else {
    iobuffer >> _iMpaEnergy; 
  }
}


void TowerType10::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer << _iMpaEnergy; 
}

//
// Printing  (incomplete)
//

void TowerType10::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType10& t)
{

  os << "MPA energy     : " <<setw(10) <<setprecision(5) << t._mpaEnergy << '\n';

  return os;
}
    
//
// Manipulation and access
//
  
float TowerType10::time(CellKey& cKey) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}

float TowerType10::time(CellKey& cKey, int ihit) const{ 
   
  //There is no time here: return -999 (a la CDF) 
  return -999;
}
	
float TowerType10::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _mpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType10 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType10::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _iMpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType10 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType10::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){

  switch ( d) {
  case CalorDetInfo::MPA:
    _mpaEnergy+= energy;
    _iMpaEnergy = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType10::addPmtEnergy(Detector det , size_t ipmt,float& energy){  

  switch ( det) {
  case MPA:
    _mpaEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType10::nTDCHits(CellKey& cKey) const{                                
  return 0;
}

int TowerType10::setTime(CalorKey& bKey,int time, int ihit){
  
  return 0;
}
int TowerType10::setTime(CellKey& cKey,int time, int ihit){
  
  return 0;
}
int TowerType10::setTDCHits(CalorKey& bKey,uint1 nhits){

  return 0;
}
int TowerType10::setTDCHits(CellKey& cKey,uint1 nhits){

  return 0;
}
TowerType10* TowerType10::clone() const{
  return new TowerType10(*this);
}

void TowerType10::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  if (CalOffLer[12] == -1.) 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * calib->calibMpa();  
  else 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * CalOffLer[12];  
  if ( _mpaEnergy < 0.)  _mpaEnergy = 0.;
}
int TowerType10::iTower() const {
  int phi = this->iPhi();
  int wedge = phi / 3;
  int phiInWedge = phi % 3;
  int tower = 2 + 14 * wedge + phiInWedge;
  return tower;
}    
//
// TowerType 11
//


  // Equality test

bool TowerType11::operator==(const TowerType11& rhs) const
{
  if (this == &rhs) return true;
  return (_mpaEnergy == rhs._mpaEnergy);
}

// EDM stuff

void TowerType11::readFromBuffer(TBuffer& iobuffer,Version_t version){
  
  if (version<8) {
    iobuffer >> _iMpaEnergy; 
    uint2 dummy;
    iobuffer >> dummy;
  } else {
    iobuffer >> _iMpaEnergy; 
  }
}


void TowerType11::writeToBuffer(TBuffer& iobuffer,Version_t version){

  iobuffer << _iMpaEnergy; 
}

//
// Printing  (incomplete)
//

void TowerType11::printTower(std::ostream& os) const{os << *this;}
ostream& operator<<(std::ostream& os, const TowerType11& t)
{
  
  os << "MPA energy     : " <<setw(10) <<setprecision(5) << t._mpaEnergy << '\n';
  
  return os;
}

//
// Manipulation and access
//

float TowerType11::time(CellKey& cKey) const{ 
  
  //There is no time here: return -999 (a la CDF) 
  return -999;
}

float TowerType11::time(CellKey& cKey, int ihit) const{ 
  
  //There is no time here: return -999 (a la CDF) 
  return -999;
}

float TowerType11::pmtEnergy(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _mpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType11 for detector " << d <<std::endl;
    return 0.0;
  }
}
uint2 TowerType11::pmtData(CellKey& cKey, size_t ipmt) const{
  
  Detector d = cKey.detector();
  switch ( d) {
  case MPA:
    return _iMpaEnergy;
  default:
    // should throw an error here...
    std::cout << "Bad Key in TowerType11 for detector " << d <<std::endl;
    return 0;
  }
}
  
int TowerType11::addPmtEnergy(CalorDetInfo::detectorTag d,size_t ipmt,float& energy, uint2 iData){
  
  switch ( d) {
  case CalorDetInfo::MPA:
    _mpaEnergy+= energy;
    _iMpaEnergy = iData;
    return 1;
  default:
    return 0;
  }
}

int TowerType11::addPmtEnergy(Detector det , size_t ipmt,float& energy){  
  
  switch ( det) {
  case MPA:
    _mpaEnergy+= energy;
    return 1;
  default:
    return 0;
  }
}

int TowerType11::nTDCHits(CellKey& cKey) const{                                
  return 0;
}

int TowerType11::setTime(CalorKey& bKey,int time, int ihit){
  
  return 0;
}
int TowerType11::setTime(CellKey& cKey,int time, int ihit){
  
  return 0;
}
int TowerType11::setTDCHits(CalorKey& bKey,uint1 nhits){

  return 0;
}
int TowerType11::setTDCHits(CellKey& cKey,uint1 nhits){

  return 0;
}
TowerType11* TowerType11::clone() const{
  return new TowerType11(*this);
}

void TowerType11::puffTower(const CalDataCalib* calib,float CalOffLer[]){
  if (CalOffLer[12] == -1.) 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * calib->calibMpa();  
  else 
    _mpaEnergy = (RawQieStorableBank::unpack_energy(_iMpaEnergy) -calib->pedMpa()) * CalOffLer[12];  
  if ( _mpaEnergy < 0.)  _mpaEnergy = 0.;
}
int TowerType11::iTower() const {
  int phi = this->iPhi();
  int wedge = phi / 2;
  int phiInWedge = phi % 2;
  int tower = 14 * wedge + phiInWedge;
  return tower;
} 
////////////////////////////////////////////////////////////////////////
//
//  The following methods calculate the energy and phi in a tower
//  given  two photo-tube energies. It can be used for CEM/CHA/WHA/.
//  Then, the TDC pulse height is corrected for CHA and WHA.
//  Finally, the CalTower member data are filled.
//  
//  This is based on code written by Anwar Bhatti. It was taken out 
//  of the routine qftowe.F.  
//
//  In subsequent iterations, this part of the code will be cleaned up 
//  so that Anwar's routines are more smoothly integrated in the new 
//  design
//
///////////////////////////////////////////////////////////////////////

void TowerType0::fillCalTower(int skcut, int paramVersion){

  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float phiCha = 0;
  int   flag   = 0; 
  int detIndex  = 0;
  float pmtEnergy0;
  float pmtEnergy1;
  bool badpmt0;
  bool badpmt1;

  CellKey cKey1(this->iEta(),this->iPhi(), CEM);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  pmtEnergy1 = this->pmtEnergy(cKey1,1);
  detIndex = 1;
 
  badpmt0 = false;
  if (this->isBad(CEMPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(CEMPMT1BAD)) badpmt1 = true;
 
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eCem,phiEm,flag,skcut,badpmt0,badpmt1,paramVersion);
 
  if(flag && skcut>= 0) this->setBad(CEMSPIKE); 
  flag =0;
    
  CellKey cKey2(this->iEta(),this->iPhi(), CHA);       
  pmtEnergy0 = this->pmtEnergy(cKey2,0);
  pmtEnergy1 = this->pmtEnergy(cKey2,1);
  detIndex = 2;
  
  badpmt0 = false;
  if (this->isBad(CHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(CHAPMT1BAD)) badpmt1 = true;
 

  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eCha,phiCha,flag,skcut,badpmt0,badpmt1,paramVersion);

  if(flag && skcut>= 0) this->setBad(CHASPIKE); 
  flag =0;

  // set rest of CalTower member data

  this->setEmEnergy(eCem);
  this->setEmPhi(phiEm);
  this->setHadEnergy(eCha);
  this->setHadPhi(phiCha);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);
}

void TowerType1::fillCalTower(int skcut,int paramVersion){

  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float phiCha = 0;
  float phiWha = 0;
  int flag =0;
  int detIndex = 0;
  float pmtEnergy0;
  float pmtEnergy1;
  bool badpmt0;
  bool badpmt1;
  

  CellKey cKey1(this->iEta(),this->iPhi(), CEM);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  pmtEnergy1 = this->pmtEnergy(cKey1,1);
  detIndex = 1;

  badpmt0 = false;
  if (this->isBad(CEMPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(CEMPMT1BAD)) badpmt1 = true;
  
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eCem,phiEm,flag,skcut,badpmt0,badpmt1,paramVersion);

  if(flag && skcut>= 0) this->setBad(CEMSPIKE);
  flag =0;

  CellKey cKey2(this->iEta(),this->iPhi(), CHA);       

  pmtEnergy0 = this->pmtEnergy(cKey2,0);
  pmtEnergy1 = this->pmtEnergy(cKey2,1);
  detIndex = 2;

  badpmt0 = false;
  if (this->isBad(CHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(CHAPMT1BAD)) badpmt1 = true;
  
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eCha,phiCha,flag,skcut,badpmt0,badpmt1,paramVersion);
  
  if(flag && skcut>= 0) this->setBad(CHASPIKE); 
  flag =0;


  CellKey cKey3(this->iEta(),this->iPhi(), WHA);       
    
  pmtEnergy0 = this->pmtEnergy(cKey3,0);
  pmtEnergy1 = this->pmtEnergy(cKey3,1);
  detIndex = 3;

  badpmt0 = false;
  if (this->isBad(WHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(WHAPMT1BAD)) badpmt1 = true;

  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eWha,phiWha,flag,skcut,badpmt0,badpmt1,paramVersion);

  if(flag && skcut>= 0) this->setBad(WHASPIKE);
  flag =0;

  // set rest of CalTower member data

  this->setEmEnergy(eCem);
  this->setEmPhi(phiEm);
  float eHad=eCha+eWha;
  this->setHadEnergy(eHad);
  if (eHad>0){
    float phi = ((eCha*phiCha+eWha*phiWha)/eHad);          
    this->setHadPhi(phi);
  } 
  else { this->setHadPhi(phiCha); }
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);
}

void TowerType2::fillCalTower(int skcut,int paramVersion){

  float eCem   = 0;
  float ePem   = 0;
  float phiEm = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float phiWha = 0;
  int flag =0;
  int detIndex = 0;
  float pmtEnergy0;
  float pmtEnergy1;
  bool badpmt0;
  bool badpmt1;

  CellKey cKey1(this->iEta(),this->iPhi(), CEM);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  pmtEnergy1 = this->pmtEnergy(cKey1,1);
  detIndex = 1;
  
  badpmt0 = false;
  if (this->isBad(CEMPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(CEMPMT1BAD)) badpmt1 = true;

  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eCem,phiEm,flag,skcut,badpmt0,badpmt1,paramVersion);

  if(flag && skcut>= 0) this->setBad(CEMSPIKE);
  flag =0;

  CellKey cKey2(this->iEta(),this->iPhi(), WHA);       

  pmtEnergy0 = this->pmtEnergy(cKey2,0);
  pmtEnergy1 = this->pmtEnergy(cKey2,1);
  detIndex = 3;
  
  badpmt0 = false;
  if (this->isBad(WHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(WHAPMT1BAD)) badpmt1 = true;
  
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eWha,phiWha,flag,skcut,badpmt0,badpmt1,paramVersion);

  if(flag && skcut>= 0) this->setBad(WHASPIKE); 
  flag =0;

  // set rest of CalTower member data

  this->setEmEnergy(eCem);
  this->setEmPhi(phiEm);
  this->setHadEnergy(eWha);
  this->setHadPhi(phiWha);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}

void TowerType3::fillCalTower(int skcut,int paramVersion){
  float pprWeight=pprWeight_0;
  switch (paramVersion) {
  case 2:
    pprWeight=pprWeight_1;
  }
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float phiWha = 0;
  int   flag   = 0;
  int   detIndex = 0;
  float pmtEnergy0;
  float pmtEnergy1;
  bool badpmt0;
  bool badpmt1;

  CellKey cKey1(this->iEta(),this->iPhi(), PEM);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  pmtEnergy1 = this->pmtEnergy(cKey1,1);
  detIndex = 1;
  ePem = pmtEnergy0+pmtEnergy1;
  if (ePem>0) { phiEm= (0.25*pmtEnergy0+0.75*pmtEnergy1)/ePem; }
  else       { phiEm= 0.5; }

  CellKey cKey2(this->iEta(),this->iPhi(), WHA);       

  pmtEnergy0 = this->pmtEnergy(cKey2,0);
  pmtEnergy1 = this->pmtEnergy(cKey2,1);
  detIndex = 3;
 
  badpmt0 = false;
  if (this->isBad(WHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(WHAPMT1BAD)) badpmt1 = true;
  
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eWha,phiWha,flag,skcut,badpmt0,badpmt1,paramVersion);
  
  if(flag && skcut>= 0) this->setBad(WHASPIKE);
  flag =0;

  CellKey cKey3(this->iEta(),this->iPhi(), PPR);       
  pmtEnergy0 = this->pmtEnergy(cKey3,0);
  pmtEnergy1 = this->pmtEnergy(cKey3,1);
  detIndex = 1;
  ePpr = pmtEnergy0+pmtEnergy1;

  // set rest of CalTower member data

  float eEm=ePem+pprWeight*ePpr;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);
  this->setHadEnergy(eWha);
  this->setHadPhi(phiWha);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);
}

void TowerType4::fillCalTower(int skcut,int paramVersion){
  float pprWeight=pprWeight_0;
  switch (paramVersion) {
  case 2:
    pprWeight=pprWeight_1;
  }
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float phiPha = 0;
  float phiWha = 0;
  int flag =0;
  int detIndex = 0;
  float pmtEnergy0;
  float pmtEnergy1;
  bool badpmt0;
  bool badpmt1;

  CellKey cKey1(this->iEta(),this->iPhi(), PEM);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  pmtEnergy1 = this->pmtEnergy(cKey1,1);
  detIndex = 1;
  ePem = pmtEnergy0+pmtEnergy1;
  if (ePem>0) { phiEm= (0.25*pmtEnergy0+0.75*pmtEnergy1)/ePem; } 
  else       { phiEm= 0.5; }

  CellKey cKey2(this->iEta(),this->iPhi(), WHA);       

  pmtEnergy0 = this->pmtEnergy(cKey2,0);
  pmtEnergy1 = this->pmtEnergy(cKey2,1);
  detIndex = 3;
  
  badpmt0 = false;
  if (this->isBad(WHAPMT0BAD)) badpmt0 = true;
  badpmt1 = false;
  if (this->isBad(WHAPMT1BAD)) badpmt1 = true;
  
  getScintillatorEnergy(detIndex,pmtEnergy0,pmtEnergy1,
			eWha,phiWha,flag,skcut,badpmt0,badpmt1,paramVersion);
  
  if(flag && skcut>= 0) this->setBad(WHASPIKE);
  flag =0;

  CellKey cKey3(this->iEta(),this->iPhi(), PHA);       

  pmtEnergy0 = this->pmtEnergy(cKey3,0);
  pmtEnergy1 = this->pmtEnergy(cKey3,1);
  ePha=pmtEnergy0+pmtEnergy1;
  if (ePha>0) { phiPha= (0.25*pmtEnergy0+0.75*pmtEnergy1)/ePha; }
  else        { phiPha= 0.5; }

  CellKey cKey4(this->iEta(),this->iPhi(), PPR);       
  pmtEnergy0 = this->pmtEnergy(cKey4,0);
  pmtEnergy1 = this->pmtEnergy(cKey4,1);
  detIndex = 1;
  ePpr = pmtEnergy0+pmtEnergy1;
	
  // set rest of CalTower member data

  float eEm=ePem+pprWeight*ePpr;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);

  float eHad=eWha+ePha;
  this->setHadEnergy(eHad);
  if (eHad>0) { 
    float phi = ((eWha*phiWha+ePha*phiPha)/eHad);         
    this->setHadPhi(phi); 
  }
  else { this->setHadPhi(phiWha); }
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}

void TowerType8::fillCalTower(int skcut,int paramVersion){
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float eMpa   = 0;
  float phiHad = 0;
  float pmtEnergy0=0;

  CellKey cKey1(this->iEta(),this->iPhi(), MPA);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  eMpa = pmtEnergy0;
  phiEm= 0.5;
    
  // set rest of CalTower member data

  float eEm=eMpa;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);

  float eHad=0;
  this->setHadEnergy(eHad);
  this->setHadPhi(phiHad);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}
void TowerType9::fillCalTower(int skcut,int paramVersion){
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float eMpa   = 0;
  float phiHad = 0;
  float pmtEnergy0=0;

  CellKey cKey1(this->iEta(),this->iPhi(), MPA);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  eMpa = pmtEnergy0;
  phiEm= 0.5;
    
  // set rest of CalTower member data

  float eEm=eMpa;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);
  float eHad=0;
  this->setHadEnergy(eHad);
  this->setHadPhi(phiHad);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}
void TowerType10::fillCalTower(int skcut,int paramVersion){
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float eMpa   = 0;
  float phiHad = 0;
  float pmtEnergy0=0;

  CellKey cKey1(this->iEta(),this->iPhi(), MPA);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  eMpa = pmtEnergy0;
  phiEm= 0.5;
    
  // set rest of CalTower member data

  float eEm=eMpa;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);
  float eHad=0;
  this->setHadEnergy(eHad);
  this->setHadPhi(phiHad);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}
void TowerType11::fillCalTower(int skcut,int paramVersion){
  float eCem   = 0;
  float ePem   = 0;
  float phiEm  = 0;
  float eCha   = 0;
  float eWha   = 0;
  float ePha   = 0;
  float ePpr   = 0;
  float eMpa   = 0;
  float phiHad = 0;
  float pmtEnergy0=0;

  CellKey cKey1(this->iEta(),this->iPhi(), MPA);       
  pmtEnergy0 = this->pmtEnergy(cKey1,0);
  eMpa = pmtEnergy0;
  phiEm= 0.5;
    
  // set rest of CalTower member data

  float eEm=eMpa;
  this->setEmEnergy(eEm);
  this->setEmPhi(phiEm);
  float eHad=0;
  this->setHadEnergy(eHad);
  this->setHadPhi(phiHad);
  this->setCemEnergy(eCem);
  this->setPemEnergy(ePem);
  this->setChaEnergy(eCha);
  this->setWhaEnergy(eWha);
  this->setPhaEnergy(ePha);
  this->setPprEnergy(ePpr);

}
//
// The other TowerTypes do not have 2 PMTs and they use 
// the CalTower base class fillTower method.
//
// } // namespace calor
#endif  // USE_CDFEDM2