weight.c

Go to the documentation of this file.

/*****************************************************************************/
#include "libtpcimgio.h"
/*****************************************************************************/
 
/*****************************************************************************/
void sifWeight(
  SIF *data,
  double halflife
) {
  int i;
  double f, d;
 
  if(SIF_TEST) printf("sifWeight(*sif, %g)\n", halflife);
  /* Calculate weights */
  for(i=0; i<data->frameNr; i++) {
    if(data->trues[i]<1.0) data->trues[i]=1.0;
    f=data->x2[i]-data->x1[i]; if(f<=0.0) f=1.0;
    if(halflife<=1.0E-8) 
      d=1.0;
    else
      d=exp( ((data->x1[i]+data->x2[i])/2.0)*0.693147/halflife );
    data->weights[i]=(f*f)/(d*data->trues[i]);
    /*printf("%3d %g %g\n", i, data->trues[i], data->weights[i]);*/
  }
 
  /* Scale weights so that average weight is 1.0 */
  sifWeightNorm(data);
 
  return;
}
/*****************************************************************************/
 
/*****************************************************************************/
void sifWeightByFrames(
  SIF *data,
  double halflife
) {
  int i;
  double f, d;
 
  if(SIF_TEST) printf("sifWeightByFrames(*sif, %g)\n", halflife);
  /* Calculate weights */
  for(i=0; i<data->frameNr; i++) {
    f=data->x2[i]-data->x1[i]; if(f<=0.0) f=1.0;
    if(halflife<=1.0E-8) 
      d=1.0;
    else
      d=exp( -((data->x1[i]+data->x2[i])/2.0)*0.693147/halflife );
    data->weights[i]=f*d;
  }
 
  /* Scale weights so that average weight is 1.0 */
  sifWeightNorm(data);
 
  return;
}
/*****************************************************************************/
 
/*****************************************************************************/
void sifWeightNorm(
  SIF *d
) {
  if(SIF_TEST) printf("sifWeightNorm(*sif)\n");
  int i;
  double f=0.0;
  for(i=0; i<d->frameNr; i++) f+=d->weights[i];
  f/=(double)d->frameNr;
  for(i=0; i<d->frameNr; i++) d->weights[i]/=f;
  return;
}
/*****************************************************************************/
 
/*****************************************************************************/
void sifModerateTrues(
  SIF *sif, 
  double limit
) {
  if(SIF_TEST) printf("sifModerateTrues(*sif, %g)\n", limit);
 
  if(sif==NULL || sif->frameNr<2) return;
  if(limit<=1.0) return;
 
  int fi;
  double w, f;
  for(w=f=sif->trues[0], fi=1; fi<sif->frameNr; fi++) {
    if(sif->trues[fi]>w) w=sif->trues[fi];
    else if(sif->trues[fi]<f) f=sif->trues[fi];
  }
  if(f*limit<w) {
    for(w/=limit, fi=0; fi<sif->frameNr; fi++)
      if(sif->trues[fi]>0.0) sif->trues[fi]+=w; else sif->trues[fi]=w;
  } else {
    for(fi=0; fi<sif->frameNr; fi++)
      if(sif->trues[fi]<0.0) sif->trues[fi]=0.0;
  }
}
/*****************************************************************************/
 
/*****************************************************************************/
void sifModerateWeights(
  SIF *sif, 
  double limit
) {
  if(SIF_TEST) printf("sifModerateWeights(*sif, %g)\n", limit);
 
  if(sif==NULL || sif->frameNr<2) return;
  if(limit<=1.0) return;
 
  int fi;
  double w, f;
  w=f=nan("");
  for(fi=0; fi<sif->frameNr; fi++) {
    if(isnan(sif->weights[fi])) continue;
    if(sif->weights[fi]<=0.0) {sif->weights[fi]=0.0; continue;}
    if(isnan(w)) w=sif->weights[fi];
    if(isnan(f)) f=sif->weights[fi];
    if(sif->weights[fi]>w) w=sif->weights[fi];
    if(sif->weights[fi]<f) f=sif->weights[fi];
  }
  if(isnan(w) || isnan(f)) return;
 
  if(f*limit<w) {
    for(w/=limit, fi=0; fi<sif->frameNr; fi++)
      if(sif->weights[fi]>0.0) sif->weights[fi]+=w;
  }
}
/*****************************************************************************/
 
/*****************************************************************************/
int sifExistentCounts(
  SIF *sif
) {
  int fi, p=0, r=0;
  double v1, v2;
  if(sif==NULL || sif->frameNr<1) return 0;
  /* If just one frame, then value > 0 is fine */
  if(sif->frameNr==1) {
    if(sif->prompts[0]>0.00000001) p=1;
    if(sif->randoms[0]>0.00000001) r=2;
    return(p+r);
  }
  /* Otherwise, check also that frames have different count level */
  for(fi=1; fi<sif->frameNr; fi++) {
    v1=sif->prompts[fi]-sif->prompts[fi-1]; if(fabs(v1)>0.001) p=1;
    v2=sif->randoms[fi]-sif->randoms[fi-1]; if(fabs(v2)>0.001) r=2;
    if((p+r)>2) break;
  }
  return(p+r);
}
/*****************************************************************************/
 
/*****************************************************************************/