{
"cells": [
{
"cell_type": "code",
"execution_count": 1,
"id": "7e7b6a8c-747b-4169-906e-dad81bde5f61",
"metadata": {},
"outputs": [],
"source": [
"# start by importing root\n",
"import ROOT"
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "81568129-ec6f-4fa0-8825-860690970500",
"metadata": {},
"outputs": [],
"source": [
"def RooConvolution():\n",
" w = ROOT.RooWorkspace(\"w\")\n",
"\n",
" w.factory(\"Landau::landau(t[-10,30],ml[5,-20,20],sl[1,0.1,10])\")\n",
" w.factory(\"Gaussian::gauss(t,mg[0],sg[2,0.1,10])\")\n",
" t = w.var(\"t\")\n",
" landau = w.pdf(\"landau\")\n",
" gauss = w.pdf(\"gauss\")\n",
"\n",
" # Construct convoluted pdf lxg(x) = landau(x) (*) gauss(x)\n",
" lxg = ROOT.RooFFTConvPdf (\"lxg\", \"landau (X) gauss\", t, landau, gauss)\n",
"\n",
" data = lxg.generate(t, 20000)\n",
"\n",
" tframe = t.frame()\n",
" data.plotOn(tframe)\n",
" lxg.plotOn(tframe)\n",
"\n",
" # Let's plot\n",
" c1 = ROOT.TCanvas(\"c1\",\"\",1000,600)\n",
" tframe.Draw()\n",
" c1.Draw()\n",
" return c1"
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "dfa40c1a-9ac3-4ff1-aa0e-6a6e6de470e7",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"[#1] INFO:Caching -- Changing internal binning of variable 't' in FFT 'lxg' from 100 to 930 to improve the precision of the numerical FFT. This can be done manually by setting an additional binning named 'cache'.\n",
"[#1] INFO:Eval -- RooRealVar::setRange(t) new range named 'refrange_fft_lxg' created with bounds [-10,30]\n",
"[#1] INFO:Caching -- RooAbsCachedPdf::getCache(lxg) creating new cache 0x55cf801d8ac0 with pdf landau_CONV_gauss_CACHE_Obs[t]_NORM_t for nset (t) with code 0\n"
]
},
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"
\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"c1 = RooConvolution()"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "fd95cfaa-c3b2-496d-a5dc-e379ae9f216f",
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"\n",
"\n",
"
\n",
"\n",
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%jsroot on\n",
"c1.Draw()"
]
},
{
"cell_type": "code",
"execution_count": null,
"id": "c650cbf6-1e75-4731-bb91-4408a396ec26",
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.13.13"
}
},
"nbformat": 4,
"nbformat_minor": 5
}