fnseriesofmissingsamples.cc

Go to the documentation of this file.

#define TF_FNSERIESOFMISSINGSAMPLES_CC_VERSION \
  "TF_FNSERIESOFMISSINGSAMPLES_CC   V1.1"

#include "gapfunctions.h"
#include "structgapanalysis.h"
#include <iostream>
#include <tfxx/misc.h>

GapSeries seriesofmissingsamples(const Gapsofstream& gos,
                                 const CompletenessBins& cb,
                                 const bool& debug)
{
  // fetch header data
  GapSeriesHeader header(cb, gos.ID);;
  GapSeries retseries(header);
  const libtime::TRelativeTime& dt=gos.ID.dt;
  // build series
  for (Gapsofstream::Tvecofgap::const_iterator G=gos.gap.begin();
       G!=gos.gap.end(); ++G)
  {
    TFXX_debug(debug, "seriesofmissingsamples",
               "fill in gap: " << (*G));
    const libtime::TAbsoluteTime& first=G->first;
    const libtime::TAbsoluteTime& last=G->last;
    libtime::TAbsoluteTime current=first;
    // do not process if this gaps actually is a break
    if (!G->isbreak())
    {
      for (unsigned int i=0; i<cb.nbins(); ++i)
      {
        TFXX_debug(debug, "seriesofmissingsamples",
                   "bin #" << i << ": "
                   << cb.bin(i).timestring() << " - " 
                   << cb.nextbin(i).timestring());
        TFXX_debug(debug, "seriesofmissingsamples",
                   " current: " << current.timestring());
        if (cb.isinbin(current, i))
        {
          TFXX_debug(debug, "seriesofmissingsamples",
                   "is in bin " << i << ": " << current.timestring());
          libtime::TAbsoluteTime nextbin=cb.nextbin(i);
          TFXX_debug(debug, "seriesofmissingsamples",
                   TFXX_value(nextbin.timestring()) << " " <<
                   TFXX_value(last.timestring()));
          if (nextbin<=last)
          {
            unsigned int n=(nextbin-current)/dt;
            while ((current+n*dt)<nextbin) 
            { 
              TFXX_debug(debug, "seriesofmissingsamples",
                   "n: " << n <<
                   " (current+n*dt): " << 
                   libtime::TAbsoluteTime(current+n*dt).timestring() <<
                   " nextbin: " << nextbin.timestring());
              ++n; 
            }
            TFXX_debug(debug, "seriesofmissingsamples",
                 "n: " << n <<
                 " (current+n*dt): " << 
                 libtime::TAbsoluteTime(current+n*dt).timestring() <<
                 " nextbin: " << nextbin.timestring());
            retseries.gapseries(i) += n;
            current += n*dt;
          }
          else
          {
            unsigned int n=(last-current)/dt;
            while ((current+n*dt)<=last) 
            { 
              TFXX_debug(debug, "seriesofmissingsamples",
                   "n: " << n <<
                   " (current+n*dt): " << 
                   libtime::TAbsoluteTime(current+n*dt).timestring() <<
                   " last: " << last.timestring());
              ++n; 
            }
            TFXX_debug(debug, "seriesofmissingsamples",
                 "n: " << n <<
                 " (current+n*dt): " << 
                 libtime::TAbsoluteTime(current+n*dt).timestring() <<
                 " last: " << last.timestring());
            retseries.gapseries(i) += n;
            // no need to cycle further
            i=cb.nbins();
          }
        } // if (cb.isinbin(current, i))
      } // end of loop over bins
    } // if (last >= first)
  } // end of loop over gaps in gos
  if (debug)
  {
    for (unsigned int j=retseries.gapseries.first();
         j<=retseries.gapseries.last(); ++j)
    {
      std::cout << TFXX_value(j) << " " 
        << TFXX_value(retseries.gapseries(j)) << "\n";
    }
  }
  return(retseries);
} // GapSeries seriesofmissingsamples()

/* ----- END OF fnseriesofmissingsamples.cc ----- */