StFmsEventClusterer.h

Go to the documentation of this file.

// $Id$
//
// $Log$
#ifndef STROOT_STFMSPOINTMAKER_STFMSEVENTCLUSTERER_H_
#define STROOT_STFMSPOINTMAKER_STFMSEVENTCLUSTERER_H_

#include <vector>

#include "TObject.h"

#include "StFmsPointMaker/StFmsClusterFinder.h"

namespace FMSCluster {  // $NMSPC
class StFmsClusterFitter;
class StFmsFittedPhoton;
class StFmsGeometry;
class StFmsEventClusterer: public TObject {
 public:
  StFmsEventClusterer(const StFmsGeometry* geometry, Int_t detectorId);
  ~StFmsEventClusterer();
  int detector() const { return mDetectorId; }
  Bool_t cluster(std::vector<FMSCluster::StFmsTower>* towers);
#ifndef __CINT__  // Hide ClusterList from CINT as it uses C++11

  ClusterList& clusters() { return mClusters; }
  const ClusterList& clusters() const { return mClusters; }
#endif  // __CINT__

 private:
#ifndef __CINT__  // Hide ClusterList from CINT as it uses C++11

  typedef ClusterList::iterator ClusterIter;
  typedef ClusterList::const_iterator ClusterConstIter;
#endif  // __CINT__

  StFmsEventClusterer(const StFmsEventClusterer&);
  StFmsEventClusterer& operator=(const StFmsEventClusterer&);
  Int_t fitEvent();
  Double_t photonEnergyInCluster(Double_t towerWidth,
                                 const StFmsTowerCluster* cluster,
                                 const StFmsFittedPhoton* photon) const;
  Double_t photonEnergyInTower(Double_t towerWidth, const StFmsTower* tower,
                               const StFmsFittedPhoton* photon) const;
  Float_t fitOnePhoton(StFmsTowerCluster* cluster);
#ifndef __CINT__  // Hide Cluster(Const)Iter from CINT as it uses C++11
  /*
   Perform a global fit of all photons in an event.

   Update the (x, y) positions end energies of the photons in each cluster based
   on a global fit including all photons.
   Only makes sense when there is more than one photon in the event.
   Arguments:
    - nPh, number of photons in the event
    - nCl, number of clusters containing those photons
    - first, iterator to the first cluster

   Returns the &chi;<sup>2</sup> of the fit.
   */
  Float_t globalFit(const Int_t, const Int_t, ClusterIter first);
  /*
   Special 2-photon fit for a single cluster.

   Cluster moments must have been calculated first

   Returns the &chi;<sup>2</sup> of the fit.
   */
  Float_t fit2PhotonClust(ClusterIter cluster);
  /*
   Run tests on the lower-energy photon in a 2-photon cluster.

   Return true if the photon passes tests, in which case it is a real photon.
   Return false if it fails, in which case it is a bogus photon due to some
   problem in reconstruction - the cluster is actually a 1-photon cluster.

   Arguments:
    - clusterIndex: index of cluster to test
    - nRealClusters: total number of clusters in the event
   */
  bool validate2ndPhoton(ClusterConstIter cluster) const;
  ClusterList mClusters;  
#endif  // __CINT__
  StFmsClusterFinder mClusterFinder;   
  const StFmsGeometry* mGeometry;   
  Int_t mDetectorId;   
  std::vector<FMSCluster::StFmsTower>* mTowers;   
  StFmsClusterFitter* mFitter;   
  std::vector<Float_t> mTowerWidthXY;   
  ClassDef(StFmsEventClusterer, 0)
};
}  // namespace FMSCluster
#endif  // STROOT_STFMSPOINTMAKER_STFMSEVENTCLUSTERER_H_