LauAbsDPDynamics.hh

Go to the documentation of this file.

// Copyright University of Warwick 2005 - 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_ABS_DP_DYNAMICS
#define LAU_ABS_DP_DYNAMICS

#include <vector>
#include <map>

#include "TString.h"

#include "LauCacheData.hh"
#include "LauParameter.hh"

class LauAbsResonance;
class LauComplex;
class LauDaughters;
class LauEffModel;
class LauFitDataTree;
class LauKinematics;
class LauResonanceMaker;


class LauAbsDPDynamics {

        public:
                typedef std::map<Int_t,LauEffModel*> LauTagCatScfFractionModelMap;


                LauAbsDPDynamics(LauDaughters* daughters, LauEffModel* effModel, LauEffModel* scfFractionModel = 0);


                LauAbsDPDynamics(LauDaughters* daughters, LauEffModel* effModel, const LauTagCatScfFractionModelMap& scfFractionModel);

                virtual ~LauAbsDPDynamics();

                enum ToyMCStatus {
                        GenOK,        
                        MaxIterError, 
                        ASqMaxError   
                };


                virtual void addResonance(const TString& resName, Int_t resPairAmpInt, const TString& resType,
                                Double_t newMass, Double_t newWidth, Int_t newSpin) = 0;


                virtual void initialise(const std::vector<LauComplex>& coeffs) = 0;


                virtual Double_t retrieveEfficiency();


                virtual Double_t retrieveScfFraction(Int_t tagCat);


                virtual void fillDataTree(const LauFitDataTree& inputFitTree) = 0;


                virtual void setDataEventNo(UInt_t iEvt);


                virtual Bool_t hasResonance(const TString& resName) const;


                virtual TString getConjResName(const TString& resName) const;


                virtual void updateCoeffs(const std::vector<LauComplex>& coeffs);


                virtual Bool_t generate() = 0;


                virtual ToyMCStatus checkToyMC(Bool_t printErrorMessages = kTRUE, Bool_t printInfoMessages = kFALSE) = 0;


                virtual void calcLikelihoodInfo(Double_t m13Sq, Double_t m23Sq) = 0;


                virtual void calcLikelihoodInfo(UInt_t iEvt) = 0;


                virtual Double_t getEvtLikelihood() const = 0;


                virtual Double_t getEvtm13Sq() const = 0;


                virtual Double_t getEvtm23Sq() const = 0;


                virtual Double_t getEvtmPrime() const = 0;


                virtual Double_t getEvtthPrime() const = 0;


                virtual Double_t getEvtEff() const = 0;


                virtual Double_t getEvtJacobian() const = 0;            //Not sure whether this should be here or only in LauIsobarDynamics.


                virtual Double_t getEvtScfFraction() const = 0;


                virtual const LauComplex& getEvtDPAmp() const = 0;


                virtual LauComplex getDynamicAmp(Int_t resID) const = 0;


                virtual void calcExtraInfo(Bool_t init = kFALSE) = 0;

                // Use this method to select events in the DP for embedded-reweighting.

                virtual Bool_t gotReweightedEvent() = 0;


                virtual Double_t getEventWeight() = 0;


                inline LauParameter getMeanEff() const {return meanDPEff_;}


                inline LauParameter getDPRate() const {return DPRate_;}


                inline const LauParArray& getFitFractions() const {return fitFrac_;}


                inline UInt_t getnDefinedResonances() const {return nResDefMax_;}


                inline UInt_t getnAmp() const {return nAmp_;}


                inline Double_t getDPNorm() const {return DPNorm_;}


                inline UInt_t nData() const {return data_.size();}


                inline const std::vector<LauCacheData*>& getCacheData() const {return data_;}


                inline LauDaughters* getDaughters() {return daughters_;}


                inline LauResonanceMaker* getResonanceMaker() {return resonanceMaker_;}


                inline LauKinematics* getKinematics() {return kinematics_;}


                inline LauEffModel* getEffModel() {return effModel_;}


                inline LauEffModel* getScfFractionModel() {return scfFractionModel_[0];}


                inline std::map <Int_t,LauEffModel*> getScfFractionModels() {return scfFractionModel_;}


                inline Bool_t usingScfModel() { return ! scfFractionModel_.empty(); }


                std::vector<LauParameter> getExtraParameters() {return extraParameters_;}

        protected:

                virtual Double_t calcSigDPNorm() = 0;


                virtual LauAbsResonance* findResonance(const TString& name) = 0;


                virtual const LauAbsResonance* findResonance(const TString& name) const = 0;

                LauDaughters* daughters_;

                LauResonanceMaker* resonanceMaker_;

                LauKinematics* kinematics_;

                LauEffModel* effModel_;


                LauTagCatScfFractionModelMap scfFractionModel_;

                UInt_t nAmp_;

                UInt_t nResDefMax_;

                std::vector<LauComplex> Amp_;

                Double_t DPNorm_;

                LauParArray fitFrac_;

                LauParameter DPRate_;

                LauParameter meanDPEff_;

                std::vector<LauCacheData*> data_;

                LauCacheData* currentEvent_;

                std::vector<LauParameter> extraParameters_;

        private:
                ClassDef(LauAbsDPDynamics,0)
};

#endif