[Color online] Estimation of PHOS nonlinearity using symmetric $\pi^0$ decays defined by $|E_{\gamma,1}-E_{\gamma,2}|< 0.05(E_{\gamma,1}+E_{\gamma,2})$. Data fit with function ($\ref{fitNL}$). The final tuned nonlinearity is shown with a dashed curve.