StFmsTowerCluster.h

Go to the documentation of this file.

// $Id$
//
// $Log$
#ifndef STROOT_STFMSPOINTMAKER_STFMSTOWERCLUSTER_H_
#define STROOT_STFMSPOINTMAKER_STFMSTOWERCLUSTER_H_

#include "Rtypes.h"  // For ClassDef macro

#include <list>
#include <memory>  // For unique_ptr

#include "StFmsPointMaker/StFmsFittedPhoton.h"

class StFmsCluster;

namespace FMSCluster {  // $NMSPC
class StFmsTower;
class StFmsTowerCluster {
 public:
  typedef std::list<StFmsTower*> Towers;  

  explicit StFmsTowerCluster(StFmsCluster* cluster);
  // Use default copy constructor and assignment operator
  ~StFmsTowerCluster();
  /* Clear photon and tower lists and reset other values to defaults */
  void Clear(const char* optionNotUsed = "");
  void calculateClusterMoments(Float_t energyCutoff);
  void findClusterAxis(Float_t Ecoff) {
    mEnergyCutoff = Ecoff;
    findClusterAxis();
  }
  Int_t index() const { return mIndex; }
  void setIndex(Int_t index) { mIndex = index; }
  float sigmaX() const { return mSigmaX; }
  float sigmaY() const { return mSigmaY; }
  float sigmaXY() const { return mSigmaXY; }
  Float_t thetaAxis() const { return mThetaAxis; }
  Float_t chiSquare() const { return mChiSquare; }
  void setChiSquare(Float_t chi2) { mChiSquare = chi2; }
  float energyCutoff() const { return mEnergyCutoff; }
  Towers& towers() { return mTowers; }
  const Towers& towers() const { return mTowers; }
  StFmsFittedPhoton* photons() { return mPhotons; }
  const StFmsFittedPhoton* photons() const { return mPhotons; }
  StFmsCluster* cluster() { return mCluster.get(); }
  const StFmsCluster* cluster() const { return mCluster.get(); }
  StFmsCluster* release() { return mCluster.release(); }

 protected:
  static const int kMaxPhotonsPerCluster = 2;  

  void findClusterAxis();
  Double_t getSigma(Double_t theta) const;
  Int_t mIndex;  
  Float_t mSigmaX;  
  Float_t mSigmaY;  
  Float_t mSigmaXY;  
  Float_t mThetaAxis;  
  Float_t mChiSquare;  
  Float_t mEnergyCutoff;  
  Towers mTowers;  
#ifndef __CINT__  // CINT won't parse unique_ptr so hide it
  std::unique_ptr<StFmsCluster> mCluster;  
#endif  // __CINT__
  StFmsFittedPhoton mPhotons[kMaxPhotonsPerCluster];  

 private:
  StFmsTowerCluster(const StFmsTowerCluster&);
  StFmsTowerCluster& operator=(const StFmsTowerCluster&);
  ClassDef(StFmsTowerCluster, 0)
};  // class StFmsTowerCluster
}  // namespace FMSCluster
#endif  // STROOT_STFMSPOINTMAKER_STFMSTOWERCLUSTER_H_