Multi-mode Gaussian fitting based on ZEMAX