JAANET::JEvtWeightCorsika Struct Reference

Implementation of event weighting for Corsika data. More...

#include <JEvtWeightCorsika.hh>

Inheritance diagram for JAANET::JEvtWeightCorsika:

Public Types
typedef JEvtWeightFactorHelper < JEvtWeightFactor_t >	JEvtWeightFactorHelper_t
typedef std::shared_ptr < JEvtWeightFactor_t >	pointer_type
typedef JClonable< JClonable_t > ::clone_type	clone_type

Public Member Functions
	JEvtWeightCorsika ()
	Default Constructor. More...
	JEvtWeightCorsika (const JHead &header)
	Constructor. More...
	JEvtWeightCorsika (const JHead &header, const JFlux &flux)
	Constructor. More...
	JEvtWeightCorsika (const JEvtWeightCorsika &object)
	Copy constructor. More...
void	configure (const JHead &header) override
	Weight configuration. More...
bool	check (const JHead &header) const override
	Check whether header is consistent with this event weighter. More...
double	getWeight (const Evt &evt) const override
	Get weight of 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...
virtual clone_type	clone () const override
	Get clone of this object. More...

Detailed Description

Implementation of event weighting for Corsika data.

Definition at line 35 of file JEvtWeightCorsika.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.

template<class JClonable_t, class JDerived_t = JNullType>

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

inherited

Definition at line 61 of file JClonable.hh.

Constructor & Destructor Documentation

JAANET::JEvtWeightCorsika::JEvtWeightCorsika ( )

inline

Default Constructor.

Definition at line 42 of file JEvtWeightCorsika.hh.

                        :
      JFluxHelper()
    {}

JAANET::JEvtWeightCorsika::JEvtWeightCorsika ( const JHead &  header )

inline

Constructor.

Parameters

header

header

Definition at line 52 of file JEvtWeightCorsika.hh.

                                           :
      JFluxHelper()
    {
      configure(header);
    }

JAANET::JEvtWeightCorsika::JEvtWeightCorsika ( const JHead &  header, const JFlux &  flux  )

inline

Constructor.

Parameters

header	header
flux	flux

Definition at line 65 of file JEvtWeightCorsika.hh.

    {
      configure(header);
      JFluxHelper::configure(flux);
    }

JAANET::JEvtWeightCorsika::JEvtWeightCorsika ( const JEvtWeightCorsika &  object )

inline

Copy constructor.

Parameters

object

original object

Definition at line 78 of file JEvtWeightCorsika.hh.

    {
      this->norm = object.norm;

      if (static_cast<const JFluxHelper&>(object)) {
        JFluxHelper::configure(*(object.get()));
      }
    }

Member Function Documentation

void JAANET::JEvtWeightCorsika::configure ( const JHead &  header )

inlineoverride

Weight configuration.

Parameters

header

header

Definition at line 93 of file JEvtWeightCorsika.hh.

    {
      using namespace JPP;
      
      if (check(header)) {

        norm = 1.0 / header.genvol.numberOfEvents;
        
      } else {

        THROW(JValueOutOfRange, "JEvtWeightCorsika::configure(): Provided header is inconsistent with Corsika.");
      }
    }

bool JAANET::JEvtWeightCorsika::check ( const JHead &  header ) const

inlineoverride

Check whether header is consistent with this event weighter.

Parameters

header

header

Returns

true if consistent; else false

Definition at line 114 of file JEvtWeightCorsika.hh.

    {
      return is_gseagen(header);
    }

double JAANET::JEvtWeightCorsika::getWeight ( const Evt &  evt ) const

inlineoverride

Get weight of given event.

Parameters

evt

event

Returns

weight [Hz]

Definition at line 126 of file JEvtWeightCorsika.hh.

    {
      using namespace JPP;
      
      const double& Nretries = evt.w2list[W2LIST_GSEAGEN_N_RETRIES];
      
      if (static_cast<const JFluxHelper&>(*this) &&
          evt.w.size() > WEIGHTLIST_DIFFERENTIAL_EVENT_RATE) {
          
        return evt.w[WEIGHTLIST_DIFFERENTIAL_EVENT_RATE] * getFactor(evt) * getNormalisation(evt) / (1 + Nretries);
        
      } else if (evt.w.size() > WEIGHTLIST_RESCALED_EVENT_RATE) {
        
        return evt.w[WEIGHTLIST_RESCALED_EVENT_RATE];

      } else if (evt.w.size() > WEIGHTLIST_EVENT_RATE) {

        return evt.w[WEIGHTLIST_EVENT_RATE] * getNormalisation(evt) / (1 + Nretries);
        
      } else {

        if (!static_cast<const JFluxHelper&>(*this)) {
          THROW(JNullPointerException, "JEvtWeightCorsika::getWeight(): Unspecified flux function.");
        } else {
          THROW(JIndexOutOfRange,      "JEvtWeightCorsika::getWeight(): Empty " << (evt.w.size() < 3 ? "w2-" : "w3-") << "weight.");
        }
      }
    }

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

template<class JClonable_t, class JDerived_t = JNullType>

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::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 >, JTOOLS::JMultiMapTransformer< N, JArgument_t >::JMultiMapDefaultTransformer, JPHYSICS::JPDFTransformer_t< JArgument_t >, JLANG::JClonable< JMatch< hit_type > >, JLANG::JClonable< JMatch< hit_type >, JMatch3D >, JLANG::JClonable< JMatch< JTRIGGER::JHitR1 > >, and JLANG::JClonable< JMatch< JTRIGGER::JHitR1 >, JMatch3G< JTRIGGER::JHitR1 > >.

Definition at line 69 of file JClonable.hh.

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

---

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

• JEvtWeightCorsika.hh