LauNRAmplitude.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 LAU_NR_AMPLITUDE
#define LAU_NR_AMPLITUDE

#include "TString.h"

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

class LauKinematics;


class LauNRAmplitude : public LauAbsResonance {

        public:

                LauNRAmplitude(LauResonanceInfo* resInfo, const Int_t resPairAmpInt, const LauDaughters* daughters);

                virtual ~LauNRAmplitude();

                virtual void initialise();


                virtual LauComplex amplitude(const LauKinematics* kinematics);


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


                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:
                virtual LauComplex resAmp(Double_t mass, Double_t spinTerm);


                Double_t f(const Double_t s, const Double_t c, const Double_t p) const;


                void setdParameter(const Double_t d);


                void setc1Parameter(const Double_t c1);


                void setc2Parameter(const Double_t c2);


                void setp1Parameter(const Double_t p1);


                void setp2Parameter(const Double_t p2);


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


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


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


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


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


                Bool_t fixdParameter() const {return (d_!=0) ? d_->fixed() : kTRUE;}


                Bool_t fixc1Parameter() const {return (c1_!=0) ? c1_->fixed() : kTRUE;}


                Bool_t fixc2Parameter() const {return (c2_!=0) ? c2_->fixed() : kTRUE;}


                Bool_t fixp1Parameter() const {return (p1_!=0) ? p1_->fixed() : kTRUE;}


                Bool_t fixp2Parameter() const {return (p2_!=0) ? p2_->fixed() : kTRUE;}

        private:
                LauNRAmplitude(const LauNRAmplitude& rhs);

                LauNRAmplitude& operator=(const LauNRAmplitude& rhs);

                LauParameter* d_; 
                LauParameter* c1_; 
                LauParameter* c2_; 
                LauParameter* p1_; 
                LauParameter* p2_;

                ClassDef(LauNRAmplitude,0) // Non-resonant amplitude model
};

#endif