JAANET::JEvtWeightFactorGSeaGen Struct Reference

Implementation of reweighting factor for simulated neutrino interactions according to a specifiable ROOT TFormula. More...

#include <JEvtWeightFactorGSeaGen.hh>

Inheritance diagram for JAANET::JEvtWeightFactorGSeaGen:

Public Types
enum	variables {   COSTH, ENERGY, BJORKEN_X, INELASTICITY,   NUMBER_OF_VARIABLES }
	Indices of reweighting variables for GSeaGen. More...
typedef JEvtWeightFactorHelper < JEvtWeightFactor_t >	JEvtWeightFactorHelper_t
typedef std::shared_ptr < JEvtWeightFactor_t >	pointer_type

Public Member Functions
	JEvtWeightFactorGSeaGen ()
	Default constructor. More...
	JEvtWeightFactorGSeaGen (const JFlux &flux, const char *name, const char *formula)
	Constructor. More...
double	operator() (const Evt &evt) const
	Get weighting factor for given event. More...
void	configure (const JEvtWeightFactor_t &factor)
	Configure event-weight factor. More...
JEvtWeightFactor_t &	getEvtWeightFactor () const
	Get reference to event-weight factor. More...
double	getFactor (const Evt &evt) const
	Get weight-factor of given event. More...

Friends
std::istream &	operator>> (std::istream &in, JEvtWeightFactorGSeaGen &object)
	Stream input. More...
std::ostream &	operator<< (std::ostream &out, const JEvtWeightFactorGSeaGen &object)
	Stream output. More...

Detailed Description

Implementation of reweighting factor for simulated neutrino interactions according to a specifiable ROOT TFormula.

Note: The ROOT TFormula may assume any number of parameters, but should be restricted to
the physical variables listed among JEvtWeightFactorGSeaGen::variables.
These variables may be specified within the ROOT TFormula as 'x[<index>]',
where <index> corresponds to the index of the desired physical variable within JEvtWeightFactorGSeaGen::variables.

Definition at line 38 of file JEvtWeightFactorGSeaGen.hh.

Member Typedef Documentation

template<class JEvtWeightFactor_t>

typedef JEvtWeightFactorHelper<JEvtWeightFactor_t> JAANET::JEvtWeightFactorHelper< JEvtWeightFactor_t >::JEvtWeightFactorHelper_t

inherited

Definition at line 32 of file JEvtWeightFactorHelper.hh.

template<class JEvtWeightFactor_t>

typedef std::shared_ptr<JEvtWeightFactor_t> JAANET::JEvtWeightFactorHelper< JEvtWeightFactor_t >::pointer_type

inherited

Definition at line 33 of file JEvtWeightFactorHelper.hh.

Member Enumeration Documentation

enum JAANET::JEvtWeightFactorGSeaGen::variables

Indices of reweighting variables for GSeaGen.

Enumerator
COSTH	Cosine zenith angle.
ENERGY	Initial state energy.
BJORKEN_X	Björken-x (= fractional momentum carried by the struck nucleon)
INELASTICITY	Inelasticity (= Björken-y)
NUMBER_OF_VARIABLES	Number of reweighting variables; N.B. This enum value needs to be specified last!

Definition at line 45 of file JEvtWeightFactorGSeaGen.hh.

                   {
      COSTH,              //!< Cosine zenith angle
      ENERGY,             //!< Initial state energy
      BJORKEN_X,          //!< Björken-x (= fractional momentum carried by the struck nucleon)
      INELASTICITY,       //!< Inelasticity (= Björken-y)
      NUMBER_OF_VARIABLES //!< Number of reweighting variables; \n
                          //!< N.B.\ This enum value needs to be specified last!
    }; 

Constructor & Destructor Documentation

JAANET::JEvtWeightFactorGSeaGen::JEvtWeightFactorGSeaGen ( )

inline

Default constructor.

Definition at line 58 of file JEvtWeightFactorGSeaGen.hh.

                              :
      JFluxHelper(),
      TFormula()
    {}

JAANET::JEvtWeightFactorGSeaGen::JEvtWeightFactorGSeaGen ( const JFlux &  flux, const char *  name, const char *  formula  )

inline

Constructor.

Parameters

flux	flux function
name	name
formula	formula

Definition at line 71 of file JEvtWeightFactorGSeaGen.hh.

                                                  :
      JFluxHelper(flux),
      TFormula(name, formula)
    {}

