gNucleonDecayEvGen.cxx File Reference

#include <cassert>
#include <cstdlib>
#include <string>
#include <vector>
#include <sstream>
#include <TSystem.h>
#include "Framework/Algorithm/AlgFactory.h"
#include "Framework/EventGen/EventRecord.h"
#include "Framework/EventGen/EventGeneratorI.h"
#include "Framework/EventGen/EventRecordVisitorI.h"
#include "Framework/EventGen/GMCJMonitor.h"
#include "Framework/Messenger/Messenger.h"
#include "Framework/Ntuple/NtpWriter.h"
#include "Physics/NucleonDecay/NucleonDecayMode.h"
#include "Physics/NucleonDecay/NucleonDecayUtils.h"
#include "Framework/Numerical/RandomGen.h"
#include "Framework/ParticleData/PDGCodes.h"
#include "Framework/ParticleData/PDGUtils.h"
#include "Framework/ParticleData/PDGLibrary.h"
#include "Framework/Utils/StringUtils.h"
#include "Framework/Utils/UnitUtils.h"
#include "Framework/Utils/PrintUtils.h"
#include "Framework/Utils/AppInit.h"
#include "Framework/Utils/RunOpt.h"
#include "Framework/Utils/CmdLnArgParser.h"

Include dependency graph for gNucleonDecayEvGen.cxx:

Go to the source code of this file.

Functions
void	GetCommandLineArgs (int argc, char **argv)
void	PrintSyntax (void)
int	SelectInitState (void)
const EventRecordVisitorI *	NucleonDecayGenerator (void)
int	main (int argc, char **argv)

Variables
string	kDefOptGeomLUnits = "mm"
string	kDefOptGeomDUnits = "g_cm3"
NtpMCFormat_t	kDefOptNtpFormat = kNFGHEP
string	kDefOptEvFilePrefix = "gntp"
Long_t	gOptRunNu = 1000
int	gOptNev = 10
NucleonDecayMode_t	gOptDecayMode = kNDNull
int	gOptDecayedNucleon = kNDNull
string	gOptEvFilePrefix = kDefOptEvFilePrefix
bool	gOptUsingRootGeom = false
map< int, double >	gOptTgtMix
string	gOptRootGeom
string	gOptRootGeomTopVol = ""
double	gOptGeomLUnits = 0
double	gOptGeomDUnits = 0
long int	gOptRanSeed = -1

Function Documentation

void GetCommandLineArgs	(	int	argc,
		char **	argv
	)

int main	(	int	argc,
		char **	argv
	)

Definition at line 176 of file gNucleonDecayEvGen.cxx.

References genie::NtpWriter::AddEventRecord(), genie::NtpWriter::CustomizeFilenamePrefix(), genie::utils::nucleon_decay::DecayedNucleonPdgCode(), GetCommandLineArgs(), gOptDecayedNucleon, gOptDecayMode, gOptEvFilePrefix, gOptNev, gOptRanSeed, gOptRunNu, genie::NtpWriter::Initialize(), genie::RunOpt::Instance(), kDefOptNtpFormat, LOG, genie::utils::app_init::MesgThresholds(), genie::Interaction::NDecay(), NucleonDecayGenerator(), pINFO, pNOTICE, genie::EventRecordVisitorI::ProcessEventRecord(), genie::utils::app_init::RandGen(), genie::NtpWriter::Save(), SelectInitState(), genie::GHepRecord::SetPrintLevel(), genie::GMCJMonitor::SetRefreshRate(), and genie::GMCJMonitor::Update().

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

  // Init messenger and random number seed
  utils::app_init::MesgThresholds(RunOpt::Instance()->MesgThresholdFiles());
  utils::app_init::RandGen(gOptRanSeed);

  // Initialize an Ntuple Writer to save GHEP records into a TTree
  NtpWriter ntpw(kDefOptNtpFormat, gOptRunNu, gOptRanSeed);
  ntpw.CustomizeFilenamePrefix(gOptEvFilePrefix);
  ntpw.Initialize();

  // Create a MC job monitor for a periodically updated status file
  GMCJMonitor mcjmonitor(gOptRunNu);
  mcjmonitor.SetRefreshRate(RunOpt::Instance()->MCJobStatusRefreshRate());

  // Set GHEP print level
  GHepRecord::SetPrintLevel(RunOpt::Instance()->EventRecordPrintLevel());

  // Get the nucleon decay generator
  const EventRecordVisitorI * mcgen = NucleonDecayGenerator();

  // Event loop
  int ievent = 0;
  int dpdg = 0;
  if (gOptDecayedNucleon > 0) {
    dpdg = gOptDecayedNucleon;
  } else {
    dpdg = utils::nucleon_decay::DecayedNucleonPdgCode(gOptDecayMode);
  }  

  while (1)
  {
     if(ievent == gOptNev) break;

     LOG("gevgen_ndcy", pNOTICE)
          << " *** Generating event............ " << ievent;

     EventRecord * event = new EventRecord;
     int target = SelectInitState();
     int decay  = (int)gOptDecayMode;
     Interaction * interaction = Interaction::NDecay(target,decay,dpdg);
     event->AttachSummary(interaction);

     // Simulate decay     
     mcgen->ProcessEventRecord(event);

     LOG("gevgen_ndcy", pINFO)
         << "Generated event: " << *event;

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

     ievent++;
  } // event loop

  // Save the generated event tree & close the output file
  ntpw.Save();

  LOG("gevgen_ndcy", pNOTICE) << "Done!";

  return 0;
}

