tacy.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include "tpcift.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpctac.h"
/*****************************************************************************/
 
/*****************************************************************************/
int tacYRange(
  TAC *d,
  int i, 
  double *ymin,
  double *ymax,
  int *smin,
  int *smax,
  int *imin,
  int *imax
) {
  if(ymin!=NULL) *ymin=nan("");
  if(ymax!=NULL) *ymax=nan("");
  if(smin!=NULL) *smin=-1;
  if(smax!=NULL) *smax=-1;
  if(imin!=NULL) *imin=-1;
  if(imax!=NULL) *imax=-1;
  /* Check the data */
  if(d==NULL || d->sampleNr<1) return(1);
  /* Check the index */
  if(i>=d->tacNr) return(2);
 
  /* Find the min and max time */
  double *xi, *xa;
  if(d->isframe) {xi=d->x1; xa=d->x2;} else xi=xa=d->x;
  double mi, ma;
  mi=ma=nan("");
  for(int j=0; j<d->sampleNr; j++) {
    if(isnan(xi[j]) || isnan(xa[j])) continue;
    for(int k=0; k<d->tacNr; k++) {
      if(i>=0 && i!=k) continue;
      if(isnan(mi) || isnan(ma)) {
        mi=d->c[k].y[j]; ma=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(smax!=NULL) *smax=j;
        if(imin!=NULL) *imin=k;
        if(imax!=NULL) *imax=k;
        continue;
      }
      if(d->c[k].y[j]<mi) {
        mi=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(imin!=NULL) *imin=k;
      } else if(d->c[k].y[j]>ma) {
        ma=d->c[k].y[j];
        if(smax!=NULL) *smax=j;
        if(imax!=NULL) *imax=k;
      }
    }
  }
  if(ymin!=NULL) *ymin=mi;
  if(ymax!=NULL) *ymax=ma;
  if(isnan(mi) || isnan(ma)) return(2);
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/
int tacYRangeInXRange(
  TAC *d,
  int i, 
  const double xmin,
  const double xmax,
  double *ymin,
  double *ymax,
  int *smin,
  int *smax,
  int *imin,
  int *imax
) {
  if(ymin!=NULL) *ymin=nan("");
  if(ymax!=NULL) *ymax=nan("");
  if(smin!=NULL) *smin=-1;
  if(smax!=NULL) *smax=-1;
  if(imin!=NULL) *imin=-1;
  if(imax!=NULL) *imax=-1;
  /* Check the data */
  if(d==NULL || d->sampleNr<1) return(1);
  /* Check the index */
  if(i>=d->tacNr) return(2);
  /* Check x range */
  if(xmin>xmax) return(3);
 
  /* Find the min and max time */
  double *xi, *xa, x;
  if(d->isframe) {xi=d->x1; xa=d->x2;} else xi=xa=d->x;
  double mi, ma;
  mi=ma=nan("");
  for(int j=0; j<d->sampleNr; j++) {
    if(isnan(xi[j]) || isnan(xa[j])) continue;
    x=0.5*(xi[j]+xa[j]); if(x<xmin || x>xmax) continue;
    for(int k=0; k<d->tacNr; k++) {
      if(i>=0 && i!=k) continue;
      if(isnan(mi) || isnan(ma)) {
        mi=d->c[k].y[j]; ma=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(smax!=NULL) *smax=j;
        if(imin!=NULL) *imin=k;
        if(imax!=NULL) *imax=k;
        continue;
      }
      if(d->c[k].y[j]<mi) {
        mi=d->c[k].y[j]; 
        if(smin!=NULL) *smin=j;
        if(imin!=NULL) *imin=k;
      } else if(d->c[k].y[j]>ma) {
        ma=d->c[k].y[j];
        if(smax!=NULL) *smax=j;
        if(imax!=NULL) *imax=k;
      }
    }
  }
  if(ymin!=NULL) *ymin=mi;
  if(ymax!=NULL) *ymax=ma;
  if(isnan(mi) || isnan(ma)) return(2);
  return(0);
}
/*****************************************************************************/
 
/*****************************************************************************/