LauLASSRes.hh
No OneTemporary
Actions

Size

5 KB

Referenced Files

None

Subscribers

None

LauLASSRes.hh
View Options


	/*
	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
	*/

	/*! \file LauLASSRes.hh
	\brief File containing declaration of LauLASSRes class.
	*/

	/*! \class LauLASSRes
	\brief Class for defining the LASS resonance model.

	Class for defining the LASS resonance model, which includes a resonant
	part, plus an effective range term. This class is the coherent sum of the two.
	*/

	#ifndef LAU_LASS_RES
	#define LAU_LASS_RES

	#include "TString.h"

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


	class LauLASSRes : public LauAbsResonance {

	public:
	//! Constructor
	/*!
	\param [in] resInfo the object containing information on the resonance name, mass, width, spin, charge, etc.
	\param [in] resPairAmpInt the number of the daughter not produced by the resonance
	\param [in] daughters the daughter particles
	*/
	LauLASSRes(LauResonanceInfo* resInfo, const Int_t resPairAmpInt, const LauDaughters* daughters);

	//! Destructor
	virtual ~LauLASSRes();

	//! Initialise the model
	virtual void initialise();

	//! Get the resonance model type
	/*!
	\return the resonance model type
	*/
	virtual LauAbsResonance::LauResonanceModel getResonanceModel() const {return LauAbsResonance::LASS;}

	//! Get the cut off parameter
	/*!
	\return the cut off parameter
	*/
	Double_t getCutOff() const {return cutOff_;}

	//! Set the cut off parameter
	/*!
	\param [in] cutOff the cut off parameter
	*/
	void setCutOff(const Double_t cutOff) {cutOff_ = cutOff;}

	//! Set value of a resonance parameter
	/*!
	\param [in] name the name of the parameter to be changed
	\param [in] value the new parameter value
	*/
	virtual void setResonanceParameter(const TString& name, const Double_t value);

	//! Allow the various parameters to float in the fit
	/*!
	\param [in] name the name of the parameter to be floated
	*/
	virtual void floatResonanceParameter(const TString& name);

	//! Access the given resonance parameter
	/*!
	\param [in] name the name of the parameter
	\return the corresponding parameter
	*/
	virtual LauParameter* getResonanceParameter(const TString& name);

	//! Retrieve the resonance parameters, e.g. so that they can be loaded into a fit
	/*!
	\return floating parameters of the resonance
	*/
	virtual const std::vector<LauParameter*>& getFloatingParameters();

	protected:
	//! Set the effective range parameter value
	/*!
	\param [in] r the new effective range parameter value
	*/
	void setEffectiveRange(const Double_t r);

	//! Set the scattering length parameter value
	/*!
	\param [in] a the new scattering length parameter value
	*/
	void setScatteringLength(const Double_t a);

	//! Get the effective range parameter
	/*!
	\return the effective range parameter
	*/
	Double_t getEffectiveRange() const {return (r_!=0) ? r_->unblindValue() : 0.0;}

	//! Get the scattering length range parameter
	/*!
	\return the scattering length parameter
	*/
	Double_t getScatteringLength() const {return (a_!=0) ? a_->unblindValue() : 0.0;}

	//! See if the effective range parameter is fixed or floating
	/*!
	\return kTRUE if the effective range parameter is fixed, kFALSE otherwise
	*/
	Bool_t fixEffectiveRange() const {return (r_!=0) ? r_->fixed() : kTRUE;}

	//! See if the scattering length parameter is fixed or floating
	/*!
	\return kTRUE if the scattering length parameter is fixed, kFALSE otherwise
	*/
	Bool_t fixScatteringLength() const {return (a_!=0) ? a_->fixed() : kTRUE;}

	//! Utility function to calculate the q0 value
	void calcQ0();

	//! Complex resonant amplitude
	/*!
	\param [in] mass appropriate invariant mass for the resonance
	\param [in] spinTerm spin term
	*/
	virtual LauComplex resAmp(Double_t mass, Double_t spinTerm);

	private:
	//! Copy constructor (not implemented)
	LauLASSRes(const LauLASSRes& rhs);

	//! Copy assignment operator (not implemented)
	LauLASSRes& operator=(const LauLASSRes& rhs);

	//! Decay momentum of either daughter in the resonance rest frame
	Double_t q0_;
	//! Sum of the daughter masses
	Double_t mDaugSum_;
	//! Square of mDaugSum
	Double_t mDaugSumSq_;
	//! Difference between the daughter masses
	Double_t mDaugDiff_;
	//! Square of mDaugDiff
	Double_t mDaugDiffSq_;
	//! The resonance mass
	Double_t resMass_;
	//! Square of the resonance mass
	Double_t resMassSq_;

	//! LASS effective range parameter
	LauParameter* r_;
	//! LASS scattering length parameter
	LauParameter* a_;

	//! LASS cut off
	Double_t cutOff_;

	ClassDef(LauLASSRes,0)
	};

	#endif

File Metadata

Mime Type: text/x-c
Expires: Tue, Sep 30, 4:40 AM (5 h, 43 m)
Storage Engine: blob
Storage Format: Raw Data
Storage Handle: 6564742
Default Alt Text: LauLASSRes.hh (5 KB)

LauLASSRes.hhNo OneTemporaryActions

LauLASSRes.hhView Options

File Metadata

Event Timeline

LauLASSRes.hh
No OneTemporary
Actions

LauLASSRes.hh
View Options