Simulation of photoelectron spectra with anharmonicity fully included: application to the X2A2 <- X1A1 band of furan