gMakeSplinesDM.cxx File Reference

#include <cassert>
#include <cstdlib>
#include <string>
#include <vector>
#include <TSystem.h>
#include "Framework/Algorithm/AlgConfigPool.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/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGLibrary.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 gMakeSplinesDM.cxx:

Go to the source code of this file.

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

Variables
string	gOptTgtPdgCodeList = ""
string	gOptGeomFilename = ""
int	gOptNKnots = -1
double	gOptMaxE = -1.
vector< double >	gOptDMMasses
vector< double >	gOptMedRatios
vector< double >	gOptZpCouplings
bool	gOptNoCopy = false
long int	gOptRanSeed = -1
string	gOptInpXSecFile = ""
string	gOptOutXSecFile = ""

Function Documentation

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

PDGCodeList* GetTargetCodes

(

void

)

int main	(	int	argc,
		char **	argv
	)

Definition at line 153 of file gMakeSplinesDM.cxx.

References genie::PDGLibrary::AddDarkMatter(), genie::RunOpt::BuildTune(), genie::AlgConfigPool::CommonList(), genie::GEVGDriver::Configure(), genie::GEVGDriver::CreateSplines(), GetCommandLineArgs(), GetTargetCodes(), gOptDMMasses, gOptInpXSecFile, gOptMaxE, gOptMedRatios, gOptNKnots, gOptNoCopy, gOptOutXSecFile, gOptRanSeed, gOptZpCouplings, genie::RunOpt::Instance(), genie::PDGLibrary::Instance(), genie::AlgConfigPool::Instance(), genie::XSecSplineList::Instance(), genie::kPdgDarkMatter, genie::Registry::Lock(), LOG, genie::utils::app_init::MesgThresholds(), pFATAL, pINFO, PrintSyntax(), genie::utils::app_init::RandGen(), genie::PDGLibrary::ReloadDBase(), genie::XSecSplineList::SaveAsXml(), genie::Registry::Set(), genie::GEVGDriver::SetEventGeneratorList(), genie::Registry::UnLock(), 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();

  for (vector<double>::iterator mass = gOptDMMasses.begin(); mass != gOptDMMasses.end(); ++mass) {
    for (vector<double>::iterator ratio = gOptMedRatios.begin(); ratio != gOptMedRatios.end(); ++ratio) {
      for (vector<double>::iterator coup = gOptZpCouplings.begin(); coup != gOptZpCouplings.end(); ++coup) {
        // Add dark matter to the table
        PDGLibrary::Instance()->ReloadDBase();
        PDGLibrary::Instance()->AddDarkMatter(*mass,*ratio);
        if (*coup > 0.) {
            Registry * r = AlgConfigPool::Instance()->CommonList("Param", "BoostedDarkMatter");
            r->UnLock();
            r->Set("ZpCoupling", *coup);
            r->Lock();
        }


        // 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(!targets || targets->size() == 0 ) {
          LOG("gmkspl_dm", pFATAL) << "Empty target PDG code list";
          PrintSyntax();
          exit(3);
        }

        LOG("gmkspl_dm", 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 dmpdgc  = kPdgDarkMatter;
          int tgtpdgc = *tgtiter;
          InitialState init_state(tgtpdgc, dmpdgc);
          GEVGDriver driver;
          driver.SetEventGeneratorList(RunOpt::Instance()->EventGeneratorList());
          driver.Configure(init_state);
          driver.CreateSplines(gOptNKnots, gOptMaxE);
        }
        delete targets;
      }
    }
  }
  //  }

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


  return 0;
}

void PrintSyntax

(

void

)

Variable Documentation

vector<double> gOptDMMasses

Definition at line 144 of file gMakeSplinesDM.cxx.

Referenced by main().

string gOptGeomFilename = ""

Definition at line 141 of file gMakeSplinesDM.cxx.

string gOptInpXSecFile = ""

Definition at line 149 of file gMakeSplinesDM.cxx.

double gOptMaxE = -1.

Definition at line 143 of file gMakeSplinesDM.cxx.

vector<double> gOptMedRatios

Definition at line 145 of file gMakeSplinesDM.cxx.

Referenced by main().

int gOptNKnots = -1

Definition at line 142 of file gMakeSplinesDM.cxx.

bool gOptNoCopy = false

Definition at line 147 of file gMakeSplinesDM.cxx.

string gOptOutXSecFile = ""

Definition at line 150 of file gMakeSplinesDM.cxx.

long int gOptRanSeed = -1

Definition at line 148 of file gMakeSplinesDM.cxx.

string gOptTgtPdgCodeList = ""

Definition at line 140 of file gMakeSplinesDM.cxx.

vector<double> gOptZpCouplings

Definition at line 146 of file gMakeSplinesDM.cxx.

Referenced by main().