Sending a Gaussian beam with aberrations through Shack Hartmann about lightpipes HOT 1 OPEN

Hello Tess,

are you referring to Fred's example in Examples/Zernike/Shack_HartmannDemo[2...3].py?
As you can see from the plots, the examples work and create the expected abberated patterns.
What might be confusing, esp. in the first demo, is that the field at the plane of the sensor is stored in F, while Forvard is applied and stored as a new variable F1 and the Zernike aberration is only applied and propagated in F2, which is never plotted. So maybe the "plot2" and "plot3" examples are more instructive.

However in all cases, the order of commands has the correct outcome:

• A large empty field is created
• A small field is created with a circular lens of the right size
• The small field is inserted in a grid in the large field
• Before propagating further, i.e. in the plane of the microlens array, the aberration is applied.
• Then the propagation to the focal plane of the microlenses (i.e. the plane of the camera) is calculated with either Forvard or Fresnel

In the same plane, e.g. the plane of the microlens array, the order of phase/intensity modulations is physically irrelevant as they are complex multiplications. So applying the phase shift due to the microlens first or the aberration first has the same outcome. You can plot the unwrapped phase after applying the microlens pattern and Zernike aberration, but before propagation.

So much for the example code. In your use case, it sounds like you would like to use an arbitrary input field and then apply the microlens array, which probably makes more sense in everday applications. This can still be done. The quickest way I can think of is to construct the Shack-Hartmann field like the example, with an unaberrated beam. Then, you can save the intensity distribution and (wrapped) phase in two variables, and use MultPhase and MultIntensity on your actual field to simulate the effect of the microlens array. A final Fresnel(F, f_microlenses) should then give you the expected distribution on the camera.

Let us know if this helps, maybe I can create an example with the method I described if you get stuck.

from lightpipes.