completeness()

CompletenessSeries completeness	(	const GapSeries &	gs,
		const bool &	debug = false
	)

convert gaps to completeness

Definition at line 46 of file fncompleteness.cc.

References CompletenessSeries::completeness, GapSeriesHeader::completenessbins, Gapid::dt, GapSeries::gapseries, GapSeries::header, CompletenessSeries::header, GapSeriesHeader::ID, and CompletenessBins::samplesinbin().

Referenced by main(), and printgaps().

{
  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)

Here is the call graph for this function:

Here is the caller graph for this function: