gEvGen.cxx File Reference

#include <cstdlib>
#include <cassert>
#include <sstream>
#include <string>
#include <vector>
#include <map>
#include <TFile.h>
#include <TTree.h>
#include <TSystem.h>
#include <TVector3.h>
#include <TH1.h>
#include <TF1.h>
#include "Framework/Conventions/XmlParserStatus.h"
#include "Framework/Conventions/GBuild.h"
#include "Framework/Conventions/Controls.h"
#include "Framework/EventGen/EventRecord.h"
#include "Framework/EventGen/GFluxI.h"
#include "Framework/EventGen/GEVGDriver.h"
#include "Framework/EventGen/GMCJDriver.h"
#include "Framework/EventGen/GMCJMonitor.h"
#include "Framework/Interaction/Interaction.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Ntuple/NtpWriter.h"
#include "Framework/Ntuple/NtpMCFormat.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/Numerical/Spline.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/RunOpt.h"
#include "Framework/Utils/XSecSplineList.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/SystemUtils.h"
#include "Framework/Utils/CmdLnArgParser.h"

Include dependency graph for gEvGen.cxx:

Go to the source code of this file.

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

Variables
int	kDefOptNevents = 0
NtpMCFormat_t	kDefOptNtpFormat = kNFGHEP
Long_t	kDefOptRunNu = 0
int	gOptNevents
double	gOptNuEnergy
double	gOptNuEnergyRange
int	gOptNuPdgCode
map< int, double >	gOptTgtMix
Long_t	gOptRunNu
string	gOptFlux
string	gOptFluxFactors
bool	gOptWeighted
bool	gOptForceInt
bool	gOptUsingFluxOrTgtMix = false
long int	gOptRanSeed
string	gOptInpXSecFile
string	gOptOutFileName
string	gOptStatFileName

Function Documentation

void GenerateEventsAtFixedInitState

(

void

)

Definition at line 293 of file gEvGen.cxx.

References genie::NtpWriter::AddEventRecord(), genie::GEVGDriver::Configure(), genie::GMCJMonitor::CustomizeFilename(), genie::NtpWriter::CustomizeFilename(), genie::GEVGDriver::GenerateEvent(), gOptNevents, gOptNuEnergy, gOptNuPdgCode, gOptOutFileName, gOptRanSeed, gOptRunNu, gOptStatFileName, gOptTgtMix, genie::NtpWriter::Initialize(), genie::RunOpt::Instance(), kDefOptNtpFormat, LOG, pNOTICE, genie::NtpWriter::Save(), genie::GEVGDriver::SetEventGeneratorList(), genie::GMCJMonitor::SetRefreshRate(), genie::GEVGDriver::SetUnphysEventMask(), and genie::GMCJMonitor::Update().

Referenced by main().

{
  int neutrino = gOptNuPdgCode;
  int target   = gOptTgtMix.begin()->first;
  double Ev    = gOptNuEnergy;
  TLorentzVector nu_p4(0.,0.,Ev,Ev); // px,py,pz,E (GeV)

  // Create init state
  InitialState init_state(target, neutrino);

  // Create/config event generation driver
  GEVGDriver evg_driver;
  evg_driver.SetEventGeneratorList(RunOpt::Instance()->EventGeneratorList());
  evg_driver.SetUnphysEventMask(*RunOpt::Instance()->UnphysEventMask());
  evg_driver.Configure(init_state);

  // Initialize an Ntuple Writer
  NtpWriter ntpw(kDefOptNtpFormat, gOptRunNu, gOptRanSeed);

  // If an output file name has been specified... use it
  if (!gOptOutFileName.empty()){
    ntpw.CustomizeFilename(gOptOutFileName);
  }
  ntpw.Initialize();


  // Create an MC Job Monitor
  GMCJMonitor mcjmonitor(gOptRunNu);
  mcjmonitor.SetRefreshRate(RunOpt::Instance()->MCJobStatusRefreshRate());

  // If a status file name has been given... use it
  if (!gOptStatFileName.empty()){
    mcjmonitor.CustomizeFilename(gOptStatFileName);
  }


  LOG("gevgen", pNOTICE)
    << "\n ** Will generate " << gOptNevents << " events for \n"
    << init_state << " at Ev = " << Ev << " GeV";

  // Generate events / print the GHEP record / add it to the ntuple
  int ievent = 0;
  while (ievent < gOptNevents) {
     LOG("gevgen", pNOTICE)
        << " *** Generating event............ " << ievent;

     // generate a single event
     EventRecord * event = evg_driver.GenerateEvent(nu_p4);

     if(!event) {
        LOG("gevgen", pNOTICE)
          << "Last attempt failed. Re-trying....";
        continue;
     }

     LOG("gevgen", pNOTICE)
        << "Generated Event GHEP Record: " << *event;

     // add event at the output ntuple, refresh the mc job monitor & clean up
     ntpw.AddEventRecord(ievent, event);
     mcjmonitor.Update(ievent,event);
     ievent++;
     delete event;
  }

  // Save the generated MC events
  ntpw.Save();
}

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

