JEvtWeightFactorDictionary.hh

Go to the documentation of this file.

#ifndef __JAANET__JEVTWEIGHTFACTORDICTIONARY__
#define __JAANET__JEVTWEIGHTFACTORDICTIONARY__
 
#include <map>
 
#include "JAAnet/JEvtWeightFactorConstant.hh"
#include "JAAnet/JEvtWeightFactorMupage.hh"
#include "JAAnet/JEvtWeightFactorGSeaGen.hh"
#include "JAAnet/JEvtWeightFactorTriggerEfficiencyRatio.hh"
 
#include "JAAnet/JEvtWeightFactorHelper.hh"
 
#include "JAAnet/JFluxDictionary.hh"
 
#include "JOscProb/JOscProbHelper.hh"
 
 
/**
 * \author bjung
 */
namespace JAANET {}
namespace JPP { using namespace JAANET; }
 
namespace JAANET {
 
  using JOSCPROB::JOscProbHelper;
  
  
  /**
   * Dictionary to map distinct event-weight factors to unique identifiers.
   * The identifiers are provided by the enum `factors`.
   */
  class JEvtWeightFactorDictionary:
    public std::map<int, JEvtWeightFactorHelper>
  {
  public:
 
    /**
     * Indices of event-weight factors in event-weight factor dictionary.
     */
    enum factors {
      CUSTOM_FACTORS_BEGIN = 100, //!< Start of custom event-weight factors
        
      CONSTANT             = 101, //!< Constant event-weight factor
      MUPAGE_TFORMULA      = 102, //!< MUPAGE  TFormula-based event-weight factor
      GSEAGEN_TFORMULA     = 103, //!< gSeaGen TFormula-based event-weight factor
      MEFF_RATIO           = 104, //!< Effective mass ratio
 
      CUSTOM_FACTORS_END   = 199  //!< End of custom event-weight factors
    };
    
 
    /**
     * Default constructor.
     *
     * \param  oscProb           oscillation probability interface   
     */
    JEvtWeightFactorDictionary(const JOscProbHelper& oscProb)
    {
      using namespace std;
      
      const JFluxDictionary fluxdict(oscProb);
 
      for (JFluxDictionary::const_iterator i = fluxdict.cbegin(); i != fluxdict.cend(); ++i) {
        this->insert(*i);
      }
 
      this->insert(make_pair(CONSTANT,         make_shared<JEvtWeightFactorConstant<JEvtWeightFactor> >()));
      this->insert(make_pair(MUPAGE_TFORMULA,  make_shared<JEvtWeightFactorMupage>()));
      this->insert(make_pair(GSEAGEN_TFORMULA, make_shared<JEvtWeightFactorGSeaGen>()));
      this->insert(make_pair(MEFF_RATIO,       make_shared<JEvtWeightFactorTriggerEfficiencyRatio>()));
    }
 
 
    /**
     * Default constructor.
     */
    JEvtWeightFactorDictionary() :
      JEvtWeightFactorDictionary(JOscProbHelper())
    {}
 
 
    /**
     * Look-up function.
     *
     * \param  identifier        flux category ID
     */
    const JEvtWeightFactor& at(const int identifier) const
    {
      using namespace std;
      using namespace JPP;
 
      const_iterator i = this->find(identifier);
 
      if (i != this->cend()) {
        return *(i->second);
      } else {
        THROW(JValueOutOfRange, "JEvtWeightFactorDictionary::at(): Invalid identifier " << identifier);
      }
    }
 
 
    /**
     * Look-up operator.
     *
     * \param  identifier        flux category ID
     */
    const JEvtWeightFactor& operator[](const int identifier) const
    {
      return this->at(identifier);
    }
  };
}
  
#endif