roclib/src_nxyter/NxDataSummary.cxx (r4864/r4112)

Go to the documentation of this file.

//============================================================================
//============================================================================

#include <iostream>

#include "boost/format.hpp"

#include "nxyter/NxDataSummary.h"

using boost::format;

//----------------------------------------------------------------------------

nxyter::NxDataSummary::NxDataSummary() :
  fNumEntries(128),
  fFull(128),
  fMedian(128),
  fWidth50(128),
  fWidth80(128),
  fMin(128),
  fMax(128)
{
}

//----------------------------------------------------------------------------

nxyter::NxDataSummary::~NxDataSummary()
{
}

//----------------------------------------------------------------------------

void nxyter::NxDataSummary::clear()
{
  for (int ch=0; ch<128; ch++) {
    fNumEntries[ch] = 0;
    fFull[ch]       = false;
    fMedian[ch]     = 0.;
    fWidth50[ch]    = 0.;
    fWidth80[ch]    = 0.;
    fMin[ch]        = 0.;
    fMax[ch]        = 0.;
  }
}

//----------------------------------------------------------------------------

void nxyter::NxDataSummary::analyse(DistFuncArray& dfa)
{
  for (int ch=0; ch<128; ch++) {
    fNumEntries[ch] = dfa[ch].numEntries();
    fFull[ch]       = dfa[ch].full();
    fMedian[ch]     = dfa[ch].getMedian();
    fWidth50[ch]    = dfa[ch].getWidth(0.25);
    fWidth80[ch]    = dfa[ch].getWidth(0.40);
    fMin[ch]        = dfa[ch].getMin();
    fMax[ch]        = dfa[ch].getMax();
  }
}

//----------------------------------------------------------------------------

void nxyter::NxDataSummary::print(std::ostream& os, int nxnum) const
{
  os << "nx cha   hits   median  w25-75  w10-90   min   max" << std::endl;
  int chaseen = 0;
  for (int cha=0; cha<128; cha++) {
    int nent = numEntries(cha);
    if (nent == 0) continue;
    chaseen += 1;
    float med = getMedian(cha);
    float w50 = getWidth50(cha);
    float w80 = getWidth80(cha);
    float min = getMin(cha);
    float max = getMax(cha);
    char cpl  = full(cha) ? '+' : ' ';
    static format fmt(" %1d %3d %6d%c  %6.1f  %6.1f  %6.1f  %4.0f  %4.0f");
    os << fmt % nxnum % cha % nent % cpl % med % w50 % w80 % min % max << std::endl;
  }
  if (chaseen == 0) {
    os << std::dec << nxnum << " ... --- no hits seen for all channels ---" 
       << std::endl;
  }
}

//----------------------------------------------------------------------------