StEEmcCluster.h

#ifndef __StEEmcCluster_h__
#define __StEEmcCluster_h__

#include <TObject.h>
#include <TVector3.h>
#include <TMath.h>
#include "StEEmcPool/StEEmcA2EMaker/StEEmcTower.h"
#include "StEEmcBaseCluster.h"

class StEmcCluster;

class StEEmcCluster : public StEEmcBaseCluster {
public:
  StEEmcCluster();
  StEEmcCluster ( const StEEmcCluster &other );
  virtual ~StEEmcCluster();

  void add(const StEEmcTower &t, Float_t weight=1.0 );
  
  Float_t energy() const {return mEnergy;}

  Float_t seedEnergy() const {return (mTowers.size())?mTowers[0].energy():0.;}

  TVector3 momentum() const {return mMomentum;}

  TVector3 position() const {return mPosition;}

  Int_t numberOfTowers() const { return (Int_t)mTowers.size(); }

  StEEmcTower &tower(Int_t t) { return mTowers[t]; }
  const StEEmcTower &tower(Int_t t) const { return mTowers[t]; }

  Float_t sigmaE() const;

  Float_t weight(Int_t t) const { return mWeights[t]; }
  StEEmcTowerVec_t &towers() { return mTowers; }
  const StEEmcTowerVec_t &towers() const { return mTowers; }

  void momentum(const TVector3 &p) {mMomentum=p;}

  StEmcCluster *stemc();

  void stemc(StEmcCluster *c){ mEmcCluster = c; }

  Bool_t operator==( const StEEmcCluster &other ) const { return this->tower(0).index() == other.tower(0).index(); }

  void print() const;

  void printLine(Bool_t Endl=false) const;

  Bool_t isNeighbor(const StEEmcTower &tower) const;

  Float_t fracEtabin() const;
  Float_t fracPhibin() const;
  
  Float_t sigmaEtabin() const { return TMath::Sqrt(mSumEta2W/mEnergy-mSumEtaW*mSumEtaW/mEnergy/mEnergy); }
  Float_t sigmaPhibin() const { return TMath::Sqrt(mSumPhi2W/mEnergy-mSumPhiW*mSumPhiW/mEnergy/mEnergy); }

  Bool_t hasTower(const StEEmcTower &tower) const;

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

  Int_t numberOfEtabins() const;
  Int_t numberOfPhibins() const;

protected:
  StEEmcTowerVec_t mTowers;
  std::vector<Float_t> mWeights;

  TVector3 mMomentum;
  TVector3 mPosition;

  StEmcCluster *mEmcCluster;

  Float_t mfEtabin; 
  Float_t mfPhibin; 
  Float_t mSumEta2W;
  Float_t mSumEtaW;
  Float_t mSumPhi2W;
  Float_t mSumPhiW;
  
  ClassDef(StEEmcCluster,1);
};

typedef std::vector<StEEmcCluster> StEEmcClusterVec_t;

#endif