sim::DumpMCTruth Class Reference

Inheritance diagram for sim::DumpMCTruth:

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

Public Types
using	Parameters = art::EDAnalyzer::Table< Config >
	Type to enable module parameters description by art. More...

Public Member Functions
	DumpMCTruth (Parameters const &config)
	Configuration-checking constructor. More...
	DumpMCTruth (DumpMCTruth const &)=delete
	DumpMCTruth (DumpMCTruth &&)=delete
DumpMCTruth &	operator= (DumpMCTruth const &)=delete
DumpMCTruth &	operator= (DumpMCTruth &&)=delete
void	analyze (art::Event const &event) override

Static Public Member Functions
template<typename Handle >
static std::string	productName (Handle const &handle)
	Returns the name of the product in the form "module_instance_process". More...

Private Attributes
std::vector< art::InputTag >	fInputTruth
	Name of MCTruth data products. More...
std::string	fOutputCategory
	Name of the stream for output. More...
bool	bAllTruth = false
	Whether to process all MCTruth collections. More...
unsigned int	fPointsPerLine
	trajectory points per output line More...

Detailed Description

Definition at line 37 of file DumpMCTruth_module.cc.

Member Typedef Documentation

using sim::DumpMCTruth::Parameters = art::EDAnalyzer::Table<Config>

Type to enable module parameters description by art.

Definition at line 65 of file DumpMCTruth_module.cc.

Constructor & Destructor Documentation

sim::DumpMCTruth::DumpMCTruth ( Parameters const &  config )

explicit

Configuration-checking constructor.

Definition at line 102 of file DumpMCTruth_module.cc.

  : EDAnalyzer(config)
  , fInputTruth()
  , fOutputCategory(config().OutputCategory())
  , bAllTruth(!config().InputTruth(fInputTruth)) // true if InputTruth omitted
  , fPointsPerLine(config().PointsPerLine())
  {}

sim::DumpMCTruth::DumpMCTruth ( DumpMCTruth const &  )

delete

sim::DumpMCTruth::DumpMCTruth ( DumpMCTruth &&  )

delete

Member Function Documentation

void sim::DumpMCTruth::analyze ( art::Event const &  event )

override

Definition at line 112 of file DumpMCTruth_module.cc.

                                                  {

  //
  // prepare the data products to be dumped
  //
  struct ProductInfo_t {
    using Thruths_t = std::vector<simb::MCTruth>;
    Thruths_t const* truths;
    std::string name;

    ProductInfo_t(art::Handle<Thruths_t> const& handle)
      : truths(handle.provenance()->isPresent()? handle.product(): nullptr)
      , name(sim::DumpMCTruth::productName(handle))
      {}
    ProductInfo_t(art::ValidHandle<Thruths_t> const& handle)
      : truths(handle.product()), name(sim::DumpMCTruth::productName(handle))
      {}

  }; // ProductInfo_t

  std::vector<ProductInfo_t> AllTruths;
  if (bAllTruth) {
    //std::vector<art::Handle<std::vector<simb::MCTruth>>> handles;
    //event.getManyByType(handles);
    auto handles = event.getMany<std::vector<simb::MCTruth>>();
    std::copy(handles.begin(), handles.end(), std::back_inserter(AllTruths));
  }
  else {
    for (auto const& inputTag: fInputTruth) {
      AllTruths.emplace_back
        (event.getValidHandle<std::vector<simb::MCTruth>>(inputTag));
    } // for
  }

  //
  // sanity check
  //
  if (AllTruths.empty()) {
    throw art::Exception(art::errors::ProductNotFound)
      << "No MC truth found to be dumped!\n";
  }

  //
  // print an introduction
  //
  unsigned int const nTruths = std::accumulate(
    AllTruths.begin(), AllTruths.end(), 0U,
    [](unsigned int total, auto const& info)
      { return total + (info.truths? info.truths->size(): 0); }
    );

  if (bAllTruth) {
    mf::LogVerbatim(fOutputCategory) << "Event " << event.id()
      << " contains " << nTruths << " MC truth blocks in "
      << AllTruths.size() << " collections";
  }
  else if (AllTruths.size() == 1) {
    mf::LogVerbatim(fOutputCategory) << "Event " << event.id();
  }
  else {
    mf::LogVerbatim(fOutputCategory) << "Dumping " << nTruths
      << " MC truth blocks from " << AllTruths.size()
      << " collections in event " << event.id();
  }

  //
  // dump data product by data product
  //
  unsigned int nParticles = 0, nNeutrinos = 0;
  for (ProductInfo_t const& truths_info: AllTruths) {

    auto const* truths = truths_info.truths;
    std::string productName = truths_info.name;

    if (!truths) {
      mf::LogVerbatim(fOutputCategory)
        << "Data product '" << productName
        << "' has been dropped. No information available.";
    }

    if (AllTruths.size() > 1) {
      mf::LogVerbatim(fOutputCategory)
        << "Data product '" << productName
        << "' contains " << truths->size() << " truth blocks:";
    }
    else if (truths->size() > 1) {
      mf::LogVerbatim(fOutputCategory)
        << truths->size() << " truth blocks:";
    }

    //
    // dump each MC truth in the data product
    //
    unsigned int iTruth = 0;
    for (auto const& truth: *truths) {

      mf::LogVerbatim log (fOutputCategory);

      if (truths->size() > 1) log << "(#" << iTruth << ") ";
      sim::dump::DumpMCTruth(log, truth, "  ", "");

      //
      // update counters
      //
      ++iTruth;
      nParticles += truth.NParticles();
      if (truth.NeutrinoSet()) ++nNeutrinos;

    } // for each truth in data product

  } // for truth data products

  //
  // all done
  //
  mf::LogVerbatim(fOutputCategory) << nNeutrinos
    << " neutrinos generated, " << nParticles
    << " generated particles to be simulated downstream.";

} // sim::DumpMCTruth::analyze()

DumpMCTruth& sim::DumpMCTruth::operator= ( DumpMCTruth const &  )

delete

DumpMCTruth& sim::DumpMCTruth::operator= ( DumpMCTruth &&  )

delete

template<typename Handle >

std::string sim::DumpMCTruth::productName ( Handle const &  handle )

static

Returns the name of the product in the form "module_instance_process".

Definition at line 236 of file DumpMCTruth_module.cc.

                                                          {
  auto const* prov = handle.provenance();
  return prov->moduleLabel()
    + '_' + prov->productInstanceName()
    + '_' + prov->processName()
    ;
} // sim::DumpMCTruth::productName()

Member Data Documentation

bool sim::DumpMCTruth::bAllTruth = false

private

Whether to process all MCTruth collections.

Definition at line 90 of file DumpMCTruth_module.cc.

std::vector<art::InputTag> sim::DumpMCTruth::fInputTruth

private

Name of MCTruth data products.

Definition at line 88 of file DumpMCTruth_module.cc.

std::string sim::DumpMCTruth::fOutputCategory

private

Name of the stream for output.

Definition at line 89 of file DumpMCTruth_module.cc.

unsigned int sim::DumpMCTruth::fPointsPerLine

private

trajectory points per output line

Definition at line 92 of file DumpMCTruth_module.cc.

---

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

• DumpMCTruth_module.cc