A method for accurate modelling of the crystal response function at a crystal sub-level applied to PET reconstruction