Jpp  18.5.0
the software that should make you happy
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
JEvtWeightFactorMupage.hh
Go to the documentation of this file.
1 #ifndef __JAANET__JEVTWEIGHTFACTORMUPAGE__
2 #define __JAANET__JEVTWEIGHTFACTORMUPAGE__
3 
6 
7 #include "JLang/JException.hh"
8 #include "JLang/JPredicate.hh"
9 
10 #include "JGeometry2D/JCircle2D.hh"
11 #include "JGeometry3D/JVector3D.hh"
13 
14 #pragma GCC diagnostic push
15 #pragma GCC diagnostic ignored "-Wall"
16 #include "TFormula.h"
17 #pragma GCC diagnostic pop
18 
19 
20 /**
21  * \author bjung
22  */
23 
24 namespace JAANET {}
25 namespace JPP { using namespace JAANET; }
26 
27 namespace JAANET {
28 
29  /**
30  * Implementation of reweighting factor for mupage events according to a specifiable ROOT TFormula.
31  *
32  * Note: The ROOT TFormula may assume any number of parameters, but should be restricted to\n
33  * the physical variables listed among `JEvtWeightFactorMupage::variables`.\n
34  * These variables may be specified within the ROOT TFormula as 'x[<index>]',\n
35  * where <index> corresponds to the index of the desired physical variable within `JEvtWeightFactorMupage::variables`.
36  */
38  public TFormula
39  {
40  /**
41  * Indices of reweighting variables for MUPAGE.
42  */
43  enum variables {
44  MUON_MULTIPLICITY, //!< Muon multiplicity
45  MEAN_ZENITH_ANGLE, //!< Average cosine of zenith angle
46  TOTAL_MUON_ENERGY, //!< Muon bundle total energy [GeV]
47  LATERAL_SPREAD, //!< Muon bundle lateral spread [m]
48 
49  NUMBER_OF_VARIABLES //!< Number of reweighting variables; \n
50  //!< N.B.\ This enum value needs to be specified last!
51  };
52 
53 
54  /**
55  * Default constructor.
56  */
58  TFormula()
59  {}
60 
61 
62  /**
63  * Constructor.
64  *
65  * \param name name
66  * \param formula formula
67  */
69  const char* formula) :
70  TFormula(name, formula)
71  {}
72 
73 
74  /**
75  * Get weighting factor for given event.
76  *
77  * \param evt event
78  * \return weighting factor
79  */
80  double operator()(const Evt& evt) const
81  {
82  using namespace std;
83  using namespace JPP;
84 
85  vector<JTrack3E> muons;
86 
87  for (vector<Trk>::const_iterator i = evt.mc_trks.cbegin(); i != evt.mc_trks.cend(); ++i) {
88  if (is_muon(*i) && is_finalstate(*i)) {
89  muons.push_back(getTrack(*i));
90  }
91  }
92 
93  if (!muons.empty()) {
94 
95  Double_t vars[NUMBER_OF_VARIABLES] = { 0.0 };
96 
97  double Etot = 0.0;
98  JVector3D dir;
99 
100  for (vector<JTrack3E>::const_iterator i = muons.begin(); i != muons.end(); ++i) {
101  Etot += i->getE();
102  dir += i->getDirection();
103  }
104 
105  const JDirection3D D(dir);
106  const JRotation3D R(D);
107 
108  for (vector<JTrack3E>::iterator i = muons.begin(); i != muons.end(); ++i) {
109  i->rotate(R);
110  }
111 
112  vector<Trk>::const_iterator iBundle = find_if(evt.mc_trks.cbegin(), evt.mc_trks.cend(),
114 
115  vars[MUON_MULTIPLICITY] = (iBundle != evt.mc_trks.cend() && iBundle->len > 0 ? iBundle->len : (double) muons.size());
116  vars[TOTAL_MUON_ENERGY] = (iBundle != evt.mc_trks.cend() ? iBundle->E : Etot);
117  vars[MEAN_ZENITH_ANGLE] = D.getDZ();
118 
119  if (this->GetNdim() > (int) LATERAL_SPREAD) {
120 
121  const JCircle2D circle = JCircle2D(muons.cbegin(), muons.cend()); // smallest enclosing circle
122 
123  vars[LATERAL_SPREAD] = circle.getRadius();
124  }
125 
126  return this->DoEval(&vars[0]);
127 
128  } else {
129 
130  THROW(JNullPointerException, "JEvtWeightFactorMupage::operator(): No muon for event " << evt.id << '.' << endl);
131  }
132  }
133  };
134 }
135 
136 #endif
Exceptions.
JPredicate< JResult_t T::*, JComparison::eq > make_predicate(JResult_t T::*member, const JResult_t value)
Helper method to create predicate for data member.
Definition: JPredicate.hh:128
Data structure for direction in three dimensions.
Definition: JDirection3D.hh:33
JTrack3E getTrack(const Trk &track)
Get track.
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.
Definition: JException.hh:712
then echo Enter input within $TIMEOUT_S seconds echo n User name
Definition: JCookie.sh:42
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.
JEvtWeightFactorMupage()
Default constructor.
static const int TRK_ST_MUONBUNDLE
initial state muon bundle (mupage)
Definition: trkmembers.hh:18
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...
int status
MC status code, see km3net-dataformat/definitions/trkmembers.csv for values.
Definition: Trk.hh:28
then JCookie sh JDataQuality D $DETECTOR_ID R
Definition: JDataQuality.sh:41
int id
offline event identifier
Definition: Evt.hh:22
do echo Generating $dir eval D
Definition: JDrawLED.sh:53
double getDZ() const
Get z direction.
Definition: JVersor3D.hh:117
std::vector< Trk > mc_trks
MC: list of MC truth tracks.
Definition: Evt.hh:49
The Evt class respresent a Monte Carlo (MC) event as well as an offline event.
Definition: Evt.hh:20