StEEmcTower.h

#ifndef __StEEmcTower_h__
#define __StEEmcTower_h__

#include "StEEmcElement.h"
#include "TVector3.h"
#include <vector>

class StEEmcTower;
typedef std::vector<StEEmcTower*> StEEmcTowerPtrVec_t;

class StEEmcTower : public StEEmcElement {
public:
  StEEmcTower();
  virtual ~StEEmcTower(){ /* nada */ };

  StEEmcTower( const StEEmcTower &other );

  TVector3 direction( Float_t zvertex=0.0 ) const; 

  void index(Int_t i);
  void layer(Int_t l) {mLayer=l;}
  void et(Float_t e) {mET=e;}

  Int_t index()const {return mIndex;}
  Int_t layer()const {return mLayer;}
  Int_t sector()const {return mSector;}
  Int_t subsector()const {return mSubsector;}
  Int_t etabin()const {return mEtabin;}
  Int_t phibin()const {return mPhibin;}
  Float_t et()const {return mET;}

  void neighbor(StEEmcTower *n) {mNeighbors.push_back(n);}
  Int_t numberOfNeighbors()const { return (Int_t)mNeighbors.size();}
  StEEmcTower &neighbor(Int_t i) { return *mNeighbors[i]; }
  const StEEmcTower &neighbor(Int_t i)const { return *mNeighbors[i]; }

  Bool_t isNeighbor( const StEEmcTower &t ) const;
  
  void print() const;
  void printLine() const;

  Bool_t operator<( const StEEmcTower &other ) const { return (this->energy() < other.energy()); }

  virtual void Clear(Option_t *opts=""){ mET=0.; StEEmcElement::Clear(opts); } 

protected:

  Int_t mIndex;      
  Int_t mLayer;      
  Int_t mSector;     
  Int_t mSubsector;  
  Int_t mEtabin;     
  Int_t mPhibin;     
  Float_t mET;

  StEEmcTowerPtrVec_t mNeighbors; 

  ClassDef(StEEmcTower,1);  

};
typedef std::vector<StEEmcTower>  StEEmcTowerVec_t;

#endif