  The event-shape and multiplicity dependence of the chemical freeze-out
temperature ($T_{\text{ch}}$), freeze-out radius ($R$), and strangeness
saturation factor ($\gamma_{s}$) are obtained by studying the particle yields
from the PYTHIA8 Monte Carlo event generator in proton-proton (pp) collisions
at the centre-of-mass $\sqrt{s}$ = 13 TeV. Spherocity is one of the transverse
event-shape techniques to distinguish jetty and isotropic events in high-energy
collisions and helps in looking into various observables in a more differential
manner. In this study, spherocity~classes are divided into three categories,
namely (i) spherocity integrated, (ii) isotropic, and~(iii) jetty. The~chemical
freeze-out parameters are extracted using a statistical thermal model as a
function of the spherocity class and charged particle multiplicity in the
canonical, strangeness canonical, and grand canonical ensembles. A clear
observation of the multiplicity and spherocity class dependence of
$T_{\text{ch}}$, $R$, and $\gamma_{s}$ is observed. A final state multiplicity,
$N_{\rm ch}\geq$ 30 in the forward multiplicity acceptance of the ALICE
detector appears to be a thermodynamic limit, where the freeze-out parameters
become almost independent of the ensembles. This~study plays an important role
in understanding the particle production mechanism in high-multiplicity pp
collisions at the Large Hadron Collider (LHC) energies in view of a finite
hadronic phase lifetime in small systems.
