1 #ifndef __JAANET__JEVTWEIGHTGSEAGEN__
2 #define __JAANET__JEVTWEIGHTGSEAGEN__
24 namespace JPP {
using namespace JAANET; }
36 public JClonable<JEvtWeight, JEvtWeightGSeaGen>
80 setNormalisation(
object.getNormalisation());
82 if (static_cast<const JFluxHelper&>(
object)) {
93 virtual const char*
const getName()
const override
114 setNormalisation(1.0 / Nevents / dt);
118 THROW(
JValueOutOfRange,
"JEvtWeightGSeaGen::configure(): Provided header is inconsistent with GSeaGen.");
145 if (static_cast<const JFluxHelper&>(*
this) &&
160 if (!static_cast<const JFluxHelper&>(*
this)) {
161 THROW(JNullPointerException,
"JEvtWeightGSeaGen::getWeight(): Unspecified flux function.");
163 THROW(
JIndexOutOfRange,
"JEvtWeightGSeaGen::getWeight(): Empty " << (evt.
w.size() < 3 ?
"w2-" :
"w3-") <<
"weight.");
static const int WEIGHTLIST_RESCALED_EVENT_RATE
Rescaled event rate [s-1].
bool is_gseagen(const JHead &header)
Check for generator.
double numberOfEvents
Number of events.
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
void configure(const JHead &header) override
Weight normalisation configuration.
std::vector< double > w
MC: Weights w[0]=w1, w[1]=w2, w[2]=w3 (see e.g. Tag list or km3net-dataformat/definitions) ...
static const int WEIGHTLIST_EVENT_RATE
Event rate [s-1].
double getFactor(const Evt &evt) const
Get weight-factor of given event.
Implementation of event weighting for GSeaGen data.
Template class for object cloning.
Helper class for event-weight factor.
void configure(const JEvtWeightFactor_t &factor)
Configure event-weight factor.
JEvtWeightGSeaGen(const JHead &header, const JFlux &flux)
Constructor.
Low-level interface for retrieving the flux corresponding to a given event.
JEvtWeightGSeaGen(const JEvtWeightGSeaGen &object)
Copy constructor.
virtual const char *const getName() const override
Get name.
bool check(const JHead &header) const override
Check whether header is consistent with this event weighter.
static const int WEIGHTLIST_DIFFERENTIAL_EVENT_RATE
Event rate per unit of flux (c.f. taglist document) [GeV m2 sr].
Exception for accessing a value in a collection that is outside of its range.
double getWeight(const Evt &evt) const override
Get weight of given event.
Exception for accessing an index in a collection that is outside of its range.
double numberOfSeconds
Time in seconds.
JEvtWeightGSeaGen(const JHead &header)
Constructor.
JEvtWeightGSeaGen()
Default Constructor.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.