Thank you for these improvements, Tom. I think this is uncontroversial.

That sounds great, thanks Tom. Let me know when you're ready for me to take a look (there is no urgency - the code above is working locally for me)

Thanks for checking through this, @tlatham !

• Remove unnecessary error message

• Add catch-all return and finish B-W updates

There are quite a few things that are backward incompatible in this first iteration of the code (e.g. demanding experiment number or axis titles). I can improve that before landing.

I've just written this up in an LHCb ANA. Can I use this, and @jback's initial documentation (from Mattermost), to take a first-pass at this?

Just to say thanks very much to @jback and @tlatham for taking time to do a really thorough review. Much appreciated!

Thanks, @jback, done!

• Correct comment

• Compilation fixes

• Factorise retrieval of the phase-space density and use to obtain barrier-factor matrix for all channels

• Reset haveCalled... variable in initialiseMatrices()

• Return propagator pointer in LauIsobarDynamics::defineKMatrixPropagator
• Change proposed by Tom to ensure that K-matrix is always treated as non-narrow for integration
• Fix: set spin type correctly to Legendre for the K-matrix

• Fix 5 compile-time syntax errors

• Remove noexcept
• Merge release notes
• Restore getFloatingParameters doxygen comment
• Fix LauKMatrixProdSVP whitespace issues
• Fix parameter naming in resAmp function declaration for ProdPole and SVP
• Further whitespace fix to SVP
• Uniform default variable
• Add const variables
• Make verbose flag static constexpr (nice, thanks Tom\!)
• Add braces for Tom
• Move setting of propagator denominator parameter to setParameters(...) routine
• Correct a initialisation and move it to a local variable in storeParameters()
• Update setParameters(...) text to describe new non-zero-spin-related input lines in configuration file for K-matrix
• Simplify default assignation for L_ and radii_ vectors
• Reinstate D39 name for pole mass sq param
• Reinstate D39 name for scatt constant param
• Reverse change to error message
• Use atoi for angular momenta not atof
• Add switch on phasespaceindex inside calcGamma(...) to catch undefined calls to that function

Do you have an example of using this generalisation of the K-matrix that could go into the examples dir?

Not yet, but very happy to add one when this implementation has been tested a bit more in my fits. Can I do this later, or shall I pause D41 landing until I've gained a bit more experience with the K-matrix in D0 K+?

• Remove noexcept
• Update release notes
• Restore getFloatingParameters doxygen comment
• Fix LauKMatrixProdSVP whitespace issues
• Fix parameter naming in resAmp function declaration for ProdPole and SVP
• Further whitespace fix to SVP
• Uniform default variable
• Add const variables
• Make verbose flag static constexpr (nice, thanks Tom\!)
• Add braces for Tom
• Move setting of propagator denominator parameter to setParameters(...) routine
• Correct a initialisation and move it to a local variable in storeParameters()
• Update setParameters(...) text to describe new non-zero-spin-related input lines in configuration file for K-matrix
• Simplify default assignation for L_ and radii_ vectors
• Reinstate D39 name for pole mass sq param
• Reinstate D39 name for scatt constant param
• Reverse change to error message
• Use atoi for angular momenta not atof
• Add switch on phasespaceindex inside calcGamma(...) to catch undefined calls to that function

• minor change cherry-picking earlier BelleNR commit

• Typo in LauKMatrixPropagator.cc
• Updates to K-matrix propagator to use refined barrier factor
• Clean up indentation
• Parallel updates to KMatrixPropagator.hh
• Cleanup indentation
• Revert "Allow slope of exponential non-resonant shape to vary positive (for low-stats fits)"
• rebase conflict
• rebase conflict
• D*K system should be in P-wave\!
• Changes after round of code-review with @jback
• Remove whitespace
• Remove whitespace
• Remove extraneous comment
• Fix compilation errors
• rebase conflict
• Simplify setting of denominator parameter in barrier factor
• Cache the a/R^2 calculation
• Neater access to D*0 mass information from LauResonanceMaker - thanks to Tom
• Add definition of angular momentum term for spin-2 (max spin) to K-matrix
• Remove done TODO comment and make const several function arguments
• rebase conflict
• Add barrier factors to numerator of K-matrix amplitude in ProdPole and ProdSVP, which need to know the spin of the K-matrix using the new getL() function
• Fix typo in ProdPole.cc and ProdSVP.cc
• - Changes to allow different orbital angular momenta for different channels inside the K-matrix
• Correctly initialise a_ vector. Remove calculation of angular factor from K-matrix propagator. Remove any need to update the propagator with the LauKinematics object
• Remove angular parts of K-matrix header, and add extra function to indicate the final-state channel of the various coupled channels
• Update ProdPole and SVP to use resAmp (and leave amplitude to LauAbsResonance)
• Update ProdPole and SVP headers
• rebase conflict
• Changes to LauResonanceMaker to allow retrieval of parent Blatt-Weisskopf factor without using resInfo that is not defined for K-matrix
• Changes to LauBlattWeisskopfFactor to accommodate tangential use by K-matrix to get parent factor without access to a resInfo
• Add K-matrix name to gCoupling parameters to avoid clashes in models with >1 K-matrix
• Resolve merge conflict
• Allow other parameters of the SVP part to be floated
• Fix merge conflict
• rebase conflict
• Remove duplicate par_mSq0_
• rebase conflict
• Protect all parameters of calls to resAmp for prodpole / svp with const
• resolve conflict
• Missing ;
• Movement of one comment and correct access to one floating parameter
• Remove unnecessary whitespace
• Allow floating of s0scatt and mSq0
• Cherry pick commit from 19th Jan that moved enum class KMatrixChannels and resolve associated conflicts
• Turn off auto-verbose-printout in propagator
• Move a comment - trivial
• Increase precision of verbose printout
• Clean whitespace

