LauDPDepBifurGaussPdf.hh

Go to the documentation of this file.

// Copyright University of Warwick 2007 - 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_BIFURGAUSS_PDF
#define LAU_DPDEP_BIFURGAUSS_PDF

#include <vector>

#include "TString.h"

#include "LauAbsPdf.hh"


class LauDPDepBifurGaussPdf : public LauAbsPdf {

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

                enum ScaleMethod {
                        poly,           
                        polyNegPower    
                };


                LauDPDepBifurGaussPdf(const TString& theVarName, const std::vector<LauParameter*>& params,
                                Double_t minAbscissa, Double_t maxAbscissa,
                                const LauDaughters* daughters,
                                const std::vector<Double_t>& meanCoeffs,
                                const std::vector<Double_t>& sigmaLCoeffs,
                                const std::vector<Double_t>& sigmaRCoeffs,
                                DPAxis dpAxis);

                virtual ~LauDPDepBifurGaussPdf();

                LauDPDepBifurGaussPdf(const LauDPDepBifurGaussPdf& other);


                virtual Bool_t isDPDependent() const {return kTRUE;}


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

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

                virtual void calcNorm();
                

                virtual void calcPDFHeight( const LauKinematics* kinematics );
                

                ScaleMethod scaleMethod() const {return scaleMethod_;}
                

                void scaleMethod(ScaleMethod method) {scaleMethod_ = method;}
                                                
        protected:

                void scalePars_poly(Double_t perEventDist);
                

                void scalePars_polyNegPower(Double_t perEventDist);
                
        private:
                const LauKinematics* kinematics_;

                LauParameter* mean_;
                LauParameter* sigmaL_;
                LauParameter* sigmaR_;

                Double_t meanVal_;
                Double_t sigmaLVal_;
                Double_t sigmaRVal_;

                const std::vector<Double_t> meanCoeffs_;
                const std::vector<Double_t> sigmaLCoeffs_;
                const std::vector<Double_t> sigmaRCoeffs_;

                DPAxis dpAxis_;

                ScaleMethod scaleMethod_;

                ClassDef(LauDPDepBifurGaussPdf,0) // Define the MVA PDF
};

#endif