Index: trunk/Mars/msimreflector/MFresnelLens.cc
===================================================================
--- trunk/Mars/msimreflector/MFresnelLens.cc	(revision 19608)
+++ trunk/Mars/msimreflector/MFresnelLens.cc	(revision 19609)
@@ -122,5 +122,5 @@
    // X^2+Y^2 = r^2
    //
-   // (dz*u.x+p.x)^2 + (dz*u.y+p.y)^2 = (z-dz)^2 * tan(theta)^2
+   // (dz*u.x+p.x)^2 + (dz*u.y+p.y)^2 = (z+dz)^2 * tan(theta)^2
    //
    // dz^2*ux^2 + px^2 + 2*dz*ux*px + dz^2*uy^2 + py^2 + 2*dz*uy*py = (z^2*tan(theta)^2 + dz^2*tan(theta)^2 + 2*z*dz*tan(theta)^2
@@ -153,6 +153,6 @@
     const double dz[2] =
     {
-        (sqrtBac-B) / a,
-        (sqrtBac+B) / a
+        (sqrtBac+B) / a,
+        (sqrtBac-B) / a
     };
 
@@ -267,5 +267,5 @@
 
     //             theta                   peak_z
-    const Groove g(TMath::Pi()/2 + theta, r0*tan(theta));
+    const Groove g(TMath::Pi()/2 + theta, -r0*tan(theta));
 
     // Calculate the insersection between the ray and the cone