LauEffModel.hh

Go to the documentation of this file.

// Copyright University of Warwick 2004 - 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 LAUEFFMODEL
#define LAUEFFMODEL

#include "LauAbsEffModel.hh"

class TH2;

class LauDaughters;
class LauKinematics;
class LauVetoes;
class Lau2DAbsDP;


class LauEffModel : public LauAbsEffModel {

        public:

                LauEffModel(const LauDaughters* daughters, const LauVetoes* vetoes);

                virtual ~LauEffModel();


                void setEffHisto(const TH2* effHisto,
                                Bool_t useInterpolation = kTRUE, Bool_t fluctuateBins = kFALSE,
                                Double_t avEff = -1.0, Double_t absError = -1.0,
                                Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void setEffHisto(const TH2* effHisto, const TH2* errorHi, const TH2* errorLo,
                                Bool_t useInterpolation = kTRUE, Bool_t fluctuateBins = kFALSE,
                                Double_t avEff = -1.0, Double_t absError = -1.0,
                                Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void setEffSpline(const TH2* effHisto,
                                 Bool_t fluctuateBins = kFALSE,
                                 Double_t avEff = -1.0, Double_t absError = -1.0,
                                 Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void setEffSpline(const TH2* effHisto, const TH2* errorHi, const TH2* errorLo,
                                 Bool_t fluctuateBins = kFALSE,
                                 Double_t avEff = -1.0, Double_t absError = -1.0,
                                 Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void addEffHisto(const TH2* effHisto,
                                Bool_t useInterpolation = kTRUE,
                                Double_t avEff = -1.0, Double_t absError = -1.0,
                                Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void addEffHisto(const TH2* effHisto, const TH2* errorHi, const TH2* errorLo,
                                Bool_t useInterpolation = kTRUE,
                                Double_t avEff = -1.0, Double_t absError = -1.0,
                                Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void addEffSpline(const TH2* effHisto,
                                 Double_t avEff = -1.0, Double_t absError = -1.0,
                                 Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                void addEffSpline(const TH2* effHisto, const TH2* errorHi, const TH2* errorLo,
                                 Double_t avEff = -1.0, Double_t absError = -1.0,
                                 Bool_t useUpperHalfOnly = kFALSE, Bool_t squareDP = kFALSE);


                Double_t calcEfficiency( const LauKinematics* kinematics ) const;


                Bool_t passVeto( const LauKinematics* kinematics ) const;

                Bool_t fluctuateEffHisto() const {return fluctuateEffHisto_;}

                /*
                    \return the LauDaughters object associated with the DP
                */
                const LauDaughters* getDaughters() const {return daughters_;}

                /*
                    \return kTRUE if the square DP is being used, kFALSE otherwise
                */
                Bool_t usingSquareDP() const {return squareDP_;}

        private:
                LauEffModel( const LauEffModel& rhs );


                Double_t getEffHistValue(Double_t xVal, Double_t yVal) const;

                const LauDaughters* daughters_;

                const LauVetoes* vetoes_;

                std::vector<Lau2DAbsDP*> effHisto_;

                Bool_t squareDP_;
                Bool_t fluctuateEffHisto_;

                mutable Bool_t lowBinWarningIssued_;

                mutable Bool_t highBinWarningIssued_;

                ClassDef(LauEffModel, 0)   // Implement the signal efficiency across the DP

};

#endif