#include "CbmEcalMaximum.h"

#include "CbmEcalCell.h"

#include <list>

using namespace std;

CbmEcalMaximum::CbmEcalMaximum(CbmEcalCell* cell)
  : fCell(cell)
{
  Double_t me=0;
  Double_t e;
  Int_t i;
  list<CbmEcalCell*> cells;
  list<CbmEcalCell*>::const_iterator p;

  fCX=cell->GetCenterX();
  fCY=cell->GetCenterY();
  fI=1;
  for(i=1;i<5;i++)
  {
    cell->GetNeighborsList(i, cells);
    e=0;
    for(p=cells.begin();p!=cells.end();++p)
      e+=(*p)->GetTotalEnergy();
    if (e>me)
    {
      fI=i;
      me=e;
    }
  }
  me+=cell->GetTotalEnergy();
  fX=cell->GetTotalEnergy()*fCX;
  fY=cell->GetTotalEnergy()*fCY;
  cell->GetNeighborsList(fI, cells);
  for(p=cells.begin();p!=cells.end();++p)
  {
    fX+=(*p)->GetTotalEnergy()*(*p)->GetCenterX();
    fY+=(*p)->GetTotalEnergy()*(*p)->GetCenterY();
  }
  fX/=me; fY/=me;
  fMark=0;
}
ClassImp(CbmEcalMaximum)