LauDPDepSumPdf.hh

Go to the documentation of this file.

/*
Copyright 2009 University of Warwick

Licensed under the Apache License, Version 2.0 (the "License");
you may not use this file except in compliance with the License.
You may obtain a copy of the License at

    http://www.apache.org/licenses/LICENSE-2.0

Unless required by applicable law or agreed to in writing, software
distributed under the License is distributed on an "AS IS" BASIS,
WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
See the License for the specific language governing permissions and
limitations under the License.
*/

/*
Laura++ package authors:
John Back
Paul Harrison
Thomas Latham
*/

#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();


                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 calcNorm();
                

                virtual void calcPDFHeight( const LauKinematics* kinematics );

        protected:

                void scaleFrac( Double_t dpPos );

        private:
                LauDPDepSumPdf(const LauDPDepSumPdf& other);

                LauDPDepSumPdf& operator=(const LauDPDepSumPdf& other);

                LauDaughters* daughters_;

                LauAbsPdf* pdf1_;
                LauAbsPdf* pdf2_;

                LauAbsRValue* 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