LauDPDepCruijffPdf.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_CRUIJFF_PDF
#define LAU_DPDEP_CRUIJFF_PDF

#include <vector>

#include "TString.h"

#include "LauAbsPdf.hh"

class LauDaughters;
class LauKinematics;
class LauParameter;


class LauDPDepCruijffPdf : public LauAbsPdf {

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


                LauDPDepCruijffPdf(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,
                                const std::vector<Double_t>& alphaLCoeffs,
                                const std::vector<Double_t>& alphaRCoeffs,
                                DPAxis dpAxis);

                virtual ~LauDPDepCruijffPdf();

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

        protected:

                void scalePars( Double_t dpPos );


                Double_t currentPDFValue(Double_t abscissa) const;


                virtual Double_t integrGaussLegendre();


                virtual Double_t integTrapezoid();

        private:
                const LauKinematics* kinematics_;

                LauParameter* mean_;
                LauParameter* sigmaL_;
                LauParameter* sigmaR_;
                LauParameter* alphaL_;
                LauParameter* alphaR_;

                Double_t meanVal_;
                Double_t sigmaLVal_;
                Double_t sigmaRVal_;
                Double_t alphaLVal_;
                Double_t alphaRVal_;

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

                DPAxis dpAxis_;

                ClassDef(LauDPDepCruijffPdf,0) // Define the Cruijff PDF
};

#endif