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 , INITIAL_STATE_ENERGY , NEUTRINO_ENERGY , FINAL_STATE_ENERGY , LEADING_LEPTON_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 JClonable< JClonable_t >::clone_type	clone_type

Public Member Functions
	JEvtWeightFactorGSeaGen ()
	Default constructor.

bool	is_valid () const override final
	Check whether this formula is valid.

double	getFactor (const Evt &evt) const override final
	Get weighting factor for given event.

virtual clone_type	clone () const override
	Get clone of this object.

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 35 of file JEvtWeightFactorGSeaGen.hh.

Member Typedef Documentation

◆ clone_type

template<class JClonable_t , class JDerived_t >

JClonable<JClonable_t>::clone_type JLANG::JClonable< JClonable_t, JDerived_t >::clone_type

inherited

Definition at line 61 of file JClonable.hh.

Member Enumeration Documentation

◆ variables

enum JAANET::JEvtWeightFactorGSeaGen::variables

Indices of reweighting variables for GSeaGen.

Enumerator
COSTH	Cosine zenith angle.
INITIAL_STATE_ENERGY	Initial state energy.
NEUTRINO_ENERGY	Initial state neutrino energy.
FINAL_STATE_ENERGY	Final state energy.
LEADING_LEPTON_ENERGY	Final state leading lepton 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 43 of file JEvtWeightFactorGSeaGen.hh.

                   {
      COSTH,                 //!< Cosine zenith angle
      INITIAL_STATE_ENERGY,  //!< Initial state energy
      NEUTRINO_ENERGY,       //!< Initial state neutrino energy
      FINAL_STATE_ENERGY,    //!< Final state energy
      LEADING_LEPTON_ENERGY, //!< Final state leading lepton 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

◆ JEvtWeightFactorGSeaGen()

JAANET::JEvtWeightFactorGSeaGen::JEvtWeightFactorGSeaGen ( )

inline

Default constructor.

Definition at line 72 of file JEvtWeightFactorGSeaGen.hh.

73 {}

Member Function Documentation

◆ is_valid()

bool JAANET::JEvtWeightFactorGSeaGen::is_valid ( ) const

inlinefinaloverride

Check whether this formula is valid.

Returns: true if valid; else false

Definition at line 81 of file JEvtWeightFactorGSeaGen.hh.

    {
      const int N = this->GetNdim();
 
      return N >= 0 && N <= (int) NUMBER_OF_VARIABLES;
    }

◆ getFactor()

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

inlinefinaloverride

Get weighting factor for given event.

Parameters

evt event

Returns: weighting factor

Definition at line 95 of file JEvtWeightFactorGSeaGen.hh.

    {
      using namespace std;
      using namespace JPP;
 
      Double_t vars[NUMBER_OF_VARIABLES] = { 0.0 };
        
      const Trk& neutrino       = get_neutrino      (evt);
      const Trk& leading_lepton = get_leading_lepton(evt);
      
      vars[COSTH]                 = -neutrino.dir.z;
      vars[INITIAL_STATE_ENERGY]  = getE0(evt);
      vars[NEUTRINO_ENERGY]       = neutrino.E;
      vars[FINAL_STATE_ENERGY]    = getE1(evt);
      vars[LEADING_LEPTON_ENERGY] = leading_lepton.E;
      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->EvalPar(&vars[0]);
    }

◆ clone()

template<class JClonable_t , class JDerived_t >

virtual clone_type JLANG::JClonable< JClonable_t, JDerived_t >::clone ( ) const

inlineoverridevirtualinherited

Get clone of this object.

Returns: pointer to newly created object

Reimplemented in JPHYSICS::JPD0Transformer_t< JArgument_t >, JPHYSICS::JPDFTransformer< 2, JArgument_t >, JPHYSICS::JPDFTransformer< 3, JArgument_t >, JPHYSICS::JPDFTransformer< 4, JArgument_t >, JPHYSICS::JPDFTransformer< 5, JArgument_t >, JPHYSICS::JPDFTransformer_t< JArgument_t >, JPHYSICS::JPDGTransformer_t< JArgument_t >, and JTOOLS::JMultiMapTransformer< N, JArgument_t >::JMultiMapDefaultTransformer.

Definition at line 69 of file JClonable.hh.

    {
      return new JDerived_t(static_cast<const JDerived_t&>(*this));
    }

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

JEvtWeightFactorGSeaGen.hh

Public Types

Public Member Functions