Pair production of scalar leptoquarks at the LHC to NLO parton shower accuracy

We present the scalar leptoquark pair production process at the LHC computed at the next-to-leading order in QCD, matched to the PYTHIA parton shower using the MC@NLO formalism. We consider the leading order decay of a leptoquark to a lepton (e, mu, tau or nu(e), nu(mu), nu(tau)) and a jet and observe the effects of parton shower on the final states. For demonstration, we display the kinametical distributions of a selection of observables along with their scale uncertainties for the 13 TeV LHC. We also present pair production cross sections and K-factors using a massless five-quark flavor scheme for different LHC center-of-mass energies. The complete stand-alone code is available online.