• Movement of one comment
• Allow floating of s0scatt and mSq0
• Cherry pick commit from 19th Jan that moved enum class KMatrixChannels and resolve associated conflicts
• Turn off auto-verbose-printout in propagator

• Allow floating of s0scatt and mSq0

• Movement of one comment and correct access to one floating parameter
• trivial reset of changes caught up in earlier commit

• Minor, aesthetic, improvements

The exact amplitude has been reproduced with another fitter. Anything left for me to do before we land?

Can we 'land' this so that we can start merging D41? Thanks!

• Protect all parameters of calls to resAmp for prodpole / svp with const
• Clean up K-matrix propagator default values. Add lots of const protection to function calls. Noexcept throughout

• Use default values at declaration in LauAbsResonance to significantly simplify the constructors

• Cleanup LauAbsResonance constructor which involved a lot of duplication for K-matrix

• Remove duplicate par_mSq0_

Thanks, @jback. A few answers to your responses, with just one that doesn't (yet) make sense to me...

• Add K-matrix name to gCoupling parameters to avoid clashes in models with >1 K-matrix
• Resolve merge conflict
• Allow other parameters of the SVP part to be floated
• Fix merge conflict

• Remove whitespace at end of file
• Fix compilation errors

I think only 1 parameter needed a changed name - @jback let me know if I missed something here though.

• Allow other parameters of the SVP to be floated

• Add K-matrix name to gCoupling parameters to avoid clashes in models with >1 K-matrix

• Correctly initialise a_ vector. Remove calculation of angular factor from K-matrix propagator. Remove any need to update the propagator with the LauKinematics object
• Remove angular parts of K-matrix header, and add extra function to indicate the final-state channel of the various coupled channels
• Update ProdPole and SVP to use resAmp (and leave amplitude to LauAbsResonance)
• Update LauAbsResonance to apply spin to constructed objects that don't use the 'normal' constructor
• Changes to LauResonanceMaker to allow retrieval of parent Blatt-Weisskopf factor without using resInfo that is not defined for K-matrix
• Changes to LauBlattWeisskopfFactor to accommodate tangential use by K-matrix to get parent factor without access to a resInfo

I've included the changes to tell the K-matrix about the orbital angular momentum of each channel, as discussed on Mattermost. I've also made 3 comments in the code for you, @jback. Thanks!

• Changes to allow different orbital angular momenta for different channels inside the K-matrix

• Changes to allow different orbital angular momenta for different channels inside the K-matrix

• Remove done TODO comment and make const several function arguments
• Add function to LauKMatrixPropagator that returns K-matrix spin, and make various function arguments const
• Add barrier factors to numerator of K-matrix amplitude in ProdPole and ProdSVP, which need to know the spin of the K-matrix using the new getL() function
• Fix typo in ProdPole.cc and ProdSVP.cc

• Correct displaced doxygen comment
• Simplify setting of denominator parameter in barrier factor from file
• Cache the a/R^2 calculation - thanks John!
• Neater access to D*0 mass information from LauResonanceMaker - thanks Tom!
• Add definition of angular momentum term for spin-2 (max spin) to K-matrix

• Correct displaced doxygen comment
• Simplify setting of denominator parameter in barrier factor from file
• Cache the a/R^2 calculation - thanks John!
• Neater access to D*0 mass information from LauResonanceMaker - thanks Tom!
• Add definition of angular momentum term for spin-2 (max spin) to K-matrix

Thanks @jback and @tlatham, I've handled all the remaining comments. There are some deeper issues to resolve (see the Laura++ mattermost channel) before we land this. D39 needs to come first, anyway, and that one is ready to go.

• Correct Doxygen comments

• D*K system should be in P-wave, not just D0K !
• Changes after round of code-review with @jback
• Remove whitespace and make minor bugfixes found during compilation tests

Thank you for the comments @jback . I've implemented them in most cases, or left a reply in the few where further input from you is needed.

• Add K-matrix name to parameter names, to allow multiple K-M's per model

• Add braces for Tom
• Allow pole mass-squared to vary

• Add doxygen lines for one function

I tested the patch and it does the job.

• Correct name of s0prod LauParameter, thus distinguishing it from s0Scatt
• Need to name scattering parameter for it to be floated

• Update release notes

• Correct name of s0prod LauParameter, thus distinguishing it from s0Scatt

• Typo in LauKMatrixPropagator.cc
• Updates to K-matrix propagator to use refined barrier factor
• Clean up indentation
• Parallel updates to KMatrixPropagator.hh
• Cleanup indentation
• Revert "Allow slope of exponential non-resonant shape to vary positive (for low-stats fits)"

Fix typo in LauKMatrixPropagator.hh

https://phab.hepforge.org/D39

• Remove more changes destined for later review

• Remove changes not intended for this review

https://phab.hepforge.org/rLAURA44c84d8721bb822b76abfd295086ac273606231e

Also extending to encompass floating of the pole-channel coupling parameters and the scattering parameters of the slowly-varying part