isotope.c File Reference

Isotope halflife functions. More...

#include "tpcclibConfig.h"
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include <time.h>
#include <string.h>
#include "tpcisotope.h"

Go to the source code of this file.

Functions
double	isotopeHalflife (int isotope_code)
double	isotopeBranching (int isotope_code)
char *	isotopeName (int isotope_code)
int	isotopeIdentifyHalflife (double halflife)
int	isotopeIdentify (const char *isotope)

Detailed Description

Isotope halflife functions.

Definition in file isotope.c.

Function Documentation

isotopeBranching()

double isotopeBranching

(

int

isotope_code

)

Get the branching ratio (fraction) of an isotope spcecified with its isotope_code.

See also

isotopeHalflife, isotopeName, lambdaFromIsotope, decayCorrectionFactorFromIsotope

Returns

branching ratio of the isotope, or NaN if not identified.

Author

Vesa Oikonen

Parameters

isotope_code

isotope_code as enum or the index of isotope in isotope table

Definition at line 81 of file isotope.c.

  {
  int n=0;
  
  while(strlen(tpc_isotope[n].name)>1) n++;
  if(isotope_code<1 || isotope_code>n-1) return nan("");
  else if(tpc_isotope[isotope_code].br<1.0E-10) return nan("");
  else return(tpc_isotope[isotope_code].br);
}

isotopeHalflife()

double isotopeHalflife

(

int

isotope_code

)

Get the half-life (in minutes) of an isotope specified with its isotope_code.

See also

isotopeBranching, isotopeName, lambdaFromIsotope, decayCorrectionFactorFromIsotope

Returns

halflife of the isotope in minutes, or NaN if not identified.

Author

Vesa Oikonen

Parameters

isotope_code

isotope_code as enum or the index of isotope in isotope table.

Definition at line 62 of file isotope.c.

  {
  int n=0;
  
  while(strlen(tpc_isotope[n].name)>1) n++;
  if(isotope_code<1 || isotope_code>n-1) return nan("");
  else return(tpc_isotope[isotope_code].hl);
}

Referenced by abssWrite(), imgWriteDICOM(), lambdaFromIsotope(), and tacDecayCorrection().

isotopeIdentify()

int isotopeIdentify

(

const char *

isotope

)

Identify the given string representation of isotope, whether it is in format like 'C-11', '11C', '^11^C', or '^11^Carbon'. Even a one-letter symbol like 'C' may be accepted for the most common PET isotopes.

See also

isotopeHalflife, isotopeName, lambdaFromIsotope, decayCorrectionFactorFromIsotope, elementIdentify

Returns

Isotope code, or 0 (enum ISOTOPE_UNKNOWN) if not identified.

Author

Vesa Oikonen

Parameters

isotope

Name of isotope to identify.

Definition at line 145 of file isotope.c.

  {
  if(isotope==NULL || strnlen(isotope, 5)<1) return ISOTOPE_UNKNOWN; 
 
  unsigned short int i, n=0;
  while(strlen(tpc_isotope[n].name)>1) n++;
 
  /* Check, if isotope can be found directly in the table */
  for(i=0; i<n; i++) {
    if(strcasecmp(tpc_isotope[i].name, isotope)==0) return(i);
  }
 
  /* Ok it was not that easy... try to figure out what it is */
  char *cptr;
  int mass_nr=0, ic_mass_nr=0;
  /* Find the element */
  unsigned short int z=elementIdentify(isotope);
  if(z==0) return(ISOTOPE_UNKNOWN);
  /* Now find the mass number */
  cptr=strpbrk(isotope, "123456789"); // mass cannot start with zero
  if(cptr!=NULL) {
    mass_nr=atoi(cptr);
  }
  //printf("fixed isotope: %s-%d\n", elementSymbol(z), mass_nr);
 
  /* Check if element and mass matches any of the isotopes in our table */
  for(i=1; i<n; i++) {
    /* check the element */
    if(elementIdentify(tpc_isotope[i].name)!=z) continue;
    /* check the mass, if found */
    cptr=strchr(tpc_isotope[i].name, '-'); if(cptr==NULL) continue;
    ic_mass_nr=atoi(cptr+1);
    if(mass_nr>0 && ic_mass_nr!=mass_nr) continue;
    /* Match was found! */
    return(i);
  }
 
  return(ISOTOPE_UNKNOWN);
}

Referenced by abssWrite(), dcmImgIsotope(), tacGetIsotope(), tacRead4DM(), and tacReadSIF().

isotopeIdentifyHalflife()

int isotopeIdentifyHalflife

(

double

halflife

)

Identify the isotope based on halflife.

Returns

Isotope code, or 0 (enum ISOTOPE_UNKNOWN) if not identified.

See also

isotopeName, lambdaFromIsotope, decayCorrectionFactorFromIsotope

Author

Vesa Oikonen

Parameters

halflife

Halflife in minutes

Definition at line 121 of file isotope.c.

  {
  if(halflife<=0.01) return ISOTOPE_UNKNOWN;
  int i, n=0;
  while(strlen(tpc_isotope[n].name)>1) n++;
  /* Check, if halflife can be found in the table */
  for(i=1; i<n; i++) {
    if(fabs(halflife/tpc_isotope[i].hl-1.0)<0.05 ) return i;
  }
  return ISOTOPE_UNKNOWN;
}

Referenced by dcmImgIsotope().

isotopeName()

char * isotopeName

(

int

isotope_code

)

Get the string representation of an isotope specified with its isotope_code.

See also

isotopeBranching, isotopeHalflife, elementName, elementSymbol

Returns

pointer to string representation of the isotope, "unknown" if not identified.

Author

Vesa Oikonen

Parameters

isotope_code

isotope_code as enum or the index of isotope in isotope table.

Definition at line 101 of file isotope.c.

  {
  int n=0;
  
  while(strlen(tpc_isotope[n].name)>1) n++;
  if(isotope_code<1 || isotope_code>n-1)
    return tpc_isotope[ISOTOPE_UNKNOWN].name;
  else
    return(tpc_isotope[isotope_code].name);
}

Referenced by dcmImgIsotope(), imgContents(), imgFillOHeader(), imgReadDICOM(), sifWeight(), tacRead4DM(), tacReadSIF(), tacSetIsotope(), tacSetWeights(), and tacWByFreq().