Abstract A measurement of the dijet production cross section is reported based on proton–proton collision data collected in 2016 at $$\sqrt{s}=13\,\text {Te}\hspace{-.08em}\text {V} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msqrt> <mml:mi>s</mml:mi> </mml:msqrt> <mml:mo>=</mml:mo> <mml:mn>13</mml:mn> <mml:mspace/> <mml:mi>Te</mml:mi> <mml:mspace/> <mml:mi>V</mml:mi> </mml:mrow> </mml:math> by the CMS experiment at the CERN LHC, corresponding to an integrated luminosity of up to 36.3 $$\,\text {fb}^{-1}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mspace/> <mml:msup> <mml:mrow> <mml:mi>fb</mml:mi> </mml:mrow> <mml:mrow> <mml:mo>-</mml:mo> <mml:mn>1</mml:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> . Jets are reconstructed with the anti- $$k_{\textrm{T}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>k</mml:mi> <mml:mtext>T</mml:mtext> </mml:msub> </mml:math> algorithm for distance parameters of $$R=0.4$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:mi>R</mml:mi> <mml:mo>=</mml:mo> <mml:mn>0.4</mml:mn> </mml:mrow> </mml:math> and 0.8. Cross sections are measured double-differentially (2D) as a function of the largest absolute rapidity $$|y |_{\text {max}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mrow> <mml:mo>|</mml:mo> <mml:mi>y</mml:mi> <mml:mo>|</mml:mo> </mml:mrow> <mml:mi>max</mml:mi> </mml:msub> </mml:math> of the two jets with the highest transverse momenta $$p_{\textrm{T}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>p</mml:mi> <mml:mtext>T</mml:mtext> </mml:msub> </mml:math> and their invariant mass $$m_{1,2} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>m</mml:mi> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn> </mml:mrow> </mml:msub> </mml:math> , and triple-differentially (3D) as a function of the rapidity separation $$y^{*} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msup> <mml:mi>y</mml:mi> <mml:mrow> <mml:mrow/> <mml:mo>∗</mml:mo> </mml:mrow> </mml:msup> </mml:math> , the total boost $$y_{\text {b}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>y</mml:mi> <mml:mi>b</mml:mi> </mml:msub> </mml:math> , and either $$m_{1,2} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>m</mml:mi> <mml:mrow> <mml:mn>1</mml:mn> <mml:mo>,</mml:mo> <mml:mn>2</mml:mn> </mml:mrow> </mml:msub> </mml:math> or the average $$p_{\textrm{T}}$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:msub> <mml:mi>p</mml:mi> <mml:mtext>T</mml:mtext> </mml:msub> </mml:math> of the two jets. The cross sections are unfolded to correct for detector effects and are compared with fixed-order calculations derived at next-to-next-to-leading order in perturbative quantum chromodynamics. The impact of the measurements on the parton distribution functions and the strong coupling constant at the mass of the $${\text {Z}} $$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mi>Z</mml:mi> </mml:math> boson is investigated, yielding a value of $$\alpha _\textrm{S} (m_{{\text {Z}}}) =0.1179\pm 0.0019$$ <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML"> <mml:mrow> <mml:msub> <mml:mi>α</mml:mi> <mml:mtext>S</mml:mtext> </mml:msub> <mml:mrow> <mml:mo>(</mml:mo> <mml:msub> <mml:mi>m</mml:mi> <mml:mi>Z</mml:mi> </mml:msub> <mml:mo>)</mml:mo> </mml:mrow> <mml:mo>=</mml:mo> <mml:mn>0.1179</mml:mn> <mml:mo>±</mml:mo> <mml:mn>0.0019</mml:mn> </mml:mrow> </mml:math> .