evgen::GENIEGen Class Reference

A module to check the results from the Monte Carlo generator. More...

Inheritance diagram for evgen::GENIEGen:

Public Types
using	ModuleType = EDProducer
template<typename UserConfig , typename KeysToIgnore = void>
using	Table = Modifier::Table< UserConfig, KeysToIgnore >

Public Member Functions
	GENIEGen (fhicl::ParameterSet const &pset)
virtual	~GENIEGen ()
void	produce (art::Event &evt)
void	beginJob ()
void	beginRun (art::Run &run)
void	beginSubRun (art::SubRun &sr)
void	endSubRun (art::SubRun &sr)
void	doBeginJob (SharedResources const &resources)
void	doEndJob ()
void	doRespondToOpenInputFile (FileBlock const &fb)
void	doRespondToCloseInputFile (FileBlock const &fb)
void	doRespondToOpenOutputFiles (FileBlock const &fb)
void	doRespondToCloseOutputFiles (FileBlock const &fb)
bool	doBeginRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doEndRun (RunPrincipal &rp, ModuleContext const &mc)
bool	doBeginSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEndSubRun (SubRunPrincipal &srp, ModuleContext const &mc)
bool	doEvent (EventPrincipal &ep, ModuleContext const &mc, std::atomic< std::size_t > &counts_run, std::atomic< std::size_t > &counts_passed, std::atomic< std::size_t > &counts_failed)
void	fillProductDescriptions ()
void	registerProducts (ProductDescriptions &productsToRegister)
ModuleDescription const &	moduleDescription () const
void	setModuleDescription (ModuleDescription const &)
std::array< std::vector< ProductInfo >, NumBranchTypes > const &	getConsumables () const
void	sortConsumables (std::string const &current_process_name)
std::unique_ptr< Worker >	makeWorker (WorkerParams const &wp)
template<typename T , BranchType BT>
ViewToken< T >	consumesView (InputTag const &tag)
template<typename T , BranchType BT>
ViewToken< T >	mayConsumeView (InputTag const &tag)

