Changeset 18525


Ignore:
Timestamp:
06/26/16 12:26:20 (5 years ago)
Author:
Daniela Dorner
Message:
handle cases where time stamp gets stuck during the night
File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/Mars/fact/processing/dust.C

    r18521 r18525  
    1818
    1919    TGraph g;
     20    Double_t timeprev;
     21    Float_t dustprev;
     22    Int_t counter=1;
     23    Int_t day=TMath::Nint((beg+end)/2);
    2024    while (file.GetNextRow())
    21         g.SetPoint(g.GetN(), time+offset, dust);
     25    {
     26        // handle doubl values
     27        if (time==timeprev && dust==dustprev)
     28            continue;
    2229
     30        //exclude day
     31        if (time+offset<day-0.25)
     32            continue;
     33        if (time+offset>day+0.35)
     34            break;
     35
     36        //handle values where time gets stuck
     37        if (time-timeprev<1/12.*counter)
     38        {
     39            g.SetPoint(g.GetN(), timeprev+1/12.*counter+offset, dust);
     40            counter++;
     41            continue;
     42        }
     43        else
     44        {
     45            g.SetPoint(g.GetN(), time+offset, dust);
     46            counter=1;
     47        }
     48
     49        timeprev=time;
     50        dustprev=dust;
     51    }
    2352    cout << "result " << MSpline3(g).Eval((beg+end)/2) << endl;
    2453}
Note: See TracChangeset for help on using the changeset viewer.