SNNueAr40CCGen_module.cc

Go to the documentation of this file.

// Framework includes
#include "art/Framework/Core/EDProducer.h"
#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Principal/Event.h"
#include "art/Framework/Principal/Run.h"
#include "art/Framework/Services/Registry/ServiceHandle.h"
#include "fhiclcpp/ParameterSet.h"

// nurandom includes
#include "nurandom/RandomUtils/NuRandomService.h"

// nusimdata includes
#include "nusimdata/SimulationBase/MCTruth.h"

// LArSoft includes
#include "larcore/Geometry/Geometry.h"
#include "larcoreobj/SummaryData/RunData.h"
#include "larsim/EventGenerator/NueAr40CCGenerator.h"

// C++ includes

namespace evgen {

  class SNNueAr40CCGen : public art::EDProducer {
  public:
    explicit SNNueAr40CCGen(fhicl::ParameterSet const& pset);

  private:
    void beginRun(art::Run& run) override;
    void produce(art::Event& event) override;

    CLHEP::HepRandomEngine& fEngine;
    evgen::NueAr40CCGenerator fGenerator;
  };

  //____________________________________________________________________________
  SNNueAr40CCGen::SNNueAr40CCGen(fhicl::ParameterSet const& pset)
    : EDProducer{pset} // Create a default random engine: obtain the random seed
    // freom NuRandomService, unless overriden in configuration with key "Seed"
    , fEngine(art::ServiceHandle<rndm::NuRandomService>()->registerAndSeedEngine(createEngine(0),
                                                                                 pset,
                                                                                 "Seed"))
    , fGenerator{pset.get<fhicl::ParameterSet>("GeneratorAlg")}
  {
    produces<std::vector<simb::MCTruth>>();
    produces<sumdata::RunData, art::InRun>();
  }

  //____________________________________________________________________________
  void SNNueAr40CCGen::beginRun(art::Run& run)
  {
    // Store information about the geometry we are using in run information
    art::ServiceHandle<geo::Geometry const> geo;
    run.put(std::make_unique<sumdata::RunData>(geo->DetectorName()), art::fullRun());
  }

  //____________________________________________________________________________
  void SNNueAr40CCGen::produce(art::Event& event)
  {
    event.put(std::make_unique<std::vector<simb::MCTruth>>(fGenerator.Generate(fEngine)));
  }

}

DEFINE_ART_MODULE(evgen::SNNueAr40CCGen)