Protected Member Functions
ConsumesCollector &	consumesCollector ()
template<typename T , BranchType = InEvent>
ProductToken< T >	consumes (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	consumesView (InputTag const &)
template<typename T , BranchType = InEvent>
void	consumesMany ()
template<typename T , BranchType = InEvent>
ProductToken< T >	mayConsume (InputTag const &)
template<typename Element , BranchType = InEvent>
ViewToken< Element >	mayConsumeView (InputTag const &)
template<typename T , BranchType = InEvent>
void	mayConsumeMany ()

Private Member Functions
std::string	ParticleStatus (int StatusCode)
std::string	ReactionChannel (int ccnc, int mode)
void	FillHistograms (simb::MCTruth mc)

Private Attributes
evgb::GENIEHelper *	fGENIEHelp
	GENIEHelper object. More...
bool	fDefinedVtxHistRange
std::vector< double >	fVtxPosHistRange
	use defined hist range; it is useful to have for asymmetric ranges like in DP FD. More...
int	fPassEmptySpills
	whether or not to kill evnets with no interactions More...
TStopwatch	fStopwatch
double	fGlobalTimeOffset
	keep track of how long it takes to run the job More...
double	fRandomTimeOffset
	The start of a simulated "beam gate". More...
::sim::BeamType_t	fBeamType
	The width of a simulated "beam gate". More...
double	fPrevTotPOT
	The type of beam. More...
double	fPrevTotGoodPOT
	Total good POT from subruns previous to current subrun. More...
TH1F *	fGenerated [6]
	Spectra as generated. More...
TH1F *	fVertexX
	vertex location of generated events in x More...
TH1F *	fVertexY
	vertex location of generated events in y More...
TH1F *	fVertexZ
	vertex location of generated events in z More...
TH2F *	fVertexXY
	vertex location in xy More...
TH2F *	fVertexXZ
	vertex location in xz More...
TH2F *	fVertexYZ
	vertex location in yz More...
TH1F *	fDCosX
	direction cosine in x More...
TH1F *	fDCosY
	direction cosine in y More...
TH1F *	fDCosZ
	direction cosine in z More...
TH1F *	fMuMomentum
	momentum of outgoing muons More...
TH1F *	fMuDCosX
	direction cosine of outgoing mu in x More...
TH1F *	fMuDCosY
	direction cosine of outgoing mu in y More...
TH1F *	fMuDCosZ
	direction cosine of outgoing mu in z More...
TH1F *	fEMomentum
	momentum of outgoing electrons More...
TH1F *	fEDCosX
	direction cosine of outgoing e in x More...
TH1F *	fEDCosY
	direction cosine of outgoing e in y More...
TH1F *	fEDCosZ
	direction cosine of outgoing e in z More...
TH1F *	fCCMode
	CC interaction mode. More...
TH1F *	fNCMode
	CC interaction mode. More...
TH1F *	fDeltaE
	difference in neutrino energy from MCTruth::Enu() vs TParticle More...
TH1F *	fECons
	histogram to determine if energy is conserved in the event More...

Detailed Description

A module to check the results from the Monte Carlo generator.

Note on random number generator

GENIE uses a TRandom generator for its purposes. Since art's RandomNumberGenerator service only provides CLHEP::HepRandomEngine, the standard LArSoft/art mechanism for handling the random stream can't be used. GENIEHelper, interface to GENIE provided by nugen, creates a TRandom that GENIE can use. It initializes it with a random seed read from RandomSeed configuration parameter. This and all the other parameters are inherited from the art module (that is, GENIEGen) configuration. LArSoft meddles with this mechanism to provide support for the standard "Seed" parameter and NuRandomService service.

Configuration parameters

• RandomSeed (integer, optional): if specified, this value is used as seed for GENIE random number generator engine
• Seed (unsigned integer, optional): if specified, this value is used as seed for GENIE random number generator engine; if RandomSeed is also specified, this value is ignored (but in the future this could turn into a configuration error)

As custom, if the random seed is not provided by the configuration, one is fetched from NuRandomService (if available), with the behaviour in lar::util::FetchRandomSeed().

Definition at line 91 of file GENIEGen_module.cc.

Member Typedef Documentation

using art::EDProducer::ModuleType = EDProducer

inherited

Definition at line 17 of file EDProducer.h.

template<typename UserConfig , typename KeysToIgnore = void>

using art::detail::Producer::Table = Modifier::Table<UserConfig, KeysToIgnore>

inherited

Definition at line 26 of file Producer.h.

Constructor & Destructor Documentation

evgen::GENIEGen::GENIEGen ( fhicl::ParameterSet const &  pset )

explicit

Definition at line 158 of file GENIEGen_module.cc.

References fBeamType, fDefinedVtxHistRange, fGENIEHelp, fGlobalTimeOffset, fPassEmptySpills, fRandomTimeOffset, fStopwatch, fVtxPosHistRange, fhicl::ParameterSet::get_if_present(), sim::kBNB, sim::kNuMI, geo::kUnknown, fhicl::ParameterSet::put(), geo::GeometryCore::ROOTFile(), geo::GeometryCore::ROOTGeoManager(), seed, and geo::GeometryCore::TotalMass().

    : EDProducer{pset}
    , fGENIEHelp(0)
    , fDefinedVtxHistRange(pset.get<bool>("DefinedVtxHistRange"))
    , fVtxPosHistRange(pset.get<std::vector<double>>("VtxPosHistRange"))
    , fPassEmptySpills(pset.get<bool>("PassEmptySpills"))
    , fGlobalTimeOffset(pset.get<double>("GlobalTimeOffset", 0))
    , fRandomTimeOffset(pset.get<double>("RandomTimeOffset", 1600.)) // BNB default value
    , fBeamType(::sim::kBNB)
  {
    fStopwatch.Start();

    produces<std::vector<simb::MCTruth>>();
    produces<std::vector<simb::MCFlux>>();
    produces<std::vector<simb::GTruth>>();
    produces<sumdata::RunData, art::InRun>();
    produces<sumdata::POTSummary, art::InSubRun>();
    produces<art::Assns<simb::MCTruth, simb::MCFlux>>();
    produces<art::Assns<simb::MCTruth, simb::GTruth>>();
    produces<std::vector<sim::BeamGateInfo>>();

    // dk2nu additions
    if (pset.get<std::string>("FluxType").find("dk2nu") != std::string::npos) {
      produces<std::vector<bsim::Dk2Nu>>();
      produces<std::vector<bsim::NuChoice>>();
      produces<art::Assns<simb::MCTruth, bsim::Dk2Nu>>();
      produces<art::Assns<simb::MCTruth, bsim::NuChoice>>();
    }

    std::string beam_type_name = pset.get<std::string>("BeamName");

    if (beam_type_name == "numi")

      fBeamType = ::sim::kNuMI;

    else if (beam_type_name == "booster")

      fBeamType = ::sim::kBNB;

    else

      fBeamType = ::sim::kUnknown;

    art::ServiceHandle<geo::Geometry const> geo;

    signed int temp_seed; // the seed read by GENIEHelper is a signed integer...
    fhicl::ParameterSet GENIEconfig(pset);
    if (!GENIEconfig.get_if_present("RandomSeed",
                                    temp_seed)) { // TODO use has_key() when it becomes available
      // no RandomSeed specified; check for the LArSoft-style "Seed" instead:
      // obtain the random seed from a service,
      // unless overridden in configuration with key "Seed"
      unsigned int seed;
      if (!GENIEconfig.get_if_present("Seed", seed))
        seed = art::ServiceHandle<rndm::NuRandomService>()->getSeed();

      // The seed is not passed to RandomNumberGenerator,
      // since GENIE uses a TRandom generator that is owned by the GENIEHelper.
      // Instead, we explicitly configure the random seed for GENIEHelper:
      GENIEconfig.put("RandomSeed", seed);
    } // if no RandomSeed present

    double detectorMass = 0;
    // detectorMass is _only_ needed by GENIEHelper in the case of
    // histogram flux and non-fixed # of event/spill (ie. POTPerSpill non-zero)
    if (pset.get<std::string>("FluxType").find("histogram") == 0 &&
        pset.get<double>("EventsPerSpill") == 0.0 && pset.get<double>("POTPerSpill") > 0.0) {
      TStopwatch timer;
      timer.Start();
      detectorMass = geo->TotalMass(pset.get<std::string>("TopVolume").c_str());
      timer.Stop();
      mf::LogInfo("GENIEProductionTime")
        << "real time to calculate TotalMass: " << timer.RealTime() << " sec";
    }

    fGENIEHelp =
      new evgb::GENIEHelper(GENIEconfig, geo->ROOTGeoManager(), geo->ROOTFile(), detectorMass);
  }

evgen::GENIEGen::~GENIEGen ( )

virtual

Definition at line 238 of file GENIEGen_module.cc.

References fGENIEHelp, and fStopwatch.

  {
    if (fGENIEHelp) delete fGENIEHelp;
    fStopwatch.Stop();
    mf::LogInfo("GENIEProductionTime") << "real time to produce file: " << fStopwatch.RealTime();
  }

Member Function Documentation

void evgen::GENIEGen::beginJob ( )

virtual

Reimplemented from art::EDProducer.

Definition at line 246 of file GENIEGen_module.cc.

References geo::GeometryCore::DetHalfHeight(), geo::GeometryCore::DetHalfWidth(), geo::GeometryCore::DetLength(), fCCMode, fDCosX, fDCosY, fDCosZ, fDefinedVtxHistRange, fDeltaE, fECons, fEDCosX, fEDCosY, fEDCosZ, fEMomentum, fGenerated, fGENIEHelp, fMuDCosX, fMuDCosY, fMuDCosZ, fMuMomentum, fNCMode, fPrevTotGoodPOT, fPrevTotPOT, fVertexX, fVertexXY, fVertexXZ, fVertexY, fVertexYZ, fVertexZ, fVtxPosHistRange, evgb::GENIEHelper::Initialize(), x, y, and z.

  {
    fGENIEHelp->Initialize();

    fPrevTotPOT = 0.;
    fPrevTotGoodPOT = 0.;

    // Get access to the TFile service.
    art::ServiceHandle<art::TFileService const> tfs;

    fGenerated[0] = tfs->make<TH1F>("fGenerated_necc", "", 100, 0.0, 20.0);
    fGenerated[1] = tfs->make<TH1F>("fGenerated_nebcc", "", 100, 0.0, 20.0);
    fGenerated[2] = tfs->make<TH1F>("fGenerated_nmcc", "", 100, 0.0, 20.0);
    fGenerated[3] = tfs->make<TH1F>("fGenerated_nmbcc", "", 100, 0.0, 20.0);
    fGenerated[4] = tfs->make<TH1F>("fGenerated_nnc", "", 100, 0.0, 20.0);
    fGenerated[5] = tfs->make<TH1F>("fGenerated_nbnc", "", 100, 0.0, 20.0);

    fDCosX = tfs->make<TH1F>("fDCosX", ";dx/ds", 200, -1., 1.);
    fDCosY = tfs->make<TH1F>("fDCosY", ";dy/ds", 200, -1., 1.);
    fDCosZ = tfs->make<TH1F>("fDCosZ", ";dz/ds", 200, -1., 1.);

    fMuMomentum = tfs->make<TH1F>("fMuMomentum", ";p_{#mu} (GeV/c)", 500, 0., 50.);
    fMuDCosX = tfs->make<TH1F>("fMuDCosX", ";dx/ds;", 200, -1., 1.);
    fMuDCosY = tfs->make<TH1F>("fMuDCosY", ";dy/ds;", 200, -1., 1.);
    fMuDCosZ = tfs->make<TH1F>("fMuDCosZ", ";dz/ds;", 200, -1., 1.);

    fEMomentum = tfs->make<TH1F>("fEMomentum", ";p_{e} (GeV/c)", 500, 0., 50.);
    fEDCosX = tfs->make<TH1F>("fEDCosX", ";dx/ds;", 200, -1., 1.);
    fEDCosY = tfs->make<TH1F>("fEDCosY", ";dy/ds;", 200, -1., 1.);
    fEDCosZ = tfs->make<TH1F>("fEDCosZ", ";dz/ds;", 200, -1., 1.);

    fCCMode = tfs->make<TH1F>("fCCMode", ";CC Interaction Mode;", 4, 0., 4.);
    fCCMode->GetXaxis()->SetBinLabel(1, "QE");
    fCCMode->GetXaxis()->SetBinLabel(2, "Res");
    fCCMode->GetXaxis()->SetBinLabel(3, "DIS");
    fCCMode->GetXaxis()->SetBinLabel(4, "Coh");
    fCCMode->GetXaxis()->CenterLabels();

    fNCMode = tfs->make<TH1F>("fNCMode", ";NC Interaction Mode;", 4, 0., 4.);
    fNCMode->GetXaxis()->SetBinLabel(1, "QE");
    fNCMode->GetXaxis()->SetBinLabel(2, "Res");
    fNCMode->GetXaxis()->SetBinLabel(3, "DIS");
    fNCMode->GetXaxis()->SetBinLabel(4, "Coh");
    fNCMode->GetXaxis()->CenterLabels();

    fDeltaE = tfs->make<TH1F>("fDeltaE", ";#Delta E_{#nu} (GeV);", 200, -1., 1.);
    fECons = tfs->make<TH1F>("fECons", ";#Delta E(#nu,lepton);", 500, -5., 5.);

    if (fDefinedVtxHistRange == false) {
      art::ServiceHandle<geo::Geometry const> geo;
      double x = 2.1 * geo->DetHalfWidth();
      double y = 2.1 * geo->DetHalfHeight();
      double z = 2. * geo->DetLength();
      int xdiv = TMath::Nint(2 * x / 5.);
      int ydiv = TMath::Nint(2 * y / 5.);
      int zdiv = TMath::Nint(2 * z / 5.);

      fVertexX = tfs->make<TH1F>("fVertexX", ";x (cm)", xdiv, -0.1 * x, x);
      fVertexY = tfs->make<TH1F>("fVertexY", ";y (cm)", ydiv, -y, y);
      fVertexZ = tfs->make<TH1F>("fVertexZ", ";z (cm)", zdiv, -0.1 * z, z);

      fVertexXY = tfs->make<TH2F>("fVertexXY", ";x (cm);y (cm)", xdiv, -0.1 * x, x, ydiv, -y, y);
      fVertexXZ =
        tfs->make<TH2F>("fVertexXZ", ";z (cm);x (cm)", zdiv, -0.2 * z, z, xdiv, -0.1 * x, x);
      fVertexYZ = tfs->make<TH2F>("fVertexYZ", ";z (cm);y (cm)", zdiv, -0.2 * z, z, ydiv, -y, y);
    }
    else {
      int xdiv = TMath::Nint((fVtxPosHistRange[1] - fVtxPosHistRange[0]) / 5.);
      int ydiv = TMath::Nint((fVtxPosHistRange[3] - fVtxPosHistRange[2]) / 5.);
      int zdiv = TMath::Nint((fVtxPosHistRange[5] - fVtxPosHistRange[4]) / 5.);

      fVertexX =
        tfs->make<TH1F>("fVertexX", ";x (cm)", xdiv, fVtxPosHistRange[0], fVtxPosHistRange[1]);
      fVertexY =
        tfs->make<TH1F>("fVertexY", ";y (cm)", ydiv, fVtxPosHistRange[2], fVtxPosHistRange[3]);
      fVertexZ =
        tfs->make<TH1F>("fVertexZ", ";z (cm)", zdiv, fVtxPosHistRange[4], fVtxPosHistRange[5]);

      fVertexXY = tfs->make<TH2F>("fVertexXY",
                                  ";x (cm);y (cm)",
                                  xdiv,
                                  fVtxPosHistRange[0],
                                  fVtxPosHistRange[1],
                                  ydiv,
                                  fVtxPosHistRange[2],
                                  fVtxPosHistRange[3]);
      fVertexXZ = tfs->make<TH2F>("fVertexXZ",
                                  ";z (cm);x (cm)",
                                  zdiv,
                                  fVtxPosHistRange[4],
                                  fVtxPosHistRange[5],
                                  xdiv,
                                  fVtxPosHistRange[0],
                                  fVtxPosHistRange[1]);
      fVertexYZ = tfs->make<TH2F>("fVertexYZ",
                                  ";z (cm);y (cm)",
                                  zdiv,
                                  fVtxPosHistRange[4],
                                  fVtxPosHistRange[5],
                                  ydiv,
                                  fVtxPosHistRange[2],
                                  fVtxPosHistRange[3]);
    }
  }

void evgen::GENIEGen::beginRun ( art::Run &  run )

virtual

Reimplemented from art::EDProducer.

Definition at line 352 of file GENIEGen_module.cc.

References geo::GeometryCore::DetectorName(), art::fullRun(), and art::Run::put().

  {
    art::ServiceHandle<geo::Geometry const> geo;
    run.put(std::make_unique<sumdata::RunData>(geo->DetectorName()), art::fullRun());
  }

void evgen::GENIEGen::beginSubRun ( art::SubRun &  sr )

virtual

Reimplemented from art::EDProducer.

Definition at line 359 of file GENIEGen_module.cc.

References fGENIEHelp, fPrevTotGoodPOT, fPrevTotPOT, and evgb::GENIEHelper::TotalExposure().

  {

    fPrevTotPOT = fGENIEHelp->TotalExposure();
    fPrevTotGoodPOT = fGENIEHelp->TotalExposure();

    return;
  }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::consumes ( InputTag const &  tag )

protectedinherited

Definition at line 61 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumes().

  {
    return collector_.consumes<T, BT>(tag);
  }

ConsumesCollector & art::ModuleBase::consumesCollector ( )

protectedinherited

Definition at line 57 of file ModuleBase.cc.

References art::ModuleBase::collector_.

  {
    return collector_;
  }

template<typename T , BranchType BT>

void art::ModuleBase::consumesMany ( )

protectedinherited

Definition at line 75 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesMany().

  {
    collector_.consumesMany<T, BT>();
  }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::consumesView ( InputTag const &  )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::consumesView ( InputTag const &  tag )

inherited

Definition at line 68 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::consumesView().

  {
    return collector_.consumesView<T, BT>(tag);
  }

void art::detail::Producer::doBeginJob ( SharedResources const &  resources )

inherited

Definition at line 22 of file Producer.cc.

References art::detail::Producer::beginJobWithFrame(), and art::detail::Producer::setupQueues().

  {
    setupQueues(resources);
    ProcessingFrame const frame{ScheduleID{}};
    beginJobWithFrame(frame);
  }

bool art::detail::Producer::doBeginRun ( RunPrincipal &  rp, ModuleContext const &  mc  )

inherited

Definition at line 65 of file Producer.cc.

References art::detail::Producer::beginRunWithFrame(), art::RangeSet::forRun(), art::RunPrincipal::makeRun(), r, art::RunPrincipal::runID(), and art::ModuleContext::scheduleID().

  {
    auto r = rp.makeRun(mc, RangeSet::forRun(rp.runID()));
    ProcessingFrame const frame{mc.scheduleID()};
    beginRunWithFrame(r, frame);
    r.commitProducts();
    return true;
  }

bool art::detail::Producer::doBeginSubRun ( SubRunPrincipal &  srp, ModuleContext const &  mc  )

inherited

Definition at line 85 of file Producer.cc.

References art::detail::Producer::beginSubRunWithFrame(), art::RangeSet::forSubRun(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::SubRunPrincipal::subRunID().

  {
    auto sr = srp.makeSubRun(mc, RangeSet::forSubRun(srp.subRunID()));
    ProcessingFrame const frame{mc.scheduleID()};
    beginSubRunWithFrame(sr, frame);
    sr.commitProducts();
    return true;
  }

void art::detail::Producer::doEndJob ( )

inherited

Definition at line 30 of file Producer.cc.

References art::detail::Producer::endJobWithFrame().

  {
    ProcessingFrame const frame{ScheduleID{}};
    endJobWithFrame(frame);
  }

bool art::detail::Producer::doEndRun ( RunPrincipal &  rp, ModuleContext const &  mc  )

inherited

Definition at line 75 of file Producer.cc.

References art::detail::Producer::endRunWithFrame(), art::RunPrincipal::makeRun(), r, art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

  {
    auto r = rp.makeRun(mc, rp.seenRanges());
    ProcessingFrame const frame{mc.scheduleID()};
    endRunWithFrame(r, frame);
    r.commitProducts();
    return true;
  }

bool art::detail::Producer::doEndSubRun ( SubRunPrincipal &  srp, ModuleContext const &  mc  )

inherited

Definition at line 95 of file Producer.cc.

References art::detail::Producer::endSubRunWithFrame(), art::SubRunPrincipal::makeSubRun(), art::ModuleContext::scheduleID(), and art::Principal::seenRanges().

  {
    auto sr = srp.makeSubRun(mc, srp.seenRanges());
    ProcessingFrame const frame{mc.scheduleID()};
    endSubRunWithFrame(sr, frame);
    sr.commitProducts();
    return true;
  }

bool art::detail::Producer::doEvent ( EventPrincipal &  ep, ModuleContext const &  mc, std::atomic< std::size_t > &  counts_run, std::atomic< std::size_t > &  counts_passed, std::atomic< std::size_t > &  counts_failed  )

inherited

Definition at line 105 of file Producer.cc.

References art::detail::Producer::checkPutProducts_, e, art::EventPrincipal::makeEvent(), art::detail::Producer::produceWithFrame(), and art::ModuleContext::scheduleID().

  {
    auto e = ep.makeEvent(mc);
    ++counts_run;
    ProcessingFrame const frame{mc.scheduleID()};
    produceWithFrame(e, frame);
    e.commitProducts(checkPutProducts_, &expectedProducts<InEvent>());
    ++counts_passed;
    return true;
  }

void art::detail::Producer::doRespondToCloseInputFile ( FileBlock const &  fb )

inherited

Definition at line 44 of file Producer.cc.

References art::detail::Producer::respondToCloseInputFileWithFrame().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToCloseInputFileWithFrame(fb, frame);
  }

void art::detail::Producer::doRespondToCloseOutputFiles ( FileBlock const &  fb )

inherited

Definition at line 58 of file Producer.cc.

References art::detail::Producer::respondToCloseOutputFilesWithFrame().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToCloseOutputFilesWithFrame(fb, frame);
  }

