GloriaClusterStatistics.h

// -*- mode: c++; mode: auto-fill; mode: flyspell-prog; -*-


#ifndef GloriaClusterStatistics_h
#define GloriaClusterStatistics_h 1

/*
 *   This source code is part of the LUIZA software package for GLORIA.
 *   You are free to use this source files for your own development as
 *   long as it stays in a public research context. You are not
 *   allowed to use it for commercial purpose. You must put this
 *   header with author names in all development based on this file.
 *
 */

namespace gloria {


  class GloriaClusterStatistics {

  public:


    GloriaClusterStatistics();

    ~GloriaClusterStatistics();

    void Clear();


    void AddPixel(int ix, int iy, double signal);


    inline int GetNpixel() { return _nPixel; };

    inline double GetSignal() { return _signal; };

    inline double GetXpos(){ return _sumX/_signal + _seedX; };

    inline double GetYpos(){ return _sumY/_signal + _seedY; };

    double GetXspread();

    double GetYspread();

    double GetXYcorr();

    double GetFWHM();

    double GetElongation();

    double GetOrientation();

    double GetLength();

    double GetWidth();


private:

    void CalculateAverages();

    // Number of pixels in the cluster
    int _nPixel;

    // Position of first (seed?) pixel
    // Used as offset for sum calculations
    double _seedX, _seedY;

    // Total signal in the cluster
    double _signal;

    // Weighted sums
    double _sumX, _sumY, _sumXX, _sumYY, _sumXY;

    //Flag for doing basic statistics
    bool _averDone;
  
    //Results of calculations
    double _cogX, _cogY, _rmsX2, _rmsY2, _rmsX, _rmsY, _corrXY;
  } ;


}

#endif