For a large NA optical system, the performance of the system needs to be analyzed through polarization ray tracing. I need to analyze transmittance apodization, diattenuation and phase retardation across the entire pupil range for each field of view, which constitute the Jones pupil. However, zemax doesn't seem to have a function for this purpose directly. I want to achieve this through Zemax Programming Language and pass
POLTRACE Hx, Hy, Px, Py, wavelength, vec, surf
Implement polarized ray tracing. But the problem is that the results returned by the keyword POLTRACE do not include transmittance apodization, diattenuation and phase delay.
Is there any way I can achieve my goal?

And, is there a way for zemax to directly implement Jones pupil? The polarized pupil function doesn't seem to serve my purpose.