castor-GATEMacToGeom.cc

Go to the documentation of this file.

#include <iostream>
#include <vector>
#include <string>
#include <cstring>
#include <fstream>
#include <sstream>
#include <math.h>

#include "gDataConversionUtilities.hh"

using namespace std;


// =====================================================================
// ---------------------------------------------------------------------
// ---------------------------------------------------------------------
// =====================================================================
void showHelp(int a_returnCode)
{
  cout << endl;
  cout << " Usage:  castor-GATEMacToGeom -m input_mac_file" << endl;
  cout << "                              -o output_geom_file" << endl;
  cout << "                              -v verbose" << endl;
  cout << endl;
  cout << "[Input settings]:" << endl;
  cout << "  -m input_mac_file     : give the path to a GATE macro file describing a cylindricalPET, ecat or SPECThead system" << endl;
  cout << "  -o output_geom_file   : give the alias of the output geom file which will be generated in the scanner repository (default location: config/scanner/)" << endl;
  //cout << "  -v verbose            : verbosity level" << endl;
  cout << endl;
  #ifdef CASTOR_VERSION
  cout << "  This program is part of the CASToR release version " << CASTOR_VERSION << "." << endl;
  cout << endl;
  #endif
  Exit(a_returnCode);
}




// =====================================================================
// ---------------------------------------------------------------------
// ---------------------------------------------------------------------
// =====================================================================
// Basic checking and processing the input

int main(int argc, char *argv[])
{
  // Path to mac file
  string path_to_mac = "";
  // Path to the output geometric file
  string path_to_geom = "";
  
  if (argc < 3)
  {
    cerr << "Invalid number of arguments. Please read the help section." << endl;
    showHelp(0);
  }
  else
  {
    for (int i=1; i<argc; i++)
    {
      string option = (string)argv[i];

      // .mac file
      if (option=="-m")
      {
        if (argv[i+1] == NULL)
        {
          Cerr("***** castor-GATEMacToGeom :: argument missing for option: " << option << endl);
          Exit(EXIT_FAILURE);
        }
        else
          path_to_mac = argv[i+1];

        // Check if the file exists
        ifstream fcheck(path_to_mac.c_str());
        if(!fcheck.good())
        {
          Cerr("***** castor-GATEMacToGeom :: cannot read provided mac file : " << path_to_mac << endl);
          Exit(EXIT_FAILURE);
        }
        fcheck.close();
        i++;
      }
      
      // output file
      else if (option=="-o")
      {
        if (argv[i+1] == NULL)
        {
          Cerr("***** castor-GATEMacToGeom :: argument missing for option: " << option << endl);
          Exit(EXIT_FAILURE);
        }
        else
          path_to_geom = argv[i+1];
        i++;
      }
      
      else
      {
        Cerr(endl << "***** castor-GATERootToCastor :: Unknown option '" << option << "' !" << endl);
        showHelp(0);
        Exit(EXIT_FAILURE);
      }
    }
  }
  

  // ============================================================================================================
  // SOutputManager object initialisation:
  // ============================================================================================================
    
  sOutputManager* p_outputManager = sOutputManager::GetInstance();  

  // Set path to the config directory
  if (p_outputManager->CheckConfigDir(""))
  {
    Cerr("***** castor-GATEMacToGeom :: A problem occured while checking for the config directory path !" << endl);
    Exit(EXIT_FAILURE);
  }
  
  string scanner_repository = sOutputManager::GetInstance()->GetPathToConfigDir() + "scanner" + OS_SEP;
  path_to_geom = scanner_repository + path_to_geom + ".geom";
  
    cout << "path mac " << path_to_mac << endl;
  // ============================================================================================================
  // Get & check system type 
  // ============================================================================================================
  switch ( GetGATESystemType(path_to_mac) ) 
  {
    case GATE_SYS_CYLINDRICAL:
      cout << endl << " --- CylindricalPET system detected. Proceeding to conversion... --- " << endl << endl;
      
      if(CreateGeomWithCylindrical(path_to_mac , path_to_geom) )
      {
        Cerr("***** castor-GATEMacToGeom :: An error occured while trying to process mac file for cylindrical system: " << path_to_mac << endl);
        Exit(EXIT_FAILURE);
      }
      break;
      
    case GATE_SYS_ECAT:
      cout << endl << " --- ECAT system detected. Proceeding to conversion... --- " << endl;
      if(CreateGeomWithECAT(path_to_mac , path_to_geom) )
      {
        Cerr("***** castor-GATEMacToGeom :: An error occured while trying to process mac file for ecat system: " << path_to_mac  << endl);
        Exit(EXIT_FAILURE);
      }
      break;

    case GATE_SYS_SPECT:
      cout << endl << " --- SPECThead system detected. Proceeding to conversion... --- " << endl;
      if(CreateGeomWithSPECT(path_to_mac , path_to_geom) )
      {
        Cerr("***** castor-GATEMacToGeom :: An error occured while trying to process mac file for SPECT system: " << path_to_mac  << endl);
        Exit(EXIT_FAILURE);
      }
      break;
      
    default: // Unknown system
      Cerr("***** castor-GATEMacToGeom :: System type not supported : " << endl);
      Cerr("                   This script only supports conversion for cylindricalPET ecat and SPECThead systems"  << endl);
      Cerr("                   The system type is recovered from the lines '/gate/systems/...'"  << endl);
      Exit(EXIT_FAILURE);  
      break;
  }
  cout << endl << " --- Conversion completed --- " << endl << endl;
  return 0;
}