1 #ifndef __JAANET__JEVTWEIGHTFACTORMUPAGE__
2 #define __JAANET__JEVTWEIGHTFACTORMUPAGE__
17 #pragma GCC diagnostic push
18 #pragma GCC diagnostic ignored "-Wall"
20 #pragma GCC diagnostic pop
28 namespace JPP {
using namespace JAANET; }
72 const char* formula) :
73 TFormula(name, formula)
105 dir +=
i->getDirection();
109 const JRotation3D
R(D);
118 vars[
MUON_MULTIPLICITY] = (iBundle != evt.
mc_trks.cend() && iBundle->len > 0 ? iBundle->len : (double) muons.size());
124 const JCircle2D circle = JCircle2D(muons.cbegin(), muons.cend());
129 return this->DoEval(&vars[0]);
133 THROW(JNullPointerException,
"JEvtWeightFactorMupage::operator(): No muon for event " << evt.
id <<
'.' << endl);
JPredicate< JResult_t T::*, JComparison::eq > make_predicate(JResult_t T::*member, const JResult_t value)
Helper method to create predicate for data member.
Data structure for direction in three dimensions.
JTrack3E getTrack(const Trk &track)
Get track.
Muon bundle lateral spread [m].
JEvtWeightFactorMupage(const char *name, const char *formula)
Constructor.
variables
Indices of reweighting variables for MUPAGE.
bool is_muon(const Trk &track)
Test whether given track is a (anti-)muon.
#define THROW(JException_t, A)
Marco for throwing exception with std::ostream compatible message.
Number of reweighting variables; N.B. This enum value needs to be specified last! ...
bool is_finalstate(const Trk &track)
Test whether given track corresponds to a final state particle.
Average cosine of zenith angle.
JEvtWeightFactorMupage()
Default constructor.
static const int TRK_ST_MUONBUNDLE
initial state muon bundle (mupage)
double operator()(const Evt &evt) const
Get weighting factor for given event.
Implementation of reweighting factor for mupage events according to a specifiable ROOT TFormula...
then fatal The output file must have the wildcard in the name
int status
MC status code, see km3net-dataformat/definitions/trkmembers.csv for values.
then JCookie sh JDataQuality D $DETECTOR_ID R
int id
offline event identifier
do echo Generating $dir eval D
Muon bundle total energy [GeV].
double getDZ() const
Get z direction.
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.