LauDPDepSumPdf.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_DPDEPSUM_PDF
#define LAU_DPDEPSUM_PDF

#include "Rtypes.h"

#include "LauAbsPdf.hh"

class LauDaughters;
class LauEffModel;
class LauParameter;
class TH2;


class LauDPDepSumPdf : public LauAbsPdf {

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


                LauDPDepSumPdf(LauAbsPdf* pdf1, LauAbsPdf* pdf2,
                                const LauDaughters* daughters,
                                const TH2* dpHisto, Bool_t upperHalf = kFALSE,
                                Bool_t useSpline = kFALSE);


                LauDPDepSumPdf(LauAbsPdf* pdf1, LauAbsPdf* pdf2,
                                LauParameter* frac,
                                const LauDaughters* daughters,
                                const std::vector<Double_t>& fracCoeffs,
                                DPAxis dpAxis);

                virtual ~LauDPDepSumPdf();

                //LauDPDepSumPdf(const LauDPDepSumPdf& 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();

                virtual void calcNorm();
                

                virtual void calcPDFHeight( const LauKinematics* kinematics );

        protected:

                void scaleFrac( Double_t dpPos );

        private:
                LauDPDepSumPdf(const LauDPDepSumPdf& other);

                LauDaughters* daughters_;

                LauAbsPdf* pdf1_;
                LauAbsPdf* pdf2_;

                LauParameter* frac_;
                Double_t fracVal_;
                LauEffModel* dpDependence_;

                const std::vector<Double_t> fracCoeffs_;

                DPAxis dpAxis_;

                std::vector<Double_t> fractions_;

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

#endif