Analytic and Monte Carlo studies of jets with heavy mesons and quarkonia

Abstract: We study jets with identified hadrons in which a family of jet-shape variables called angularities are measured, extending the concept of fragmenting jet functions (FJFs) to these observables. FJFs determine the fraction of energy, z, carried by an identified hadron in a jet with angularity, τ_a. The FJFs are convolutions of fragmentation functions (FFs), evolved to the jet energy scale, with perturbatively calculable matching coefficients. Renormalization group equations are used to provide resummed calculations with next-to-leading logarithm prime (NLL’) accuracy. We apply this formalism to two-jet events in e⁺e⁻ collisions with B mesons in the jets, and three-jet events in which a J/ψ is produced in the gluon jet. In the case of B mesons, we use a phenomenological FF extracted from e⁺e⁻ collisions at the Z⁰ pole evaluated at the scale μ = m_b. For events with J/ψ, the FF can be evaluated in terms of Non-Relativistic QCD (NRQCD) matrix elements at the scale μ = 2m_c. The z and τ_a distributions from our NLL’ calculations are compared with predictions from monte carlo event generators. While we find consistency between the predictions for B mesons and the J/ψ distributions in τ_a, we find the z distributions for J/ψ differ significantly. We describe an attempt to merge PYTHIA showers with NRQCD FFs that gives good agreement with NLL’ calculations of the z distributions.