LauEffModel.hh

Go to the documentation of this file.

/*
Copyright 2004 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 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_;}

        private:
                LauEffModel( const LauEffModel& rhs );

                LauEffModel& operator=( const LauEffModel& rhs );


                Double_t getEffHistValue( const LauKinematics* kinematics ) const;

                const LauDaughters* daughters_;

                const LauVetoes* vetoes_;

                std::vector<Lau2DAbsDP*> effHisto_;

                Bool_t fluctuateEffHisto_;

                mutable Bool_t lowBinWarningIssued_;

                mutable Bool_t highBinWarningIssued_;

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

};

#endif