evwgh::EventWeight Class Reference

Inheritance diagram for evwgh::EventWeight:

Public Member Functions
	EventWeight (fhicl::ParameterSet const &p)
	EventWeight (EventWeight const &)=delete
	EventWeight (EventWeight &&)=delete
EventWeight &	operator= (EventWeight const &)=delete
EventWeight &	operator= (EventWeight &&)=delete

Private Member Functions
void	produce (art::Event &e) override
void	endJob () override

Private Attributes
WeightManager	_wgt_manager
std::string	fGenieModuleLabel

Detailed Description

Definition at line 32 of file EventWeight_module.cc.

Constructor & Destructor Documentation

evwgh::EventWeight::EventWeight ( fhicl::ParameterSet const &  p )

explicit

Definition at line 58 of file EventWeight_module.cc.

    : EDProducer{p}
    , fGenieModuleLabel{p.get<std::string>("genie_module_label", "generator")}
  {
    // Configure the appropriate GENIE tune if needed (important for v3+ only)
    // NOTE: In all normal use cases, relying on the ${GENIE_XSEC_TUNE}
    // environment variable set by the genie_xsec package should be sufficient.
    // Only include the "TuneName" FHiCL parameter for EventWeight if you
    // really know what you're doing! The same goes for the
    // "EventGeneratorList" parameter.
    std::string genie_tune_name = p.get<std::string>("TuneName",
                                                     "${GENIE_XSEC_TUNE}");

    // The default empty string used here will cause the subsequent call to
    // evgb::SetEventGeneratorListAndTune() to leave GENIE's current event
    // generator list name (probably "Default") alone
    std::string evgen_list_name = p.get<std::string>("EventGeneratorList", "");

    // Tell GENIE about the event generator list and tune
    evgb::SetEventGeneratorListAndTune( evgen_list_name, genie_tune_name );

    auto const n_func = _wgt_manager.Configure(p, *this);
    if ( n_func > 0 )
      produces<std::vector<MCEventWeight> >();
  }

evwgh::EventWeight::EventWeight ( EventWeight const &  )

delete

evwgh::EventWeight::EventWeight ( EventWeight &&  )

delete

Member Function Documentation

void evwgh::EventWeight::endJob ( )

overrideprivate

Definition at line 106 of file EventWeight_module.cc.

  {
    // Get the map from sting to Weight_t from the manager
    std::map<std::string, Weight_t*> weightCalcMap = _wgt_manager.GetWeightCalcMap();

    std::stringstream job_summary;
    job_summary  <<  std::setprecision(2);
    for (int i=1; i <= 110 ;i++) job_summary << "=";
    job_summary << std::endl;
    job_summary << std::setw(20) << "WeightCalc"
                << std::setw(15) << "Type"
                << std::setw(15) << "#RW neutrinos"
                << std::setw(15) << "#Multisims"
                << std::setw(15) << "Min"
                << std::setw(15) << "Max"
                << std::setw(15) << "Avg"
                << std::endl;
    for (int i=1; i <= 110; i++) job_summary << "=";
    job_summary << std::endl;
    for (auto it = weightCalcMap.begin(); it!=weightCalcMap.end(); it++) {
      job_summary << std::setw(20) << it->first
                  << std::setw(15) << (it->second->fWeightCalcType)
                  << std::setw(15) << (it->second->fNcalls)
                  << std::setw(15) << (it->second->fNmultisims)
                  << std::setw(15) << (it->second->fMinWeight)
                  << std::setw(15) << (it->second->fMaxWeight)
                  << std::setw(15) << (it->second->fAvgWeight)
                  << std::endl;
    }
    for (int i=1; i<=110; i++) job_summary << "=";
    job_summary << std::endl;
    mf::LogInfo("") << job_summary.str();
  }

EventWeight& evwgh::EventWeight::operator= ( EventWeight const &  )

delete

EventWeight& evwgh::EventWeight::operator= ( EventWeight &&  )

delete

void evwgh::EventWeight::produce ( art::Event &  e )

overrideprivate

Definition at line 84 of file EventWeight_module.cc.

  {
    // Implementation of required member function here.
    auto mcwghvec = std::make_unique<std::vector<MCEventWeight>>();

    // Get the MC generator information out of the event
    // these are all handles to mc information.
    std::vector<art::Ptr<simb::MCTruth> > mclist;

      // Actually go and get the stuff
    auto const mcTruthHandle = e.getValidHandle<std::vector<simb::MCTruth>>(fGenieModuleLabel);
      art::fill_ptr_vector(mclist, mcTruthHandle);

    // Loop over all neutrinos in this event
    for (unsigned int inu = 0; inu < mclist.size(); ++inu) {
      auto const mcwgh = _wgt_manager.Run(e, inu);
      mcwghvec->push_back(mcwgh);
    }

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

Member Data Documentation

WeightManager evwgh::EventWeight::_wgt_manager

private

Definition at line 54 of file EventWeight_module.cc.

std::string evwgh::EventWeight::fGenieModuleLabel

private

Definition at line 55 of file EventWeight_module.cc.

---

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

• EventWeight_module.cc