search

ComPWA / gluex-nstar

Investigation of N* (and K*) in GlueX photoproduction experiment
https://compwa.github.io/gluex-nstar
MIT License
2 stars 0 forks source link

Get the intensity result of the model from `ampform` and put it in weights of Dalitz Plot #32

Closed shenvitor closed 1 week ago

shenvitor commented 1 week ago

Problem description

Find a way to get the model.intensity to be used in weighted of hist2d in https://github.com/ComPWA/gluex-nstar/pull/28 after phase space and data sample Dalitz plots

Proposed solution

maybe something more than model.intensity

What should the interface look like?

No response

Additional context

No response

redeboer commented 1 week ago

You mean you should only plot the model by histogramming over phase space with the model (intensity function) as phase space, right? Data generation is not needed for now, as we only want to investigate model behavior (see also #31).

shenvitor commented 1 week ago

The data sample is extra I know. Just thinking about how to make the model use the phase space sample to output the intensity then I feed it to the weights in hist2d

redeboer commented 1 week ago

Essentially just e.g.

plt.hist2d(
    phsp["m12"],
    phsp["m23"],
    weights=intensity_func(phsp),
)

with phsp a DataSample.

shenvitor commented 1 week ago

Closes by https://github.com/ComPWA/gluex-nstar/pull/28/commits/7585cc5d3c5a1be4b5f427c0f934029bfb8cc356