iScannerPET.hh

Go to the documentation of this file.

/*
This file is part of CASToR.

    CASToR is free software: you can redistribute it and/or modify it under the
    terms of the GNU General Public License as published by the Free Software
    Foundation, either version 3 of the License, or (at your option) any later
    version.

    CASToR is distributed in the hope that it will be useful, but WITHOUT ANY
    WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
    FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
    details.

    You should have received a copy of the GNU General Public License along with
    CASToR (in file GNU_GPL.TXT). If not, see <http://www.gnu.org/licenses/>.

Copyright 2017-2019 all CASToR contributors listed below:

    --> Didier BENOIT, Claude COMTAT, Marina FILIPOVIC, Thibaut MERLIN, Mael MILLARDET, Simon STUTE, Valentin VIELZEUF

This is CASToR version 3.0.
*/

#ifndef ISCANNERPET_HH
#define ISCANNERPET_HH 1

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

class iScannerPET : public vScanner
{
  // Constructor & Destructor
  public:
    iScannerPET();

    ~iScannerPET();


  // -------------------------------------------------------------------
  // 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(iScannerPET)
    void DescribeSpecific();
    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 GetEdgesCenterPositions( int a_index1, int a_index2,
      FLTNB ap_pos_line_point1[3], FLTNB ap_pos_line_point2[3],
      FLTNB ap_pos_point1_x[4], FLTNB ap_pos_point1_y[4], FLTNB ap_pos_point1_z[4],
      FLTNB ap_pos_point2_x[4], FLTNB ap_pos_point2_y[4], FLTNB ap_pos_point2_z[4] );
    inline int GetSystemNbElts()
           {return m_nbCrystals;};

    int IsAvailableLOR(int a_elt1, int a_elt2);
    int GetGeometricInfoFromDataFile(string a_pathToDataFilename);
    inline int SetPETMaxAxialDiffmm(FLTNB a_maxAxialDiffmm)
           {m_maxAxialDiffmm = a_maxAxialDiffmm; return 0;}
    void ShowHelp();
    inline FLTNB GetDetectionElementSizeTrans()
           {return mp_sizeCrystalTrans[0];}
    inline FLTNB GetDetectionElementSizeAxial()
           {return mp_sizeCrystalAxial[0];}

  // -------------------------------------------------------------------
  // Functions dedicated to Analytic Projection
    int PROJ_GetPETSpecificParameters(FLTNB* ap_maxAxialDiffmm);


  // -------------------------------------------------------------------
  // Private member functions
  private:
    int LoadLUT();
    int ComputeLUT();
    int GetLayer(int a_idx);


  // -------------------------------------------------------------------
  // Data members
  private:
    int m_nbLayers;                    
    int m_nbCrystals;                  
    int* mp_nbCrystalsInLayer;         
    FLTNB* mp_crystalCentralPositionX; 
    FLTNB* mp_crystalCentralPositionY; 
    FLTNB* mp_crystalCentralPositionZ; 
    FLTNB* mp_crystalOrientationX;     
    FLTNB* mp_crystalOrientationY;     
    FLTNB* mp_crystalOrientationZ;     
    FLTNB* mp_sizeCrystalTrans;        
    FLTNB* mp_sizeCrystalAxial;        
    FLTNB* mp_sizeCrystalDepth;        
    FLTNB* mp_meanDepthOfInteraction;  
    FLTNB m_minAngleDifference;        
    FLTNB m_maxAxialDiffmm;            
};


// 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(PET,iScannerPET)

#endif