Blame |
Last modification |
View Log
| RSS feed
//
// ********************************************************************
// * License and Disclaimer *
// * *
// * The GAMOS software is copyright of the Copyright Holders of *
// * the GAMOS Collaboration. It is provided under the terms and *
// * conditions of the GAMOS Software License, included in the file *
// * LICENSE and available at http://fismed.ciemat.es/GAMOS/license .*
// * These include a list of copyright holders. *
// * *
// * Neither the authors of this software system, nor their employing *
// * institutes,nor the agencies providing financial support for this *
// * work make any representation or warranty, express or implied, *
// * regarding this software system or assume any liability for its *
// * use. Please see the license in the file LICENSE and URL above *
// * for the full disclaimer and the limitation of liability. *
// * *
// * This code implementation is the result of the scientific and *
// * technical work of the GAMOS collaboration. *
// * By using, copying, modifying or distributing the software (or *
// * any work based on the software) you agree to acknowledge its *
// * use in resulting scientific publications, and indicate your *
// * acceptance of all terms of the GAMOS Software license. *
// ********************************************************************
//
#ifndef PETProjDataMgr_h
#define PETProjDataMgr_h
/*---------------------------------------------------------------------------
ClassName: PETProjDataMgr
Author: M. Canadas, P. Arce
Changes: 01/11: creation
-------------------------------------------------------------------------
// Description
PET output (List-mode X Y Z) into projection data (sinograms).
Output data: Sinograms for PET image reconstruction. Interfile format, STIR compatible (.hs/.s),
STIR Software for Tomographic Image Reconstruction: http://stir.sourceforge.net/main.htm
-----------------------------------------------------------------------*/
#include <stdlib.h>
#include <stdio.h>
#include <stdint.h>
#include "H1D.h"
#include "H2D.h"
#include "H3D.h"
#include "daq.h"
typedef unsigned short SINO_TYPE; //!!NOTE: Check "Write_sino3D" (.hv file) if SINO_TYPE changes !!
struct PETProjDataMgr {
int m_Debug;
double m_AxialDistance;
double m_RingDiameter;
int m_NOfPlanes;
int m_NOfBins;
int m_NOfAngles;
int m_MaxRingDifference;
int m_OutFormat;
int m_nch;
int m_TotalAxialPlanes;
char m_Filename;
//G4bool bDumpCout;
SINO_TYPE ***m_projections;
SINO_TYPE *m_Buffer;
unsigned long int m_TotalCoincidences;
unsigned long int m_TotalProjectionCoincidences;
};
//-----------------------------------------------------------------------
int PETProjDataMgrInit();
void PETProjDataMgrFree();
void AddEvent(const HVector3 pos1, const HVector3 pos2);
//void ReadFile();
//PETOutput ReadEvent( G4bool& bEof );
void SetProjection(int axialplane, int h2d);
void WriteInterfile(const char *fname);
void SetNPlanes(int n);
void SetNBin(int n);
void SetNAng(int n);
void SetDebug(int n);
void SetNumberOfChannels(int n);
void SetRingDiameter( double x);
void SetAxialDistance( double x);
int Phantom(int kaj);
int FwdProject2d(int img2d, int h2d);
int FwdProject3d(int img3d, int h3d);
int FwdProject(double x,double y, double z, int nmax, int h1d);
HVector3 Hits2Digits(const HVector3 r);
struct PETProjDataMgr* GetInstance();
#endif