StFmsTowerCluster.h

Go to the documentation of this file.

// $Id: StFmsTowerCluster.h,v 1.5 2016/01/26 14:42:48 akio Exp $
//
// $Log: StFmsTowerCluster.h,v $
// Revision 1.5  2016/01/26 14:42:48  akio
// better chi2 handling
//
// Revision 1.4  2015/10/21 15:58:05  akio
// Code speed up (~x2) by optimizing minimization fuctions and showershape function
// Add option to merge small cells to large, so that it finds cluster at border
// Add option to perform 1photon fit when 2photon fit faield
// Add option to turn on/off global refit
// Moment analysis done without ECUTOFF when no tower in cluster exceed ECUTOFF=0.5GeV
//
// Revision 1.3  2015/10/01 19:55:48  akio
// *** empty log message ***
//
// Revision 1.2  2015/09/02 15:01:32  akio
// Removing StFmsGeometry class, and now it uses StFmsDbMaker to get appropriate parameters.
//
// Revision 1.1  2015/03/10 14:38:54  jeromel
// First version of FmsUtil from Yuxi Pan - reviewd 2015/02
//
#ifndef STROOT_STFMSPOINTMAKER_STFMSTOWERCLUSTER_H_
#define STROOT_STFMSPOINTMAKER_STFMSTOWERCLUSTER_H_

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

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

#include "StFmsUtil/StFmsFittedPhoton.h"

class StFmsCluster;

namespace FMSCluster {  // $NMSPC
class StFmsTower;
class StFmsTowerCluster {
 public:
  typedef std::list<StFmsTower*> Towers;  
  typedef std::vector<StFmsFittedPhoton> Photons;
  explicit StFmsTowerCluster(StFmsCluster* cluster, Int_t detectorId);
  // Use default copy constructor and assignment operator
  virtual ~StFmsTowerCluster();
  /* Clear photon and tower lists and reset other values to defaults */
  void Clear(const char* optionNotUsed = "");
  void calculateClusterMoments(Double_t energyCutoff);
  void findClusterAxis(Double_t Ecoff, Double_t xwidth, Double_t ywidth) {
    mEnergyCutoff = Ecoff;
    findClusterAxis(xwidth, ywidth);
  }
  Int_t index() const { return mIndex; }
  void setIndex(Int_t index) { mIndex = index; }
  Int_t detectorId() const { return mDetectorId; }
  double etot() const { return mEtot; }
  double sigmaX() const { return mSigmaX; }
  double sigmaY() const { return mSigmaY; }
  double sigmaXY() const { return mSigmaXY; }
  Double_t thetaAxis() const { return mThetaAxis; }
  Double_t chiSquare()  const { return mChiSquare; }
  Double_t chiSquare1() const { return mChiSquare1; }
  Double_t chiSquare2() const { return mChiSquare2; }
  void setChiSquare(Double_t chi2)  { mChiSquare = chi2; }
  void setChiSquare1(Double_t chi2) { mChiSquare1 = chi2; }
  void setChiSquare2(Double_t chi2) { mChiSquare2 = chi2; }
  double energyCutoff() const { return mEnergyCutoff; }
  Towers& towers() { return mTowers; }
  const Towers& towers() const { return mTowers; }
  Photons& photons() { return mPhotons; }
  const Photons& photons() const { return mPhotons; }
  StFmsCluster* cluster() { return mCluster.get(); }
  const StFmsCluster* cluster() const { return mCluster.get(); }
  StFmsCluster* release(); 

 protected:
  void findClusterAxis(Double_t xwidth, Double_t ywidth);
  Double_t getSigma(Double_t theta, Double_t xwidth, Double_t ywidth);
  Int_t mIndex;  
  Int_t mDetectorId;
  Float_t mEtot;
  Double_t mSigmaX;  
  Double_t mSigmaY;  
  Double_t mSigmaXY;  
  Double_t mThetaAxis;  
  Double_t mChiSquare;  
  Double_t mChiSquare1;  
  Double_t mChiSquare2;  
  Double_t mEnergyCutoff;  
  Towers mTowers;  
#ifndef __CINT__  // CINT won't parse unique_ptr so hide it
  std::unique_ptr<StFmsCluster> mCluster;  
#endif  // __CINT__
  Photons mPhotons;  

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