LauDPDepGaussPdf.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_DPDEP_GAUSS_PDF
#define LAU_DPDEP_GAUSS_PDF

#include <vector>

#include "TString.h"

#include "LauAbsPdf.hh"

class LauDaughters;
class LauKinematics;
class LauParameter;


class LauDPDepGaussPdf : public LauAbsPdf {

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


                LauDPDepGaussPdf(const TString& theVarName, const std::vector<LauAbsRValue*>& params,
                                Double_t minAbscissa, Double_t maxAbscissa,
                                const LauDaughters* daughters,
                                const std::vector<Double_t>& meanCoeffs,
                                const std::vector<Double_t>& sigmaCoeffs,
                                DPAxis dpAxis);

                virtual ~LauDPDepGaussPdf();

                LauDPDepGaussPdf(const LauDPDepGaussPdf& other);


                virtual Bool_t isDPDependent() const {return kTRUE;}


                virtual void calcLikelihoodInfo(const LauAbscissas& abscissas);
                
                using LauAbsPdf::calcLikelihoodInfo;

                virtual void calcNorm();
                

                virtual void calcPDFHeight( const LauKinematics* kinematics );

        protected:

                void scalePars( Double_t dpPos );

        private:
                const LauKinematics* kinematics_;

                LauAbsRValue* mean_;
                LauAbsRValue* sigma_;

                Double_t meanVal_;
                Double_t sigmaVal_;

                const std::vector<Double_t> meanCoeffs_;
                const std::vector<Double_t> sigmaCoeffs_;

                DPAxis dpAxis_;

                ClassDef(LauDPDepGaussPdf,0) // Define the Gaussian PDF
};

#endif