Event shape variables measured using multijet final states in proton-proton collisions at $\sqrt{s} =$ 13 TeV