gMakeSplines.cxx File Reference

#include <cassert>
#include <cstdlib>
#include <string>
#include <vector>
#include <TSystem.h>
#include "Framework/Conventions/GBuild.h"
#include "Framework/EventGen/GEVGDriver.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/ParticleData/PDGCodeList.h"
#include "Framework/Utils/RunOpt.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/XSecSplineList.h"
#include "Framework/Utils/CmdLnArgParser.h"

Include dependency graph for gMakeSplines.cxx:

Go to the source code of this file.

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

Variables
string	gOptNuPdgCodeList = ""
string	gOptTgtPdgCodeList = ""
string	gOptGeomFilename = ""
int	gOptNKnots = -1
double	gOptMaxE = -1.
bool	gOptNoCopy = false
long int	gOptRanSeed = -1
string	gOptInpXSecFile = ""
string	gOptOutXSecFile = ""

Function Documentation

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

PDGCodeList * GetNeutrinoCodes

(

void

)

Definition at line 366 of file gMakeSplines.cxx.

References gOptNuPdgCodeList, genie::PDGCodeList::push_back(), and genie::utils::str::Split().

Referenced by main().

{
  // split the comma separated list
  vector<string> nuvec = utils::str::Split(gOptNuPdgCodeList,  ",");

  // fill in the PDG code list
  PDGCodeList * list = new PDGCodeList;
  vector<string>::const_iterator iter;
  for(iter = nuvec.begin(); iter != nuvec.end(); ++iter) {
    list->push_back( atoi(iter->c_str()) );
  }
  return list;
}

PDGCodeList * GetTargetCodes

(

void

)

Definition at line 380 of file gMakeSplines.cxx.

References genie::gAbortingInErr, geom, gOptGeomFilename, gOptTgtPdgCodeList, genie::geometry::ROOTGeomAnalyzer::ListOfTargetNuclei(), LOG, pFATAL, pINFO, genie::PDGCodeList::push_back(), and genie::utils::str::Split().

Referenced by main().

{
  bool from_geom_file = ( gOptGeomFilename.size()   > 0 );
  bool from_cmd_line  = ( gOptTgtPdgCodeList.size() > 0 );

  if (from_cmd_line) {
     // split the comma separated list
     vector<string> tgtvec = utils::str::Split(gOptTgtPdgCodeList, ",");

     // fill in the PDG code list
     PDGCodeList * list = new PDGCodeList;
     vector<string>::const_iterator iter;
     for(iter = tgtvec.begin(); iter != tgtvec.end(); ++iter) {
        list->push_back( atoi(iter->c_str()) );
     }
     return list;
  }

  if (from_geom_file) {
#ifdef __GENIE_GEOM_DRIVERS_ENABLED__
     // create/configure a geometry driver
     LOG("gmkspl", pINFO) << "Creating/configuring a ROOT geom. driver";
     ROOTGeomAnalyzer * geom = new ROOTGeomAnalyzer(gOptGeomFilename);

     PDGCodeList * list = new PDGCodeList(geom->ListOfTargetNuclei());

     delete geom;
     return list;
#else
     LOG("gmkspl", pFATAL)
      << "To read-in a ROOT geometry you need to enable the geometry drivers!";
     gAbortingInErr = true;
     exit(1);
     return 0;
#endif

  }
  return 0;
}

int main	(	int	argc,
		char **	argv
	)

Definition at line 145 of file gMakeSplines.cxx.

References genie::RunOpt::BuildTune(), genie::GEVGDriver::Configure(), genie::GEVGDriver::CreateSplines(), GetCommandLineArgs(), GetNeutrinoCodes(), GetTargetCodes(), gOptInpXSecFile, gOptMaxE, gOptNKnots, gOptNoCopy, gOptOutXSecFile, gOptRanSeed, genie::RunOpt::Instance(), genie::XSecSplineList::Instance(), LOG, genie::utils::app_init::MesgThresholds(), pFATAL, pINFO, PrintSyntax(), genie::utils::app_init::RandGen(), genie::XSecSplineList::SaveAsXml(), genie::GEVGDriver::SetEventGeneratorList(), and genie::utils::app_init::XSecTable().

