% Petzval Lens parameters:
theta_x1    0.0[deg]                                             Field angle 1, x-component
theta_y1    0.0[deg]                                             Field angle 1, y-component
theta_x2    0.0[deg]                                             Field angle 2, x-component
theta_y2    6.0[deg]                                             Field angle 2, y-component
theta_x3    0.0[deg]                                             Field angle 3, x-component
theta_y3    9.0[deg]                                             Field angle 3, y-component

N_ring    12                                                     Number of hexapolar rings
P_nom     41.5[mm]                                               Nominal entrance pupil diameter

vx1       tan(theta_x1)                                          Ray direction vector 1, x-component
vy1       tan(theta_y1)                                          Ray direction vector 1, y-component
vx2       tan(theta_x2)                                          Ray direction vector 2, x-component
vy2       tan(theta_y2)                                          Ray direction vector 2, y-component
vx3       tan(theta_x3)                                          Ray direction vector 3, x-component
vy3       tan(theta_y3)                                          Ray direction vector 3, y-component
vz        1                                                      Ray direction vector, z-component

z_stop    Tc_1+T_1+Tc_2+T_2                                      Stop z-coordinate
z_image   Tc_1+T_1+Tc_2+T_2+Tc_3+T_3+Tc_4+T_4+Tc_5+T_5+Tc_6+T_6  Image plane nominal z-coordinate

P1_fac    -1.142                                                 Pupil shift constant 1
P2_fac    -0.080                                                 Pupil shift constant 2
P_fac1    P1_fac+P2_fac*sin(sqrt(theta_x1^2+theta_y1^2))         Pupil shift factor, field 1
P_fac2    P1_fac+P2_fac*sin(sqrt(theta_x2^2+theta_y2^2))         Pupil shift factor, field 2
P_fac3    P1_fac+P2_fac*sin(sqrt(theta_x3^2+theta_y3^2))         Pupil shift factor, field 3
dx1       (dz+P_fac1*z_stop)*tan(theta_x1)                       Pupil center, field 1, x-coordinate
dy1       (dz+P_fac1*z_stop)*tan(theta_y1)                       Pupil center, field 1, y-coordinate
dx2       (dz+P_fac2*z_stop)*tan(theta_x2)                       Pupil center, field 2, x-coordinate
dy2       (dz+P_fac2*z_stop)*tan(theta_y2)                       Pupil center, field 2, y-coordinate
dx3       (dz+P_fac3*z_stop)*tan(theta_x3)                       Pupil center, field 3, x-coordinate
dy3       (dz+P_fac3*z_stop)*tan(theta_y3)                       Pupil center, field 3, y-coordinate
dz        -5[mm]                                                 Pupil center, z-component