Index: /trunk/FACT++/src/makedata.cc
===================================================================
--- /trunk/FACT++/src/makedata.cc	(revision 16376)
+++ /trunk/FACT++/src/makedata.cc	(revision 16377)
@@ -188,6 +188,7 @@
 
         // Relative  energy threshold prediction
-        const double cs = cos((90+hrz.alt)*M_PI/180);
-        const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
+        //const double cs = cos((90+hrz.alt)*M_PI/180);
+        //const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
+        const double ratio = pow(cos((90-hrz.alt)*M_PI/180, -2.664);
 
         // Add points to curve
Index: /trunk/FACT++/src/makeplots.cc
===================================================================
--- /trunk/FACT++/src/makeplots.cc	(revision 16376)
+++ /trunk/FACT++/src/makeplots.cc	(revision 16377)
@@ -309,6 +309,7 @@
 
             // Relative  energy threshold prediction
-            const double cs = cos((90+hrz.alt)*M_PI/180);
-            const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
+            //const double cs = cos((90+hrz.alt)*M_PI/180);
+            //const double ratio = (10.*sqrt(409600.*cs*cs+9009.) + 6400.*cs - 60.)/10.;
+            const double ratio = pow(cos((90-hrz.alt)*M_PI/180, -2.664);
 
             // Add points to curve