Index: trunk/Mars/hawc/fresnellens_traceray.C
===================================================================
--- trunk/Mars/hawc/fresnellens_traceray.C	(revision 19648)
+++ trunk/Mars/hawc/fresnellens_traceray.C	(revision 19649)
@@ -32,24 +32,23 @@
     const double D = 54.92;  // [cm] Diameter of the refractive surface
     const double F = 50.21;  // [cm] Nominal focal length of the lens
-    const double Z = F+H/2;  // [cm] camera position w.r.t. lens exit (see also MGeomCamFAMOUS!)
 
-    const double R = D/2;    // [cm] radius of lens
+    const double H = 0.25;   // [cm] Width of a single groove
     const double w = 0.01;   // [cm] Width of a single groove
-    const double H = 0.25;   // [cm] Width of a single groove
     //const double w = 1;        // [cm] Width of a single groove
     //const double H = 2.5;    // [cm] Thickness of lens
 
-    const double lambda = 546; // [nm] Wavelength for lens surface definition
+    const double Z = F+H;    // [cm] camera position w.r.t. lens exit (see also MGeomCamFAMOUS!)
+    const double R = D/2;    // [cm] radius of lens
 
-    double angle = 0;  // Check ix=0 at angle=-55 and N=10
+    const double lambda0 = 546; // [nm] Wavelength for lens surface definition
+    const double lambda  = 546; // [nm] Wavelength for the simulated photons
 
-    double Z0 = 3;
+    double angle = 6; // Angle of incidence of the simulated rays
 
-    const double n0 = MFresnelLens::RefractiveIndex(lambda);
+    double Z0 = 3;    // Starting Z-position of rays
 
     MFresnelLens lens;
-    lens.DefineLens(F, D, w, H, lambda);
+    lens.DefineLens(F, D, w, H, lambda0);
     //lens.SetPSF(psf);
-    //lens.DefineLens(F, D, w, H, lambda);
     //lens.EnableSlopeAbsorption();
     //lens.EnableDraftAbsorption();