void art::detail::Producer::doRespondToOpenInputFile ( FileBlock const &  fb )

inherited

Definition at line 37 of file Producer.cc.

References art::detail::Producer::respondToOpenInputFileWithFrame().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToOpenInputFileWithFrame(fb, frame);
  }

void art::detail::Producer::doRespondToOpenOutputFiles ( FileBlock const &  fb )

inherited

Definition at line 51 of file Producer.cc.

References art::detail::Producer::respondToOpenOutputFilesWithFrame().

  {
    ProcessingFrame const frame{ScheduleID{}};
    respondToOpenOutputFilesWithFrame(fb, frame);
  }

void evgen::GENIEGen::endSubRun ( art::SubRun &  sr )

virtual

Reimplemented from art::EDProducer.

Definition at line 369 of file GENIEGen_module.cc.

References fGENIEHelp, fPrevTotGoodPOT, fPrevTotPOT, art::SubRun::put(), art::subRunFragment(), and evgb::GENIEHelper::TotalExposure().

  {

    auto p = std::make_unique<sumdata::POTSummary>();

    p->totpot = fGENIEHelp->TotalExposure() - fPrevTotPOT;
    p->totgoodpot = fGENIEHelp->TotalExposure() - fPrevTotGoodPOT;

    sr.put(std::move(p), art::subRunFragment());

    return;
  }

