Index: /trunk/FACT++/src/makedata.cc
===================================================================
--- /trunk/FACT++/src/makedata.cc	(revision 17365)
+++ /trunk/FACT++/src/makedata.cc	(revision 17366)
@@ -112,5 +112,5 @@
 
     const mysqlpp::StoreQueryResult res =
-        Database(fDatabase).query("SELECT fRightAscension, fDeclination FROM source WHERE fSourceName='"+source_name+"'").store();
+        Database(fDatabase).query("SELECT fRightAscension, fDeclination FROM Source WHERE fSourceName='"+source_name+"'").store();
 
     // ------------- Create canvases and frames ---------------------
@@ -151,6 +151,9 @@
         const double calt = sin(hrzm.alt*M_PI/180);
 
-        const double lc = cang==0 ? -1 : calt*pow(disk, 2.2)*pow(dist, -2);
-        const double cur = lc>0 ? 4+103*lc : 4;
+        //const double lc = cang==0 ? -1 : calt*pow(disk, 2.2)*pow(dist, -2);
+        //const double cur = lc>0 ? 4+103*lc : 4;
+
+        const double lc  = sqrt(calt)*pow(disk, 2.3)*pow(cang+1, 1)*pow(dist, -2);
+        const double cur = lc>0 ? 7.2 + 69*lc : 7.2;
 
         // Relative  energy threshold prediction
Index: /trunk/FACT++/src/makeplots.cc
===================================================================
--- /trunk/FACT++/src/makeplots.cc	(revision 17365)
+++ /trunk/FACT++/src/makeplots.cc	(revision 17366)
@@ -140,5 +140,5 @@
 
     const mysqlpp::StoreQueryResult res =
-        Database(fDatabase).query("SELECT fSourceName, fRightAscension, fDeclination FROM source").store();
+        Database(fDatabase).query("SELECT fSourceName, fRightAscension, fDeclination FROM Source").store();
 
     // ------------- Create canvases and frames ---------------------
@@ -284,6 +284,9 @@
             const double calt = sin(hrzm.alt*M_PI/180);
 
-            const double lc = cang==0 ? -1 : calt*pow(disk, 2.2)*pow(dist, -2);
-            const double cur = lc>0 ? 4+103*lc : 4;
+            //const double lc = cang==0 ? -1 : calt*pow(disk, 2.2)*pow(dist, -2);
+            //const double cur = lc>0 ? 4+103*lc : 4;
+
+            const double lc  = sqrt(calt)*pow(disk, 2.3)*pow(cang+1, 1)*pow(dist, -2);
+            const double cur = lc>0 ? 7.2 + 69*lc : 7.2;
 
             // Relative  energy threshold prediction