evwgh::WeightManager Class Reference

#include <WeightManager.h>

Public Member Functions
	WeightManager (const std::string name="WeightManager")
	Default constructor. More...
	~WeightManager ()
	Default destructor. More...
const std::string &	Name () const
	Name getter. More...
template<typename Module >
size_t	Configure (fhicl::ParameterSet const &cfg, Module &module)
	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"

)

Default constructor.

Definition at line 6 of file WeightManager.cxx.

    : _name(name)
  {
    _configured = false;
  }

evwgh::WeightManager::~WeightManager ( )

inline

Default destructor.

Definition at line 36 of file WeightManager.h.

{}

Member Function Documentation

template<typename Module >

size_t evwgh::WeightManager::Configure	(	fhicl::ParameterSet const &	cfg,
		Module &	module
	)

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 85 of file WeightManager.h.

  {

    ::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->createEngine(module, "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 69 of file WeightManager.h.

{ return fWeightCalcMap; }

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

(

void

)

const

Name getter.

Definition at line 12 of file WeightManager.cxx.

  { return _name; }

void evwgh::WeightManager::PrintConfig

(

void

)

Definition at line 51 of file WeightManager.cxx.

                                  {

    return;
  }

void evwgh::WeightManager::Reset ( void  )

inline

Reset.

Definition at line 72 of file WeightManager.h.

    { _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 19 of file WeightManager.cxx.

  {

    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 80 of file WeightManager.h.

std::string evwgh::WeightManager::_name

private

Name.

Definition at line 81 of file WeightManager.h.

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

private

A set of custom weight calculators.

Definition at line 79 of file WeightManager.h.

---

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

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