fncompleteness.cc

Go to the documentation of this file.

#define TF_FNCOMPLETENESS_CC_VERSION \
  "TF_FNCOMPLETENESS_CC   V1.0   "

#include "gapfunctions.h"
#include <tfxx/misc.h>
#include <aff/seriesoperators.h>

CompletenessSeries completeness(const GapSeries& gs, 
                                const bool& debug)
{
  TFXX_debug(debug, "completeness",
             "just entered function" << "\n" << gs.header.ID);
  CompletenessSeries retseries(gs.header);
  retseries.completeness.copyin(gs.gapseries);
  retseries.header=gs.header;
  double samplesperbin
    =static_cast<double>(gs.header.completenessbins.samplesinbin(gs.header.ID.dt));
  retseries.completeness *= -100./samplesperbin;
  retseries.completeness += 100.;
  // handle first and last bin separately
  for (unsigned int j=0; j<2; ++j)
  {
    unsigned int i=0;
    if (j==1) { i=gs.gapseries.last(); }
    samplesperbin
      =static_cast<double>(gs.header.completenessbins.samplesinbin(i, gs.header.ID.dt));
    TFXX_debug(debug,  "completeness",
               TFXX_value(j) << " " << TFXX_value(samplesperbin));
    retseries.completeness(i)=
      100.*(1.-static_cast<double>(gs.gapseries(i))/samplesperbin);
  }
  if (debug)
  {
    for (unsigned int j=retseries.completeness.first();
         j<=retseries.completeness.last(); ++j)
    {
      std::cout << TFXX_value(j) << " " 
        << TFXX_value(retseries.completeness(j)) << "\n";
    }
  }
  TFXX_debug(debug, "completeness",
             "about to leave function");
  return retseries;
} // CompletenessSeries completeness(const GapSeries& gs)

/* ----- END OF fncompleteness.cc ----- */