Changeset 18525

Timestamp:

06/26/16 12:26:20 (8 years ago)

Author:

Daniela Dorner

Message:

handle cases where time stamp gets stuck during the night

File:

• trunk/Mars/fact/processing/dust.C (modified) (1 diff)

trunk/Mars/fact/processing/dust.C

@@ -18 +18 @@
 
     TGraph g;
+    Double_t timeprev;
+    Float_t dustprev;
+    Int_t counter=1;
+    Int_t day=TMath::Nint((beg+end)/2);
     while (file.GetNextRow())
-        g.SetPoint(g.GetN(), time+offset, dust);
+    {
+        // handle doubl values
+        if (time==timeprev && dust==dustprev)
+            continue;
 
+        //exclude day
+        if (time+offset<day-0.25)
+            continue;
+        if (time+offset>day+0.35)
+            break;
+
+        //handle values where time gets stuck
+        if (time-timeprev<1/12.*counter)
+        {
+            g.SetPoint(g.GetN(), timeprev+1/12.*counter+offset, dust);
+            counter++;
+            continue;
+        }
+        else
+        {
+            g.SetPoint(g.GetN(), time+offset, dust);
+            counter=1;
+        }
+
+        timeprev=time;
+        dustprev=dust;
+    }
     cout << "result " << MSpline3(g).Eval((beg+end)/2) << endl;
 }