#ifndef RICH_MIRRORS_CBMRICHRONCHIANA_H_
#define RICH_MIRRORS_CBMRICHRONCHIANA_H_

#include "TObject.h"
#include <string>
#include <vector>

#include <cmath>

using namespace std;

class TH2;

class CbmRichRonchiAna : public TObject 
{
public:
    CbmRichRonchiAna();
    
    virtual ~CbmRichRonchiAna();

    void Run();

    void SetTiffFileNameV(const string& fileName) {fTiffFileNameV = fileName;}
    void SetTiffFileNameH(const string& fileName) {fTiffFileNameH = fileName;}

private:

    string fTiffFileNameV;
    string fTiffFileNameH;
    
    // constant values
    double fPi;
    int fRadiusMirror;                 // in microns
    int fEdgeLengthCCD;                // in microns
    int fEdgeLengthPixel;              // in microns
    int fPitchGrid;                    // in microns; distance of lines
    int fImageWidth;                   // in pixels    
    
    // values to be measured first
    int fDistRulingCCD;                // in microns 
    int fDistMirrorRuling;             // in microns
    int fEffectiveEdgeLengthGridPlate; // in microns
    int fEffectiveEdgeLengthMirror;    // in microns; edge length of part of mirror that is projected onto CCD
    int fEffectiveEdgeLengthCCD;       // in microns; edge length of mirror image on CCD (NumberOfLines*LineDistance)
    int fLineDistance;                 // in pixels    
    int fCenterX;                      // in pixels;  approximate center of image
    int fCenterY;                      // in pixels
            
    vector<vector<int> > ReadTiffFile(const string& fileName);

    void DoMeanIntensityY(vector<vector<int> >& data);

    void DoPeakFinderY(vector<vector<int> >& data);

    void DoRotation(vector<vector<int> >& data);

    void DoMeanY(vector<vector<int> >& data);

    void FillH2WithVector(TH2* hist, const vector<vector<int> >& data);
    
    void DoEjectSingle(vector<vector<int> >& data);
    
    void DoFillGaps(vector<vector<int> >& data);
    
    void DoShadowCenter(vector<vector<int> >& data);
    
    vector<vector<int> > DoSuperpose(const vector<vector<int> >& dataH, const vector<vector<int> >& dataV);
    
    vector<vector<int> > DoIntersection(vector<vector<int> >& dataH, const vector<vector<int> >& dataV);
    
    void DoIntersectionFinder(vector<vector<int> >& dataSup);
    
    vector<vector<int> > DoNumberIntersections(vector<vector<int> >& intersectionXY);
    
    void DoScanLine(vector<vector<int> >& intersectionXY, vector<vector<int> >& matrix);     
        
    void DoSearchNextLine(vector<vector<int> >& matrix);
    
    void DoSearchAbove(vector<vector<int>>& intersectionXY, vector<vector<int> >& matrix);
    
    void DoSearchBelow(vector<vector<int>>& intersectionXY, vector<vector<int> >& matrix);
    
    vector<vector<double>> DoLocalNormal(vector<vector<int> >& matrix);
    
    vector<vector<double>> DoHeight(vector<vector<double> >& normal);
    
    void DoSphere(vector<vector<double> >& sphere);
    
    void DoScanLineHeight(vector<vector<double> >& sphere);
    
    void DoCalculateRemaining(vector<vector<double> >& sphere);
    
    void DoIntegrate(vector<vector<double> >& sphere);
    
    void DoAverageSurroundings(vector<vector<double> >& sphere);
    
    /**
    * \brief Copy constructor.
    */
    CbmRichRonchiAna(const CbmRichRonchiAna&);

    /**
    * \brief Assignment operator.
    */
    CbmRichRonchiAna& operator=(const CbmRichRonchiAna&);

    ClassDef(CbmRichRonchiAna,1)
};

#endif /* RICH_MIRRORS_CBMRICHRONCHIANA_H_ */