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,   INTERACTION_TYPE, CURRENT_TYPE, XSEC, XSEC_MEAN,   XSEC_DIFFERENTIAL, XSEC_WATER, INT_LENGTH_WATER, COLUMN_DEPTH,   P_EARTH, P_SCALE, TARGET_A, TARGET_Z,   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)
INTERACTION_TYPE	Interaction channel type.
CURRENT_TYPE	Weak current type (CC or NC)
XSEC	Exclusive total cross-section of the interaction.
XSEC_MEAN	Average interaction cross-section per nucleon along neutrino path [m2].
XSEC_DIFFERENTIAL	Differential cross-section of the interaction.
XSEC_WATER	Inclusive cross-section in water.
INT_LENGTH_WATER	Interaction length in pure water.
COLUMN_DEPTH	Column density [m.w.e].
P_EARTH	Earth transmission probability.
P_SCALE	GENIE ineraction probability scale.
TARGET_A	Number of nucleons in the target.
TARGET_Z	Number of protons in the target.
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)
      INTERACTION_TYPE,    //!< Interaction channel type
      CURRENT_TYPE,        //!< Weak current type (CC or NC)
      XSEC,                //!< Exclusive total cross-section of the interaction
      XSEC_MEAN,           //!< Average interaction cross-section per nucleon along neutrino path [m2]
      XSEC_DIFFERENTIAL,   //!< Differential cross-section of the interaction
      XSEC_WATER,          //!< Inclusive cross-section in water
      INT_LENGTH_WATER,    //!< Interaction length in pure water
      COLUMN_DEPTH,        //!< Column density [m.w.e]
      P_EARTH,             //!< Earth transmission probability
      P_SCALE,             //!< GENIE ineraction probability scale
      TARGET_A,            //!< Number of nucleons in the target
      TARGET_Z,            //!< Number of protons  in the target
      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 70 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 83 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 97 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];
        vars[INTERACTION_TYPE]  = evt.w2list[W2LIST_GSEAGEN_ICHAN];
        vars[CURRENT_TYPE]      = evt.w2list[W2LIST_GSEAGEN_CC];
        vars[XSEC]              = evt.w2list[W2LIST_GSEAGEN_XSEC];
        vars[XSEC_MEAN]         = evt.w2list[W2LIST_GSEAGEN_XSEC_MEAN];
        vars[XSEC_DIFFERENTIAL] = evt.w2list[W2LIST_GSEAGEN_DXSEC];
        vars[XSEC_WATER]        = evt.w2list[W2LIST_GSEAGEN_WATERXSEC];
        vars[INT_LENGTH_WATER]  = evt.w2list[W2LIST_GSEAGEN_WATER_INT_LEN];
        vars[COLUMN_DEPTH]      = evt.w2list[W2LIST_GSEAGEN_COLUMN_DEPTH];
        vars[P_EARTH]           = evt.w2list[W2LIST_GSEAGEN_P_EARTH];
        vars[P_SCALE]           = evt.w2list[W2LIST_GSEAGEN_P_SCALE];
        vars[TARGET_A]          = evt.w2list[W2LIST_GSEAGEN_TARGETA];
        vars[TARGET_Z]          = evt.w2list[W2LIST_GSEAGEN_TARGETZ];

        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 142 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 177 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