La dependencia de la forma del evento y de la multiplicidad de la temperatura de congelación química, el radio de congelación y el factor de saturación de extrañeza se obtienen al estudiar los rendimientos de partículas del generador de eventos Monte Carlo PYTHIA8 en colisiones de protones-protones (pp) a una energía en el centro de masa de 13 TeV. La esfericidad es una de las técnicas transversales de forma de evento para distinguir eventos jetty e isotrópicos en colisiones de alta energía y ayuda a analizar varios observables de manera más diferencial. En este estudio, las clases de esfericidad se dividen en tres categorías, a saber, (i) esfericidad integrada, (ii) isotrópica y (iii) jetty. Los parámetros de congelación química se extraen utilizando un modelo térmico estadístico como función de la clase de esfericidad y la multiplicidad de partículas cargadas en los conjuntos canónicos, canónicos de extrañeza y gran canónicos. Se observa claramente la dependencia de la multiplicidad y la clase de esfericidad de los parámetros mencionados. Una multiplicidad en el estado final, 30 en la aceptación de multiplicidad hacia adelante del detector ALICE, parece ser un límite termodinámico, donde los parámetros de congelación se vuelven casi independientes de los conjuntos. Este estudio juega un papel importante en la comprensión del mecanismo de producción de partículas en colisiones pp de alta multiplicidad a las energías del Gran Colisionador de Hadrones (LHC) en vista de una vida útil de fase hadrónica finita en sistemas pequeños.