-
Notifications
You must be signed in to change notification settings - Fork 15
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Issue in boost() #244
Labels
🐛bug
Something isn't working
Comments
Hi @ajueid, could you explain what the problem is? |
In the code, there is gamma^2 instead of gamma.
Cheers,
Benj
… On 13 Feb 2024, at 18:31, Jack Y. Araz ***@***.***> wrote:
Hi @ajueid <https://github.com/ajueid>, could you explain what the problem is?
—
Reply to this email directly, view it on GitHub <#244 (comment)>, or unsubscribe <https://github.com/notifications/unsubscribe-auth/AIWJNJ3XUW6ZN2CCKLSNTADYTOPOLAVCNFSM6AAAAABDEHK662VHI2DSMVQWIX3LMV43OSLTON2WKQ3PNVWWK3TUHMYTSNBSGA3DKNBTGE>.
You are receiving this because you are subscribed to this thread.
|
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
System Settings
All operating systems
Describe the bug
There was a problem with void boost(MALorentzVector &p) {} (Line 146 in MABoost.h). The old function does not properly boost the Lorentz vector. Below you can find a fix
void boost(MALorentzVector &p) const
{
MAdouble64 bp = bx_ * p.X() + by_ * p.Y() + bz_ * p.Z();
p.SetX(p.X() - gamma_ * bx_ * p.T() + (gamma_ - 1) * bx_ * bp / b2_);
p.SetY(p.Y() - gamma_ * by_ * p.T() + (gamma_ - 1) * by_ * bp / b2_);
p.SetZ(p.Z() - gamma_ * bz_ * p.T() + (gamma_ - 1) * bz_ * bp / b2_);
p.SetT(gamma_ * (p.T() - bp));
}
To Reproduce
All modes
Expected behaviour
No response
Log files
No response
Additional information
No response
The text was updated successfully, but these errors were encountered: