Forward & Far-Forward Heavy Hadrons with J ETHAD : A High-energy Viewpoint

Inspired by recent findings that semi-inclusive detections of heavy hadrons exhibit fair stabilization patterns in high-energy resummed distributions against (missing) higher-order corrections, we review and extend our studies on the hadroproduction of light and heavy hadrons tagged in forward and far-forward rapidity ranges. We analyze the NLL∕NLO + behavior of rapidity rates and angular multiplicities via the J ETHAD method, where the resummation of next-to-leading energy logarithms and beyond is consistently embodied in the collinear picture. We explore kinematic regions that are within LHC typical acceptances, as well as novel sectors accessible thanks the combined tagging of a far-forward light or heavy hadron at future Forward Physics Facilities and a of central particle at LHC experiments via a precise timing-coincidence setup.


Introduction
New paths in the exploration of fundamental interactions by next-generation colliders  mark the turn of a new era in particle physics.Venturing into uncharted kinematic territories permits stringent analyses of the Standard Model (SM) and direct or indirect quests for deviations from its predictions.The strong-interaction sector presents significant challenges within the SM.Here, the interplay between perturbative and nonperturbative aspects of Quantum ChromoDynamics (QCD) gives rise to unresolved puzzles regarding fundamental questions, including the origin of hadron mass and spin, as well as the behavior of QCD observables in critical kinematic regimes.
Precise examinations of the dynamics underlying strong interactions rely on two essential components.Firstly, the ability to perform increasingly accurate calculations of high-energy parton scatterings through higher-order perturbative QCD techniques.Secondly, the understanding of proton structure, dictated by the motion and spin interactions among constituent partons.A series of successes in describing data for hadron, lepton, and lepton-hadron reactions have been achieved through collinear factorization [34,35], where partonic cross sections, computed within pure perturbative QCD, are convoluted with collinear Parton Distribution Functions (PDFs).
These PDFs encode information about the likelihood of finding a parton inside the struck hadron with a specific longitudinal momentum fraction, .They evolve according to the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [36][37][38][39][40]. Collinear PDFs are well-suited for describing inclusive or semi-inclusive observables weakly sensitive to low-transverse momentum regimes.Similarly, collinear Fragmentation Functions (FFs) portray the production mechanism of identified hadrons, detailing the probability of generating a specific final-state hadron with momentum fraction  from an outgoing collinear parton with longitudinal fraction,  ≡ ∕.
However, by relying on collinear PDFs only, one overlooks information about the transverse-space distribution and motion of partons.Thus, the description provided by collinear factorization can be seen as a one-dimensional mapping of hadron properties.Achieving an accurate portrayal of low-transverse-momentum observables necessitates embracing a three-dimensional perspective, which permits the capture of intrinsic effects stemming from the transverse motion and spin of partons, and their interaction with the polarization state of the parent hadron.Such a tomographic representation of hadrons is naturally provided by the Transverse-Momentum-Dependent (TMD) factorization (see Refs. [41,42] and reference therein).Furthermore, given the nonperturbative nature of parton densities and fragmentation functions, they must be extracted from data through global fits encompassing various hadronic processes.Despite the significant successes of the collinear approach, certain kinematic regions demand a departure from the fixed-order, DGLAP-driven description.Enhanced logarithmic contributions, which enter the perturbative expansion of the strong running coupling   with increasingly higher powers, must be included to all orders to restore the convergence of the perturbative series.These logarithms, depending on the kinematic regimes, necessitate specific all-order resummation techniques.
Another crucial kinematic regime, sensitive to logarithmic enhancements, is the semi-hard (or Regge-Gribov) sector [132][133][134][135], characterized by the scale hierarchy √  ≫ {} ≪ Λ QCD , with √  being the center-of-mass energy, {} one or a set of process-characteristic energy scales, and Λ QCD the QCD hadronization scale.Here, large energy logarithms of the form ln(∕ 2 ) compensate for the narrowness of   , potentially spoiling the perturbative-series convergence.
Like in previous cases, these logarithms must be resummed to all orders.The Balitsky-Fadin-Kuraev-Lipatov (BFKL) formalism is the most powerful tool for such resummations [136][137][138][139][140]. BFKL resummation accounts for contributions proportional to    ln(∕ 2 )  up to the Leading Logarithmic (LL) level, and terms proportional to  +1  ln(∕ 2 )  at the Next-to-Leading Logarithmic (NLL) level.In the BFKL framework, the imaginary part of amplitudes (and also cross sections of inclusive processes, via to the optical theorem [141]) takes the form of a high-energy factorization, where TM-dependent functions play a key role [142][143][144].Specifically, the BFKL amplitude is a convolution between two singly off-shell emission functions (also known in the BFKL jargon as forward-production impact factors), describing the transition from each parent particle to the outgoing objects in its fragmentation region, and a Green's function, which evolves according to the BFKL integral equation.The kernel of this equation has been computed with Next-to-Leading Order (NLO) accuracy for forward scatterings [145][146][147][148][149][150][151] and partially at next-to-NLO [152][153][154][155][156][157].
A primary category of processes that act as probe channels for high-energy resummation involves single forward emissions.In these cases, the cross section for an inclusive process follows the typical BFKL-factorized structure.Specifically, when at least one hadron participates in the initial state, the impact factor governing the production of the forward identified particle is convoluted with the BFKL Green's function and a nonperturbative quantity known as the hadron impact factor.The sub-convolution of the last two components provides an operational definition of the BFKL Unintegrated Gluon Distribution (UGD).
The hadron impact factor serves as the initial-scale UGD, while the Green's function governs its low- evolution.Given the forward kinematics, the struck parton is predominantly a gluon with a small longitudinal-momentum fraction.Hence, in this context, high-energy resummation effectively amounts to low- resummation.A similar high-energy factorization formula applies to the imaginary part of the amplitude for exclusive single forward processes.This is feasible because, in the forward limit, skewness effects are suppressed, allowing for the use of the same UGD.In more general off-forward configurations, one would consider low- enhanced GPDs [187][188][189][190].
Besides single forward emissions, the low- QCD sector can also be probed via gluon-induced single central productions.In these processes, the cross sections are formulated in a pure high-energy factorized form, involving a convolution between two BFKL UGDs and a central-production impact factor, embodying a doubly off-shell coefficient function.Due to its doubly off-shell nature (involving virtual gluons  *  * ), with gluon virtualities determined by their transverse momenta, computing the coefficient function is considerably more intricate compared to the forward-case scenarios.To our knowledge, the coefficient function depicting the inclusive emission of a central, light-flavored jet is the only one known at NLO [222].
A very powerful formalism aimed at enhancing standard fixed-order computation of central processes via the low- resummation is the Altarelli-Ball-Forte (ABF) scheme [223][224][225][226][227][228][229], where   -factorization theorems [142-144, 230, 231] permits to combine DGLAP and BFKL inputs.The high-energy series is stabilized by imposing consistency conditions based on duality aspects, symmetrizing the BFKL kernel in (anti-)collinear regions of the phase space, and encoding such contributions to the running coupling.
Significant progress has been achieved in small- studies within the ABF formalism, particularly in the context of inclusive central emissions.Notable advancements include investigations into the inclusive central emissions of Higgs bosons in gluon-gluon fusion [232][233][234][235] and higher-order corrections to top-quark pair emissions [236].Applications ABF to resummed inclusive or differential distributions for Higgs-boson and heavy-flavor hadroproductions were done by means of the High Energy Large Logarithms (HELL) method [85,237,238].Moreover, this framework has been instrumental in extracting low- enhanced collinear PDFs for the first time [239][240][241].Subsequently, the information from these collinear PDFs was utilized to constrain the parameters of initial-scale unpolarized and helicity gluon TMD PDFs [242].
Another class of processes serving as promising channel whereby unveiling the onset of high-energy QCD dynamics, is represented by inclusive hadroproductions of two particles emitted with transverse momenta well above Λ QCD and strongly separated in rapidity. 1Contrarily to single forward and central processes, forward-plus-backward two-particle hadroproduction rates are sensitive to enhanced energy logarithms even at moderate values of , due to the peculiar kinematic ranges in transverse momentum and rapidity currently covered by acceptances of LHC detectors.
Thus, on the one hand, a collinear treatment still holds here.On the other hand, however, high rapidity intervals (Δ ) correspond to large TM exchanges in the -channel, leading to the emergence of energy logarithms.In such cases, a high-energy factorization treatment is needed, and it is inherently provided by the BFKL formalism.Consequently, another form of hybrid high-energy and collinear factorization is established [249][250][251][252][253][254][255] (see Refs. [256][257][258][259][260][261][262] for a close-in-spirit formalism), wherein high-energy resummed partonic cross sections are derived directly from BFKL and subsequently convoluted with collinear PDFs.This approach allows for a comprehensive description of processes occurring in kinematic regimes characterized by both large rapidity intervals and transverse momenta.
Analyses on angular correlations for light-jet and/or light-hadron detections have played a crucial role in distinguishing between high-energy resummed and fixed-order calculations.By employing asymmetric TM ranges, we have been able to decisively discriminate between these approaches, shedding light on the underlying dynamics of highenergy QCD processes [250,270,271].However, these studies have also revealed significant challenges associated with higher-order BFKL corrections.Specifically, NLL contributions have been found to be comparable in magnitude to LL terms but with opposite signs, leading to instability in the high-energy series.This sensitivity to the variations of renormalization (  ) and factorization (  ) scales, aimed at gauging the size of Missing Higher-Order Uncertainties (MHOUs), poses a significant obstacle to achieving reliable theoretical predictions.
Efforts to address these challenges have included the adoption of scale-optimization methods such as the Brodsky-Lepage-Mackenzie (BLM) optimization [243,[324][325][326] in its semi-hard oriented version [272].While this approach has shown some success in partially mitigating instabilities in azimuthal correlations, it has proven ineffective for cross sections.In particular, the optimal scales obtained using this method have often been much larger than the natural scales dictated by kinematics [250], resulting in a substantial and unphysical reduction in statistical precision.As a result, achieving precision in the study of inclusive forward-backward light-flavored objects has remained elusive despite these efforts.
A key observation is a clear signature of stabilization, which is directly linked to the distinctive behavior exhibited by Variable-Flavor Number-Scheme (VFNS) [334][335][336] collinear FFs governing the production of these singly heavyflavored particles at high transverse momentum.These findings mark a significant advancement in our understanding of high-energy QCD processes and offer promising prospects for future precision studies in this area.Subsequent analyses on vector quarkonia [254,255,337,338],   mesons [317,320], and   Q q tetraquarks [321], confirmed that this remarkable feature, known as natural stability of the BFKL resummation [315], comes out as basic property connected to the inclusive emission of any given heavy-flavored particle.
In this review, we will provide predictions for rapidity-interval rates and angular distributions for a novel selection of forward-plus-backward two-particle semi-hard reactions.These processes involve final states characterized by identified hadrons only (see Fig. 1).Specifically, the first hadron can be a charged pion or a charged  * meson.The second hadron singly -flavored particle, i.e. and inclusive state consisting of the sum of fragmentation channels to noncharmed  mesons and Λ 0  baryons.Comparing predictions from the BFKL-driven approach with those from a high-energy fixed-order treatment will help gauge the impact of high-energy resummation on top of the DGLAP approach.To this end, a numerical tool for calculating NLO cross sections for the inclusive production of two identified hadrons widely separated in rapidity in proton collisions would be needed and crucial for a systematic high-energy versus DGLAP analysis.While the LO limit for such reactions can be extracted from higher-order works [339][340][341][342][343][344], it currently cannot be compared with calculations based on our hybrid factorization due to kinematic constraints.Two-particle LO computations without resummation typically result in a back-to-back final state, which is incompatible with the asymmetric windows for observed transverse momenta in our calculations.
Therefore, to assess the weight of the NLL resummation on top of the DGLAP approach, we will compare BFKL predictions with the corresponding ones obtained from a high-energy fixed-order treatment, firstly developed to address light-flavored di-jet [270,271] and hadron-plus-jet [250] studies.It builds upon truncating of the high-energy series up to the NLO level, so that we can reconstruct the high-energy signal of a pure NLO computation.
We will consider two different rapidity configurations.The first one is a standard LHC tagging, where both particles ( ± or  * ± and the -hadron) are detected by a current LHC detector, say CMS or ATLAS.This selection provides symmetric rapidity windows for both particles, making it an ideal channel for further testing high-energy QCD dynamics, akin to previous investigations into two-particle semi-hard processes.In the second scenario, we consider the tag of the pion or the  meson in a far-forward rapidity range accessible at future Forward Physics Facilities (FPFs), while the bottomed hadron is simultaneously detected in the barrel of a current LHC detector.Our exploration fairly takes inspiration for a prospective FPF + ATLAS study [253,318] proposed in the context of the FPF program [1,2] (see Refs. [345][346][347][348][349][350][351][352] for related work).
The simultaneous tag of a far-forward object and a central one results in an asymmetric configuration between the longitudinal fractions  of the two struck partons.One parton has a large , while the other assumes more moderate  values.Thus, a coincidence between the FPF and a LHC detector give us a faultless chance to unravel not only the highenergy dynamics arising from to the very large rapidity intervals accessed, but also the connection between BFKL and threshold resummations.This interplay can shed light on the dynamics of parton interactions in this kinematic regime, providing valuable insights into the underlying physics.
The impact of double BFKL-plus-threshold logarithms for central Higgs-bosons inclusive rates is small at ongoing LHC energies, but it becomes quite relevant at the nominal ones of the Future Circular Collider (FCC) [29][30][31][32][33]. the other side, the high-energy resummed TM rates for the Higgs-plus-jet hadroproduction already deviate from the fixed-order background at current LHC energies [251].Thus, our two-particle observables are expected to be very sensitive to the co-action of the two resummations.
Future analyses at FPFs will be crucial for deepening our understanding of perturbative QCD and the structure of protons and nuclei in previously unexplored regimes.The sensitivity of FPF detectors to far-forward production of light hadrons and charmed mesons will allow us to investigate BFKL effects and gluon-recombination dynamics.Additionally, TeV-scale neutrino-induced DIS experiments at FPFs will serve as valuable probes of proton structure and the production mechanisms of heavy or light decaying hadrons.Our work on light mesons at the FPF via hybrid factorization can provide a common framework for describing the production and decays of these particles.
QCD studies are fundamental components of the multi-frontier research program at FPFs.Searches for long-lived particles, indirect detection of dark matter, sterile neutrinos, as well as investigations into the muon puzzle, lepton universality, and the connection between high-energy particle physics and modern astroparticle physics, all depend on a deep understanding of the Standard Model.Progress towards precision QCD in the kinematic sectors reachable at FPFs will be essential for driving scientific interest towards new and compelling directions.This review is organized as follows: highlights on the NLL∕NLO + hybrid factorization for our reference processes are given in Section 2; the phenomenological analysis is discussed in Section 3; conclusions and outlook are drawn in Section 4.

Hybrid factorization at work
In this Section we presents basic features of the NLL∕NLO + hybrid factorization well-adapted to the description of our reactions.After a brief introduction of the process kinematics (Section 2.1), we give details of the NLL resummed cross section (Section 2.2).Choices for collinear PDFs and FFs are explained later (Section 2.3).

Kinematics of the process
We consider the two following processes (Fig. 1) ( The upper equation of ( 1) says that a  ± meson with mass   = 139.57MeV, four-momentum   , rapidity   , and azimuthal angle   is emitted in association with a singly bottomed hadron   with four-momentum   , rapidity   , and azimuthal angle   .In the lower equation the pion is replaced by a  * ± meson with mass   = 1.968GeV, fourmomentum   , rapidity   , and azimuthal angle   .The   particle is given by an inclusive sum of all noncharmed  mesons and Λ  baryons.Initial-state partons posses four-momentum    1 and    2 , where  1,2 are the momenta of the incoming protons.The  system in Eq. ( 1) depicts all the final-state inclusively irradiated gluons.One can take  1,2 as Sudakov vectors satisfying  2 1 =  2 2 = 0 and 2 1 ⋅  2 = , to decompose the final-state four momenta as where the  subscript inclusively refer to  or .The outgoing-object longitudinal momentum fractions,  , , can be calculated by inverting the relation so that The semi-hard nature of our reactions follows () from the size of final-state transverse momenta, which are taken to respect the hierarchy Λ QCD ≪ | ⃗   , | ≪ √ , and () from imposing a large rapidity distance Δ =   −   between the  meson and the -hadron.Furthermore, to ensure the validity of a VFNS description for heavy-hadron production [334][335][336], | ⃗   , | ranges must remain sufficiently above the thresholds for DGLAP evolution determined by the charm and bottom masses.

NLO cross section resummed at NLL and beyond
Following a pure QCD collinear vision for the LO cross section of our processes (Eq.( 1)), one would take a one-dimensional convolution between the partonic hard subprocess, incoming proton PDFs and tagged hadron FFs ( Here the ,  indices run over all the parton species except for the top quark which does not hadronize, with standing for  meson and -hadron FFs.Then,  1,2 are longitudinal-momentum fractions for the partons initiation the hard scattering and  1,2 the ones for partons fragmenting to detected hadrons.The hard factor σ, ( ŝ,   ,   ) depends on the squared center-of-mass energy of the partonic collision, ŝ ≡     , and on factorization (  ) and renormalization (  ) scales.
Vice versa, to construct differential observables in our NLL∕NLO + hybrid formalism one first performs the highenergy resummation of logarithms connected to transverse-momentum exchanges in the t-channel.Subsequently, one adds collinear ingredients, namely PDFs and FFs.The NLL∕NLO + differential cross section can be rewritten as a Fourier sum of angular coefficients where  =   −   −  is the distance between the azimuthal angles of the light and the heavy hadron.The angular coefficients   are calculated within the BFKL framework and they embody the resummation of energy logarithms.A NLL-consistent formula obtained in the Modified Minimal Subtraction () renormalization scheme [353] is cast as follows (for details on the derivation see, e.g., Ref. [163]) ] , (7) with ᾱ (  ) ≡   (  )  ∕,   the number of colors and  0 the QCD -function leading coefficient.The  NLL (, ) stands for the NLL high energy kernel: where is its eigenvalue at LL and () = Γ ′ ()∕Γ().The χ(, ) NLO term was calculated in Ref. [354,355] and can be found in the Appendix A. The  NLO ℎ (, , | ⃗   |, ) emission function describe the forward inclusive emission of a given hadron, labeled as ℎ.It was obtained in Ref. [166] in the light-quark cases and embodies collinear inputs.It also can be employed for heavy hadrons within a VFNS scheme, provided that observed transverse momenta | ⃗    | are sufficiently higher than the charm (for a  meson) or bottom (for a  hadron) masses.One has with its LO limit and ĉℎ (, , | ⃗   |, ) its NLO correction (see Appendix B).The  () function in Eq. ( 7) embodies the logarithmic derivative of the two LO emission functions From Eqs. ( 7)-( 14), we gather the implementation of our hybrid factorization approach.The cross section is expressed as a factorized formula reminiscent of the BFKL formalism.In this formulation, the Green's function read as a highenergy convolution between the emission functions of the two tagged hadrons.These functions are further expressed as collinear convolutions between PDFs and FFs, along with the hard-scattering term.The '+' label indicates that our representation for angular coefficients in Eq. ( 7) encode terms beyond the NLL accuracy generated by the cross product of the NLO emission functions, ĉ (, , Expanding and truncating to the ( 3  ) order the NLL∕NLO + angular coefficients in Eq. ( 7), one gets the highenergy fixed-order (HE-NLO + ) counterpart of the BFKL-resummed cross section [250,270,271,289,312]).It is sensitive to the leading-power asymptotic dynamics of a pure NLO DGLAP computation and discards factors suppressed by inverse powers of ŝ.Our HE-NLO + expression for angular coefficients in the MS scheme [353] reads In our approach, an expansion up to contributions proportional to   (  ) takes over the BFKL exponentiated kernel.Analogous to Eq. ( 7), our fixed-order formula at high energies is denoted as HE-NLO + , meaning contributions the beyond NLL accuracy arise from the cross product of NLO emission functions.
Then, we can obtain a genuine LL∕LO formula by simply discarding all the NLO contributions, thus having In our phenomenological analysis (see Section 3) we compare observables built in terms of NLL∕NLO + , HE-NLO + , and LL∕LO angular coefficients.We fix  , scales at the natural energies given by kinematics, thus having being the transverse mass of the specific hadron ℎ.To assess the impact of MHOUs, scales will be varied from 1∕2 to 2 times their natural values, as specified by the   parameter (see Section 3.2).
We employ a two-loop QCD-coupling setup with   (   ) = 0.11707 and five quark flavors active.In the MS renormalization scheme [353], we write with We emphasize that in our approach, energy scales are inherently tied to the transverse masses of observed particles.Consequently, they consistently fall within the perturbative regime, obviating the need for any infrared enhancement of the running coupling (see, for instance, Ref. [356]).Moreover, the utilization of large-scale values shields us from a regime where the significance of the diffusion pattern (see, for instance, Refs.[357][358][359]) becomes pronounced.

Choice of collinear PDFs and FFs
As we already mentioned, due to the moderate parton  values, we build our framework upon standard collinear inputs, at NLO while the high-energy resummation is performed via BFKL at NLL.As for proton PDFs, make us of the novel NNPDF4.0NLO determination [360,361].
When considering charged pions, there is a broad array of NLO FFs available for potential use.The NNFF1.0 FFs [362,363], derived from SIA data using a neural-network approach, feature NLO gluon FFs.On the other hand, DEHSS14 sets [364] were obtained from a combination of SIA, Semi-Inclusive Deep-Inelastic Scattering (SIDIS), and proton-proton collision data.They assume a partial  (2) isospin symmetry, resulting in specific relations among the different flavor components.Meanwhile, JAM20 determinations [365] incorporate datasets from both SIA and SIDIS and are determined concurrently with collinear PDFs in DIS and fixed-target Drell-Yan measurements.These FFs adhere to a full  (2) isospin symmetry.
Recently, MAPFF1.0 functions [366,367] have been derived from SIA and SIDIS data using neural-network techniques.Notably, these FFs allow for a deviation from isospin symmetry, with separate parametrizations for   + and   + d , which varies with the momentum fraction .Gluon FFs are generated at NLO, although the data are collected at lower energies, where the gluon distribution has a more pronounced impact.Moreover, the methodology employed for extracting MAPFF1.0FFs has been extended to study the fragmentation of the Ξ − ∕ Ξ+ octet baryon [368] and to establish a new FF set for describing unidentified charged light hadrons, both from SIA and SIDIS data [369].In our phenomenological analysis, we make use of two pion NLO FF determinations: NNFF1.0 and MAPFF1.0.
As regards heavy hadrons, we employ the KKKS08 NLO set [370][371][372][373] to describe parton fragmentation into Λ  baryons.These FFs were extracted from OPAL and Belle data for SIA and mainly rely on a Bowler-like description [374] for charm and bottom flavors.We depict emissions of  flavored hadrons in terms of the KKSS07 parametrization [375] based on data of the inclusive -meson production in SIA events at CERN LEP1 and SLAC SLC and portrayed by a simple, three-parameter power-like function [376] for heavy-quark species.The KKSS19 and KKSS07 determinations use the VFNS.We remark that the employment of given VFNS PDFs or FFs is admitted in our approach, provided that typical energy scales are much larger than thresholds for the DGLAP evolution of charm and bottom quarks.As highlighted in Section 3.3, this requirement is always fulfilled. 2e note that KKSS08 and KKSS07 FF sets do not carry any quantitative information regarding the extraction uncertainty.Future investigations incorporating potential new parametrizations of Λ  and   fragmentation functions, including uncertainties, are essential to supplement our analysis of systematic errors in high-energy distributions.

Heavy hadrons with JETHAD
All predictions presented in this review were generated using JETHAD, a hybrid code that consistently integrates both PYTHON-and FORTRAN-based modules.JETHAD is specifically designed for computing, managing, and processing physical distributions defined within various formalisms [250,253,255,322].Numeric calculations of differential distributions were primarily performed using FORTRAN 2008 modular routines within JETHAD, while the built-in PYTHON 3.0 analyzer was utilized for final data analysis and interpretation.All computation of observables are conducted within the MS scheme [353].Section 3.1 briefly introduces core elements of the v0.5.2 version of JETHAD.Our strategy to gauge systematic uncertainties is discussed in Section 3.2.Final-state kinematic ranges can be found in Section 3.3.Predictions rapidityinterval rates and angular multiplicities are discussed in Sections 3.4 and 3.5, respectively.

Highlights of JETHAD v0.5.2
The inception of the JETHAD project dates back to late 2017, driven by the necessity for precise predictions of semi-hard hadron [281,283] and jet [270,273,286] sensitive final states at the LHC.Phenomenological analyses of such reactions, proposed as probe channels for high-energy resummation in QCD, required the development of a reference numerical framework dedicated to computing and analyzing high-energy related distributions.
From the fundamental core to service modules and routines, JETHAD has been designed to dynamically achieve high levels of computational performance.The multidimensional integrators within JETHAD support extensive parallel computing to actively choose the most suitable integration algorithm based on the shape of the integrand.
Any process implemented in JETHAD can be dynamically selected through an intuitive, structure based smartmanagement interface.Physical final-state particles are represented by object prototypes within this interface, where particle objects encapsulate all pertinent information about their physical counterparts, ranging from mass and charge to kinematic ranges and rapidity tags.These particle objects are initially loaded from a master database using a dedicated particle generation routine, and custom particle generation is also supported.Then, these objects are cloned into a final-state vector and injected from the integrand routine to the corresponding, process-specific module by a dedicated controller.
The flexibility in generating the physical final states is accompanied by a range of options for selecting the initial state.A unique particle-ascendancy structure attribute enables JETHAD to rapidly learn whether an object is hadroproduced, electroproduced, photoproduced, etc.This dynamic feature ensures that only relevant modules are initialized, optimizing computing-time efficiency.
JETHAD is structured as an object-based interface that is entirely independent of the specific reaction under investigation.While originally inspired by high-energy QCD and TMD factorization phenomenology, the code's design allows for easy encoding of different approaches by simply implementing novel, dedicated (super)modules.These can be straightforwardly linked to the core structure of the code by means of a natively-equipped point-to-routine system, making JETHAD a versatile, particle-physics oriented environment.
With the aim of providing the Scientific Community with a standard computation technology tailored for the management of diverse processes (described by distinct formalisms), we envision releasing the first public version of JETHAD in the medium-term future.

Uncertainty estimation
A commonly adopted approach to assess the impact of MHOUs involves examining the sensitivity of our observables to variations in the renormalization scale and the factorization scale around their natural values.It is widely acknowledged that MHOUs significantly contribute to the overall uncertainty [253].To gauge their influence, we simultaneously vary   and   around   ∕2 and 2  , with the   parameter in the figures of Section 3 denoted as   ≡   ∕  =   ∕  .
Another potential source of uncertainty arises from proton PDFs.Recent analyses of high-energy production rates suggest that selecting different PDF parametrizations and members within the same set has minimal effect [250,253,286,313].Hence, our observables will be computed using only the central member of the NNPDF4.0parametrization.Additional uncertainties may stem from a collinear improvement of the NLO kernel, which entails incorporating renormalization-group (RG) terms to align the BFKL equation with the DGLAP one in the collinear limit, or from changes in the renormalization scheme [392][393][394][395][396][397][398].The impact of collinear-improvement techniques on semi-hard rapidity-differential rates is observed to be encapsulated within the error bands generated by MHOUs [253].
Furthermore, the MS [353] to MOMentum (MOM) [399,400] renormalization-scheme transition was assessed in Ref. [253], resulting in systematically higher MOM results for rapidity distributions.However, these outcomes remain within the MHOUs bands.Notably, a proper MOM analysis should be grounded on MOM-evolved PDFs and FFs, which are presently unavailable.
To establish uncertainty bands for our distributions, we combine MHOUs with the numerical errors arising from multidimensional integration (see Section 3.3).The latter is consistently kept below 1% owing to the integrators employed in JETHAD.

Final-state kinematic ranges
By making use of formulae of Eqs. ( 7), (13), and ( 14), we construct physical observables as functions of angular coefficients integrated over the final-state phase-space variables, while the rapidity separation Δ =   −   between  the two observed particles is kept fixed.Thus, we have The   terms inclusively represents all the NLL∕NLO + , HE-NLO + , and LL∕LO integrated angular coefficients.This approach allows us to impose and investigate various windows in transverse momenta and rapidities, based on realistic kinematic configurations employed in current and future experimental studies at the LHC.In particular, we will focus on the following two kinematic ranges.

Standard LHC tagging
We delve into the emission dynamics of both hadrons within the standard detection acceptances of CMS or ATLAS barrel calorimeters.Unlike the Mueller-Navelet channel, which extends into the end-cap region (as illustrated in panel () of Fig. 2), allowing for jet tagging with | jet | < 4.7, hadron detection is confined to the barrel regions.A realistic proxy for the rapidity range of hadron tags at the LHC can be derived from recent CMS analyses on Λ  particles [401], yielding | Λ  | < 2. For our study, we adopt a slightly broader range, matching the coverage of the CMS barrel detector, Following precedent studies [283,289,312], we set a 10 < | ⃗    |∕GeV < 20 window for the transverse momentum of the  hadron.In contrast, we select a wider and disjoint transverse momentum range for the -hadron, specifically 20 < | ⃗    |∕GeV < 60, as recently suggested in studies of high-energy -flavored emissions [313,320].This selection ensures the robustness of our treatment within the VFNS framework, as energy scales consistently exceed the thresholds for DGLAP evolution of heavy-quark FFs (we refer to Section 2.3 for further details).
On the one hand, the symmetric tagging in rapidity of both light and heavy particles enables us to rigorously apply our formalism in a well-defined regime, allowing for stringent tests of the NLL∕NLO + hybrid factorization.On the other hand, as highlighted in Ref. [250], the use of disjoint transverse momentum intervals amplifies effects of additional, undetected gluon radiation by suppressing the Born contribution.This accentuates the signatures of the high-energy resummation over the fixed-order background.Furthermore, asymmetric | ⃗   |-windows mitigate potential instabilities arising in NLO calculations [402,403], as well as violations of energy-momentum at NLL [404].However, the achievable Δ values with standard LHC tagging may not be sufficiently large to unambiguously discriminate between BFKL and fixed-order signals.This challenge can be addressed by exploring far-forward emissions of one of the two particles, as proposed in Section 3.3.2.

FPF + LHC coincidence
Expanding upon the standard detection regime outlined in Section 3.3.1,we advocate for the simultaneous identification of an far-forward particle (our choice falls into the  hadron) alongside a more centrally located one (the   hadron, as depicted in panel () of Fig. 2).Once the planned FPF [1,2] becomes operational, this novel configuration might be achievable by coordinating FPF detectors with ATLAS.Integrating information from both ATLAS and the FPF hinges on the capability to trigger ATLAS events using far-forward signatures, necessitating precise timing protocols and influencing FPF detector design.Technical specifics regarding a possible FPF + ATLAS tight timing synchronization are detailed in Section VI E of Ref. [1].
As a preliminary investigation, we examine the emission of a  meson within the far-forward range given by 5 <   < 7.While greater rapidities could be explored, we adopt a conservative approach, opting for a rapidity interval disjoint from and more forward than that accessible by end-caps of standard LHC detector.Investigations into larger rapidity ranges are deferred to future endeavors.The meson is paired with a -hadron detected by the LHC barrel calorimeter within the standard rapidity spectrum, |  | < 2.4.The transverse momenta of both final-state hadrons remain consistent with those specified in Section 3.3.1.
The combined FPF + LHC tagging strategy yields a hybrid and markedly asymmetric range selection, presenting an excellent avenue for disentangling clear high-energy signals from collinear backgrounds.However, as highlighted in Refs.[252,289] and mentioned previously, the joint detection of an far-forward particle alongside a central one results in an asymmetric configuration between the longitudinal momentum fractions of the corresponding incoming partons.This configuration significantly suppresses the contribution of undetected gluon radiation at LO and has a notable impact at NLO.This kinematic constraint leads to an incomplete cancellation between virtual and real contributions from gluon emissions, resulting in the emergence of large threshold logarithms in the perturbative series.Given that the BFKL framework encompasses the resummation of energy-type single logarithms while systematically neglecting threshold ones, we anticipate a partial degradation in the convergence of our resummed calculation in this FPF + LHC coincidence setup compared to the LHC standard scenario.Despite these challenges, these features motivate future developments aimed at incorporating the threshold resummation into our framework [72][73][74][75][76][77][78][79][80][81][82][83][84][85][86][87][88][89][90].
In conclusion, explorations via a FPF + LHC coincidence method offer an unparalleled opportunity for rigorous and in-depth examinations of strong interaction dynamics in the high-energy regime.In this regard, the advent of the planned FPF [1,2] may complement the capabilities of the ATLAS detector, enabling us to (i) assess the feasibility of precision analyses using hybrid high-energy and collinear factorization, and (ii) explore potential commonalities among distinct resummation techniques.

Rapidity-interval rates
We delve into the analysis of Δ -rates for our studied processes, as illustrated in Figure 1.These distributions represent -summed, Δ -differential cross sections.We note that clues of a natural stabilization of these observables are strongly expected when a standard LHC tagging is considered, while they are awaited in a FPF + LHC coincidence setup.Observing such stability would not only be a further reliability test for the NLL∕NLO + factorization, but it also would mark a milestone toward future precision analyses.
Stabilization effects can manifest at various levels.Firstly, the involve the ability to study Δ -distributions for all considered final states around the natural energy scales prescribed by kinematics.This prerequisite is essential for claiming evidence of stability, a feat unattainable when considering light hadrons and/or jets [250,266,268,283].In such cases, large NLL contributions, of similar magnitude but opposite sign to their LL counterparts, can lead to unphysical Δ -distributions, sometimes even resulting in negative values for large Δ .Another sign of instability arises in the analysis of mean values of cosines of multiples of the azimuthal-angle distance, ⟨cos()⟩, which can exceed one.Secondly, achieving a substantial reduction in the discrepancy between pure LL calculations and NLLresummed ones would push natural stability to a higher level.
To assess the impact of our resummed calculations on fixed-order predictions, we contrast LL∕LO and NLL∕NLO + results with the corresponding ones derived using our HE-NLO + formula (Eq.13), which effectively mimics the highenergy signal of a pure NLO computation.Given a dedicated numerical tool for higher-order perturbative calculations of cross sections for semi-hard hadroproduction of two identified hadrons is not yet available, our HE-NLO + approach remains the most valid and efficient method for a BFKL versus fixed-order comparison.
Plots of Fig. 3 depict the Δ -rate for inclusive  ± +   productions at the LHC (left panels) and at FPF + LHC (right panel).Analogously, plots of Fig. 4 refer to the inclusive  * ± +   detections.Uncertainty bands are constructed considering the combined effect of MHOUs from energy-scale variation and numerical multi-dimensional integration over the final-state phase space, with the former being significantly dominant.Upon inspection, our Δ -distributions generally exhibit favorable statistics, lying in the range 10 −1 to 5 × 10 2 nb.The observed reflects the typical dynamics of our hybrid factorization.While the BFKL resummation predicts an increase in the partonic hard-scattering cross section with energy, its convolution with collinear PDFs and FFs leads to a decrease with Δ for LL∕LO, NLL∕NLO + , and NLL∕NLO + results.This decrease is more pronounced for LHC kinematic configurations and appears smoother in the FPF + LHC case.The noncontiguous rapidity ranges covered by the FPF and current LHC detectors may contribute to this smoother shape, compensating for the increment with Δ of the available phase space by the absence of detected events in the interval between  min FPF and  max LHC .A similar pattern was also noted in a complementary configuration, namely the CMS + CASTOR setup (see Fig. 10 of Ref. [250]).Results for rapidity-interval rates obtained with different pion FF parametrizations (upper versus lower panels of Fig. 3) exhibit qualitatively similar trends, with their differences remaining significant, beyond a factor of three.We observe the emergence of clear and natural stabilization effects in the high-energy series when standard LHC cuts are applied.The Δ -distributions for all the channels in Fig. 1 feature NLL∕NLO + bands that are partially or entirely nested within LL∕LO bands for small and moderate values of Δ .However, as the rapidity interval increases, NLL∕NLO corrections become increasingly negative, causing NLL∕NLO + predictions to become smaller than pure LL∕LO ones.NLL∕NLO and HE-NLO + uncertainty bands are generally narrower than NLL∕NLO and HE-NLO + ones.Furthermore, their width generally diminishes in the large Δ range, where high-energy effects dominate over pure DGLAP ones.These observations are consistent with previous analyses on heavy-flavored emissions at CMS [254,312,313], highlighting the convergence of the energy-resummed series thanks to the natural stabilizing effect of VFNS FFs, which accurately describe the hadronization mechanisms of the detected heavy-flavor species.
The high-energy stabilizing pattern is also evident in the FPF + LHC coincidence setup, although its effects are less pronounced (see right panels of Figs. 5 to 7).While our required condition for asserting evidence of stability is met, allowing for precision studies of cross sections around the natural energy scales provided by kinematics, it is notable that FPF + LHC NLL∕NLO + predictions consistently remain below LL∕LO results.Additionally, NLL∕NLO + uncertainty bands are narrower than LL∕LO ones, but slightly wider than NLL∕NLO + bands for the same channels investigated in the standard LHC configurations.Furthermore, LL∕LO results consistently exceed HE-NLO + ones, while NLL∕NLO + results are smaller.Although further dedicated studies are required to determine if the observed natural-stability signals degrade when the FPF rapidity acceptances are expanded beyond those imposed in our analysis, an explanation for the increased sensitivity of Δ -rates to the resummation accuracy can be provided based on our current understanding of the dynamics behind other resummation mechanisms.
We stress that the semi-hard nature of the final states considered leads to high energies but not necessarily to small- dynamics.This is particularly evident in the FPF + LHC coincidence setup, where the strongly asymmetric final-state rapidity ranges result in one of the two parton longitudinal fractions being consistently large, while the other takes more moderate values.As highlighted in Section 3.3.2,our approach does not capture large- logarithms, whose inclusion would be performed via an appropriate resummation mechanism, the aforementioned threshold resummation.A significant finding from Ref. [343] on inclusive di-hadron detections in hadronic collisions is that incorporating NLL threshold resummation on top of pure NLO calculations leads to a notable increase in cross sections.Remarkably, this increase is comparable to the gap between our LL∕LO and NLL∕NLO + high-energy predictions for Δ -distributions in FPF + LHC configurations.
Recent works [89,405] have demonstrated how NLL instabilities arising from forward hadron hadroproduction, described within the saturation framework [406], can be substantially mitigated when threshold logarithms are incorporated into these calculations.Considering these findings, we contend that the natural stability observed in our high-energy studies is not compromised by the adoption of FPF + LHC coincidence setups.It remains robust and provides a satisfactory description of Δ -rates at natural scales.The discrepancy between LL∕LO and NLL∕NLO + predictions could come from those large-, threshold logarithms, which are currently not accounted for by our hybrid factorization.Incorporating the large- resummation represents a crucial step toward enhancing the description of the considered observables.Therefore, it should be pursued as the next logical step to evaluate the feasibility of precision studies of Δ -distributions at FPF + LHC.
The analysis presented in this section underscores that light-meson plus heavy-flavor production processes offer a promising avenue for stabilizing the high-energy resummation, as anticipated.Δ -distributions emerge as particularly promising observables for detecting signals of high-energy dynamics and potentially discriminating between BFKLdriven and fixed-order computations.Further exploration of these distributions holds significant potential for delving into the interplay between high-energy QCD dynamics and other resummations, particularly the large- threshold effects.By probing these aspects in greater detail, we can gain deeper insights into the underlying mechanisms governing semi-hard phenomenology and refine our understanding of QCD dynamics in the high-energy regime.

Angular multiplicities
Semi-hard phenomenology delves into observables that become increasingly sensitive to final-state rapidity intervals.When these observables are also differential in azimuthal angles, it exposes a core aspect of high-energy QCD.In the context of two-particle hadroproduction reactions, significant insights into the onset of BFKL dynamics emerge when large rapidity distances (Δ ) enhance the weight of undetected gluons strongly ordered in rapidity.These gluon emissions, resummed as energy logarithms, induce a growing-with-Δ decorrelation on the azimuthal plane of the outgoing particles.
The angular decorrelation was initially observed through the Δ dependence of cross-section azimuthal moment [249,265,266,268,407], defined as the ratios  0 ≡   ∕ 0 between a specific azimuthal coefficient ( ≥1 ) and the -summed  0 coefficient, which genuinely corresponds to the Δ rate of Section 3.4.The  10 ratio effectively measures the azimuthal decorrelation between the two outgoing particles, akin to the mean value of ⟨cos ⟩.The  0 ratios represent the higher moments ⟨cos()⟩, while further probes of BFKL were proposed through ratios between azimuthal moments,   ≡   ∕  = ⟨cos()⟩∕⟨cos()⟩, in earlier studies [408,409].NLL-resummed predictions for angular correlations of Mueller-Navelet jets exhibited satisfactory agreement with LHC data at √  = 7 TeV, particularly for symmetric | ⃗   |-windows [265,266,268].However, due to instabilities affecting the NLL series for processes involving the emission of two light-flavored particles, comparisons between theory and experiment had to be conducted at energy scales optimized via different procedures [265,268,272].
Recent investigations into inclusive hadroproductions of heavy-flavored hadrons have revealed that the stabilizing effects associated with the use of heavy-flavor VFNS FFs are less pronounced when a heavy hadron is emitted alongside a light jet, as opposed to another heavy object [312,313].Consequently, only a partial reduction in instabilities for azimuthal moments is observed.
As recently shown [254,277], starting from the angular coefficients, we can construct a more stable observable that contains signals of high-energy dynamics coming from all azimuthal modes.We refer to the angular multiplicity with  ≡   −   − .This observable, initially proposed in Ref. [264] within the framework of Mueller-Navelet analyses, was further examined with NLL precision in Ref. [265] for the same process. 3Exploring its characteristics holds significant advantages from both theoretical and experimental standpoints.On the one hand, it collects signals of high-energy dynamics from all azimuthal modes, rendering it one of the most robust observables for spotting BFKL effects.As a normalized distribution, i.e., a multiplicity, its sensitivity to uncertainties arising from uncertainties propagating from PDFs and/or FFs, as well as those from different replicas within the same set, is notably diminished.This allows for a focused examination of uncertainties intrinsic to high-energy resummation, facilitating stringent BFKL tests.On the other hand, experimental detector acceptances typically do not cover the entire (2) azimuthal-angle range.Thus, comparing -dependent observables, like our azimuthal distribution, with experimental data is considerably more straightforward than for   correlation moments.
From a numerical standpoint, ensuring reliable predictions for our -distribution necessitates computing a large number of   coefficients.We assessed the numerical stability of our calculations by progressively increasing the effective upper bound of the  sum in Eq. ( 18), achieving satisfactory convergence at  [bound] ≃ 20.The prominent feature shared among all the depicted multiplicities is the emergence of a peak centered at  = 0, corresponding of the physical configuration where the  meson and the   hadron are emitted (almost) back-to-back.With increasing Δ , a characteristic trend emerges: the NLL∕NLO + peak height decreases while its width expands.This trend stems from the heightened number of secondary gluons emitted with substantial rapidity separation, as predicted by the BFKL equation.Consequently, the correlation in the azimuthal plane between the two tagged particles diminishes, leading to a decrease in the number of nearly back-to-back events.

Results for angular multiplicities as functions of
Conversely, LL∕LO distributions exhibit an opposing pattern: the peak grows while its width decreases with increasing Δ .This behavior generates from the connection between the strongly asymmetric transverse-momentum windows at which the two hadrons are tagged (see Section 3.3) and the corresponding longitudinal-momentum fractions.This bring to a reduction in the combinations of these fractions for the given Δ .However, this results in an unphysical recorrelation pattern in the azimuthal plane for LL∕LO distributions, contradicting the expected loss of correlation due to the weight of rapidity-ordered gluons forming the inclusive system .The correct behavior is reinstated when full NLL corrections are considered.
Additionally, the size of uncertainty bands due to scale variations sensibly diminishes as Δ increases.This becomes strongly manifest in FPF + LHC selections.Multiplicities taken at lower reference values of the rapidity interval exhibit two symmetric minima at || ≳ ∕2, which extend to unphysical values beyond zero, whereas no negative values are observed for larger Δ values.This indicates a natural stabilization of the high-energy series in the large rapidity-interval regime, as as expected.Upon qualitative comparison of our predictions with corresponding results studied in other semi-hard channels, novel features emerge.As an example, angular distributions for our reactions appear less peaked compared to vector quarkonium + jet distributions (Fig. 6 of Ref. [254]).They exhibit similarities with light-hadron or jet rates (Figs.17 and 18 of Ref. [250]).
In summary, multiplicities for -meson plus -hadron productions within the hybrid factorization at NLL∕NLO + allow for stringent tests of high-energy QCD dynamics.The natural stabilization of the high-energy resummation becomes significant in its expected applicability domain, particularly in the large Δ sector.This makes azimuthal distributions easily measurable in current LHC experimental configurations and future analyses facilitated by the FPF + LHC coincidence methods, this enhances our ability to conduct rigorous analyses of high-energy QCD dynamics.

Toward precision studies of high-energy QCD
By employing the hybrid high-energy and collinear factorization approach at NLL∕NLO + , we have investigated the inclusive detection of pion or a  meson, in association with a singly-bottomed hadron, at current LHC energies and kinematic configurations.Additionally, we explored configurations accessible via a FPF + LHC tight timingcoincidence setup.
Our analysis of distributions differential in the observed rapidity interval (Δ ) or in the azimuthal-angle distance () between the two tagged objects has corroborated the remarkable property of natural stability of the high-energy resummation.This feature, recently uncovered in the context of heavy-flavor studies in forward directions [254,312,313,317,321], enables a reliable description of the considered observables around the natural values of energy scales dictated by process kinematics.This stability is a prerequisite and an initial stride toward precision investigations of high-energy QCD through inclusive di-hadron system emissions in proton collisions.
We have demonstrated that Δ -rates effectively discriminate between the high-energy signal from the fixed-order background.On the other hand, -distributions exhibit robust stability in the large Δ regime, offering a means to identify new and distinctive high-energy features.The promising statistical outcomes of our observables in the FPF + LHC configuration underscore the interest of the FPF Community [1,2] in exploring the intriguing prospect of enabling FPF and LHC detectors to operate in coincidence.Achieving this will necessitate extremely precise timing procedures, the technical feasibility of which should be actively pursued and complemented by positive feedback from the theoretical domain.
A striking conclusion from our recent investigation into the interplay between BFKL and DGLAP in inclusive semi-hard emissions of light jets and hadrons at the LHC [250] is the imperative for a multi-lateral formalism.Such a framework would entail the simultaneous and consistent incorporation of several distinct resummation mechanisms, serving as a fundamental element for conducting precision studies of high-energy QCD.The sensitivity of FPF + LHC results presented in this study to both high-energy and threshold resummation underscores the urgency of developing such a unified description as a top priority in the medium-term future.
As a first prospect, we will complement our investigation on  mesons emitted in FPF-like kinematic configurations, where heavy particles are detected by LHC detectors, by examining the opposite configuration.In this setup, a far-forward heavy-flavored hadron is tagged at the FPF while another central object remains within LHC cuts.Subsequently, we will explore the high-energy behavior of observables sensitive to single inclusive emissions of heavy hadrons reconstructed by FPF detectors.Our aim is to gain access to the proton content in far-forward (very low-) regimes provided by FPF cuts.Here, our hybrid factorization framework could serve as a theoretical common basis for exploring production mechanisms and decays of heavy-flavored particles.
Then, by making use of estimates for differential cross sections of heavy-quark spectra at the particle-generated level, such as the ones provided via the FONLL method [418][419][420][421], it will be possible to make JETHAD work with parameters tuned for current high-energy hadron-hadron and lepton-hadron colliders which can produce heavy hadrons in forward directions of rapidity.
Mapping the proton structure in the very low- regime will hinge on a comprehensive exploration of the connections among different approaches.Specifically, we aim to investigate the interplay between our hybrid factorization, which permits the description of cross sections for single forward emissions in terms of a   -factorization between off-shell matrix elements and the UGD and the ABF formalism which, as mentioned before, allowed us to obtain low- enhanced collinear PDFs (see also Section 6.1.2 of Ref. [2]).
We view the analyses presented in this work as a significant step forward toward conducting precision studies of high-energy QCD.Our hybrid factorization framework, potentially enhanced via the integration of additional resummation techniques, offers a systematic approach to reducing uncertainties stemming from both perturbative calculations of high-energy scatterings and from collinear inputs.This serves a dual purpose: it provides a benchmark for SM measurements and establishes a shared foundation for the exploration of Beyond-the-Standard-Model (BSM) physics.

Appendix B: Forward-hadron emission function at NLO
The NLO correction to the forward-hadron singly off-shell emission function can be written as [166] ℎ and with 2  1 the Gauss hypergeometric function.The plus-prescription in Eqs. ( 24) and ( 27) act as on any function  ( ) regular at  = 1.

Figure 1 :
Figure 1: Diagrammatic representation of the hybrid factorization at work for the  ± +   (, left) and  * ± +   (, right) channels.The big yellow ovals denote collinear PDFs.Maroon (navy blue) blobs stand for  ± ( * ± meson collinear FFs, while green rectangles depict   hadron collinear FFs.The high-energy Green's function (red oval blob) is connected to singly off-shell coefficient functions by Reggeon zigzag lines.

Figure 2 :
Figure 2: Left: A typical forward-backward hadron-jet tag in standard LHC rapidity ranges.Right: concurrent detection of a far-forward hadron at FPFs [1, 2] and a central one at an LHC experiment thanks a narrow timing coincidence.
and for two different sets of values for the rapidity interval, Δ = 3, 4 (LHC) or Δ = 6, 7 (FPF + LHC), are shown in Figs. 6, 5, and 7. Panels in these figures are structured as follows.Upper (lowers) plots refer to NLL∕NLO + (LL∕LO) rates.Standard LHC tagging and the FPF + LHC coincidence setups are used in left and right plots, respectively.Figures 5 and 6 refer to pion channels described via NNFF1.0 and MAPFF1.0pion collinear FFs, whereas Fig. 7 show results for -meson channels depicted by means of KKKS08 FF determinations.