Lau2DHistPdf.hh

Go to the documentation of this file.

/*
Copyright 2008 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_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();

                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:
                Lau2DHistPdf(const Lau2DHistPdf& other);

                Lau2DHistPdf& operator=(const Lau2DHistPdf& other);

                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