void evgen::GENIEGen::FillHistograms ( simb::MCTruth  mc )

private

< fill the vertex histograms from the neutrino - that is always particle 0 in the list

look for the outgoing lepton in the particle stack just interested in the first one

Definition at line 545 of file GENIEGen_module.cc.

References util::abs(), simb::MCNeutrino::CCNC(), DEFINE_ART_MODULE, simb::MCParticle::E(), energy, fCCMode, fDCosX, fDCosY, fDCosZ, fECons, fEDCosX, fEDCosY, fEDCosZ, fEMomentum, fGenerated, fMuDCosX, fMuDCosY, fMuDCosZ, fMuMomentum, fNCMode, fVertexX, fVertexXY, fVertexXZ, fVertexY, fVertexYZ, fVertexZ, simb::MCTruth::GetNeutrino(), simb::MCTruth::GetParticle(), simb::kCC, art::left(), MF_LOG_DEBUG, simb::MCNeutrino::Mode(), simb::MCTruth::NParticles(), simb::MCNeutrino::Nu(), simb::MCParticle::P(), part, ParticleStatus(), simb::MCParticle::PdgCode(), simb::MCParticle::Px(), simb::MCParticle::Py(), simb::MCParticle::Pz(), ReactionChannel(), simb::MCParticle::Vx(), simb::MCParticle::Vy(), and simb::MCParticle::Vz().

