simplifyGDML.cc File Reference

Reprocesses a GDML file via GEANT4. More...

#include "Geant4/G4GDMLParser.hh"
#include <getopt.h>
#include <unistd.h>
#include <string>
#include <iostream>
#include <cstdio>
#include <cstdlib>

Go to the source code of this file.

Classes
struct	ConfigurationParameters
struct	ConfigurationParser
struct	ConfigurationParser::opt

Functions
std::string	addNameSuffix (std::string name, std::string suffix)
bool	exists (std::string path)
G4VPhysicalVolume *	readWorldVolume (G4GDMLParser &parser, ConfigurationParameters const &params)
int	writeWorld (G4GDMLParser &parser, ConfigurationParameters const &params)
int	main (int argc, char **argv)

Detailed Description

Reprocesses a GDML file via GEANT4.

Author

Gianluca Petrillo (petri.nosp@m.llo@.nosp@m.fnal..nosp@m.gov)

Date

June 13, 2017

Run with --help argument for usage instructions.

This program needs to be linked to:

• GEANT4
• CLHEP
• XERCES (C)

Example of build command in UPS environment:

g++ -Wall -pedantic -std=c++14 \
  -I"${GEANT4_FQ_DIR}/include" -I"$CLHEP_INC" -I"$XERCES_C_INC" \
  -L"${GEANT4_FQ_DIR}/lib64" -lG4persistency -lG4geometry \
  -L"$XERCES_C_LIB" -lxerces-c \
  -o simplifyGDML.exe simplifyGDML.cc

Note that this program is a glorified 3-line code: parse XML file, get the world volume, write it.

Definition in file simplifyGDML.cc.

Function Documentation

std::string addNameSuffix	(	std::string	name,
		std::string	suffix
	)

Definition at line 283 of file simplifyGDML.cc.

                                                          {
  auto const iExt = name.find(".gdml");
  if (iExt == std::string::npos) name.append(suffix);
  else name.insert(iExt, suffix);
  return name;
} // addNameSuffix()

bool exists

(

std::string

path

)

Definition at line 290 of file simplifyGDML.cc.

                            {
  return access(path.c_str(), F_OK);
} // exists()

int main	(	int	argc,
		char **	argv
	)

Definition at line 345 of file simplifyGDML.cc.

                                {

  //
  // argument parsing
  //
  ConfigurationParameters params;

  auto parseRes = ConfigurationParser::parse(params, argc, argv);

  if (params.destPath.empty() && !params.dontWrite)
    params.destPath = addNameSuffix(params.sourcePath, "-simplified");

  if (params.debugLevel > 0) params.print();

  // print usage instructions before exiting
  if (params.help) {
    ConfigurationParser::printHelp(0, argv[0]);
    return 0;
  }
  if (parseRes != ConfigurationParser::Success) return (int) parseRes;

  //
  // the magic
  //
  G4GDMLParser parser;

  std::cout << std::string(80, '-')
    << "\nReading '" << params.sourcePath << "'..." << std::endl;
  parser.Read(params.sourcePath, params.validate);

  std::cout << std::string(80, '-')
    << "\nFetching the world volume for setup '" << params.setupName << "'..."
    << std::endl;
  // for c2: worldVolume is unused here
  //decltype(auto) worldVolume = parser.GetWorldVolume(params.setupName);
  parser.GetWorldVolume(params.setupName);

  if (params.dontWrite) {
    if (!params.destPath.empty()) {
      std::cerr << "Output path ignored since we don't write anything."
        << std::endl;
    }
  }
  else {
    int res = writeWorld(parser, params);
    if (res != 0) return res;
  }

  //
  // the missing error code
  //
  return 0;
} // main()

G4VPhysicalVolume* readWorldVolume	(	G4GDMLParser &	parser,
		ConfigurationParameters const &	params
	)

Definition at line 296 of file simplifyGDML.cc.

{
  decltype(auto) worldVolume = parser.GetWorldVolume(params.setupName);
  if (!worldVolume) {
    std::cerr << "Failed to find the world volume for setup '"
      << params.setupName << "'" << std::endl;
  }
  return worldVolume;
} // readWorldVolume()

int writeWorld	(	G4GDMLParser &	parser,
		ConfigurationParameters const &	params
	)

Definition at line 306 of file simplifyGDML.cc.

                                                                            {

  std::cout << std::string(80, '-')
    << "\nFetching the world volume for setup '" << params.setupName << "'..."
    << std::endl;
  decltype(auto) worldVolume = readWorldVolume(parser, params);
  if (!worldVolume) {
    std::cerr << "Failed to find the world volume for setup '"
      << params.setupName << "'" << std::endl;
    return 1;
  }

  if (exists(params.destPath)) {
    if (params.overwrite) {
      std::cout << "Overwriting the destination file..." << std::endl;
      std::remove(params.destPath.c_str());
    }
    else {
      std::cerr << "Destination file '" << params.destPath
        << "' already exists." << std::endl;
      return 2;
    }
  }

  std::cout  << std::string(80, '-')
    << "\nWriting '" << params.destPath << "'..." << std::endl;
  parser.Write
    (params.destPath, worldVolume, false /* do not add references to names*/);

  std::cout  << std::string(80, '-')
    << "\n'" << params.sourcePath << "' => [" << params.setupName << "] => '"
    << params.destPath << "': done."
    << std::endl;

  return 0;
} // writeVolume()