Lau2DCubicSpline.hh

Go to the documentation of this file.

/*
Copyright 2013 University of Warwick

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
*/

/*
Laura++ package authors:
John Back
Paul Harrison
Thomas Latham
*/

#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