{
  // Parse command line arguments
  GetCommandLineArgs(argc,argv);

  if ( ! RunOpt::Instance()->Tune() ) {
    LOG("gmkspl", pFATAL) << " No TuneId in RunOption";
    exit(-1);
  }
  RunOpt::Instance()->BuildTune();

  // throw on NaNs and Infs...
#if defined(HAVE_FENV_H) && defined(HAVE_FEENABLEEXCEPT)
  feenableexcept(FE_DIVBYZERO | FE_INVALID | FE_OVERFLOW);
#endif

  // Init
  utils::app_init::MesgThresholds(RunOpt::Instance()->MesgThresholdFiles());
  utils::app_init::RandGen(gOptRanSeed);
  utils::app_init::XSecTable(gOptInpXSecFile, false);

  // Get list of neutrinos and nuclear targets

  PDGCodeList * neutrinos = GetNeutrinoCodes();
  PDGCodeList * targets   = GetTargetCodes();

  if(!neutrinos || neutrinos->size() == 0 ) {
     LOG("gmkspl", pFATAL) << "Empty neutrino PDG code list";
     PrintSyntax();
     exit(2);
  }
  if(!targets || targets->size() == 0 ) {
     LOG("gmkspl", pFATAL) << "Empty target PDG code list";
     PrintSyntax();
     exit(3);
  }

  LOG("gmkspl", pINFO) << "Neutrinos: " << *neutrinos;
  LOG("gmkspl", pINFO) << "Targets: "   << *targets;

  // Loop over all possible input init states and ask the GEVGDriver
  // to build splines for all the interactions that its loaded list
  // of event generators can generate.

  PDGCodeList::const_iterator nuiter;
  PDGCodeList::const_iterator tgtiter;
  for(nuiter = neutrinos->begin(); nuiter != neutrinos->end(); ++nuiter) {
    for(tgtiter = targets->begin(); tgtiter != targets->end(); ++tgtiter) {
      int nupdgc  = *nuiter;
      int tgtpdgc = *tgtiter;
      InitialState init_state(tgtpdgc, nupdgc);
      GEVGDriver driver;
      driver.SetEventGeneratorList(RunOpt::Instance()->EventGeneratorList());
      driver.Configure(init_state);
      driver.CreateSplines(gOptNKnots, gOptMaxE);
    }
  }

  // Save the splines at the requested XML file
  XSecSplineList * xspl = XSecSplineList::Instance();
  bool save_init = !gOptNoCopy;
  xspl->SaveAsXml(gOptOutXSecFile, save_init);

  delete neutrinos;
  delete targets;

  return 0;
}

void PrintSyntax

(

void

)

Variable Documentation

string gOptGeomFilename = ""

Definition at line 136 of file gMakeSplines.cxx.

Referenced by GetTargetCodes(), and main().

string gOptInpXSecFile = ""

Definition at line 141 of file gMakeSplines.cxx.

double gOptMaxE = -1.

Definition at line 138 of file gMakeSplines.cxx.

Referenced by main().

int gOptNKnots = -1

Definition at line 137 of file gMakeSplines.cxx.

Referenced by main().

bool gOptNoCopy = false

Definition at line 139 of file gMakeSplines.cxx.

Referenced by main().

string gOptNuPdgCodeList = ""

Definition at line 134 of file gMakeSplines.cxx.

Referenced by GetNeutrinoCodes().

string gOptOutXSecFile = ""

Definition at line 142 of file gMakeSplines.cxx.

Referenced by main().

long int gOptRanSeed = -1

Definition at line 140 of file gMakeSplines.cxx.

string gOptTgtPdgCodeList = ""

Definition at line 135 of file gMakeSplines.cxx.

Referenced by GetTargetCodes().