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. More...

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

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

virtual clone_type	clone () const override
	Get clone of this object. 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 35 of file JEvtWeightFactorGSeaGen.hh.

Member Typedef Documentation

◆ clone_type

template<class JClonable_t , class JDerived_t >

typedef 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 = getFormula().GetNpar();
  
       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 getFormula().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 JTOOLS::JMultiMapTransformer< N, JArgument_t >::JMultiMapDefaultTransformer, JPHYSICS::JPDFTransformer< 5, JArgument_t >, JPHYSICS::JPDFTransformer< 4, JArgument_t >, JPHYSICS::JPDFTransformer< 3, JArgument_t >, JPHYSICS::JPDFTransformer< 2, JArgument_t >, JPHYSICS::JPDGTransformer_t< JArgument_t >, JPHYSICS::JPD0Transformer_t< JArgument_t >, and JPHYSICS::JPDFTransformer_t< JArgument_t >.

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