StFmsClusterFinder.h

Go to the documentation of this file.

// $Id: StFmsClusterFinder.h,v 1.5 2015/10/30 21:33:55 akio Exp $
//
// $Log: StFmsClusterFinder.h,v $
// Revision 1.5  2015/10/30 21:33:55  akio
// fix parameter initialization
// adding new cluster categorization method
//
// Revision 1.4  2015/10/29 21:14:55  akio
// increase max number of clusters
// a bug fixes in valley tower association
// removing some debug comments
//
// Revision 1.3  2015/10/21 15:58:04  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.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:53  jeromel
// First version of FmsUtil from Yuxi Pan - reviewd 2015/02
//
#ifndef STROOT_STFMSPOINTMAKER_STFMSCLUSTERFINDER_H_
#define STROOT_STFMSPOINTMAKER_STFMSCLUSTERFINDER_H_

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

#include "Rtypes.h"  // Provides ROOT ClassDef macro
class TObjArray;

namespace FMSCluster {  // $NMSPC
class StFmsTowerCluster;
class StFmsTower;
// Typedef a list of cluster unique_ptrs for convenience.
// Hide it from CINT as it can't handle parsing the header :(
#ifndef __CINT__
typedef std::list<std::unique_ptr<StFmsTowerCluster>> ClusterList;
#endif  // __CINT__

class StFmsClusterFinder {
 public:
  // Typedef StFmsTower pointer list for convenience
  typedef std::list<FMSCluster::StFmsTower*> TowerList;
  StFmsClusterFinder(double energyCutoff = 0.5);
  // Use default copy constructor and assignment operator.
  virtual ~StFmsClusterFinder();
  void calculateClusterMoments(StFmsTowerCluster* cluster) const;
  int categorise(StFmsTowerCluster* cluster);
  int categorise2(StFmsTowerCluster* cluster);

  double momentEnergyCutoff() const { return mEnergyCutoff; }
#ifndef __CINT__  // Hide ClusterList from CINT

   int findClusters(TowerList* towers, ClusterList* clusters, int detetorId);
#endif  // __CINT__

 private:
  static const unsigned kMaxNClusters = 25;  
#ifndef __CINT__  // Hide ClusterList from CINT

  unsigned locateClusterSeeds(TowerList* towers, TowerList* neighbors,
                              ClusterList* clusters) const;
  unsigned associateTowersWithClusters(TowerList* neighbors,
                                       ClusterList* clusters,
                                       TObjArray* valleys) const;
  unsigned associateValleyTowersWithClusters(TowerList* neighbors,
                                             ClusterList* clusters,
                                             TObjArray* valleys) const;
  unsigned associateResidualTowersWithClusters(TowerList* neighbors,
                                               ClusterList* clusters) const;
  void associateSubThresholdTowersWithClusters(TowerList* towers,
                                               ClusterList* clusters) const;
  
#endif  // __CINT__
  Double_t mEnergyCutoff;  
  Int_t mNClusts;  
  Int_t mDetectorId; //current working detectorId

  ClassDef(StFmsClusterFinder, 0)
};  // class StFmsClusterFinder
}  // namespace FMSCluster
#endif  // STROOT_STFMSPOINTMAKER_STFMSCLUSTERFINDER_H_