Skip to content

Commit

Permalink
Merge branch 'documentation' of https://github.com/RUB-EP1/AmplitudeM…
Browse files Browse the repository at this point in the history
…odelSerialization into documentation
  • Loading branch information
mmikhasenko committed May 10, 2024
2 parents 649d047 + c87411a commit 5b3db52
Showing 1 changed file with 5 additions and 5 deletions.
10 changes: 5 additions & 5 deletions docs/julia/lc2pkpi.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -108,7 +108,7 @@ let
# compute, compare
_parameters = array2dict(parameter_point["parameters"];
key="name", apply=v -> v["value"])
#
#
# @assert value ≈ dist(_parameters) "Check-point validation failed with $distribution 🥕"
status = (value ≈ dist(_parameters)) ? "🟢" : "🔴"
(; Distribution=distribution,
Expand Down Expand Up @@ -156,10 +156,10 @@ let k = 2
i, j = ij_from_k(k)
xlab = "m($i$j) [GeV]"
model = model_dist.model
#
#
mlims = sqrt.(lims(k, masses(model)))
mv = range(mlims..., 150) |> shift_by_half
#
#
plot()
plot!(mv, lab="Total") do m
I = Base.Fix1(unpolarized_intensity, model)
Expand Down Expand Up @@ -188,7 +188,7 @@ The contribution of different resonances to the overall decay process is quantit
let
model = model_dist.model
ms = masses(model)
#
#
x2 = rand(10000, 2)
data = map(eachslice(x2; dims=1)) do (x, y)
σ1 = lims1(ms)[1] + x * diff(lims1(ms) |> collect)[1]
Expand All @@ -198,7 +198,7 @@ let
filter!(data) do σs
Kibble(σs, ms^2) < 0
end
#
#
chain_names = Set(model.names) |> collect |> sort
_int_i = map(chain_names) do name
_calv = unpolarized_intensity.(model[model.names.==name] |> Ref, data)
Expand Down

0 comments on commit 5b3db52

Please sign in to comment.