Parton probability density functions are derived via a DGLAP-type equation involving QCD and electroweak terms. The Parton-Branching (PB) method is employed in order to derive the solution and the evolution is obtained at leading order (LO) for electroweak bosons and next-to-leading order (NLO) for QCD partons.
This novel solution is performed through an initial fit of QCD parton density functions to HERA deep inelastic scattering (DIS) data, along a perturbative generation of electroweak boson densities.
Finally, the obtained densities in both collinear (PDF) and transverse-momentum dependent (TMD) form are provided in a parametrised form for phenomenological studies and its implementation in standard HEP libraries

