Lau2DHistPdf.hh

Go to the documentation of this file.

// Copyright University of Warwick 2008 - 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_2DHIST_PDF
#define LAU_2DHIST_PDF

#include "LauAbsPdf.hh"

class TH2;
class TH1;
class Lau1DHistPdf;

class Lau2DHistPdf : public LauAbsPdf {

        public:

                Lau2DHistPdf(const std::vector<TString>& theVarNames, const TH2* hist,
                                const LauFitData& minVals, const LauFitData& maxVals,
                                Bool_t useInterpolation = kTRUE, Bool_t fluctuateBins = kFALSE);
                
                virtual ~Lau2DHistPdf();

                Lau2DHistPdf(const Lau2DHistPdf& other);

                Double_t getMinX() const {return minX_;}

                Double_t getMaxX() const {return maxX_;}

                Double_t getRangeX() const {return rangeX_;}

                Double_t getMinY() const {return minY_;}

                Double_t getMaxY() const {return maxY_;}

                Double_t getRangeY() const {return rangeY_;}


                virtual void cacheInfo(const LauFitDataTree& inputData);


                virtual void calcLikelihoodInfo(const LauAbscissas& abscissas);
                

                virtual void calcLikelihoodInfo(UInt_t iEvt);



                virtual Double_t getLikelihood( const TString& theVarName ) const;
                using LauAbsPdf::getLikelihood;

                virtual void calcNorm();


                virtual void calcPDFHeight( const LauKinematics* kinematics );


                virtual LauFitData generate(const LauKinematics* kinematics);

        protected:
                void doBinFluctuation();

                void checkNormalisation();
                

                Double_t getBinHistValue(Int_t i, Int_t j) const;
                

                Double_t interpolateXY(Double_t x, Double_t y) const;


                Double_t interpolateXYNorm(Double_t x, Double_t y) const;

        private:
                TH2* hist_;
                TH1* xProj_;
                TH1* yProj_;

                Lau1DHistPdf* xVarPdf_;
                Lau1DHistPdf* yVarPdf_;

                TString xName_;
                TString yName_;

                Int_t nBinsX_;
                Int_t nBinsY_;
                Double_t minX_;
                Double_t maxX_;
                Double_t minY_;
                Double_t maxY_;
                Double_t rangeX_;
                Double_t rangeY_;
                
                Double_t binXWidth_;
                Double_t binYWidth_; 

                Double_t invBinXWidth_;
                Double_t invBinYWidth_;

                Bool_t useInterpolation_;
                Bool_t fluctuateBins_;

                ClassDef(Lau2DHistPdf,0) // 2D histogram pdf class
};

#endif