Referenced by produce().

  {
    // Decide which histograms to put the spectrum in
    int id = -1;
    if (mc.GetNeutrino().CCNC() == simb::kCC) {
      fCCMode->Fill(mc.GetNeutrino().Mode());
      if (mc.GetNeutrino().Nu().PdgCode() == 12)
        id = 0;
      else if (mc.GetNeutrino().Nu().PdgCode() == -12)
        id = 1;
      else if (mc.GetNeutrino().Nu().PdgCode() == 14)
        id = 2;
      else if (mc.GetNeutrino().Nu().PdgCode() == -14)
        id = 3;
      else
        return;
    }
    else {
      fNCMode->Fill(mc.GetNeutrino().Mode());
      if (mc.GetNeutrino().Nu().PdgCode() > 0)
        id = 4;
      else
        id = 5;
    }
    if (id == -1) abort();

    // Fill the specta histograms
    fGenerated[id]->Fill(mc.GetNeutrino().Nu().E());

    simb::MCNeutrino mcnu = mc.GetNeutrino();
    const simb::MCParticle nu = mcnu.Nu();

    fVertexX->Fill(nu.Vx());
    fVertexY->Fill(nu.Vy());
    fVertexZ->Fill(nu.Vz());

    fVertexXY->Fill(nu.Vx(), nu.Vy());
    fVertexXZ->Fill(nu.Vz(), nu.Vx());
    fVertexYZ->Fill(nu.Vz(), nu.Vy());

    double mom = nu.P();
    if (std::abs(mom) > 0.) {
      fDCosX->Fill(nu.Px() / mom);
      fDCosY->Fill(nu.Py() / mom);
      fDCosZ->Fill(nu.Pz() / mom);
    }

    MF_LOG_DEBUG("GENIEInteractionInformation")
      << std::endl
      << "REACTION:  " << ReactionChannel(mc.GetNeutrino().CCNC(), mc.GetNeutrino().Mode())
      << std::endl
      << "-----------> Particles in the Stack = " << mc.NParticles() << std::endl
      << std::setiosflags(std::ios::left) << std::setw(20) << "PARTICLE"
      << std::setiosflags(std::ios::left) << std::setw(32) << "STATUS" << std::setw(18) << "E (GeV)"
      << std::setw(18) << "m (GeV/c2)" << std::setw(18) << "Ek (GeV)" << std::endl
      << std::endl;

    const TDatabasePDG* databasePDG = TDatabasePDG::Instance();

    // Loop over the particle stack for this event
    for (int i = 0; i < mc.NParticles(); ++i) {
      simb::MCParticle part(mc.GetParticle(i));
      std::string name = databasePDG->GetParticle(part.PdgCode())->GetName();
      int code = part.StatusCode();
      std::string status = ParticleStatus(code);
      double mass = part.Mass();
      double energy = part.E();
      double Ek = (energy - mass); // Kinetic Energy (GeV)
      if (status == "kIStStableFinalState" || status == "kIStHadronInTheNucleus")
        MF_LOG_DEBUG("GENIEFinalState")
          << std::setiosflags(std::ios::left) << std::setw(20) << name
          << std::setiosflags(std::ios::left) << std::setw(32) << status << std::setw(18) << energy
          << std::setw(18) << mass << std::setw(18) << Ek << std::endl;
      else
        MF_LOG_DEBUG("GENIEFinalState")
          << std::setiosflags(std::ios::left) << std::setw(20) << name
          << std::setiosflags(std::ios::left) << std::setw(32) << status << std::setw(18) << energy
          << std::setw(18) << mass << std::endl;
    }

    if (mc.GetNeutrino().CCNC() == simb::kCC) {

      for (int i = 0; i < mc.NParticles(); ++i) {
        simb::MCParticle part(mc.GetParticle(i));
        if (abs(part.PdgCode()) == 11) {
          fEMomentum->Fill(part.P());
          fEDCosX->Fill(part.Px() / part.P());
          fEDCosY->Fill(part.Py() / part.P());
          fEDCosZ->Fill(part.Pz() / part.P());
          fECons->Fill(nu.E() - part.E());
          break;
        }
        else if (abs(part.PdgCode()) == 13) {
          fMuMomentum->Fill(part.P());
          fMuDCosX->Fill(part.Px() / part.P());
          fMuDCosY->Fill(part.Py() / part.P());
          fMuDCosZ->Fill(part.Pz() / part.P());
          fECons->Fill(nu.E() - part.E());
          break;
        }
      } // end loop over particles
    }   //end if CC interaction

    return;
  }

