mftac.c

Go to the documentation of this file.

/*****************************************************************************/
#include "tpcclibConfig.h"
/*****************************************************************************/
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
/*****************************************************************************/
#include "tpcextensions.h"
#include "tpcmodels.h"
#include "tpctac.h"
/*****************************************************************************/
#include "tpctacmod.h"
/*****************************************************************************/
 
/*****************************************************************************/ 
int mfCreateTAC(
  PAR *par,
  double endx,
  double dx,
  TAC *tac,
  TPCSTATUS *status
) {
  int verbose=0; if(status!=NULL) verbose=status->verbose;
  if(verbose>2) {
    printf("%s(par, %g, %g, tac, status)\n", __func__, endx, dx);
    fflush(stdout);
  }
 
  /* Check provided data */
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_FAIL);
  if(par==NULL || tac==NULL || !(endx>0.0) || !(dx>0.0)) return(status->error);
  if(par->tacNr<1) {
    statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_NO_DATA);
    return(status->error);
  }
 
  /* Make sure that TAC is empty */
  tacFree(tac);
 
  /* Allocate memory for TAC(s) */
  int xNr=1+endx/dx;  
  int ret=tacAllocate(tac, xNr, par->tacNr); 
  if(ret!=TPCERROR_OK) {
    statusSet(status, __func__, __FILE__, __LINE__, ret);
    return(status->error);
  }
  tac->sampleNr=xNr;
  tac->tacNr=par->tacNr;
 
  /* Fill x (times) */
  tac->isframe=0;
  for(int i=0; i<xNr; i++) tac->x[i]=(double)i*dx;
 
  /* Calculate curves */
  for(int ci=0; ci<par->tacNr; ci++) {
    strcpy(tac->c[ci].name, par->r[ci].name);
    if(verbose>6)
      printf("'%s' model=%u code=%s\n", par->r[ci].name, par->r[ci].model, modelCode(par->r[ci].model));
    if(mfEvalY(modelCode(par->r[ci].model), modelParNr(par->r[ci].model), par->r[ci].p, 
               tac->sampleNr, tac->x, tac->c[ci].y, verbose-3))
    {
      statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_INCOMPATIBLE_DATA);
      tacFree(tac); return(status->error);
    }
  }
 
  /* Try to find units */
  int i;
  i=iftFindKey(&par->h, "unit", 0); if(i<0) i=iftFindKey(&par->h, "calibration_unit", 0);
  if(i>=0) tac->cunit=unitIdentify(par->h.item[i].value);
  i=iftFindKey(&par->h, "timeunit", 0); if(i<0) i=iftFindKey(&par->h, "time_unit", 0);
  if(i>=0) tac->tunit=unitIdentify(par->h.item[i].value);
 
  statusSet(status, __func__, __FILE__, __LINE__, TPCERROR_OK);
  return(TPCERROR_OK);
}
/*****************************************************************************/
 
/*****************************************************************************/