void Initialize

(

void

)

int main	(	int	argc,
		char **	argv
	)

Definition at line 245 of file gEvGen.cxx.

References GenerateEventsAtFixedInitState(), GetCommandLineArgs(), gOptUsingFluxOrTgtMix, Initialize(), LOG, and pERROR.

{
  GetCommandLineArgs(argc,argv);
  Initialize();

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

  if(gOptUsingFluxOrTgtMix) {
#ifdef __CAN_GENERATE_EVENTS_USING_A_FLUX_OR_TGTMIX__
        GenerateEventsUsingFluxOrTgtMix();
#else
  LOG("gevgen", pERROR)
    << "\n   To be able to generate neutrino events from a flux and/or a target mix"
    << "\n   you need to add the following config options at your GENIE installation:"
    << "\n   --enable-flux-drivers  --enable-geom-drivers \n" ;
#endif
  } else {
     GenerateEventsAtFixedInitState();
  }
  return 0;
}

void PrintSyntax

(

void

)

Variable Documentation

string gOptFlux

Definition at line 234 of file gEvGen.cxx.

Referenced by BuildSpectrum().

string gOptFluxFactors

Definition at line 235 of file gEvGen.cxx.

bool gOptForceInt

Definition at line 237 of file gEvGen.cxx.

string gOptInpXSecFile

Definition at line 240 of file gEvGen.cxx.

int gOptNevents

Definition at line 228 of file gEvGen.cxx.

Referenced by GenerateEventsAtFixedInitState(), and main().

double gOptNuEnergy

Definition at line 229 of file gEvGen.cxx.

Referenced by GenerateEventsAtFixedInitState().

double gOptNuEnergyRange

Definition at line 230 of file gEvGen.cxx.

int gOptNuPdgCode

Definition at line 231 of file gEvGen.cxx.

Referenced by GenerateEventsAtFixedInitState().

string gOptOutFileName

Definition at line 241 of file gEvGen.cxx.

Referenced by ConvertToGHad(), ConvertToGHepMock(), ConvertToGINuke(), ConvertToGRooTracker(), ConvertToGST(), ConvertToGTracker(), ConvertToGXML(), GenerateEventsAtFixedInitState(), and RunCherryPicker().

long int gOptRanSeed

Definition at line 239 of file gEvGen.cxx.

Long_t gOptRunNu

Definition at line 233 of file gEvGen.cxx.

string gOptStatFileName

Definition at line 242 of file gEvGen.cxx.

Referenced by GenerateEventsAtFixedInitState().

map<int,double> gOptTgtMix

Definition at line 232 of file gEvGen.cxx.

bool gOptUsingFluxOrTgtMix = false

Definition at line 238 of file gEvGen.cxx.

Referenced by main().

bool gOptWeighted

Definition at line 236 of file gEvGen.cxx.

int kDefOptNevents = 0

Definition at line 223 of file gEvGen.cxx.

NtpMCFormat_t kDefOptNtpFormat = kNFGHEP

Definition at line 224 of file gEvGen.cxx.

Long_t kDefOptRunNu = 0

Definition at line 225 of file gEvGen.cxx.