You have a couple of things to think about here
"How many" quarks (and anti-quarks) of initial momentum $\vec{p}$ are there?
Given a pair with momenta $p_1$ and $p_2$ what is cross-section for scattering with invariant mass $\sqrt{s}$.
The first part is answered by the PDF. At LHC energies you can reasonably take the momentum of a parton to first order as $x p_p$ where $x$ is the Bjorken variable and $p_p$ is the momentum of the proton.
The second part you have already got modulo boosting from the parton CoM frame to the lab frame. But you can ditch the boost if you can express the cross-sections in Lorentz invariant terms (i.e. in terms of the Mandelstam variables.
In total you have an integration which is notionally of the form
$$
\sigma(Q) = \int_{x_1} dx_1 F_q(x_1,Q) \int_{x_2} dx_2 F_\bar{q}(x_2,Q) \sigma_{q\bar{q}}\left(\sqrt{(\mathbf{p}_1 + \mathbf{p}_2)^2} \right) \quad ,
$$
where $\mathbf{p}_1$ and $\mathbf{p}_2$ are the 4-vectos of the two participating partons and $Q = -\sqrt{t}$. You can, of course, also integrate over $Q$ to get a total cross-section. For comparison against data you will probably need to do a little more work to get into detector coordinates.
For more precise results you need generalized PDFs so that you can take both transverse momenta and spins into account.