  The study of global event shape variables can provide sensitive tests of
predictions for multijet production in proton-proton collisions. This paper
presents a study of several event shape variables calculated using jet four
momenta in proton-proton collisions at a centre-of-mass energy of 13 TeV and
uses data recorded with the CMS detector at the LHC corresponding to an
integrated luminosity of 2.2 fb$^{-1}$. After correcting for detector effects,
the resulting distributions are compared with several theoretical predictions.
The agreement generally improves as the energy, represented by the average
transverse momentum of the two leading jets, increases.
