evgen::MarleyGen Class Reference

Inheritance diagram for evgen::MarleyGen:

Classes
struct	Config
	Collection of configuration parameters for the module. More...
struct	KeysToIgnore

Public Types
using	Name = fhicl::Name
using	Comment = fhicl::Comment
using	Parameters = art::EDProducer::Table< Config, KeysToIgnore >

Public Member Functions
	MarleyGen (const Parameters &p)
virtual void	produce (art::Event &e) override
virtual void	beginRun (art::Run &run) override
virtual void	reconfigure (const Parameters &p)

Private Attributes
std::unique_ptr < evgen::MARLEYHelper >	fMarleyHelper
std::unique_ptr < evgen::ActiveVolumeVertexSampler >	fVertexSampler
std::unique_ptr< marley::Event >	fEvent
TTree *	fEventTree
uint_fast32_t	fRunNumber
uint_fast32_t	fSubRunNumber
uint_fast32_t	fEventNumber

Detailed Description

Definition at line 40 of file MARLEYGen_module.cc.

Member Typedef Documentation

using evgen::MarleyGen::Comment = fhicl::Comment

Definition at line 45 of file MARLEYGen_module.cc.

using evgen::MarleyGen::Name = fhicl::Name

Definition at line 44 of file MARLEYGen_module.cc.

using evgen::MarleyGen::Parameters = art::EDProducer::Table<Config, KeysToIgnore>

Definition at line 73 of file MARLEYGen_module.cc.

Constructor & Destructor Documentation

evgen::MarleyGen::MarleyGen ( const Parameters &  p )

explicit

Definition at line 105 of file MARLEYGen_module.cc.

  : EDProducer{ p.get_PSet() },
    fEvent(new marley::Event), fRunNumber(0), fSubRunNumber(0), fEventNumber(0)
{
  // Configure the module (including MARLEY itself) using the FHiCL parameters
  this->reconfigure( p );

  // Create a ROOT TTree using the TFileService that will store the MARLEY
  // event objects (useful for debugging purposes)
  art::ServiceHandle<art::TFileService const> tfs;
  fEventTree = tfs->make<TTree>("MARLEY_event_tree",
    "Neutrino events generated by MARLEY");
  fEventTree->Branch("event", "marley::Event", fEvent.get());

  // Add branches that give the art::Event run, subrun, and event numbers for
  // easy match-ups between the MARLEY and art TTrees. All three are recorded
  // as 32-bit unsigned integers.
  fEventTree->Branch("run_number", &fRunNumber, "run_number/i");
  fEventTree->Branch("subrun_number", &fSubRunNumber, "subrun_number/i");
  fEventTree->Branch("event_number", &fEventNumber, "event_number/i");

  produces< std::vector<simb::MCTruth>   >();
  produces< sumdata::RunData, art::InRun >();
}

Member Function Documentation

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

overridevirtual

Definition at line 131 of file MARLEYGen_module.cc.

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

void evgen::MarleyGen::produce ( art::Event &  e )

overridevirtual

Definition at line 138 of file MARLEYGen_module.cc.

{
  // Get the run, subrun, and event numbers from the current art::Event
  fRunNumber = e.run();
  fSubRunNumber = e.subRun();
  fEventNumber = e.event();

  std::unique_ptr< std::vector<simb::MCTruth> >
    truthcol(new std::vector<simb::MCTruth>);

  // Get the primary vertex location for this event
  art::ServiceHandle<geo::Geometry const> geo;
  TLorentzVector vertex_pos = fVertexSampler->sample_vertex_pos(*geo);

  // Create the MCTruth object, and retrieve the marley::Event object
  // that was generated as it was created
  simb::MCTruth truth = fMarleyHelper->create_MCTruth(vertex_pos,
    fEvent.get());

  // Write the marley::Event object to the event tree
  fEventTree->Fill();

  truthcol->push_back(truth);

  e.put(std::move(truthcol));
}

void evgen::MarleyGen::reconfigure ( const Parameters &  p )

virtual

Definition at line 166 of file MARLEYGen_module.cc.

{
  const auto& seed_service = art::ServiceHandle<rndm::NuRandomService>();
  const auto& geom_service = art::ServiceHandle<geo::Geometry const>();

  // Create a new evgen::ActiveVolumeVertexSampler object based on the current
  // configuration
  fVertexSampler = std::make_unique<evgen::ActiveVolumeVertexSampler>(
    p().vertex_, *seed_service, *geom_service, "MARLEY_Vertex_Sampler");

  // Create a new marley::Generator object based on the current configuration
  fhicl::ParameterSet marley_pset = p.get_PSet().get< fhicl::ParameterSet >(
    "marley_parameters" );
  fMarleyHelper = std::make_unique<MARLEYHelper>( marley_pset,
    *seed_service, "MARLEY" );
}

Member Data Documentation

std::unique_ptr<marley::Event> evgen::MarleyGen::fEvent

private

Definition at line 93 of file MARLEYGen_module.cc.

uint_fast32_t evgen::MarleyGen::fEventNumber

private

Definition at line 101 of file MARLEYGen_module.cc.

TTree* evgen::MarleyGen::fEventTree

private

Definition at line 96 of file MARLEYGen_module.cc.

std::unique_ptr<evgen::MARLEYHelper> evgen::MarleyGen::fMarleyHelper

private

Definition at line 86 of file MARLEYGen_module.cc.

uint_fast32_t evgen::MarleyGen::fRunNumber

private

Definition at line 99 of file MARLEYGen_module.cc.

uint_fast32_t evgen::MarleyGen::fSubRunNumber

private

Definition at line 100 of file MARLEYGen_module.cc.

std::unique_ptr<evgen::ActiveVolumeVertexSampler> evgen::MarleyGen::fVertexSampler

private

Definition at line 90 of file MARLEYGen_module.cc.

---

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

• MARLEYGen_module.cc