function [ Z ] = gaussShape( xPitch,yPitch,XI,ETA )
% Generates a 2D Gaussian function profile

Z = exp(-4*(XI.*xPitch).^2/2).*exp(-4*(ETA.*yPitch).^2);
end

