StFmsClusterFitter.h

Go to the documentation of this file.

// $Id$
//
// $Log$
#ifndef STROOT_STFMSPOINTMAKER_STFMSCLUSTERFITTER_H_
#define STROOT_STFMSPOINTMAKER_STFMSCLUSTERFITTER_H_

#include <list>

#include "TMinuit.h"
#include "TObject.h"

#include "StFmsPointMaker/StFmsFittedPhoton.h"
#include "StFmsPointMaker/StFmsTowerCluster.h"

class TF2;

namespace FMSCluster {  // $NMSPC
typedef std::list<StFmsFittedPhoton> PhotonList;
class StFmsGeometry;
class StFmsTower;
class StFmsClusterFitter : public TObject {
 public:
  StFmsClusterFitter(const StFmsGeometry* geometry, Int_t detectorId);
  StFmsClusterFitter() { }
  ~StFmsClusterFitter();
  TF2* showerShapeFunction();
  void setTowers(StFmsTowerCluster::Towers* towers) { mTowers = towers; }
  Double_t fit(const Double_t* par, const Double_t* step,
               const Double_t* low, const Double_t* up, PhotonList* photons);
  Int_t fit2PhotonCluster(const Double_t* para, const Double_t* step,
                          const Double_t* low, const Double_t* up,
                          PhotonList* photons);
  static Double_t energyDepositionInTower(Double_t* x, Double_t* par);
  static int maxNFittedPhotons() { return kMaxNPhotons; }

 private:
  static const Int_t kMaxNPhotons = 7;  

  StFmsClusterFitter(const StFmsClusterFitter&);
  StFmsClusterFitter& operator=(const StFmsClusterFitter&);
  static Double_t energyDepositionDistribution(Double_t* x, Double_t* par);
  static void minimizationFunctionNPhoton(Int_t& npar, Double_t* grad,
                                          Double_t& fval, Double_t* par,
                                          Int_t iflag);
  static void minimizationFunction2Photon(Int_t& nparam, Double_t* grad,
                                          Double_t& fval, Double_t* param,
                                          Int_t iflag);
  Double_t mSteps[3 * kMaxNPhotons + 1];  
  Double_t mTowerWidth;  
  TMinuit mMinuit;  
  static StFmsTowerCluster::Towers* mTowers;  
  static Float_t mTowerWidthXY[2];  
  ClassDef(StFmsClusterFitter, 0)
};  // class StFmsClusterFitter
}  // namespace FMSCluster
#endif  // STROOT_STFMSPOINTMAKER_STFMSCLUSTERFITTER_H_