void art::Modifier::fillProductDescriptions ( )

inherited

Definition at line 10 of file Modifier.cc.

References art::ProductRegistryHelper::fillDescriptions(), and art::ModuleBase::moduleDescription().

  {
    ProductRegistryHelper::fillDescriptions(moduleDescription());
  }

std::array< std::vector< ProductInfo >, NumBranchTypes > const & art::ModuleBase::getConsumables ( ) const

inherited

Definition at line 43 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::getConsumables().

  {
    return collector_.getConsumables();
  }

std::unique_ptr< Worker > art::ModuleBase::makeWorker ( WorkerParams const &  wp )

inherited

Definition at line 37 of file ModuleBase.cc.

References art::ModuleBase::doMakeWorker(), and art::NumBranchTypes.

  {
    return doMakeWorker(wp);
  }

template<typename T , BranchType BT>

ProductToken< T > art::ModuleBase::mayConsume ( InputTag const &  tag )

protectedinherited

Definition at line 82 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsume().

  {
    return collector_.mayConsume<T, BT>(tag);
  }

template<typename T , BranchType BT>

void art::ModuleBase::mayConsumeMany ( )

protectedinherited

Definition at line 96 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeMany().

  {
    collector_.mayConsumeMany<T, BT>();
  }

template<typename Element , BranchType = InEvent>

ViewToken<Element> art::ModuleBase::mayConsumeView ( InputTag const &  )

protectedinherited

template<typename T , BranchType BT>

ViewToken<T> art::ModuleBase::mayConsumeView ( InputTag const &  tag )

inherited

Definition at line 89 of file ModuleBase.h.

References art::ModuleBase::collector_, and art::ConsumesCollector::mayConsumeView().

  {
    return collector_.mayConsumeView<T, BT>(tag);
  }

ModuleDescription const & art::ModuleBase::moduleDescription ( ) const

inherited

Definition at line 13 of file ModuleBase.cc.

References art::errors::LogicError.

Referenced by art::OutputModule::doRespondToOpenInputFile(), art::OutputModule::doWriteEvent(), art::Modifier::fillProductDescriptions(), art::OutputModule::makePlugins_(), art::OutputWorker::OutputWorker(), reco::shower::LArPandoraModularShowerCreation::produce(), art::Modifier::registerProducts(), and art::OutputModule::registerProducts().

  {
    if (md_.has_value()) {
      return *md_;
    }

    throw Exception{errors::LogicError,
                    "There was an error while calling moduleDescription().\n"}
      << "The moduleDescription() base-class member function cannot be called\n"
         "during module construction.  To determine which module is "
         "responsible\n"
         "for calling it, find the '<module type>:<module "
         "label>@Construction'\n"
         "tag in the message prefix above.  Please contact artists@fnal.gov\n"
         "for guidance.\n";
  }

std::string evgen::GENIEGen::ParticleStatus ( int  StatusCode )

private

Definition at line 496 of file GENIEGen_module.cc.

References sim::ParticleStatusName().

Referenced by FillHistograms().

  {
    int code = StatusCode;
    std::string ParticleStatusName;

    switch (code) {
    case -1: ParticleStatusName = "kIStUndefined"; break;
    case 0: ParticleStatusName = "kIStInitialState"; break;
    case 1: ParticleStatusName = "kIStStableFinalState"; break;
    case 2: ParticleStatusName = "kIStIntermediateState"; break;
    case 3: ParticleStatusName = "kIStDecayedState"; break;
    case 11: ParticleStatusName = "kIStNucleonTarget"; break;
    case 12: ParticleStatusName = "kIStDISPreFragmHadronicState"; break;
    case 13: ParticleStatusName = "kIStPreDecayResonantState"; break;
    case 14: ParticleStatusName = "kIStHadronInTheNucleus"; break;
    case 15: ParticleStatusName = "kIStFinalStateNuclearRemnant"; break;
    case 16: ParticleStatusName = "kIStNucleonClusterTarget"; break;
    default: ParticleStatusName = "Status Unknown";
    }
    return ParticleStatusName;
  }

void evgen::GENIEGen::produce ( art::Event &  evt )

virtual

Implements art::EDProducer.

Definition at line 383 of file GENIEGen_module.cc.

