We study a sample of 28 S0 galaxies extracted from the integral-field spectroscopic (IFS) survey CALIFA. We combine an accurate two-dimensional (2D) multi-component photometric decomposition with the IFS kinematic properties of their bulges to understand their formation scenario. Our final sample is representative of S0s with high stellar masses (Mstar/Msun>1010). They lay mainly on the red sequence and live in relatively isolated environments similar to that of the field and loose groups. We use our 2D photometric decomposition to define the size and photometric properties of the bulges, as well as their location within the galaxies. We perform mock spectroscopic simulations mimicking our observed galaxies to quantify the impact of the und...