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

Add checksums for dynamics functions #44

Closed
redeboer opened this issue May 26, 2024 · 3 comments · Fixed by #55
Closed

Add checksums for dynamics functions #44

redeboer opened this issue May 26, 2024 · 3 comments · Fixed by #55
Assignees
Labels
⚙️ Enhancement Improvements and optimizations of existing features

Comments

@redeboer
Copy link
Collaborator

ComPWA/ampform-dpd#132 implements deserialization for Lc2ppiK.json. I can't get the model to match yet, so it would be good to have an intermediate check on the dynamics functions.

@mmikhasenko could you add some checksums for the validation point in Lb2pKg.json and Lc2ppiK.json?

@mmikhasenko
Copy link
Collaborator

mmikhasenko commented May 28, 2024

D1232_BW: 0.9851192189816919 + 7.368771683700546im
L1520_BW: 0.9987860691365062 + 4.275745929347674e-16im
L1600_BW: 2.5634585626407667 - 2.498511047573852e-15im
K892_BW: -1.3907780720661012 + 0.24088888081536192im
BlattWeisskopf_resonance_l1: 0.9138115486202572
BlattWeisskopf_resonance_l2: 0.7172636782885041
L2000_BW: 0.36358007778802165 + 3.975360709286902e-17im
D1600_BW: 0.8379622518372177 + 0.04718049804791146im
L1670_BW: 0.7429878125956567 + 1.308995543600444e-16im
BlattWeisskopf_b_decay_l1: 0.9912279006826347
L1405_Flatte: 1.3676695116556137 + 5.453150898559425e-16im
L1690_BW: 0.7027464265893806 + 1.4456550031031727e-16im
D1700_BW: 0.7373735704444893 + 0.005092865592925969im
Notebook that does it
### A Pluto.jl notebook ###
# v0.19.42

using Markdown
using InteractiveUtils

# ╔═╡ 245bddfe-1d04-11ef-2ae2-ef1e63a8a12e
# ╠═╡ show_logs = false
begin
	using Pkg: Pkg
	Pkg.activate(tempdir())
	Pkg.add([
		Pkg.PackageSpec(url="https://github.com/mmikhasenko/HadronicLineshapes.jl"),
		Pkg.PackageSpec(url="https://github.com/mmikhasenko/ThreeBodyDecaysIO.jl")
	])
	# 
	using ThreeBodyDecaysIO
	using ThreeBodyDecaysIO.ThreeBodyDecays
	using ThreeBodyDecaysIO.HadronicLineshapes
	using ThreeBodyDecaysIO.JSON
	using ThreeBodyDecaysIO.Parameters
end

# ╔═╡ 29deb2a2-f9e0-43d9-91b3-1ffaa6240063
filename = "/path_to/AmplitudeModelSerialization/models/Lc2ppiK.json"

# ╔═╡ 08bae6c1-6040-4390-8615-bc18a9a84d75
input = open(filename) do io
	JSON.parse(io)
end;

# ╔═╡ 9ae5fe51-7aaf-43b8-ad3b-6494cceb8d16
begin
	workspace = Dict{String,Any}()
	
	@unpack functions = input
	for fn in functions
	    @unpack name, type = fn
		type == "BreitWignerWidthExpLikeBugg" && continue
	    instance_type = eval(Symbol(type))
	    workspace[name] = dict2instance(instance_type, fn)
	end
end

# ╔═╡ 8d3aa585-d69b-40d4-9ea9-38c290b29faa
for k in keys(workspace)
	println(k, ": ",  workspace[k](1.5))
end

# ╔═╡ Cell order:
# ╠═245bddfe-1d04-11ef-2ae2-ef1e63a8a12e
# ╠═29deb2a2-f9e0-43d9-91b3-1ffaa6240063
# ╠═08bae6c1-6040-4390-8615-bc18a9a84d75
# ╠═9ae5fe51-7aaf-43b8-ad3b-6494cceb8d16
# ╠═8d3aa585-d69b-40d4-9ea9-38c290b29faa

@mmikhasenko
Copy link
Collaborator

for pKg, at sigma=3.4

❔
L1890_BW: 2.417204697341433 + 2.8639956760324514im
BlattWeisskopf_b_decay_l5: 0.9731465844805017
L1810_BW: -2.167205026003352 + 2.5231003086528014im
L1800_BW: -0.9632501958568631 + 2.25664857350167im
MomentumPower1: 3.4
L1520_BW: -0.46355244347919783 + 0.4575690662529205im
L2350_BW: 0.47114186161792027 + 0.0005572072559064473im
MomentumPower3: 39.303999999999995
BlattWeisskopf_b_decay_l2: 0.9947967049304253
L1600_BW: -0.5504262623492924 + 0.5935488940734557im
BlattWeisskopf_resonance_l1: 0.9813137496771572
L2110_BW: 1.0229596072487914 + 0.10109864421641898im
BlattWeisskopf_resonance_l2: 0.9412859645745634
L1820_BW: -2.1684546597475265 + 4.477929545657705im
BlattWeisskopf_b_decay_l4: 0.9823366286390378
BlattWeisskopf_resonance_l4: 0.7675036602366174
L1680_BW: -1.6555196910221934 + 0.17017771710551707im
L2100_BW: 0.9885812994879498 + 0.038734740638986874im
BlattWeisskopf_b_decay_l1: 0.9982743731749959
MomentumPower2: 11.559999999999999
BlattWeisskopf_resonance_l5: 0.6094653239092164
L1405_Flatte: -0.6471951617416996 + 0.18708267991827043im
L1690_BW: -1.2737676644936642 + 0.8481928515271395im
L1830_BW: -1.8043003526171761 + 4.769951073988454im
BlattWeisskopf_b_decay_l6: 0.9617756605290125
BlattWeisskopf_b_decay_l3: 0.9895126377303913
BlattWeisskopf_resonance_l3: 0.8737093369245293

@IlyaSegal
Copy link
Collaborator

@assert _fρ(1.1)  -1.8472929896027317 + 0.6744244890043742im
@assert _ff2(1.1)  0.30305342103185806 + 0.11181942166047641im
@assert _fρ3(1.1)  0.1589433409235323 + 0.02906252876860443im
@assert _ff0_1500(1.1)  0.13756331374474612 + 0.019296000940740514im
@assert _ff0_980(1.1)  -0.9212576583634419 + 2.1470398931994494im
@assert _ff0_980(0.4)  0.7246075861113888 + 0.07865591658208379im
@assert _ff0_980(0.13)  0.38814688954687493 + 0.015760515277007826im
@assert _fσ(1.1)  0.10172436035046228 + 1.0273440332286132im
@assert _fσ(0.135)  (0.5840174783557608 + 0.26875840850408017im)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
⚙️ Enhancement Improvements and optimizations of existing features
Projects
None yet
Development

Successfully merging a pull request may close this issue.

3 participants