gMaxPathLengths.cxx File Reference

#include <cassert>
#include <string>
#include <TMath.h>
#include "Framework/EventGen/PathLengthList.h"
#include "Tools/Geometry/ROOTGeomAnalyzer.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/Utils/CmdLnArgParser.h"
#include "Framework/Utils/UnitUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/RunOpt.h"

Go to the source code of this file.

Functions
void	GetCommandLineArgs (int argc, char **argv)
void	PrintSyntax (void)
int	main (int argc, char **argv)

Variables
string	kDefOptXMLFilename = "maxpl.xml"
string	kDefOptGeomLUnits = "mm"
string	kDefOptGeomDUnits = "g_cm3"
string	gOptGeomFilename = ""
string	gOptXMLFilename = ""
string	gOptRootGeomTopVol = ""
double	gOptGeomLUnits = 0
double	gOptGeomDUnits = 0
int	gOptNPoints = -1
int	gOptNRays = -1
long int	gOptRanSeed = -1

Function Documentation

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

Definition at line 150 of file gMaxPathLengths.cxx.

{
  LOG("gmxpl", pINFO) << "Parsing command line arguments";

  // Common run options. 
  RunOpt::Instance()->ReadFromCommandLine(argc,argv);

  // Parse run options for this app

  CmdLnArgParser parser(argc,argv);

  // output XML file name
  if( parser.OptionExists('o') ) {
    LOG("gmxpl", pDEBUG) << "Reading output filename";
    gOptXMLFilename = parser.ArgAsString('o');
  } else {
    LOG("gmxpl", pDEBUG) 
       << "Unspecified output filename - Using default";
    gOptXMLFilename = kDefOptXMLFilename;
  } // -o

  // legth & density units
  string lunits, dunits;
  if( parser.OptionExists('L') ) {
    LOG("gmxpl", pDEBUG) << "Checking for input geometry length units";
    lunits = parser.ArgAsString('L');
  } else {
    LOG("gmxpl", pDEBUG) << "Using default geometry length units";
    lunits = kDefOptGeomLUnits;
  } // -L
  if( parser.OptionExists('D') ) {
    LOG("gmxpl", pDEBUG) << "Checking for input geometry density units";
    dunits = parser.ArgAsString('D');
  } else {
    LOG("gmxpl", pDEBUG) << "Using default geometry density units";
    dunits = kDefOptGeomDUnits;
  } // -D
  gOptGeomLUnits = genie::utils::units::UnitFromString(lunits);
  gOptGeomDUnits = genie::utils::units::UnitFromString(dunits);

  // root geometry top volume name
  if( parser.OptionExists('t') ) {
    LOG("gmxpl", pDEBUG) 
       << "Reading root geometry top volume name";
    gOptRootGeomTopVol = parser.ArgAsString('t');
  } else {
    LOG("gmxpl", pDEBUG) 
       << "Unspecified geometry top volume - Using default";
    gOptRootGeomTopVol = "";
  } // -o
  
  // number of scanning points / surface
  if( parser.OptionExists('n') ) {
    LOG("gmxpl", pDEBUG) 
       << "Reading input number of scanning points/surface";
    gOptNPoints = parser.ArgAsInt('n');
  } else {
    LOG("gmxpl", pDEBUG)
      << "Unspecified number of points - Using driver's default";
  } //-n

  // number of scanning rays / point
  if( parser.OptionExists('r') ) {
    LOG("gmxpl", pDEBUG) 
       << "Reading input number of scanning rays/point";
    gOptNRays = parser.ArgAsInt('r');
  } else {
    LOG("gmxpl", pDEBUG)
      << "Unspecified number of rays - Using driver's default";
  } //-r

  // input geometry file
  if( parser.OptionExists('f') ) {
    LOG("gmxpl", pDEBUG) 
       << "Reading ROOT/GEANT geometry filename";
    gOptGeomFilename = parser.ArgAsString('f');
  } else {
    LOG("gmxpl", pFATAL) 
       << "No geometry file was specified - Exiting";
    PrintSyntax();
    exit(1);
  } //-f

  // random number seed
  if( parser.OptionExists("seed") ) {
    LOG("gmxpl", pINFO) << "Reading random number seed";
    gOptRanSeed = parser.ArgAsLong("seed");
  } else {
    LOG("gmxpl", pINFO) << "Unspecified random number seed - Using default";
    gOptRanSeed = -1;
  }

  // print the command line arguments
  LOG("gmxpl", pNOTICE)
     << "\n"
     << utils::print::PrintFramedMesg("gmxpl job inputs");
  LOG("gmxpl", pNOTICE) << "Command line arguments";
  LOG("gmxpl", pNOTICE) << "Input ROOT geometry     : " << gOptGeomFilename;
  LOG("gmxpl", pNOTICE) << "Output XML file         : " << gOptXMLFilename;
  LOG("gmxpl", pNOTICE) << "Geometry length units   : " << gOptGeomLUnits;
  LOG("gmxpl", pNOTICE) << "Geometry density units  : " << gOptGeomDUnits;
  LOG("gmxpl", pNOTICE) << "Scanner points/surface  : " << gOptNPoints;
  LOG("gmxpl", pNOTICE) << "Scanner rays/point      : " << gOptNRays;
  LOG("gmxpl", pNOTICE) << "Random number seed      : " << gOptRanSeed;

  LOG("gmxpl", pNOTICE) << "\n";
  LOG("gmxpl", pNOTICE) << *RunOpt::Instance();
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 112 of file gMaxPathLengths.cxx.

{
  GetCommandLineArgs(argc,argv);

  utils::app_init::RandGen(gOptRanSeed);
  utils::app_init::MesgThresholds(RunOpt::Instance()->MesgThresholdFiles());

  // Create the geometry driver
  LOG("gmxpl", pINFO)
     << "Creating/configuring a ROOT geom. driver";

  ROOTGeomAnalyzer * geom = new ROOTGeomAnalyzer(gOptGeomFilename);
  geom -> SetLengthUnits       (gOptGeomLUnits);
  geom -> SetDensityUnits      (gOptGeomDUnits);
  geom -> SetWeightWithDensity (true);

  // Set the top volume name
  geom -> SetTopVolName        (gOptRootGeomTopVol);
  geom -> SetWeightWithDensity (true);

  if(gOptNPoints > 0) geom->SetScannerNPoints(gOptNPoints);
  if(gOptNRays   > 0) geom->SetScannerNRays  (gOptNRays);

  // Compute the maximum path lengths
  LOG("gmxpl", pINFO)
      << "Asking input GeomAnalyzerI for the max path-lengths";
  const PathLengthList & plmax = geom->ComputeMaxPathLengths();

  // Print & save the maximum path lengths in XML format
  LOG("gmxpl", pINFO)
      << "Maximum path lengths: " << plmax;
  plmax.SaveAsXml(gOptXMLFilename);

  delete geom;

  return 0;
}

void PrintSyntax

(

void

)

Definition at line 259 of file gMaxPathLengths.cxx.

{
  LOG("gmxpl", pNOTICE)
      << "\n\n" << "Syntax:" << "\n"
      << "   gmxpl"
      << " -f geom_file"
      << " [-L length_units]"
      << " [-D density_units]" 
      << " [-t top_volume_name]"
      << " [-o output_xml_file]"
      << " [-seed random_number_seed]"
      << " [--message-thresholds xml_file]\n";

}

Variable Documentation

double gOptGeomDUnits = 0

Definition at line 106 of file gMaxPathLengths.cxx.

string gOptGeomFilename = ""

Definition at line 102 of file gMaxPathLengths.cxx.

double gOptGeomLUnits = 0

Definition at line 105 of file gMaxPathLengths.cxx.

int gOptNPoints = -1

Definition at line 107 of file gMaxPathLengths.cxx.

int gOptNRays = -1

Definition at line 108 of file gMaxPathLengths.cxx.

long int gOptRanSeed = -1

Definition at line 109 of file gMaxPathLengths.cxx.

string gOptRootGeomTopVol = ""

Definition at line 104 of file gMaxPathLengths.cxx.

string gOptXMLFilename = ""

Definition at line 103 of file gMaxPathLengths.cxx.

string kDefOptGeomDUnits = "g_cm3"

Definition at line 99 of file gMaxPathLengths.cxx.

string kDefOptGeomLUnits = "mm"

Definition at line 98 of file gMaxPathLengths.cxx.

string kDefOptXMLFilename = "maxpl.xml"

Definition at line 97 of file gMaxPathLengths.cxx.