const EventRecordVisitorI * NucleonDecayGenerator

(

void

)

Definition at line 298 of file gNucleonDecayEvGen.cxx.

References genie::gAbortingInErr, genie::AlgFactory::GetAlgorithm(), genie::AlgFactory::Instance(), LOG, and pFATAL.

Referenced by main().

{
  string sname   = "genie::EventGenerator";
  string sconfig = "NucleonDecay";
  AlgFactory * algf = AlgFactory::Instance();
  const EventRecordVisitorI * mcgen =
     dynamic_cast<const EventRecordVisitorI *> (algf->GetAlgorithm(sname,sconfig));
  if(!mcgen) {
     LOG("gevgen_ndcy", pFATAL) << "Couldn't instantiate the nucleon decay generator";
     gAbortingInErr = true;
     exit(1);
  }
  return mcgen;
}

void PrintSyntax

(

void

)

int SelectInitState

(

void

)

Variable Documentation

int gOptDecayedNucleon = kNDNull

Definition at line 165 of file gNucleonDecayEvGen.cxx.

Referenced by main().

NucleonDecayMode_t gOptDecayMode = kNDNull

Definition at line 164 of file gNucleonDecayEvGen.cxx.

string gOptEvFilePrefix = kDefOptEvFilePrefix

Definition at line 166 of file gNucleonDecayEvGen.cxx.

double gOptGeomDUnits = 0

Definition at line 172 of file gNucleonDecayEvGen.cxx.

double gOptGeomLUnits = 0

Definition at line 171 of file gNucleonDecayEvGen.cxx.

int gOptNev = 10

Definition at line 163 of file gNucleonDecayEvGen.cxx.

long int gOptRanSeed = -1

Definition at line 173 of file gNucleonDecayEvGen.cxx.

string gOptRootGeom

Definition at line 169 of file gNucleonDecayEvGen.cxx.

string gOptRootGeomTopVol = ""

Definition at line 170 of file gNucleonDecayEvGen.cxx.

Long_t gOptRunNu = 1000

Definition at line 162 of file gNucleonDecayEvGen.cxx.

map<int,double> gOptTgtMix

Definition at line 168 of file gNucleonDecayEvGen.cxx.

bool gOptUsingRootGeom = false

Definition at line 167 of file gNucleonDecayEvGen.cxx.

string kDefOptEvFilePrefix = "gntp"

Definition at line 159 of file gNucleonDecayEvGen.cxx.

string kDefOptGeomDUnits = "g_cm3"

Definition at line 157 of file gNucleonDecayEvGen.cxx.

string kDefOptGeomLUnits = "mm"

Definition at line 156 of file gNucleonDecayEvGen.cxx.

NtpMCFormat_t kDefOptNtpFormat = kNFGHEP

Definition at line 158 of file gNucleonDecayEvGen.cxx.