When a system deviates from equilibrium, it is possible to manipulate and control it to drive it towards equilibrium within a finite time tf, even reducing its natural relaxation time scale τrelax. Although numerous theoretical and experimental studies have explored these shortcut protocols, few have yielded analytical results for the probability distribution of work, heat and entropy production. In this paper, we propose a two-step protocol that captures the essential characteristics of more general protocols and has analytical solution for the relevant thermodynamic probability distributions. Additionally, we demonstrate that for very short protocol duration, all protocols exhibit universal behavior in their shape and for the ratio of probability distribution functions of positive and negative work, heat and the entropy production.