Rivet 3.1.9
DISKinematics.hh
1// -*- C++ -*-
2#ifndef RIVET_DISKinematics_HH
3#define RIVET_DISKinematics_HH
4
5#include "Rivet/Particle.hh"
6#include "Rivet/Event.hh"
7#include "Rivet/Projection.hh"
8#include "Rivet/Projections/DISLepton.hh"
9#include "Rivet/Projections/Beam.hh"
10
11namespace Rivet {
12
13
15 class DISKinematics : public Projection {
16 public:
17
19 DISKinematics(const DISLepton & lepton = DISLepton(),
20 const std::map<std::string,std::string> & opts =
21 std::map<std::string,std::string>())
22 : _theQ2(-1.0), _theW2(-1.0), _theX(-1.0), _theY(-1.0), _theS(-1.0), _theGH(-1.0)
23 {
24 setName("DISKinematics");
25 //addPdgIdPair(ANY, hadid);
26 declare(Beam(), "Beam");
27 declare(lepton, "Lepton");
28 }
29
32
34 using Projection::operator =;
35
36
37 protected:
38
40 virtual void project(const Event& e);
41
43 virtual CmpState compare(const Projection& p) const;
44
45
46 public:
47
52 double Q2() const { return _theQ2; }
53
55 double W2() const { return _theW2; }
56
58 double x() const { return _theX; }
59
61 double y() const { return _theY; }
62
64 double s() const { return _theS; }
65
67 double gammahad() const { return _theGH; }
68
69
70
72 const LorentzTransform& boostHCM() const {
73 return _hcm;
74 }
75
77 const LorentzTransform& boostBreit() const {
78 return _breit;
79 }
80
82 const Particle& beamHadron() const {
83 return _inHadron;
84 }
85
87 const Particle& beamLepton() const {
88 return _inLepton;
89 }
90
92 const Particle& scatteredLepton() const {
93 return _outLepton;
94 }
95
99 int orientation() const {
100 return sign(_inHadron.pz());
101 }
102
103
104 protected:
105
107 double _theQ2;
108
110 double _theW2;
111
113 double _theX;
114
116 double _theY;
117
119 double _theS;
121 double _theGH;
122
123
125 Particle _inHadron, _inLepton, _outLepton;
126
128 LorentzTransform _hcm;
129
131 LorentzTransform _breit;
132
133 };
134
135
136}
137
138#endif
Project out the incoming beams.
Definition Beam.hh:129
Get the DIS kinematic variables and relevant boosts for an event.
Definition DISKinematics.hh:15
virtual CmpState compare(const Projection &p) const
Compare with other projections.
DEFAULT_RIVET_PROJ_CLONE(DISKinematics)
Clone on the heap.
double y() const
The inelasticity .
Definition DISKinematics.hh:61
double Q2() const
The , i.e. the virtuality of the DIS photon.
Definition DISKinematics.hh:52
double W2() const
The , i.e. the squared CoM energy of the DIS-photon+hadron system.
Definition DISKinematics.hh:55
double x() const
The Bjorken .
Definition DISKinematics.hh:58
DISKinematics(const DISLepton &lepton=DISLepton(), const std::map< std::string, std::string > &opts=std::map< std::string, std::string >())
The default constructor.
Definition DISKinematics.hh:19
virtual void project(const Event &e)
Perform the projection operation on the supplied event.
double gammahad() const
The angle $f\gamma_{had}.
Definition DISKinematics.hh:67
double s() const
The centre of mass energy .
Definition DISKinematics.hh:64
Get the incoming and outgoing leptons in a DIS event.
Definition DISLepton.hh:18
Representation of a HepMC event, and enabler of Projection caching.
Definition Event.hh:22
Object implementing Lorentz transform calculations and boosts.
Definition LorentzTrans.hh:21
const PROJ & declare(const PROJ &proj, const std::string &name)
Register a contained projection (user-facing version)
Definition ProjectionApplier.hh:170
Base class for all Rivet projections.
Definition Projection.hh:29
void setName(const std::string &name)
Used by derived classes to set their name.
Definition Projection.hh:142
double p(const ParticleBase &p)
Unbound function access to p.
Definition ParticleBaseUtils.hh:684
Definition MC_Cent_pPb.hh:10
constexpr std::enable_if< std::is_arithmetic< NUM >::value, int >::type sign(NUM val)
Find the sign of a number.
Definition MathUtils.hh:266