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 JPointer< JClass_t >	pointer_type

typedef JStorage< JClass_t, JMemory_t >	storage_type

typedef JMemory_t< JClass_t >	memory_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...

double	getFactor (const Evt &evt) const
	Get weight-factor of given event. More...

const JSharedPointer &	getSharedPointer () const
	Get shared pointer. More...

JSharedPointer &	getSharedPointer ()
	Get shared pointer. More...

void	setSharedPointer (const JSharedPointer &object)
	Set shared pointer. More...

virtual void	reset () override
	Reset pointer. More...

template<class T >
void	reset (const JPointer< T > &pointer)
	Reset pointer. More...

void	reset (JClass_t *p)
	Reset pointer. More...

template<class T >
void	set (const JPointer< T > &pointer)
	Set pointer. More...

void	initialise ()
	Initialise counter. More...

void	attach (const JSharedCounter &object)
	Attach this counter to given shared counter object. More...

bool	detach ()
	Detach. More...

const int	getCount ()
	Get count. More...

void	recreate ()
	Recreate object in memory. More...

void	create ()
	Create object in memory. More...

void	create (const unsigned int size)
	Create array of objects in memory. More...

virtual JClass_t *	get () const override
	Get pointer. More...

JClass_t *const &	getReference () const
	Get rereference to internal pointer. More...

JClass_t *&	getReference ()
	Get rereference to internal pointer. More...

virtual bool	equals (const JAbstractPointer &object) const
	Equals. More...

bool	is_valid () const
	Check validity of pointer. More...

JClass_t *	operator-> () const
	Smart pointer operator. More...

	operator JClass_t * () const
	Type conversion operator. More...

Protected Member Functions
void	set (const JSharedPointer &object)
	Set pointer. More...

virtual void	set (JClass_t *p) override
	Set pointer. More...

void	release ()
	Release memory. More...

Protected Attributes
int *	counter

JClass_t *	__p
	pointer to object 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 JClass_t, template< class > class JMemory_t = JNew>

typedef JPointer<JClass_t> JLANG::JSharedPointer< JClass_t, JMemory_t >::pointer_type

inherited

Definition at line 34 of file JSharedPointer.hh.

template<class JClass_t, template< class > class JMemory_t = JNew>

typedef JStorage<JClass_t, JMemory_t> JLANG::JSharedPointer< JClass_t, JMemory_t >::storage_type

inherited

Definition at line 35 of file JSharedPointer.hh.

template<class JClass_t, template< class > class JMemory_t = JNew>

typedef JMemory_t<JClass_t> JLANG::JStorage< JClass_t, JMemory_t >::memory_type

inherited

