LauDPDepBifurGaussPdf.hh

Go to the documentation of this file.

 
/*
Copyright 2007 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_DPDEP_BIFURGAUSS_PDF
#define LAU_DPDEP_BIFURGAUSS_PDF
 
#include "LauAbsPdf.hh"
 
#include "TString.h"
 
#include <vector>
 
class LauDaughters;
 
class LauDPDepBifurGaussPdf : public LauAbsPdf {
 
  public:
    enum DPAxis {
        M12,        
        M13,        
        M23,        
        CentreDist, 
        MMIN,       
        MMAX        
    };
 
    enum ScaleMethod {
        poly,        
        polyNegPower 
    };
 
 
    LauDPDepBifurGaussPdf( 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>& sigmaLCoeffs,
                           const std::vector<Double_t>& sigmaRCoeffs,
                           DPAxis dpAxis );
 
    virtual ~LauDPDepBifurGaussPdf();
 
 
    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 );
 
 
    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:
    LauDPDepBifurGaussPdf( const LauDPDepBifurGaussPdf& other );
 
    LauDPDepBifurGaussPdf& operator=( const LauDPDepBifurGaussPdf& other );
 
    const LauKinematics* kinematics_;
 
    LauAbsRValue* mean_;
    LauAbsRValue* sigmaL_;
    LauAbsRValue* 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