Как использовать FFTConvPDFV1 в zfit?Python

Программы на Python
Ответить
Anonymous
 Как использовать FFTConvPDFV1 в zfit?

Сообщение Anonymous »

Я пытаюсь свернуть два PDF-файла (CrystalBall и Gauss), используя FFTConvPDFV1 в zfit.
Я использовал следующий код:

Код: Выделить всё

mass_obs = zfit.Space('mass', 0, 1000)

mu_sig = zfit.Parameter('mu_sig', 400, 100, 600)
sigma_sig = zfit.Parameter('sigma_sig', 90, 1, 100)
alpha_sig = zfit.Parameter('alpha_sig', 200, 100, 400, floating=False)  # won't be used in the fit
npar_sig = zfit.Parameter('n sig', 4, 0.1, 30, floating=False)
signal = zfit.pdf.CrystalBall(obs=mass_obs, mu=mu_sig, sigma=sigma_sig, alpha=alpha_sig, n=npar_sig)

mu = zfit.Parameter('mu_res', 0, -1, 1, floating=False)
sigma = zfit.Parameter('sigma_res', 2, 0, 10, floating=False)
gauss_res = zfit.pdf.Gauss(obs=mass_obs, mu=mu, sigma=sigma)
sig_conv = zfit.pdf.FFTConvPDFV1(signal, gauss_res)
Я рисую «сигнал» и «sig_conv», как показано ниже:

Код: Выделить всё

x = znp.linspace(0, 1000, num=10000)
y = sig_conv.pdf(x)
y1 = signal.pdf(x)

plt.figure(figsize=(8,6))
plt.plot(x, y, label='model1')
plt.plot(x, y1, label='model2')

но произошло что-то странное, фигура выглядит так:
введите описание изображения здесь
Похоже, что синяя кривая неправильная.< /p>
Если вы дадите какую-либо подсказку, мы будем очень признательны.

Подробнее здесь: https://stackoverflow.com/questions/791 ... v1-in-zfit
Ответить

Быстрый ответ

Изменение регистра текста: 
Смайлики
:) :( :oops: :roll: :wink: :muza: :clever: :sorry: :angel: :read: *x)
Ещё смайлики…
   
К этому ответу прикреплено по крайней мере одно вложение.

Если вы не хотите добавлять вложения, оставьте поля пустыми.

Максимально разрешённый размер вложения: 15 МБ.

Вернуться в «Python»