Definition at line 33 of file JStorage.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 (is_valid()) {
       
         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 63 of file JEvtWeightFactorHelper.hh.

     {
       using namespace JPP;
       
       JEvtWeightFactor_t* p = dynamic_cast<JEvtWeightFactor_t*>(factor.clone());
 
       if (p != NULL) {
         reset(p);
       } else {
         THROW(JNullPointerException, "JEvtWeightFactorHelper::configure(): Invalid downcast to class derived from JEvtWeightFactor.");
       }
     }

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 83 of file JEvtWeightFactorHelper.hh.

     {
       using namespace JPP;
       
       if (is_valid()) {
         return get()->getFactor(evt);
       } else {
         THROW(JNullPointerException, "JEvtWeightFactorHelper::getFactor(): Unspecified event-weighting factor.");
       }
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

const JSharedPointer& JLANG::JSharedPointer< JClass_t, JMemory_t >::getSharedPointer ( ) const

inlineinherited

Get shared pointer.

Returns: this shared pointer

Definition at line 95 of file JSharedPointer.hh.

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

template<class JClass_t, template< class > class JMemory_t = JNew>

JSharedPointer& JLANG::JSharedPointer< JClass_t, JMemory_t >::getSharedPointer ( )

inlineinherited

Get shared pointer.

Returns: this shared pointer

Definition at line 106 of file JSharedPointer.hh.

     {
       return static_cast<JSharedPointer&>(*this);
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JSharedPointer< JClass_t, JMemory_t >::setSharedPointer ( const JSharedPointer< JClass_t, JMemory_t > & object )

inlineinherited

Set shared pointer.

Parameters

object shared pointer

Definition at line 117 of file JSharedPointer.hh.

     {
       if (this->get() != object.get()) {
 
         this->reset();
 
         if (object.is_valid()) {
           this->set(object);
         }
       }
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

virtual void JLANG::JSharedPointer< JClass_t, JMemory_t >::reset ( )

inlineoverridevirtualinherited

Reset pointer.

The reference counter is decremented by one and the object pointed to previously is deleted when its reference counter is zero.

Reimplemented from JLANG::JStorage< JClass_t, JMemory_t >.

Definition at line 171 of file JSharedPointer.hh.

     {
       if (this->detach()) {
         storage_type::reset();
       }
 
       pointer_type::reset();
     }

template<class JClass_t>

template<class T >

void JLANG::JPointer< JClass_t >::reset ( const JPointer< T > & pointer )

inlineinherited

Reset pointer.

Parameters

pointer pointer to object

Definition at line 108 of file JPointer.hh.

     {
       this->reset(pointer.get());
     }

template<class JClass_t>

void JLANG::JAbstractPointer< JClass_t >::reset ( JClass_t * p )

inlineinherited

Reset pointer.

Parameters

p	pointer to object

Definition at line 94 of file JAbstractPointer.hh.

     {
       if (this->get() != p) {
 
         this->reset();
 
         if (p != NULL) {
           this->set(p);
         }
       }
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JSharedPointer< JClass_t, JMemory_t >::set ( const JSharedPointer< JClass_t, JMemory_t > & object )

inlineprotectedinherited

Set pointer.

The reference counter of the shared object pointed to is incremented by one.

Parameters

object shared pointer

Definition at line 188 of file JSharedPointer.hh.

     {
       pointer_type::set(object.get());
 
       this->attach(object);
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

virtual void JLANG::JSharedPointer< JClass_t, JMemory_t >::set ( JClass_t * p )

inlineoverrideprotectedvirtualinherited

Set pointer.

The reference counter of the shared object pointed to is initialised to one.

Parameters

p	pointer to derived class object

Reimplemented from JLANG::JPointer< JClass_t >.

Definition at line 202 of file JSharedPointer.hh.

     {
       pointer_type::set(p);
 
       this->initialise();
     }

template<class JClass_t>

template<class T >

void JLANG::JPointer< JClass_t >::set ( const JPointer< T > & pointer )

inlineinherited

Set pointer.

Parameters

pointer pointer to object

Definition at line 96 of file JPointer.hh.

     {
       this->set(pointer.get());
     }

void JLANG::JSharedCounter::initialise ( )

inlineinherited

Initialise counter.

Definition at line 33 of file JSharedCounter.hh.

     {
       detach();
 
       counter = new int(1);
     }

void JLANG::JSharedCounter::attach ( const JSharedCounter & object )

inlineinherited

Attach this counter to given shared counter object.

Parameters

object shared counter

Definition at line 46 of file JSharedCounter.hh.

     {
       detach();
 
       counter = object.counter;
 
       if (counter != NULL) {
         ++(*counter);
       }
     }

bool JLANG::JSharedCounter::detach ( )

inlineinherited

Detach.

Returns: true if counter at zero; else false

Definition at line 63 of file JSharedCounter.hh.

     {
       if (counter != NULL) {
 
         if (--(*counter) == 0) {
 
           delete counter;
 
           counter = NULL;
 
           return true;
         }
 
         counter = NULL;
       }
 
       return false;
     }

const int JLANG::JSharedCounter::getCount ( )

inlineinherited

Get count.

Returns: count

Definition at line 88 of file JSharedCounter.hh.

     {
       return (counter != NULL ? *counter : 0);
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JStorage< JClass_t, JMemory_t >::recreate ( )

inlineinherited

Recreate object in memory.

A new object is created if no memory is allocated yet, else the previously created object is maintained.

Definition at line 57 of file JStorage.hh.

     {
       if (!this->is_valid()) {
         this->set(memory_type::create());
       }
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JStorage< JClass_t, JMemory_t >::create ( )

inlineinherited

Create object in memory.

The memory allocated by a previously created object will be released.

Definition at line 69 of file JStorage.hh.

     {
       this->reset(memory_type::create());
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JStorage< JClass_t, JMemory_t >::create ( const unsigned int size )

inlineinherited

Create array of objects in memory.

The memory allocated by previously created objects will be released.

Parameters

size	number of elements

Definition at line 81 of file JStorage.hh.

     {
       this->reset(memory_type::create(size));
     }

template<class JClass_t, template< class > class JMemory_t = JNew>

void JLANG::JStorage< JClass_t, JMemory_t >::release ( )

inlineprotectedinherited

Release memory.

Definition at line 91 of file JStorage.hh.

     { 
       memory_type::release(this->get());
     }

template<class JClass_t>

virtual JClass_t* JLANG::JPointer< JClass_t >::get ( ) const

inlineoverridevirtualinherited

Get pointer.

Returns: pointer to object

Implements JLANG::JAbstractPointer< JClass_t >.

Reimplemented in JLANG::JMultiPointer< JTypeList< JHead_t, JTail_t > >, JLANG::JMultiPointer< JClass_t >, JLANG::JMultiPointer< JTail_t >, JLANG::JMultiPointer< JHead_t >, JLANG::JMultiPointer< T >, and JLANG::JMultiPointer< JTypeList< JDAQEvent, JTypelist_t > >.

Definition at line 64 of file JPointer.hh.

     {
       return this->__p;
     }

template<class JClass_t>

JClass_t* const& JLANG::JPointer< JClass_t >::getReference ( ) const

inlineinherited

Get rereference to internal pointer.

Returns: reference to internal pointer

Definition at line 119 of file JPointer.hh.

     {
       return __p;
     }

template<class JClass_t>

JClass_t* & JLANG::JPointer< JClass_t >::getReference ( )

inlineinherited

Get rereference to internal pointer.

Returns: reference to internal pointer

Definition at line 130 of file JPointer.hh.

     {
       return __p;
     }

template<class JClass_t>

virtual bool JLANG::JAbstractPointer< JClass_t >::equals ( const JAbstractPointer< JClass_t > & object ) const

inlinevirtualinherited

Equals.

The equality is evaluated by comparison of the internal pointers.

Parameters

object abstract pointer

Returns: true if equals; else false

Definition at line 50 of file JAbstractPointer.hh.

     {
       return this->get() == object.get();
     }

template<class JClass_t>

bool JLANG::JAbstractPointer< JClass_t >::is_valid ( ) const

inlineinherited

Check validity of pointer.

Returns: true if pointer not null; else false

Definition at line 83 of file JAbstractPointer.hh.

     { 
       return this->get() != NULL; 
     }

template<class JClass_t>

JClass_t* JLANG::JAbstractPointer< JClass_t >::operator-> ( ) const

inlineinherited

Smart pointer operator.

Returns: pointer to object

Definition at line 112 of file JAbstractPointer.hh.

     { 
       if (!is_valid())
         THROW(JNullPointerException, "JAbstractPointer::operator->()");
       else
         return this->get();
     }

template<class JClass_t>

JLANG::JAbstractPointer< JClass_t >::operator JClass_t * ( ) const

inlineinherited

Type conversion operator.

Returns: pointer to object

Definition at line 126 of file JAbstractPointer.hh.

     {
       return this->get();
     }

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;
       }

Member Data Documentation

int* JLANG::JSharedCounter::counter

protectedinherited

Definition at line 94 of file JSharedCounter.hh.

template<class JClass_t>

JClass_t* JLANG::JPointer< JClass_t >::__p

protectedinherited

pointer to object

Definition at line 136 of file JPointer.hh.

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

JEvtWeightFactorGSeaGen.hh

Public Types

Public Member Functions

Protected Member Functions

Protected Attributes

Friends

Detailed Description

Member Typedef Documentation

Member Enumeration Documentation

Constructor & Destructor Documentation

Member Function Documentation

Friends And Related Function Documentation

Member Data Documentation