Abstract
The Nearest Neighbors estimation of likelihood is implemented for the multivariate amplitude analysis of the cascades of particle decays. In this approach Monte Carlo simulated events that are used to describe data are assigned weights dependent on the parameters of the theoretical fitting model. Each of such weighted events is considered as contributing to N-dimensional probability density function (p.d.f.) in its neighborhood in the parameter space. No analytic p.d.f. or cpu-intensive re-generation of N-dimensional p.d.f. at each fit step are required. The method allows accurate accounting for reconstruction effects, demands modest cpu resources and can be effectively multi-threaded. General setup of the fit as well as toy example of the amplitude analysis of B-meson decays are demonstrated.
References
[1] LHCb Collaboration, Observation of j/ψp resonances consistent with pentaquark states in Λ0b → j/ψK−p decays, Phys. Rev. Lett. 115 (2015) 072001. doi:10.1103/PhysRevLett.115.072001.
[2] W. Verkerke, D. Kirkby, The roofit toolkit for data modeling (07 2003). doi:10.1142/9781860948985_0039.
[3] BESIII Collaboration, Study of the decay ψ(3686) → Σ0Σ¯0ϕ, Phys. Rev. D 111 (2025) 052012. doi:10.1103/PhysRevD.111.052012.
[4] Belle Collaboration, Observation of a new charged charmonium-like state in B¯0→J/ψK−π+ decays, Phys. Rev. D 90 (2014) 112009. doi:10.1103/PhysRevD.90.112009.
[5] M. F. Schilling, Multivariate two-sample tests based on nearest neighbors, J. Am. Stat. Assoc.81 (395) (1986) 799–806. doi:10.1080/01621459.1986.10478337.
[6] N. Henze, A multivariate two-sample test based on the number of nearest neighbor type coincidences, Ann. Stat. 16 (2) (1988) 772–783. doi:10.1214/aos/1176350835.
[7] T. Sjöstrand, S. Mrenna, P. Skands, A brief introduction to pythia 8.1, Comp. Phys. Commun. 178 (11) (2008) 852–867. doi:https://doi.org/10.1016/j.cpc.2008.01.036.