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...
size_t	Configure (fhicl::ParameterSet const &cfg, art::EngineCreator &ec)
	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 39 of file WeightManager.h.

Constructor & Destructor Documentation

evwgh::WeightManager::WeightManager

(

const std::string

name = "WeightManager"

)

Default constructor.

Definition at line 12 of file WeightManager.cxx.

References _configured.

    : _name(name)
  {
    _configured = false;
  }

evwgh::WeightManager::~WeightManager ( )

inline

Default destructor.

Definition at line 47 of file WeightManager.h.

References Configure(), e, Name(), and Run().

{}

Member Function Documentation

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

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 21 of file WeightManager.cxx.

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

Referenced by evwgh::EventWeight::EventWeight(), and ~WeightManager().

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

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

    // Loop over all the functions and register them
    for (auto ifunc=rw_func.begin(); ifunc!=rw_func.end(); ifunc++) 
    {
      fhicl::ParameterSet const &ps_func = p.get<fhicl::ParameterSet> (*ifunc);
      std::string func_type = ps_func.get<std::string>("type");
  
      WeightCalc* wcalc=WeightCalcFactory::Create(func_type+"WeightCalc");
      if ( wcalc == NULL )
        throw cet::exception(__FUNCTION__) << "Function " << *ifunc << " requested in fcl file has not been registered!" << std::endl;
      if ( fWeightCalcMap.find(*ifunc)!= fWeightCalcMap.end() )
        throw cet::exception(__FUNCTION__) << "Function " << *ifunc << " has been requested multiple times in fcl file!" << std::endl;
  
      //mf::LogInfo("") << "Configuring weight calculator " << *ifunc;
  
      // Create random engine for each rw function (name=*ifunc) (and seed it with random_seed set in the fcl)
      seedservice->createEngine(engine, "HepJamesRandom", *ifunc, ps_func, "random_seed");
  
      wcalc->SetName(*ifunc);
      wcalc->Configure(p);
      Weight_t* winfo=new Weight_t();
      winfo->fWeightCalcType=func_type;
      winfo->fWeightCalc=wcalc;
      winfo->fNmultisims=ps_func.get<int>("number_of_multisims");
  
      std::pair<std::string, Weight_t*> pwc(*ifunc,winfo);
      fWeightCalcMap.insert(pwc);
    }

    _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 79 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 18 of file WeightManager.cxx.

References _name.

Referenced by ~WeightManager().

  { return _name; }

void evwgh::WeightManager::PrintConfig

(

)

Definition at line 99 of file WeightManager.cxx.

Referenced by Reset().

                                  {
    
    return; 
  }

void evwgh::WeightManager::Reset ( )

inline

Reset.

Definition at line 82 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 67 of file WeightManager.cxx.

References _configured, evwgh::MCEventWeight::fWeight, and fWeightCalcMap.

Referenced by evwgh::EventWeight::produce(), and ~WeightManager().

  {     

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

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

std::string evwgh::WeightManager::_name

private

Name.

Definition at line 94 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 90 of file WeightManager.h.

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

---

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

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