LauDPDepMapPdf.hh

Go to the documentation of this file.

// Copyright University of Warwick 2009 - 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_DPDEPMAP_PDF
#define LAU_DPDEPMAP_PDF

#include <vector>

#include "Rtypes.h"

#include "LauAbsPdf.hh"

class LauDaughters;
class Lau2DHistDP;
class LauParameter;
class TH2;
class TH1;


class LauDPDepMapPdf : public LauAbsPdf {

        public:
                enum DPAxis {
                        M12,            
                        M13,            
                        M23,            
                        CentreDist,     
                        MMIN,           
                        MMAX            
                };


                LauDPDepMapPdf( const std::vector<LauAbsPdf*> &pdfs,
                                const LauDaughters* daughters,
                                const TH2* dpHisto, Bool_t upperHalf = kFALSE);


                LauDPDepMapPdf( const std::vector<LauAbsPdf*> &pdfs,
                                const LauDaughters* daughters,
                                const TH1* dpAxisHisto,
                                DPAxis dpAxis);

                virtual ~LauDPDepMapPdf();

                LauDPDepMapPdf(const LauDPDepMapPdf& other);


                virtual Bool_t isDPDependent() const {return kTRUE;}


                virtual void cacheInfo(const LauFitDataTree& inputData);


                virtual void calcLikelihoodInfo(const LauAbscissas& abscissas);
                

                virtual void calcLikelihoodInfo(UInt_t iEvt);

                virtual void checkPositiveness(); // Nothing to check here.

                virtual void calcNorm();
                

                virtual void calcPDFHeight( const LauKinematics* kinematics );

        protected:

                UInt_t determineDPRegion( Double_t m13Sq, Double_t m23Sq ) const;

        private:
                LauDaughters* daughters_;

                std::vector<LauAbsPdf*> pdfs_;

                Lau2DHistDP* dpDependence_;
                TH1* dpAxisDependence_;

                DPAxis dpAxis_;

                std::vector<UInt_t> indices_;

                ClassDef(LauDPDepMapPdf,0) // Define the sum PDF
};

#endif