Changeset 14883 for trunk/FACT++

Timestamp:

02/07/13 13:37:13 (12 years ago)

Author:

tbretz

Message:

Fixed a mistake in the new formula.

File:

• trunk/FACT++/src/makeplots.cc (modified) (2 diffs)

trunk/FACT++/src/makeplots.cc

@@ -312 +312 @@
             const double calt = sin(hrzm.alt*M_PI/180);
 
-            const double lc = calt>0 ? cang*sqrt(calt)*pow(disk, 3) : -1;
-            const double cur = lc>0 ? 8+104.5*lc : 8;
+            const double lc = cang==0 ? -1 : calt*pow(disk, 3)/sqrt(cang);
+            const double cur = lc>0 ? 8.1+94.6*lc : 8.1;
 
             // Relative  energy threshold prediction
@@ -337 +337 @@
 
             if (no_limits || (cur<max_current && 90-hrz.alt<max_zd))
-                g3.SetPoint(g3.GetN(), axis, ratio*cur/8);
+                g3.SetPoint(g3.GetN(), axis, ratio*cur/8.1);
 
             if (no_limits || (cur<max_current && 90-hrz.alt<max_zd))