/home/epics/devel/tomoRecon/tomoReconApp/src/tomoRecon.h

Go to the documentation of this file.

/*
 * tomoRecon.h
 * 
 * C++ class for doing computed tomography reconstruction using Gridrec.
 *
 * This runs the reconstruction using multiple threads, each thread reconstructing a set of slices.
 *
 * It uses the EPICS libCom library for OS-independent functions for threads, mutexes, message queues, etc.
 *
 * Author: Mark Rivers
 *
 * Created: July 9, 2012
 */

#include <epicsMessageQueue.h>
#include <epicsThread.h>

#include "grid.h"

typedef struct {
  int sliceNumber;  
  float center;     
  float *pIn1;      
  float *pIn2;      
  float *pOut1;     
  float *pOut2;     
} toDoMessage_t;

typedef struct {
  int sliceNumber;      
  int numSlices;        
  double sinogramTime;  
  double reconTime;     
} doneMessage_t;

typedef struct {
  int numThreads;           
  int numPixels;            
  int maxSlices;            
  int numProjections;       
  int paddedSinogramWidth;  
  int airPixels;            
  int ringWidth;            
  int fluorescence;         
  int debug;                
  char debugFileName[256];  
  // These are gridRec parameters
  int geom;                 
  float pswfParam;          
  float sampl;              
  float MaxPixSiz;          
  float R;                  
  float X0;                 
  float Y0;                 
  char fname[16];           
  int ltbl;                 
} tomoParams_t;

#ifdef __cplusplus

typedef struct {
  class tomoRecon *pTomoRecon; 
  int taskNum;                 
} workerCreateStruct;

class tomoRecon {
public:
  tomoRecon(tomoParams_t *pTomoParams, float *pAngles);
  virtual ~tomoRecon();
  virtual int reconstruct(int numSlices, float *center, float *pInput, float *pOutput);
  virtual void supervisorTask();
  virtual void workerTask(int taskNum);
  virtual void sinogram(float *pIn, float *pOut);
  virtual void poll(int *pReconComplete, int *pSlicesRemaining);
  virtual void logMsg(const char *pFormat, ...);

private:
  void shutDown();
  tomoParams_t *pTomoParams_;
  int numPixels_;
  int numSlices_;
  int numProjections_;
  int paddedWidth_;
  int numThreads_;
  float *pAngles_;
  float *pInput_;
  float *pOutput_;
  int queueElements_;
  int debug_;
  FILE *debugFile_;
  int reconComplete_;
  int slicesRemaining_;
  int shutDown_;
  epicsMessageQueueId toDoQueue_;
  epicsMessageQueueId doneQueue_;
  epicsEventId supervisorWakeEvent_;
  epicsEventId supervisorDoneEvent_;
  epicsEventId *workerWakeEvents_;
  epicsEventId *workerDoneEvents_;
  epicsMutexId fftwMutex_;
};
#endif