ENH: Add relativistic Breit-Wigner Distribution
Hi team, gh-17505 proposes to add a new continuous distribution to stats, the Relativistic Breit-Wigner distribution, which is used in high energy physics for modeling unstable particles. The addition of this distribution was first suggested in issue gh-6414 by contributor Andrew Fowlie (andrewfowlie), who is a physicist working in the area. Discussion faltered for several years but has picked up again recently. It has been determined that this is not a reparametrization or transformation of an existing distribution. Particle physicist and developer of relevant statistical software Hans Dembinski (HDembinski) recently provided context and voiced support for its inclusion in gh-6414. Hdembinski also pointed out that there is precedent for including distributions from high energy physics due to scipy.stats.crystalball. This seemed like a fun project to work on and a good excuse to try to refresh some long-faded physics knowledge so I pieced this together. Feedback is welcome at: https://github.com/scipy/scipy/pull/17505 Thanks, Albert
On Tue, Nov 29, 2022 at 3:07 AM Albert Steppi <albert.steppi@gmail.com> wrote:
Hi team,
gh-17505 proposes to add a new continuous distribution to stats, the Relativistic Breit-Wigner distribution, which is used in high energy physics for modeling unstable particles. The addition of this distribution was first suggested in issue gh-6414 by contributor Andrew Fowlie (andrewfowlie), who is a physicist working in the area. Discussion faltered for several years but has picked up again recently.
It has been determined that this is not a reparametrization or transformation of an existing distribution. Particle physicist and developer of relevant statistical software Hans Dembinski (HDembinski) recently provided context and voiced support for its inclusion in gh-6414. Hdembinski also pointed out that there is precedent for including distributions from high energy physics due to scipy.stats.crystalball. This seemed like a fun project to work on and a good excuse to try to refresh some long-faded physics knowledge so I pieced this together.
Feedback is welcome at: https://github.com/scipy/scipy/pull/17505
Thanks Albert. The discussion on the issue is quite in-depth and interesting. I think a good case is made for inclusion, and given that multiple maintainers are willing to work on this, that's good enough I'd say. The point that it's only used in some physics subfields is relevant, but should not be the sole determining factor for rejecting a feature, especially if the maintenance cost is fairly low. If it's that useful and there are people willing to implement and review, and most maintainers who weighed in are in favor, that seems good enough to me. Cheers, Ralf
Thanks Ralf! On Tue, Nov 29, 2022 at 03:47 Ralf Gommers <ralf.gommers@gmail.com> wrote:
On Tue, Nov 29, 2022 at 3:07 AM Albert Steppi <albert.steppi@gmail.com> wrote:
Hi team,
gh-17505 proposes to add a new continuous distribution to stats, the Relativistic Breit-Wigner distribution, which is used in high energy physics for modeling unstable particles. The addition of this distribution was first suggested in issue gh-6414 by contributor Andrew Fowlie (andrewfowlie), who is a physicist working in the area. Discussion faltered for several years but has picked up again recently.
It has been determined that this is not a reparametrization or transformation of an existing distribution. Particle physicist and developer of relevant statistical software Hans Dembinski (HDembinski) recently provided context and voiced support for its inclusion in gh-6414. Hdembinski also pointed out that there is precedent for including distributions from high energy physics due to scipy.stats.crystalball. This seemed like a fun project to work on and a good excuse to try to refresh some long-faded physics knowledge so I pieced this together.
Feedback is welcome at: https://github.com/scipy/scipy/pull/17505
Thanks Albert. The discussion on the issue is quite in-depth and interesting. I think a good case is made for inclusion, and given that multiple maintainers are willing to work on this, that's good enough I'd say. The point that it's only used in some physics subfields is relevant, but should not be the sole determining factor for rejecting a feature, especially if the maintenance cost is fairly low. If it's that useful and there are people willing to implement and review, and most maintainers who weighed in are in favor, that seems good enough to me.
Cheers, Ralf
_______________________________________________ SciPy-Dev mailing list -- scipy-dev@python.org To unsubscribe send an email to scipy-dev-leave@python.org https://mail.python.org/mailman3/lists/scipy-dev.python.org/ Member address: albert.steppi@gmail.com
participants (2)
-
Albert Steppi -
Ralf Gommers