Changeset View
Changeset View
Standalone View
Standalone View
inc/LauKMatrixPropagator.hh
| Show First 20 Lines • Show All 43 Lines • ▼ Show 20 Lines | |||||
| #include "TMatrixD.h" | #include "TMatrixD.h" | ||||
| #include "TString.h" | #include "TString.h" | ||||
| #include <map> | #include <map> | ||||
| #include <vector> | #include <vector> | ||||
| class LauKMatrixPropagator { | class LauKMatrixPropagator { | ||||
| public: | public: | ||||
| //! Constructor | //! Constructor | ||||
| /*! | /*! | ||||
| \param [in] name name of the propagator | \param [in] name name of the propagator | ||||
| \param [in] paramFileName the parameter file name | \param [in] paramFileName the parameter file name | ||||
| \param [in] resPairAmpInt the number of the daughter not produced by the resonance | \param [in] resPairAmpInt the number of the daughter not produced by the resonance | ||||
| \param [in] nChannels the number of channels | \param [in] nChannels the number of channels | ||||
| \param [in] nPoles the number of poles | \param [in] nPoles the number of poles | ||||
| \param [in] rowIndex this specifies which row of the propagator should be used when summing over the amplitude channels | \param [in] rowIndex this specifies which row of the propagator should be used when summing over the amplitude channels | ||||
| */ | */ | ||||
| LauKMatrixPropagator(const TString& name, const TString& paramFileName, | LauKMatrixPropagator( const TString& name, const TString& paramFileName, | ||||
| Int_t resPairAmpInt, Int_t nChannels, Int_t nPoles, | Int_t resPairAmpInt, Int_t nChannels, Int_t nPoles, | ||||
| Int_t rowIndex = 1); | Int_t rowIndex = 1 ); | ||||
| //! Destructor | //! Destructor | ||||
| virtual ~LauKMatrixPropagator(); | virtual ~LauKMatrixPropagator(); | ||||
| //! Calculate the invariant mass squared s | //! Calculate the invariant mass squared s | ||||
| /*! | /*! | ||||
| \param [in] kinematics the kinematics of the current event | \param [in] kinematics the kinematics of the current event | ||||
| */ | */ | ||||
| void updatePropagator(const LauKinematics* kinematics); | void updatePropagator(const LauKinematics* kinematics); | ||||
| //! Calculate the K-matrix propagator for the given s value | //! Calculate the K-matrix propagator for the given s value | ||||
| /*! | /*! | ||||
| \param [in] s the invariant mass squared | \param [in] s the invariant mass squared | ||||
| */ | */ | ||||
| void updatePropagator(Double_t s); | void updatePropagator(Double_t s); | ||||
| //! Read an input file to set parameters | //! Read an input file to set parameters | ||||
| /*! | /*! | ||||
| \param [in] inputFile name of the input file | \param [in] inputFile name of the input file | ||||
| */ | */ | ||||
| void setParameters(const TString& inputFile); | void setParameters(const TString& inputFile); | ||||
| //! Get the scattering K matrix | //! Get the scattering K matrix | ||||
| /*! | /*! | ||||
| \return the real, symmetric scattering K matrix | \return the real, symmetric scattering K matrix | ||||
| */ | */ | ||||
| TMatrixD getKMatrix() const {return ScattKMatrix_;} | TMatrixD getKMatrix() const {return ScattKMatrix_;} | ||||
| //! Get the real part of the propagator full matrix | //! Get the real part of the propagator full matrix | ||||
| /*! | /*! | ||||
| \return the real part of the propagator full matrix | \return the real part of the propagator full matrix | ||||
| */ | */ | ||||
| TMatrixD getRealPropMatrix() const {return realProp_;} | TMatrixD getRealPropMatrix() const {return realProp_;} | ||||
| //! Get the negative imaginary part of the full propagator matrix | //! Get the negative imaginary part of the full propagator matrix | ||||
| /*! | /*! | ||||
| \return the negative imaginary part of the full propagator matrix | \return the negative imaginary part of the full propagator matrix | ||||
| */ | */ | ||||
| TMatrixD getNegImagPropMatrix() const {return negImagProp_;} | TMatrixD getNegImagPropMatrix() const {return negImagProp_;} | ||||
| //! Get the real part of the term of the propagator | //! Get the real part of the term of the propagator | ||||
| /*! | /*! | ||||
| \param [in] channelIndex the channel number | \param [in] channelIndex the channel number | ||||
| \return the real part of the propagator term | \return the real part of the propagator term | ||||
| */ | */ | ||||
| Double_t getRealPropTerm(Int_t channelIndex) const; | Double_t getRealPropTerm(Int_t channelIndex) const; | ||||
| Show All 14 Lines | public: | ||||
| //! Get coupling constants that were loaded from the input file | //! Get coupling constants that were loaded from the input file | ||||
| /*! | /*! | ||||
| \param [in] poleIndex number of the required pole | \param [in] poleIndex number of the required pole | ||||
| \param [in] channelIndex number of the required channel | \param [in] channelIndex number of the required channel | ||||
| \return the value of the coupling constant | \return the value of the coupling constant | ||||
| */ | */ | ||||
| Double_t getCouplingConstant(Int_t poleIndex, Int_t channelIndex) const; | Double_t getCouplingConstant(Int_t poleIndex, Int_t channelIndex) const; | ||||
| //! Get coupling parameters, set according to the input file | |||||
| /*! | |||||
| \param [in] poleIndex number of the required pole | |||||
| \param [in] channelIndex number of the required channel | |||||
| \return the parameter of the coupling constant | |||||
| */ | |||||
| LauParameter& getCouplingParameter(Int_t poleIndex, Int_t channelIndex); | |||||
| //! Get scattering constants that were loaded from the input file | //! Get scattering constants that were loaded from the input file | ||||
| /*! | /*! | ||||
| \param [in] channel1Index number of the first channel index | \param [in] channel1Index number of the first channel index | ||||
| \param [in] channel2Index number of the second channel index | \param [in] channel2Index number of the second channel index | ||||
| \return the value of the scattering constant | \return the value of the scattering constant | ||||
| */ | */ | ||||
| Double_t getScatteringConstant(Int_t channel1Index, Int_t channel2Index) const; | Double_t getScatteringConstant(Int_t channel1Index, Int_t channel2Index) const; | ||||
| //! Get the "slowly-varying part" term of the amplitude | //! Get the "slowly-varying part" term of the amplitude | ||||
| /*! | /*! | ||||
| \return the svp term | \return the svp term | ||||
| */ | */ | ||||
| //! Get scattering parameters, set according to the input file | |||||
| /*! | |||||
| \param [in] channel1Index number of the first channel index | |||||
| \param [in] channel2Index number of the second channel index | |||||
| \return the parameter of the scattering constant | |||||
| */ | |||||
| LauParameter& getScatteringParameter(Int_t channel1Index, Int_t channel2Index); | |||||
| //! Get s0 production parameter | |||||
| /*! | |||||
| \return the s0Prod parameter | |||||
| */ | |||||
| LauParameter& gets0Prod() {return s0Prod_;} | |||||
tlatham: These Doxygen lines should be attached to what is now line 169. | |||||
| Double_t getProdSVPTerm() const {return prodSVP_;} | Double_t getProdSVPTerm() const {return prodSVP_;} | ||||
| //! Get the full complex propagator term for a given channel | //! Get the full complex propagator term for a given channel | ||||
| /*! | /*! | ||||
| \param [in] channelIndex the number of the required channel | \param [in] channelIndex the number of the required channel | ||||
| \return the complex propagator term | \return the complex propagator term | ||||
| */ | */ | ||||
| LauComplex getPropTerm(Int_t channelIndex) const; | LauComplex getPropTerm(Int_t channelIndex) const; | ||||
| //! Get the DP axis identifier | //! Get the DP axis identifier | ||||
| /*! | /*! | ||||
| /return the value to identify the DP axis in question | \return the value to identify the DP axis in question | ||||
| */ | */ | ||||
| Int_t getResPairAmpInt() const {return resPairAmpInt_;} | Int_t getResPairAmpInt() const {return resPairAmpInt_;} | ||||
| //! Get the number of channels | //! Get the number of channels | ||||
| /*! | /*! | ||||
| /return the number of channels | \return the number of channels | ||||
| */ | */ | ||||
| Int_t getNChannels() const {return nChannels_;} | Int_t getNChannels() const {return nChannels_;} | ||||
| //! Get the number of poles | //! Get the number of poles | ||||
| /*! | /*! | ||||
| /return the number of poles | \return the number of poles | ||||
| */ | */ | ||||
| Int_t getNPoles() const {return nPoles_;} | Int_t getNPoles() const {return nPoles_;} | ||||
| //! Get the propagator name | //! Get the propagator name | ||||
| /*! | /*! | ||||
| /return the name of the propagator | \return the name of the propagator | ||||
| */ | */ | ||||
| TString getName() const {return name_;} | TString getName() const {return name_;} | ||||
| //! Get the unitary transition amplitude for the given channel | //! Get the unitary transition amplitude for the given channel | ||||
| /*! | /*! | ||||
| \param [in] s The invariant mass squared | \param [in] s The invariant mass squared | ||||
| \param [in] channel The index number of the channel process | \param [in] channel The index number of the channel process | ||||
| \return the complex amplitude T | \return the complex amplitude T | ||||
| */ | */ | ||||
| LauComplex getTransitionAmp(Double_t s, Int_t channel); | LauComplex getTransitionAmp(Double_t s, Int_t channel); | ||||
| //! Get the complex phase space term for the given channel and invariant mass squared | //! Get the complex phase space term for the given channel and invariant mass squared | ||||
| /*! | /*! | ||||
| \param [in] s The invariant mass squared | \param [in] s The invariant mass squared | ||||
| \param [in] channel The index number of the channel process | \param [in] channel The index number of the channel process | ||||
| \return the complex phase space term rho(channel, channel) | \return the complex phase space term rho(channel, channel) | ||||
| */ | */ | ||||
| LauComplex getPhaseSpaceTerm(Double_t s, Int_t channel); | LauComplex getPhaseSpaceTerm(Double_t s, Int_t channel); | ||||
| //! Get the Adler zero factor, which is set when updatePropagator is called | //! Get the Adler zero factor, which is set when updatePropagator is called | ||||
| /*! | /*! | ||||
| \return the Adler zero factor | \return the Adler zero factor | ||||
| */ | */ | ||||
| Double_t getAdlerZero() const {return adlerZeroFactor_;} | Double_t getAdlerZero() const {return adlerZeroFactor_;} | ||||
| //! Get the THat amplitude for the given s and channel number | //! Get the THat amplitude for the given s and channel number | ||||
| /*! | /*! | ||||
| \param [in] s The invariant mass squared | \param [in] s The invariant mass squared | ||||
| \param [in] channel The index number of the channel process | \param [in] channel The index number of the channel process | ||||
| \return the complex THat amplitude | \return the complex THat amplitude | ||||
| */ | */ | ||||
| LauComplex getTHat(Double_t s, Int_t channel); | LauComplex getTHat(Double_t s, Int_t channel); | ||||
| protected: | protected: | ||||
| //! Calculate the scattering K-matrix for the given value of s | //! Calculate the scattering K-matrix for the given value of s | ||||
| /*! | /*! | ||||
| \param [in] s the invariant mass squared | \param [in] s the invariant mass squared | ||||
| */ | */ | ||||
| void calcScattKMatrix(Double_t s); | void calcScattKMatrix(Double_t s); | ||||
| //! Calculate the real and imaginary part of the phase space density diagonal matrix | //! Calculate the real and imaginary part of the phase space density diagonal matrix | ||||
| /*! | /*! | ||||
| \param [in] s the invariant mass squared | \param [in] s the invariant mass squared | ||||
| */ | */ | ||||
| void calcRhoMatrix(Double_t s); | void calcRhoMatrix(Double_t s); | ||||
| //! Calculate the real and imaginary part of the gamma angular-momentum barrier matrix | |||||
| /*! | |||||
| \param [in] s the invariant mass squared | |||||
| */ | |||||
| void calcGammaMatrix(Double_t s); | |||||
| //! Calculate the DK P-wave gamma angular-momentum barrier | |||||
| /*! | |||||
| \param [in] s the invariant mass squared | |||||
| */ | |||||
| Double_t calcD0KPwaveGamma(Double_t s) const; | |||||
| //! Calulate the term 1/(m_pole^2 - s) for the scattering and production K-matrix formulae | //! Calulate the term 1/(m_pole^2 - s) for the scattering and production K-matrix formulae | ||||
| /*! | /*! | ||||
| \param [in] s the invariant mass squared | \param [in] s the invariant mass squared | ||||
| */ | */ | ||||
| void calcPoleDenomVect(Double_t s); | void calcPoleDenomVect(Double_t s); | ||||
| //! Calculate the D0K+ phase space factor | |||||
| /*! | |||||
| \param [in] s the invariant mass squared | |||||
| \return the complex phase space factor | |||||
| */ | |||||
| LauComplex calcD0KRho(Double_t s) const; | |||||
| //! Calculate the D*0K+ phase space factor | |||||
| /*! | |||||
| \param [in] s the invariant mass squared | |||||
| \return the complex phase space factor | |||||
| */ | |||||
| LauComplex calcDstar0KRho(Double_t s) const; | |||||
| //! Calculate the pipi phase space factor | //! Calculate the pipi phase space factor | ||||
| /*! | /*! | ||||
| \param [in] s the invariant mass squared | \param [in] s the invariant mass squared | ||||
| \return the complex phase space factor | \return the complex phase space factor | ||||
| */ | */ | ||||
| LauComplex calcPiPiRho(Double_t s) const; | LauComplex calcPiPiRho(Double_t s) const; | ||||
| //! Calculate the KK phase space factor | //! Calculate the KK phase space factor | ||||
| ▲ Show 20 Lines • Show All 74 Lines • ▼ Show 20 Lines | protected: | ||||
| //! Check the phase space factors that need to be used | //! Check the phase space factors that need to be used | ||||
| /*! | /*! | ||||
| \param [in] phaseSpaceInt phase space types | \param [in] phaseSpaceInt phase space types | ||||
| \return true of false | \return true of false | ||||
| */ | */ | ||||
| Bool_t checkPhaseSpaceType(Int_t phaseSpaceInt) const; | Bool_t checkPhaseSpaceType(Int_t phaseSpaceInt) const; | ||||
| //! Get the unitary transition amplitude matrix for the given kinematics | //! Get the unitary transition amplitude matrix for the given kinematics | ||||
| /*! | /*! | ||||
| \param [in] kinematics The pointer to the constant kinematics | \param [in] kinematics The pointer to the constant kinematics | ||||
| */ | */ | ||||
| void getTMatrix(const LauKinematics* kinematics); | void getTMatrix(const LauKinematics* kinematics); | ||||
| //! Get the unitary transition amplitude matrix for the given kinematics | //! Get the unitary transition amplitude matrix for the given kinematics | ||||
| /*! | /*! | ||||
| \param [in] s The invariant mass squared of the system | \param [in] s The invariant mass squared of the system | ||||
| */ | */ | ||||
| void getTMatrix(Double_t s); | void getTMatrix(Double_t s); | ||||
| //! Get the square root of the phase space matrix | //! Get the square root of the phase space matrix | ||||
| void getSqrtRhoMatrix(); | void getSqrtRhoMatrix(); | ||||
| private: | private: | ||||
| //! Copy constructor (not implemented) | //! Copy constructor (not implemented) | ||||
| LauKMatrixPropagator(const LauKMatrixPropagator& rhs); | LauKMatrixPropagator(const LauKMatrixPropagator& rhs); | ||||
| //! Copy assignment operator (not implemented) | //! Copy assignment operator (not implemented) | ||||
| LauKMatrixPropagator& operator=(const LauKMatrixPropagator& rhs); | LauKMatrixPropagator& operator=(const LauKMatrixPropagator& rhs); | ||||
| //! Create a map for the K-matrix parameters | //! Create a map for the K-matrix parameters | ||||
| typedef std::map<int, std::vector<LauParameter> > KMatrixParamMap; | typedef std::map<int, std::vector<LauParameter> > KMatrixParamMap; | ||||
| //! Initialise and set the dimensions for the internal matrices and parameter arrays | //! Initialise and set the dimensions for the internal matrices and parameter arrays | ||||
| void initialiseMatrices(); | void initialiseMatrices(); | ||||
| //! Store the (phase space) channel indices from a line in the parameter file | //! Store the (phase space) channel indices from a line in the parameter file | ||||
| /*! | /*! | ||||
| \param [in] theLine Vector of strings corresponding to the line from the parameter file | \param [in] theLine Vector of strings corresponding to the line from the parameter file | ||||
| */ | */ | ||||
| void storeChannels(const std::vector<std::string>& theLine); | void storeChannels(const std::vector<std::string>& theLine); | ||||
| //! Store the pole mass and couplings from a line in the parameter file | //! Store the pole mass and couplings from a line in the parameter file | ||||
| /*! | /*! | ||||
| \param [in] theLine Vector of strings corresponding to the line from the parameter file | \param [in] theLine Vector of strings corresponding to the line from the parameter file | ||||
| */ | */ | ||||
| void storePole(const std::vector<std::string>& theLine); | void storePole(const std::vector<std::string>& theLine); | ||||
| //! Store the scattering coefficients from a line in the parameter file | //! Store the scattering coefficients from a line in the parameter file | ||||
| /*! | /*! | ||||
| \param [in] theLine Vector of strings corresponding to the line from the parameter file | \param [in] theLine Vector of strings corresponding to the line from the parameter file | ||||
| */ | */ | ||||
| void storeScattering(const std::vector<std::string>& theLine); | void storeScattering(const std::vector<std::string>& theLine); | ||||
| //! Store miscelleanous parameters from a line in the parameter file | //! Store miscelleanous parameters from a line in the parameter file | ||||
| /*! | /*! | ||||
| \param [in] keyword the name of the parameter to be set | \param [in] keyword the name of the parameter to be set | ||||
| \param [in] parString the string containing the value of the parameter | \param [in] parString the string containing the value of the parameter | ||||
| */ | */ | ||||
| void storeParameter(const TString& keyword, const TString& parString); | void storeParameter(const TString& keyword, const TString& parString); | ||||
| //! String to store the propagator name | //! String to store the propagator name | ||||
| TString name_; | TString name_; | ||||
| //! Name of the input parameter file | //! Name of the input parameter file | ||||
| TString paramFileName_; | TString paramFileName_; | ||||
| //! Number to identify the DP axis in question | //! Number to identify the DP axis in question | ||||
| Int_t resPairAmpInt_; | Int_t resPairAmpInt_; | ||||
| //! Row index - 1 | //! Row index - 1 | ||||
| Int_t index_; | Int_t index_; | ||||
| //! s value of the previous pole | //! s value of the previous pole | ||||
| Double_t previousS_; | Double_t previousS_; | ||||
| //! "slowly-varying part" for the scattering K-matrix | //! "slowly-varying part" for the scattering K-matrix | ||||
| Double_t scattSVP_; | Double_t scattSVP_; | ||||
| //! "slowly-varying part" for the production K-matrix | //! "slowly-varying part" for the production K-matrix | ||||
| Double_t prodSVP_; | Double_t prodSVP_; | ||||
| //! Real part of the propagator matrix | //! Real part of the propagator matrix | ||||
| TMatrixD realProp_; | TMatrixD realProp_; | ||||
| //! Imaginary part of the propagator matrix | //! Imaginary part of the propagator matrix | ||||
| TMatrixD negImagProp_; | TMatrixD negImagProp_; | ||||
| // Integers to specify the allowed channels for the phase space calculations. | // Integers to specify the allowed channels for the phase space calculations. | ||||
| // Please keep Zero at the start and leave TotChannels at the end | // Please keep Zero at the start and leave TotChannels at the end | ||||
| // whenever more channels are added to this. | // whenever more channels are added to this. | ||||
| //! Integers to specify the allowed channels for the phase space calculations | //! Integers to specify the allowed channels for the phase space calculations | ||||
| enum KMatrixChannels {Zero, PiPi, KK, FourPi, EtaEta, EtaEtaP, | enum KMatrixChannels {Zero, PiPi, KK, FourPi, EtaEta, EtaEtaP, | ||||
| KPi, KEtaP, KThreePi, TotChannels}; | KPi, KEtaP, KThreePi, D0K_Pwave, Dstar0K_Swave, TotChannels}; | ||||
| //! Scattering K-matrix | //! Scattering K-matrix | ||||
| TMatrixD ScattKMatrix_; | TMatrixD ScattKMatrix_; | ||||
| //! Real part of the phase space density diagonal matrix | //! Real part of the phase space density diagonal matrix | ||||
| TMatrixD ReRhoMatrix_; | TMatrixD ReRhoMatrix_; | ||||
| //! Imaginary part of the phase space density diagonal matrix | //! Imaginary part of the phase space density diagonal matrix | ||||
| TMatrixD ImRhoMatrix_; | TMatrixD ImRhoMatrix_; | ||||
| //! Gmma angular-momentum barrier matrix | |||||
| TMatrixD GammaMatrix_; | |||||
| //! Identity matrix | //! Identity matrix | ||||
| TMatrixD IMatrix_; | TMatrixD IMatrix_; | ||||
| //! Null matrix | //! Null matrix | ||||
| TMatrixD zeroMatrix_; | TMatrixD zeroMatrix_; | ||||
| //! Real part of the square root of the phase space density diagonal matrix | //! Real part of the square root of the phase space density diagonal matrix | ||||
| TMatrixD ReSqrtRhoMatrix_; | TMatrixD ReSqrtRhoMatrix_; | ||||
| //! Imaginary part of the square root of the phase space density diagonal matrix | //! Imaginary part of the square root of the phase space density diagonal matrix | ||||
| TMatrixD ImSqrtRhoMatrix_; | TMatrixD ImSqrtRhoMatrix_; | ||||
| //! Real part of the unitary T matrix | //! Real part of the unitary T matrix | ||||
| TMatrixD ReTMatrix_; | TMatrixD ReTMatrix_; | ||||
| //! Imaginary part of the unitary T matrix | //! Imaginary part of the unitary T matrix | ||||
| TMatrixD ImTMatrix_; | TMatrixD ImTMatrix_; | ||||
| //! Number of channels | //! Number of channels | ||||
| Int_t nChannels_; | Int_t nChannels_; | ||||
| //! Number of poles | //! Number of poles | ||||
| Int_t nPoles_; | Int_t nPoles_; | ||||
| //! Vector of squared pole masses | //! Vector of squared pole masses | ||||
| std::vector<LauParameter> mSqPoles_; | std::vector<LauParameter> mSqPoles_; | ||||
| //! All parameters of the propagator (copuling/scattering constants) | |||||
| std::vector<LauParameter*> propParameters_; | |||||
| //! Array of coupling constants | //! Array of coupling constants | ||||
| LauParArray gCouplings_; | LauParArray gCouplings_; | ||||
| //! Array of scattering SVP values | //! Array of scattering SVP values | ||||
| LauParArray fScattering_; | LauParArray fScattering_; | ||||
| //! Vector of phase space types | //! Vector of phase space types | ||||
| std::vector<Int_t> phaseSpaceTypes_; | std::vector<Int_t> phaseSpaceTypes_; | ||||
| //! Vector of angular momentum states | |||||
| std::vector<Int_t> angMomStates_; | |||||
| //! Vector of squared masses | //! Vector of squared masses | ||||
| std::vector<Double_t> mSumSq_; | std::vector<Double_t> mSumSq_; | ||||
| //! Vector of mass differences | //! Vector of mass differences | ||||
| std::vector<Double_t> mDiffSq_; | std::vector<Double_t> mDiffSq_; | ||||
| //! Vector of 1/(m_pole^2 - s) terms for scattering and production K-matrix formulae | //! Vector of 1/(m_pole^2 - s) terms for scattering and production K-matrix formulae | ||||
| std::vector<Double_t> poleDenomVect_; | std::vector<Double_t> poleDenomVect_; | ||||
| //! Constant from input file | //! Constant from input file | ||||
| Show All 30 Lines | protected: | ||||
| //! Defined as (mK-3*mPi)^2 | //! Defined as (mK-3*mPi)^2 | ||||
| Double_t mK3piDiffSq_; | Double_t mK3piDiffSq_; | ||||
| //! Factor used to calculate the Kpipipi phase space term | //! Factor used to calculate the Kpipipi phase space term | ||||
| Double_t k3piFactor_; | Double_t k3piFactor_; | ||||
| //! Factor used to calculate the pipipipi phase space term | //! Factor used to calculate the pipipipi phase space term | ||||
| Double_t fourPiFactor1_; | Double_t fourPiFactor1_; | ||||
| //! Factor used to calculate the pipipipi phase space term | //! Factor used to calculate the pipipipi phase space term | ||||
| Double_t fourPiFactor2_; | Double_t fourPiFactor2_; | ||||
| //! Defined as (mD0+mK)^2 | |||||
| Double_t mD0KSumSq_; | |||||
| //! Defined as (mD0-mK)^2 | |||||
| Double_t mD0KDiffSq_; | |||||
| //! Defined as (mD*0+mK)^2 | |||||
| Double_t mDstar0KSumSq_; | |||||
| //! Defined as (mD*0-mK)^2 | |||||
| Double_t mDstar0KDiffSq_; | |||||
| //! Multiplicative factor containing various Adler zero constants | //! Multiplicative factor containing various Adler zero constants | ||||
| Double_t adlerZeroFactor_; | Double_t adlerZeroFactor_; | ||||
| //! Tracks if all params have been set | //! Tracks if all params have been set | ||||
| Bool_t parametersSet_; | Bool_t parametersSet_; | ||||
| //! Control the output of the functions | //! Control the output of the functions | ||||
| Bool_t verbose_; | Bool_t verbose_; | ||||
| //! Control if scattering constants are channel symmetric: f_ji = f_ij | //! Control if scattering constants are channel symmetric: f_ji = f_ij | ||||
| Bool_t scattSymmetry_; | Bool_t scattSymmetry_; | ||||
| ClassDef(LauKMatrixPropagator,0) // K-matrix amplitude model | ClassDef(LauKMatrixPropagator,0) // K-matrix amplitude model | ||||
| }; | }; | ||||
| #endif | #endif | ||||
These Doxygen lines should be attached to what is now line 169.