iScannerSPECTConv.hh

Go to the documentation of this file.

#ifndef ISCANNERSPECTConv_HH
#define ISCANNERSPECTConv_HH 1

#include "gVariables.hh"
#include "vScanner.hh"
#include "sAddonManager.hh"


class iScannerSPECTConv : public vScanner
{
  // -------------------------------------------------------------------
  // Constructor & Destructor
  public:
    iScannerSPECTConv();
    ~iScannerSPECTConv();


  // -------------------------------------------------------------------
  // Public member functions
  public:
    // Function for automatic insertion (put the class name as the parameters and do not add semi-colon at the end of the line)
    FUNCTION_SCANNER(iScannerSPECTConv)
    int Instantiate(bool a_scannerFileIsLUT);
    int CheckParameters();
    int Initialize();
    int BuildLUT( bool a_scannerFileIsLUT );
    int GetPositionsAndOrientations( int a_index1, int a_index2,
                                     FLTNB ap_Position1[3], FLTNB ap_Position2[3],
                                     FLTNB ap_Orientation1[3], FLTNB ap_Orientation2[3],
                                     FLTNB* ap_POI1 = NULL, FLTNB* ap_POI2 = NULL );
    int GetRdmPositionsAndOrientations( int a_index1, int a_index2,
                                        FLTNB ap_Position1[3], FLTNB ap_Position2[3],
                                        FLTNB ap_Orientation1[3], FLTNB ap_Orientation2[3] );
    int GetPositionWithRandomDepth( int a_index1, int a_index2, FLTNB ap_Position1[3], FLTNB ap_Position2[3] );
    int GetTwoCorners( int a_index1, int a_index2,
                       FLTNB ap_CornerInf1[3], FLTNB ap_CornerSup1[3],
                       FLTNB ap_CornerInf2[3], FLTNB ap_CornerSup2[3] );
    int GetGeometricInfoFromDatafile( string a_pathToDF );
    inline int GetSystemNbElts()
           {return m_nbCrystals;}
    void ShowHelp();


  // -------------------------------------------------------------------
  // Functions dedicated to Analytic Projection
  public:
    inline int IsAvailableLOR(int a_elt1, int a_elt2)
           {return 1;};
    int GetSPECTSpecificParameters( uint16_t* ap_nbOfProjections, 
                                    uint16_t* ap_nbHeads, 
                                    uint16_t* ap_nbOfBins,
                                    FLTNB*  ap_pixSizeXY, 
                                    FLTNB*& ap_angles, 
                                    FLTNB*& ap_CORtoDetectorDistance,
                                    int* ap_headRotDirection);
    int PROJ_SetSPECTAngles(FLTNB* ap_projectionAngles);
    int PROJ_SetSPECTCORtoDetectorDistance(FLTNB a_distance);
    inline uint16_t PROJ_GetSPECTNbProjections()
           {return m_nbOfProjections;}
    inline uint16_t PROJ_GetSPECTNbPixels()
           {return m_vNbPixelsTrans*m_vNbPixelsAxial;}
    int PROJ_GetSPECTNbBins(uint16_t* ap_nbOfBins);
    int PROJ_SetSPECTNbBins(uint16_t* ap_nbOfBins);
    int PROJ_SetSPECTNbProjections(uint32_t a_nbOfProjections);

  // -------------------------------------------------------------------
  // Private member functions
  private:
    int LoadLUT();
    int ComputeLUT();
    int ComputeFocalPositions( FLTNB a_posX, FLTNB a_posY, FLTNB a_posZ, int a_headID, int a_cryID );
  
  
  // -------------------------------------------------------------------
  // Data members
  private:
    int m_nbCrystals;                  
    int m_nbHeads;                     
    uint16_t mp_nbOfBins[2];           
    uint16_t m_nbOfProjections;        
    FLTNB* mp_projectionAngles;        
    FLTNB* mp_CORtoDetectorDistance;   
    FLTNB* mp_radius;                  
    uint32_t m_nbPixelsTrans;          
    FLTNB m_pixelsSizeTrans;           
    FLTNB m_gapSizeTrans;              
    uint32_t m_nbPixelsAxial;          
    FLTNB m_pixelsSizeAxial;           
    FLTNB m_gapSizeAxial;              
    uint32_t m_vNbPixelsTrans;         
    uint32_t m_vNbPixelsAxial;         
    FLTNB m_vPixelsSizeTrans;          
    FLTNB m_vPixelsSizeAxial;          
    FLTNB m_crystalDepth;              
    string* mp_focalModelTrans;        
    uint8_t* mp_nbCoefModelTrans;      
    FLTNB** m2p_transFocalParameters;  
    string* mp_focalModelAxial;        
    uint8_t* mp_nbCoefModelAxial;      
    FLTNB** m2p_axialFocalParameters;  
    FLTNB* mp_crystalCentralPositionX; 
    FLTNB* mp_crystalCentralPositionY; 
    FLTNB* mp_crystalCentralPositionZ; 
    FLTNB* mp_crystalOrientationX;     
    FLTNB* mp_crystalOrientationY;     
    FLTNB* mp_crystalOrientationZ;     
    FLTNB* mp_crystalFocalPositionX;   
    FLTNB* mp_crystalFocalPositionY;   
    FLTNB* mp_crystalFocalPositionZ;   
};


// Class for automatic insertion (set here the visible scanner type name, put the class name as the parameters and do not add semi-colon at the end of the line)
CLASS_SCANNER(SPECT_CONVERGENT,iScannerSPECTConv)

#endif