evwgh::WeightManager Class Reference

#include "WeightManager.h"

Public Member Functions
	WeightManager (const std::string name="WeightManager")
const std::string &	Name () const
	Name getter. More...
template<typename EngineCreator >
size_t	Configure (fhicl::ParameterSet const &cfg, EngineCreator engineCreator)
	Configuration function. More...
MCEventWeight	Run (art::Event &e, const int inu)
	Core function (previous call to Configure is needed) More...
std::map< std::string, Weight_t * >	GetWeightCalcMap ()
	Returns the map between calculator name and Weight_t product. More...
void	Reset ()
	Reset. More...
void	PrintConfig ()

Private Attributes
std::map< std::string, Weight_t * >	fWeightCalcMap
	A set of custom weight calculators. More...
bool	_configured {false}
	Readiness flag. More...
std::string	_name
	Name. More...

Detailed Description

Definition at line 28 of file WeightManager.h.

Constructor & Destructor Documentation

evwgh::WeightManager::WeightManager

(

const std::string

name = "WeightManager"

)

Definition at line 5 of file WeightManager.cxx.

References _configured.

                                                   : _name(name)
  {
    _configured = false;
  }

Member Function Documentation

template<typename EngineCreator >

size_t evwgh::WeightManager::Configure	(	fhicl::ParameterSet const &	cfg,
		EngineCreator	engineCreator
	)

Configuration function.

Parameters

cfg	the input parameters for settings
the	enging creator for the random seed (usually passed with *this) CONFIGURE FUNCTION: created the weights algorithms in the following way: 0) Looks at the weight_functions fcl parameter to get the name of the calculators 1) Creates the Calculators requested in step 0, and assigne a different random seed to each one 3) The future call WeightManager::Run will run the calculators

Definition at line 77 of file WeightManager.h.

References _configured, evwgh::WeightCalcFactory::Create(), evwgh::Weight_t::fNmultisims, evwgh::Weight_t::fWeightCalc, fWeightCalcMap, evwgh::Weight_t::fWeightCalcType, fhicl::ParameterSet::get(), and rndm::NuRandomService::registerAndSeedEngine().

Referenced by evwgh::EventWeight::EventWeight().

  {

    ::art::ServiceHandle<rndm::NuRandomService> seedservice;

    // Get list of weight functions
    auto const rw_func = p.get<std::vector<std::string>>("weight_functions");

    // Loop over all the functions and register them
    auto const module_label = p.get<std::string>("module_label");
    for (auto const& func : rw_func) {
      auto const ps_func = p.get<fhicl::ParameterSet>(func);
      std::string func_type = ps_func.get<std::string>("type");

      WeightCalc* wcalc = WeightCalcFactory::Create(func_type + "WeightCalc");
      if (wcalc == nullptr)
        throw cet::exception(__FUNCTION__)
          << "Function " << func << " requested in fcl file has not been registered!" << std::endl;
      if (fWeightCalcMap.find(func) != fWeightCalcMap.end())
        throw cet::exception(__FUNCTION__)
          << "Function " << func << " has been requested multiple times in fcl file!" << std::endl;

      // Create random engine for each rw function (name=func) (and seed it with random_seed set in the fcl)
      CLHEP::HepRandomEngine& engine = seedservice->registerAndSeedEngine(
        createEngine("HepJamesRandom", func), "HepJamesRandom", func, ps_func, "random_seed");
      wcalc->SetName(func);
      wcalc->Configure(p, engine);
      Weight_t* winfo = new Weight_t();
      winfo->fWeightCalcType = func_type;
      winfo->fWeightCalc = wcalc;
      winfo->fNmultisims = ps_func.get<int>("number_of_multisims", 0);

      fWeightCalcMap.emplace(func, winfo);
    }

    _configured = true;
    return fWeightCalcMap.size();
  }

std::map<std::string, Weight_t*> evwgh::WeightManager::GetWeightCalcMap ( )

inline

Returns the map between calculator name and Weight_t product.

Definition at line 63 of file WeightManager.h.

References fWeightCalcMap.

Referenced by evwgh::EventWeight::endJob().

{ return fWeightCalcMap; }

const std::string & evwgh::WeightManager::Name

(

)

const

Name getter.

Definition at line 10 of file WeightManager.cxx.

References _name.

  {
    return _name;
  }

void evwgh::WeightManager::PrintConfig

(

)

Definition at line 47 of file WeightManager.cxx.

Referenced by Reset().

  {
    return;
  }

void evwgh::WeightManager::Reset ( )

inline

Reset.

Definition at line 66 of file WeightManager.h.

References _configured, and PrintConfig().

{ _configured = false; }

MCEventWeight evwgh::WeightManager::Run	(	art::Event &	e,
		const int	inu
	)

Core function (previous call to Configure is needed)

Parameters

e

the art event

inu

the index of the simulated neutrino in the event CORE FUNCTION: executes algorithms to assign a weight to the event as requested users. WeightManager::Configure needs to be called first The execution takes following steps: 0) Loos over all the previously emplaced calculators 1) For each of them calculates the weights (more weight can be requested per calculator) 3) Returns a map from "calculator name" to vector of weights calculated which is available inside MCEventWeight

Definition at line 18 of file WeightManager.cxx.

References _configured, util::empty(), evwgh::MCEventWeight::fWeight, and fWeightCalcMap.

Referenced by evwgh::EventWeight::produce().

  {

    if (!_configured)
      throw cet::exception(__PRETTY_FUNCTION__) << "Have not configured yet!" << std::endl;

    //
    // Loop over all functions ang calculate weights
    //
    MCEventWeight mcwgh;
    for (auto it = fWeightCalcMap.begin(); it != fWeightCalcMap.end(); it++) {

      auto const& weights = it->second->GetWeight(e);

      if (weights.size() == 0) {
        std::vector<double> empty;
        std::pair<std::string, std::vector<double>> p("empty", empty);
        mcwgh.fWeight.insert(p);
      }
      else {
        std::pair<std::string, std::vector<double>> p(it->first + "_" + it->second->fWeightCalcType,
                                                      weights[inu]);
        mcwgh.fWeight.insert(p);
      }
    }

    return mcwgh;
  }

Member Data Documentation

bool evwgh::WeightManager::_configured {false}

private

Readiness flag.

Definition at line 72 of file WeightManager.h.

Referenced by Configure(), Reset(), Run(), and WeightManager().

std::string evwgh::WeightManager::_name

private

Name.

Definition at line 73 of file WeightManager.h.

Referenced by Name().

std::map<std::string, Weight_t*> evwgh::WeightManager::fWeightCalcMap

private

A set of custom weight calculators.

Definition at line 71 of file WeightManager.h.

Referenced by Configure(), GetWeightCalcMap(), and Run().

---

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

• larsim/v09_43_00/source/larsim/EventWeight/Base/WeightManager.h
• larsim/v09_43_00/source/larsim/EventWeight/Base/WeightManager.cxx