No Arabic abstract
Because of properties of QED, the bremsstrahlung corrections to decays of particles or resonances can be calculated, with a good precision, separately from other effects. Thanks to the widespread use of event records such calculations can be embodied into a separate module of Monte Carlo simulation chains, as used in High Energy Experiments of today. The PHOTOS Monte Carlo program is used for this purpose since nearly 20 years now. In the following talk let us review the main ideas and constraints which shaped the program version of today and enabled it widespread use. We will concentrate specially on conflicting requirements originating from the properties of QED matrix elements on one side and degrading (evolving) with time standards of event record(s). These issues, quite common in other modular software applications, become more and more difficult to handle as precision requirements become higher.
For high precision measurements of K decays, the presence of radiated photons cannot be neglected. The Monte Carlo simulations must include the radiative corrections in order to compute the correct event counting and efficiency calculations. In this paper we briefly describe a method for simulating such decays.
Because of properties of QED, the bremsstrahlung corrections to decays of particles or resonances can be calculated, with a good precision, separately from other effects. Thanks to the widespread use of event records such calculations can be embodied into a separate module of Monte Carlo simulation chains, as used in High Energy Experiments of today. The PHOTOS Monte Carlo program is used for this purpose since nearly 20 years now. In the following talk let us review the main ideas and constraints which shaped the program version of today and enabled it widespread use. Finally, we will underline importance of aspects related to reliability of program results: event record contents and implementation of channel specific matrix elements.
With the approaching start-up of the experiments at LHC, the urgency to quantify systematic uncertainties of the generators, used in the interpretation of the data, is becoming pressing. The PHOTOS Monte Carlo program is often used for the simulationof experimental, selection-sensitive, QED radiative corrections in decays of Z bosons and other heavy resonances and particles. Thanks to its complete phase-space coverage it is possible, with no approximations for any decay channel, to implement the matrix-element. The present paper will be devoted to those parts of the next-to-leading order corrections for Z decays which are normally missing in PHOTOS. The analytical form of the exact and truncated (standard) kernel used in PHOTOS will be explicitly given. The correction, being the ratio of the exact to the approximate kernel, can be activated as an optional contribution to the internal weight of PHOTOS. To calculate the weight, the information on the effective Born-level Z/gamma* couplings and even directions of the incoming beams, is needed. A universal implementation would have made the PHOTOS solution less modular and less convenient for the users. That is why, for the time being, we will keep the correcting weight as an extra option, available for special tests only. We will quantify the numerical effect of the approximation with the help of a multitude of distributions. The numerical size of the effect is in general below 0.1%; however, in some corners of the phase-space (well defined and contributing less than 0.5% to the total rate), it may reach up to about 20% of their relative size.
The algorithm for Monte Carlo simulation of parton-level events based on an Artificial Neural Network (ANN) proposed in arXiv:1810.11509 is used to perform a simulation of $Hto 4ell$ decay. Improvements in the training algorithm have been implemented to avoid numerical instabilities. The integrated decay width evaluated by the ANN is within 0.7% of the true value and unweighting efficiency of 26% is reached. While the ANN is not automatically bijective between input and output spaces, which can lead to issues with simulation quality, we argue that the training procedure naturally prefers bijective maps, and demonstrate that the trained ANN is bijective to a very good approximation.
Monte Carlo simulations are widely used in many areas including particle accelerators. In this lecture, after a short introduction and reviewing of some statistical backgrounds, we will discuss methods such as direct inversion, rejection method, and Markov chain Monte Carlo to sample a probability distribution function, and methods for variance reduction to evaluate numerical integrals using the Monte Carlo simulation. We will also briefly introduce the quasi-Monte Carlo sampling at the end of this lecture.