Changeset 16725

Timestamp:

06/05/13 16:14:53 (12 years ago)

Author:

mknoetig

Message:

edit currents.C to work with fCurrentsDiffToPrediction

File:

• trunk/Mars/fact/processing/currents.C (modified) (6 diffs)

trunk/Mars/fact/processing/currents.C

@@ +1 @@
+
+
 void currents(const char *fname, double beg=0, double end=100000)
 {
@@ -7 +9 @@
 
     Double_t time;
+    Double_t jd;
+
     Float_t med, dev;
     file.SetPtrAddress("Time",  &time);
@@ -32 +36 @@
     double diff = -1;
 
+    double calc_dev_avg  = 0;
+
     while (file.GetNextRow())
     {
         time += offset;
+        jd    = time + 2400000.5;
 
         if (time>end)
@@ -63 +70 @@
         }
 
+        //calculate the moon
+        Nova::EquPosn moon = Nova::GetLunarEquCoords(jd, 0.01);
+        // get local position of moon
+        Nova::HrzPosn hrzm = Nova::GetHrzFromEqu(moon, jd);
+
+        //calculate the light condition
+        double calt = sin(hrzm.alt*TMath::DegToRad());
+        double disk = Nova::GetLunarDisk(jd);
+        double lc = calt*pow(disk, 2.5);
+
+        double Ical = lc>0 ? 4.5+103.0*lc : 4.5;
+        calc_dev_avg += med - Ical;
+
         cnt++;
     }
@@ -71 +91 @@
             return;
 
-        cout << "result " << med_last << " 0 " << dev_last << " 0 0 0" << endl;
+        cout << "result " << med_last << " 0 " << dev_last << " 0 0 0 0" << endl;
         return;
     }
@@ -85 +105 @@
     dev_rms /= cnt;
     dev_rms = sqrt(dev_rms-dev_avg*dev_avg);
+    calc_dev_avg = calc_dev_avg / cnt;
 
-    cout << "result " << med_avg << " " << med_rms << " " << dev_avg << " " << dev_rms << " " << med_beg << " " << med_end << endl;
+    cout << "result " << med_avg << " " << med_rms << " " << dev_avg << " " << dev_rms << " " << med_beg << " " << med_end << " " << calc_dev_avg << endl;
 }