We develop a novel Monte Carlo parton branching algorithm based on the Gribov-Levin-Ryskin (GLR) equation. The formulations of both forward evolution and backward evolution for the GLR equation are presented. The results from the Monte Carlo implementation of the GLR equation are in full agreement with its numerical solutions. Our work thus paves the way for developing an event generator that embodies the saturation effect.

• Received 30 November 2022
• Accepted 11 January 2023

DOI:https://doi.org/10.1103/PhysRevD.107.016017

Published by the American Physical Society under the terms of the Creative Commons Attribution 4.0 International license. Further distribution of this work must maintain attribution to the author(s) and the published article’s title, journal citation, and DOI. Funded by SCOAP³.

Published by the American Physical Society