Member Function Documentation

double JAANET::JEvtWeightFactorGSeaGen::operator() ( const Evt &  evt ) const

inline

Get weighting factor for given event.

Parameters

evt

event

Returns

weighting factor

Definition at line 85 of file JEvtWeightFactorGSeaGen.hh.

    {
      using namespace std;
      using namespace JPP;

      if (static_cast<const JFluxHelper&>(*this)) {
      
        Double_t vars[NUMBER_OF_VARIABLES] = { 0.0 };
        
        const Trk&   neutrino = get_neutrino(evt);
        const double flux     = getFactor(evt);
      
        vars[COSTH]        = neutrino.dir.z;
        vars[ENERGY]       = getE0(evt);
        vars[BJORKEN_X]    = evt.w2list[W2LIST_GSEAGEN_BX];
        vars[INELASTICITY] = evt.w2list[W2LIST_GSEAGEN_BY];
      
        return this->DoEval(&vars[0]) * flux;
        
      } else {

        THROW(JNullPointerException, "JEvtWeightFactorGSeaGen::operator(): Unspecified flux function.");
      }
    }

template<class JEvtWeightFactor_t>

void JAANET::JEvtWeightFactorHelper< JEvtWeightFactor_t >::configure ( const JEvtWeightFactor_t &  factor )

inlineinherited

Configure event-weight factor.

Parameters

factor

event-weight factor

Definition at line 59 of file JEvtWeightFactorHelper.hh.

    {
      using namespace JPP;
      
      JEvtWeightFactor_t* p = dynamic_cast<JEvtWeightFactor_t*>(factor.clone());

      if (p != NULL) {
        this->reset(p);
      } else {
        THROW(JNullPointerException, "JEvtWeightFactorHelper::configure(): Invalid downcast to class derived from JEvtWeightFactor.");
      }
    }

template<class JEvtWeightFactor_t>

JEvtWeightFactor_t& JAANET::JEvtWeightFactorHelper< JEvtWeightFactor_t >::getEvtWeightFactor ( ) const

inlineinherited

Get reference to event-weight factor.

Returns

reference to event-weight factor

Definition at line 78 of file JEvtWeightFactorHelper.hh.

    {
      using namespace JPP;

      if (static_cast<const JEvtWeightFactorHelper_t&>(*this)) {
        return *(this->get());
      } else {
        THROW(JNullPointerException, "JEvtWeightFactorHelper::getEvtWeightFactor(): Event-weight factor is not set.");
      }
    }

template<class JEvtWeightFactor_t>

double JAANET::JEvtWeightFactorHelper< JEvtWeightFactor_t >::getFactor ( const Evt &  evt ) const

inlineinherited

Get weight-factor of given event.

Parameters

evt

event

Returns

event-weight factor

Definition at line 96 of file JEvtWeightFactorHelper.hh.

    {
      const JEvtWeightFactor_t& weightFactor = getEvtWeightFactor();

      return weightFactor.getFactor(evt);
    }

Friends And Related Function Documentation

std::istream& operator>> ( std::istream &  in, JEvtWeightFactorGSeaGen &  object  )

friend

Stream input.

Parameters

in	input stream
object	gSeaGen event-weight factor

Returns

input stream

Definition at line 118 of file JEvtWeightFactorGSeaGen.hh.

    {
      using namespace std;
      using namespace JPP;

      typedef JToken<';'>   JToken_t;

      JToken_t formula;
      JToken_t equation;
      
      in >> formula;

      object.Clear();
      object.Compile(formula.c_str());
      
      for (int count = 0; count < object.GetNpar() && in >> equation; ++count) {

        const int    index = getParameter(equation);
        const double value = getValue(equation, 0);

        object.SetParameter(index, value);
      }
      
      return in;
    }

std::ostream& operator<< ( std::ostream &  out, const JEvtWeightFactorGSeaGen &  object  )

friend

Stream output.

Parameters

out	output stream
object	gSeaGen event-weight factor

Returns

output stream

Definition at line 153 of file JEvtWeightFactorGSeaGen.hh.

      {
        using namespace std;
        using namespace JPP;

        out << object.GetExpFormula() << ';';
        
        for (int i = 0; i < object.GetNpar(); ++i) {
          out << " p" << i << " = " << SCIENTIFIC(6,3) << object.GetParameter(i) << ';';
        }
        
        return out;
      }

---

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

• JEvtWeightFactorGSeaGen.hh