Lau2DCubicSpline.hh

Go to the documentation of this file.

// Copyright University of Warwick 2006 - 2013.
// Distributed under the Boost Software License, Version 1.0.
// (See accompanying file LICENSE_1_0.txt or copy at http://www.boost.org/LICENSE_1_0.txt)

// Authors:
// Thomas Latham
// John Back
// Paul Harrison

#ifndef LAU_2DCUBICSPLINE
#define LAU_2DCUBICSPLINE

class TH2;

class Lau2DCubicSpline
{
    private:
        enum { NCoeff = 16, CoeffRecLen = 17 };

    public:

        Lau2DCubicSpline(const TH2& h);

        virtual ~Lau2DCubicSpline();


        virtual Double_t evaluate(Double_t x, Double_t y) const;


        virtual Double_t analyticalIntegral(Double_t x1, Double_t x2, Double_t y1, Double_t y2) const;
        
        virtual Double_t analyticalIntegral() const;


    private:
        Lau2DCubicSpline( const Lau2DCubicSpline& rhs );

        Lau2DCubicSpline& operator=(const Lau2DCubicSpline& rhs);

        Int_t nBinsX;

        Int_t nBinsY;

        Double_t binSizeX;

        Double_t binSizeY;

        Double_t xmin;

        Double_t xmax;

        Double_t ymin;

        Double_t ymax;

        std::vector<Double_t> coeffs;


        Double_t histcont(const TH2& h, Int_t xbin, Int_t ybin) const;


        Double_t dhistdx(const TH2& h, Int_t xbin, Int_t ybin) const;


        Double_t dhistdy(const TH2& h, Int_t xbin, Int_t ybin) const;


        Double_t d2histdxdy(const TH2& h, Int_t xbin, Int_t ybin) const;


        inline const Double_t& coeff(Int_t binx, Int_t biny, Int_t theCoeff) const
        { return coeffs[theCoeff + CoeffRecLen * (binx + nBinsX * biny)]; }


        inline Double_t& coeff(Int_t binx, Int_t biny, Int_t theCoeff)
        { return coeffs[theCoeff + CoeffRecLen * (binx + nBinsX * biny)]; }


        Double_t evalX(Double_t x1, Double_t x2, Double_t y) const;
        

        Double_t evalY(Double_t x, Double_t y1, Double_t y2) const;
        

        Double_t evalXY(Double_t x1, Double_t x2, Double_t y1, Double_t y2) const;

        ClassDef(Lau2DCubicSpline, 0); // Class for defining a 2D cubic spline based on an input histogram
};

#endif