References util::CreateAssn(), sim::dump::DumpGTruth(), sim::dump::DumpMCTruth(), fBeamType, fGENIEHelp, fGlobalTimeOffset, FillHistograms(), fPassEmptySpills, fRandomTimeOffset, evgb::GENIEHelper::GetFluxDriver(), simb::MCTruth::GetNeutrino(), simb::MCNeutrino::InteractionType(), simb::kNuanceOffset, MF_LOG_DEBUG, simb::MCTruth::NeutrinoSet(), art::Event::put(), evgb::GENIEHelper::Sample(), and evgb::GENIEHelper::Stop().

  {
    std::unique_ptr<std::vector<simb::MCTruth>> truthcol(new std::vector<simb::MCTruth>);
    std::unique_ptr<std::vector<simb::MCFlux>> fluxcol(new std::vector<simb::MCFlux>);
    std::unique_ptr<std::vector<simb::GTruth>> gtruthcol(new std::vector<simb::GTruth>);
    std::unique_ptr<art::Assns<simb::MCTruth, simb::MCFlux>> tfassn(
      new art::Assns<simb::MCTruth, simb::MCFlux>);
    std::unique_ptr<art::Assns<simb::MCTruth, simb::GTruth>> tgtassn(
      new art::Assns<simb::MCTruth, simb::GTruth>);
    std::unique_ptr<std::vector<sim::BeamGateInfo>> gateCollection(
      new std::vector<sim::BeamGateInfo>);

    std::unique_ptr<std::vector<bsim::Dk2Nu>> dk2nucol(new std::vector<bsim::Dk2Nu>);
    std::unique_ptr<std::vector<bsim::NuChoice>> nuchoicecol(new std::vector<bsim::NuChoice>);
    std::unique_ptr<art::Assns<simb::MCTruth, bsim::Dk2Nu>> dk2nuassn(
      new art::Assns<simb::MCTruth, bsim::Dk2Nu>);
    std::unique_ptr<art::Assns<simb::MCTruth, bsim::NuChoice>> nuchoiceassn(
      new art::Assns<simb::MCTruth, bsim::NuChoice>);

    genie::flux::GDk2NuFlux* dk2nuDriver =
      dynamic_cast<genie::flux::GDk2NuFlux*>(fGENIEHelp->GetFluxDriver(true));
    while (truthcol->size() < 1) {
      while (!fGENIEHelp->Stop()) {

        simb::MCTruth truth;
        simb::MCFlux flux;
        simb::GTruth gTruth;

        // GENIEHelper returns a false in the sample method if
        // either no neutrino was generated, or the interaction
        // occurred beyond the detector's z extent - ie something we
        // would never see anyway.
        if (fGENIEHelp->Sample(truth, flux, gTruth)) {

          truthcol->push_back(truth);
          fluxcol->push_back(flux);
          gtruthcol->push_back(gTruth);
          auto const truthPtr = art::PtrMaker<simb::MCTruth>{evt}(truthcol->size() - 1);
          tfassn->addSingle(truthPtr, art::PtrMaker<simb::MCFlux>{evt}(fluxcol->size() - 1));
          tgtassn->addSingle(truthPtr, art::PtrMaker<simb::GTruth>{evt}(gtruthcol->size() - 1));

          FillHistograms(truth);

          if (dk2nuDriver) {
            const bsim::Dk2Nu& dk2nuObj = dk2nuDriver->GetDk2Nu();
            dk2nucol->push_back(dk2nuObj);
            const bsim::NuChoice& nuchoiceObj = dk2nuDriver->GetNuChoice();
            nuchoicecol->push_back(nuchoiceObj);
            util::CreateAssn(
              evt, *truthcol, *dk2nucol, *dk2nuassn, dk2nucol->size() - 1, dk2nucol->size());
            util::CreateAssn(evt,
                             *truthcol,
                             *nuchoicecol,
                             *nuchoiceassn,
                             nuchoicecol->size() - 1,
                             nuchoicecol->size());
          }

          // check that the process code is not unsupported by GENIE
          // (see issue #18025 for reference);
          // if it is, print all the information we can about this truth record
          if (truth.NeutrinoSet() &&
              (truth.GetNeutrino().InteractionType() == simb::kNuanceOffset)) {
            mf::LogWarning log("GENIEmissingProcessMapping");
            log << "Found an interaction that is not represented by the interaction type code in "
                   "GENIE:"
                   "\nMCTruth record:"
                   "\n";
            sim::dump::DumpMCTruth(log, truth, 2U); // 2 trajectory points per line
            log << "\nGENIE truth record:"
                   "\n";
            sim::dump::DumpGTruth(log, gTruth);
          } // if

        } // end if genie was able to make an event

      } // end event generation loop

      // check to see if we are to pass empty spills
      if (truthcol->size() < 1 && fPassEmptySpills) {
        MF_LOG_DEBUG("GENIEGen") << "no events made for this spill but "
                                 << "passing it on and ending the event anyway";
        break;
      }

    } // end loop while no interactions are made

    // Create a simulated "beam gate" for these neutrino events.
    // We're creating a vector of these because, in a
    // distant-but-possible future, we may be generating more than one
    // beam gate within a simulated time window.
    gateCollection->push_back(sim::BeamGateInfo(fGlobalTimeOffset, fRandomTimeOffset, fBeamType));

    // put the collections in the event
    evt.put(std::move(truthcol));
    evt.put(std::move(fluxcol));
    evt.put(std::move(gtruthcol));
    evt.put(std::move(tfassn));
    evt.put(std::move(tgtassn));
    evt.put(std::move(gateCollection));

    // dk2nu additions
    if (dk2nuDriver) {
      evt.put(std::move(dk2nucol));
      evt.put(std::move(nuchoicecol));
      evt.put(std::move(dk2nuassn));
      evt.put(std::move(nuchoiceassn));
    }

    return;
  }

std::string evgen::GENIEGen::ReactionChannel ( int  ccnc, int  mode  )

private

Definition at line 519 of file GENIEGen_module.cc.

