Skip to content
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

PhaseSpace3 and PhaseSpace3s do not agree #19

Open
VanyaBelyaev opened this issue Dec 30, 2020 · 0 comments
Open

PhaseSpace3 and PhaseSpace3s do not agree #19

VanyaBelyaev opened this issue Dec 30, 2020 · 0 comments

Comments

@VanyaBelyaev
Copy link
Member

Ostap::Math::PhaseSpace3 and Ostap::Math::PhaeSpace3s do not agree

import ROOT, math
from   ostap.core.pyrouts  import Ostap
from   ostap.logger.logger import getLogger

logger = getLogger ("ChicPi")

PS3 = Ostap.Math.PhaseSpace3  ## <-- HERE 
PS3 = Ostap.Math.PhaseSpace3s ## <-- HERE  
GeV = 1.0
MeV = 0.001 * GeV


m_chic1 = 3510.67 * MeV
m_K     = 493.677 * MeV
m_pi    = 139.570 * MeV
m_p     = 938.272 * MeV
m_Lb    = 5619.60 * MeV


ps_pi = PS3 ( m_chic1 , m_pi , m_p )
ps_K  = PS3 ( m_chic1 , m_K  , m_p )

logger.info ( "PS(pi)                       : %s "  % ps_pi ( m_Lb ) )
logger.info ( "PS(K)                        : %s "  % ps_K  ( m_Lb ) )
r = ps_pi ( m_Lb ) / ps_K ( m_Lb )

theta_C = 13.02 * math.pi / 180.0
tc2     = math.tan (  theta_C ) **2
logger.info ( "PS(pi)/PS(K)                 : %s "  % r    )
logger.info ( "tan^2 theta_C                : %s "  % tc2  )
logger.info ( "PS(pi)/PS(K) * tan^2 theta_C : %s "  % ( r  * tc2 ) )
@VanyaBelyaev VanyaBelyaev changed the title PhaseSpace3 andPhaeSpace3s do not agree PhaseSpace3 and PhaseSpace3s do not agree Dec 30, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Projects
None yet
Development

No branches or pull requests

1 participant