LauRhoOmegaMix.hh

Go to the documentation of this file.

// Copyright University of Warwick 2004 - 2016.
// 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_RHO_OMEGA_MIX
#define LAU_RHO_OMEGA_MIX

#include "TString.h"

#include "LauAbsResonance.hh"
#include "LauComplex.hh"

#include <vector>

class LauResonanceInfo;
class LauDaughters;
class LauParameter;

class LauRhoOmegaMix : public LauAbsResonance {

public:

    LauRhoOmegaMix(LauResonanceInfo* resInfo, const LauAbsResonance::LauResonanceModel resType,
                   const Int_t resPairAmpInt, const LauDaughters* daughters);

    virtual ~LauRhoOmegaMix();

    virtual void initialise();

    void initialiseRho();

    void initialiseOmega();


    virtual LauComplex amplitude(const LauKinematics* kinematics);


    virtual LauAbsResonance::LauResonanceModel getResonanceModel() const {return model_;}


    virtual void setResonanceParameter(const TString& name, const Double_t value);


    virtual void floatResonanceParameter(const TString& name);
    

    virtual LauParameter* getResonanceParameter(const TString& name);
    

    virtual const std::vector<LauParameter*>& getFloatingParameters();
    
protected:


    void setmOmegaValue(const Double_t mOmega);    
    

    Double_t getmOmegaValue() const { return (mOmega_!=0) ? mOmega_->unblindValue() : 0.0; }
    

    Bool_t fixmOmegaValue() const { return (mOmega_!=0) ? mOmega_->fixed() : 0.0; }



    void setwOmegaValue(const Double_t wOmega);
        

    Double_t getwOmegaValue() const { return (wOmega_!=0) ? wOmega_->unblindValue() : 0.0; }
    

    Bool_t fixwOmegaValue() const { return (wOmega_!=0) ? wOmega_->fixed() : 0.0; }

    

    void setmagBValue(const Double_t magB);


    Double_t getmagBValue() const { return (magB_!=0) ? magB_->unblindValue() : 0.0; }
    

    Bool_t fixmagBValue() const { return (magB_!=0) ? magB_->fixed() : 0.0; }



    void setphiBValue(const Double_t phiB);


    Double_t getphiBValue() const { return (phiB_!=0) ? phiB_->unblindValue() : 0.0; }
    

    Bool_t fixphiBValue() const { return (phiB_!=0) ? phiB_->fixed() : 0.0; }



    void setdeltaValue(const Double_t delta);


    Double_t getdeltaValue() const { return (delta_!=0) ? delta_->unblindValue() : 0.0; }
    

    Bool_t fixdeltaValue() const { return (delta_!=0) ? delta_->fixed() : 0.0; }
   
    virtual LauComplex resAmp(Double_t mass, Double_t spinTerm);
    
    void checkDaughterTypes() const;

private:

    LauRhoOmegaMix(const LauRhoOmegaMix& rhs);
    
    LauRhoOmegaMix& operator=(const LauRhoOmegaMix& rhs);
    
    LauAbsResonance::LauResonanceModel model_;

    Double_t rhoMass_;

    Double_t rhoResRadius_;

    Double_t rhoParRadius_;

    LauParameter* mOmega_;

    Double_t mOmega0_;

    Double_t mOmegaCur_;

    LauParameter* wOmega_;

    Double_t wOmega0_;

    Double_t wOmegaCur_;

    LauParameter* magB_;

    LauParameter* phiB_;

    LauParameter* delta_;

    Bool_t useDenom_;

    Bool_t doneFirstInit_;

    LauAbsResonance* rhoRes_;

    LauAbsResonance* omegaRes_;

    ClassDef(LauRhoOmegaMix,0) // Rho-omega mixing model

};

#endif