
F
C++ OE code for STS
class STS
{
public:
// voxel number (I)
static const int VOX = 3;
// detector elements (K)
static const int DE = 3;
enum PRIORTYPE {FLAT=0, GAMMA, JEFFREYS, ENTROPY};
// These are statistical functions used / code is not provided
ArekStat *stat;
private:
PRIORTYPE prior;
// prior parameters for gamma prior
double g1[VOX],g2[VOX];
// prior parameter for entropy prior
double beta;
// lookups
// Size if equal to number of counts
// detector bin in which count was detected
int *lookupBins;
// voxel in which the origin is located
int *lookupVoxels;
// data
int g[DE];
// Number of events emitted and detected per voxel
int c[VOX];
// Number of event emitted ...