mftac.c File Reference

Make mathematical function based TAC. More...

#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"

Go to the source code of this file.

Functions
int	mfCreateTAC (PAR *par, double endx, double dx, TAC *tac, TPCSTATUS *status)
	Make TAC(s) based on mathematical functions in PAR format.

Detailed Description

Make mathematical function based TAC.

Definition in file mftac.c.

Function Documentation

mfCreateTAC()

int mfCreateTAC	(	PAR *	par,
		double	endx,
		double	dx,
		TAC *	tac,
		TPCSTATUS *	status )

Make TAC(s) based on mathematical functions in PAR format.

See also

mfEvalY, parExampleTTACs

Returns

enum tpcerror (TPCERROR_OK when successful).

Author

Vesa Oikonen

Parameters

par	Pointer to PAR structure containing functions for TACs.
endx	End time (x) of created TAC data. TAC will always start at zero.
dx	Step size (delta x).
tac	Pointer to initiated TAC structure. Any previous contents are deleted.
status	Pointer to status data; enter NULL if not needed.

Definition at line 26 of file mftac.c.

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