bpmprocess/basic_stats.c

Go to the documentation of this file.

#include <math.h>

#include <bpm/bpm_messages.h>
#include <bpm/bpm_process.h>

int basic_stats( int *wf, int ns, int range, int nbits, double *offset, double *rms , 
                 int *max, int *min, int *unsat_sample) {
  
  double avg = 0., sqr_sum = 0.;
  int i;

  if ( ! wf || ! offset ) {
    bpm_error( "Invalid pointer argument in get_pedestal(...)",
               __FILE__, __LINE__ );
    return BPM_FAILURE;
  }
  
  if ( range == 0 || ns == 0 ) {
    bpm_error("Maximum range or number of samples equal zero in get_pedestal(...)",
              __FILE__, __LINE__ );
    return BPM_FAILURE;
  }

  if ( ns < range ) {
    bpm_warning("Number of samples < max range in get_pedestal(...)", 
                __FILE__, __LINE__ );
    range = ns;
  }
  
  // find the mean and keep track of sum of squares
  *min = 1e6;
  *max = 0;
  for ( i = 0; i < range; i++) {
    avg     += (double) wf[i];
    sqr_sum += (double) wf[i]*wf[i];
    if( wf[i] > *max) *max = wf[i];
    if( wf[i] < *min) *min = wf[i];
  }

  *offset = avg / (double) range;
  if ( rms ) {
    *rms = sqrt( sqr_sum / (double) range - SQR( *offset ) );
  } else {
    bpm_warning( "NULL pointer given for rms in get_pedestal(...)",
                 __FILE__, __LINE__ );
  }

  handle_saturation( wf, ns, ns - 20, nbits, 64, unsat_sample );

  return BPM_SUCCESS;
}

---