source: trunk/MagicSoft/Mars/mimage/MNewImagePar.h@ 6893

#ifndef MARS_MNewImagePar
#define MARS_MNewImagePar

#ifndef MARS_MParContainer
#include "MParContainer.h"
#endif

class MHillas;
class MGeomCam;
class MSignalCam;

class MNewImagePar : public MParContainer
{
private:
    // FIXME: MOVE ALL PARAMETRS WHICH DOES'T DEPEND ON THE IMAGE
    //        CLEANING TO A NEW CLASS!
    Float_t fLeakage1;             // (photons in most outer ring of pixels) over fSize
    Float_t fLeakage2;             // (photons in the 2 outer rings of pixels) over fSize
    Float_t fInnerLeakage1;        // (photons in most outer rings of inner pixels) over fInnerSize
    Float_t fInnerLeakage2;        // (photons in the 2 outer rings of inner pixels) over fInnerSize
    Float_t fInnerSize;            //

    Float_t fConc;                 // [ratio] concentration ratio: sum of the two highest pixels / fSize
    Float_t fConc1;                // [ratio] concentration ratio: sum of the highest pixel / fSize

    Float_t fUsedArea;             // Area of pixels which survived the image cleaning
    Float_t fCoreArea;             // Area of core pixels
    Short_t fNumUsedPixels;        // Number of pixels which survived the image cleaning
    Short_t fNumCorePixels;        // number of core pixels

public:
    MNewImagePar(const char *name=NULL, const char *title=NULL);

    void Reset();

    Float_t GetLeakage1() const { return fLeakage1; }
    Float_t GetLeakage2() const { return fLeakage2; }

    Float_t GetInnerLeakage1() const { return fInnerLeakage1; }
    Float_t GetInnerLeakage2() const { return fInnerLeakage2; }
    Float_t GetInnerSize()     const { return fInnerSize; }

    Float_t GetConc() const  { return fConc;  }
    Float_t GetConc1() const { return fConc1; }

    Short_t GetNumUsedPixels() const { return fNumUsedPixels; }
    Short_t GetNumCorePixels() const { return fNumCorePixels; }

    Float_t GetUsedArea() const { return fUsedArea; }
    Float_t GetCoreArea() const { return fCoreArea; }

    void Print(Option_t *opt=NULL) const;
    void Print(const MGeomCam &geom) const;

    void Calc(const MGeomCam &geom, const MSignalCam &evt,
              const MHillas &hillas, Int_t island=-1);

    ClassDef(MNewImagePar, 4) // Container to hold new image parameters
};

#endif