Referenced by FillHistograms().

  {
    std::string ReactionChannelName = " ";

    if (ccnc == 0)
      ReactionChannelName = "kCC";
    else if (ccnc == 1)
      ReactionChannelName = "kNC";
    else
      std::cout << "Current mode unknown!! " << std::endl;

    if (mode == 0)
      ReactionChannelName += "_kQE";
    else if (mode == 1)
      ReactionChannelName += "_kRes";
    else if (mode == 2)
      ReactionChannelName += "_kDIS";
    else if (mode == 3)
      ReactionChannelName += "_kCoh";
    else
      std::cout << "interaction mode unknown!! " << std::endl;

    return ReactionChannelName;
  }

void art::Modifier::registerProducts ( ProductDescriptions &  productsToRegister )

inherited

Definition at line 16 of file Modifier.cc.

References art::ModuleBase::moduleDescription(), and art::ProductRegistryHelper::registerProducts().

  {
    ProductRegistryHelper::registerProducts(productsToRegister,
                                            moduleDescription());
  }

void art::ModuleBase::setModuleDescription ( ModuleDescription const &  md )

inherited

Definition at line 31 of file ModuleBase.cc.

References art::ModuleBase::md_.

  {
    md_ = md;
  }

void art::ModuleBase::sortConsumables ( std::string const &  current_process_name )

inherited

Definition at line 49 of file ModuleBase.cc.

References art::ModuleBase::collector_, and art::ConsumesCollector::sortConsumables().

  {
    // Now that we know we have seen all the consumes declarations,
    // sort the results for fast lookup later.
    collector_.sortConsumables(current_process_name);
  }

Member Data Documentation

::sim::BeamType_t evgen::GENIEGen::fBeamType

private

The width of a simulated "beam gate".

Definition at line 118 of file GENIEGen_module.cc.

Referenced by GENIEGen(), and produce().

TH1F* evgen::GENIEGen::fCCMode

private

CC interaction mode.

Definition at line 147 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fDCosX

private

direction cosine in x

Definition at line 133 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fDCosY

private

direction cosine in y

Definition at line 134 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fDCosZ

private

direction cosine in z

Definition at line 135 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

bool evgen::GENIEGen::fDefinedVtxHistRange

private

Definition at line 110 of file GENIEGen_module.cc.

Referenced by beginJob(), and GENIEGen().

TH1F* evgen::GENIEGen::fDeltaE

private

difference in neutrino energy from MCTruth::Enu() vs TParticle

Definition at line 150 of file GENIEGen_module.cc.

Referenced by beginJob().

TH1F* evgen::GENIEGen::fECons

private

histogram to determine if energy is conserved in the event

Definition at line 151 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fEDCosX

private

direction cosine of outgoing e in x

Definition at line 143 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fEDCosY

private

direction cosine of outgoing e in y

Definition at line 144 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fEDCosZ

private

direction cosine of outgoing e in z

Definition at line 145 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fEMomentum

private

momentum of outgoing electrons

Definition at line 142 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fGenerated[6]

private

Spectra as generated.

Definition at line 123 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

evgb::GENIEHelper* evgen::GENIEGen::fGENIEHelp

private

GENIEHelper object.

Definition at line 108 of file GENIEGen_module.cc.

Referenced by beginJob(), beginSubRun(), endSubRun(), GENIEGen(), produce(), and ~GENIEGen().

double evgen::GENIEGen::fGlobalTimeOffset

private

keep track of how long it takes to run the job

Definition at line 116 of file GENIEGen_module.cc.

Referenced by GENIEGen(), and produce().

TH1F* evgen::GENIEGen::fMuDCosX

private

direction cosine of outgoing mu in x

Definition at line 138 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fMuDCosY

private

direction cosine of outgoing mu in y

Definition at line 139 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fMuDCosZ

private

direction cosine of outgoing mu in z

Definition at line 140 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fMuMomentum

private

momentum of outgoing muons

Definition at line 137 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fNCMode

private

CC interaction mode.

Definition at line 148 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

int evgen::GENIEGen::fPassEmptySpills

private

whether or not to kill evnets with no interactions

Definition at line 113 of file GENIEGen_module.cc.

Referenced by GENIEGen(), and produce().

double evgen::GENIEGen::fPrevTotGoodPOT

private

Total good POT from subruns previous to current subrun.

Definition at line 121 of file GENIEGen_module.cc.

Referenced by beginJob(), beginSubRun(), and endSubRun().

double evgen::GENIEGen::fPrevTotPOT

private

The type of beam.

Total POT from subruns previous to current subrun

Definition at line 120 of file GENIEGen_module.cc.

Referenced by beginJob(), beginSubRun(), and endSubRun().

double evgen::GENIEGen::fRandomTimeOffset

private

The start of a simulated "beam gate".

Definition at line 117 of file GENIEGen_module.cc.

Referenced by GENIEGen(), and produce().

TStopwatch evgen::GENIEGen::fStopwatch

private

Definition at line 114 of file GENIEGen_module.cc.

Referenced by GENIEGen(), and ~GENIEGen().

TH1F* evgen::GENIEGen::fVertexX

private

vertex location of generated events in x

Definition at line 125 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH2F* evgen::GENIEGen::fVertexXY

private

vertex location in xy

Definition at line 129 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH2F* evgen::GENIEGen::fVertexXZ

private

vertex location in xz

Definition at line 130 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fVertexY

private

vertex location of generated events in y

Definition at line 126 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH2F* evgen::GENIEGen::fVertexYZ

private

vertex location in yz

Definition at line 131 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

TH1F* evgen::GENIEGen::fVertexZ

private

vertex location of generated events in z

Definition at line 127 of file GENIEGen_module.cc.

Referenced by beginJob(), and FillHistograms().

std::vector<double> evgen::GENIEGen::fVtxPosHistRange

private

use defined hist range; it is useful to have for asymmetric ranges like in DP FD.

Definition at line 111 of file GENIEGen_module.cc.

Referenced by beginJob(), and GENIEGen().

---

The documentation for this class was generated from the following file:

• larsim/v09_43_00/source/larsim/EventGenerator/GENIE/GENIEGen_module.cc