Rivet  3.1.5
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 
11 namespace 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 
33 
34  protected:
35 
37  virtual void project(const Event& e);
38 
40  virtual CmpState compare(const Projection& p) const;
41 
42 
43  public:
44 
46  double Q2() const { return _theQ2; }
47 
49  double W2() const { return _theW2; }
50 
52  double x() const { return _theX; }
53 
55  double y() const { return _theY; }
56 
58  double s() const { return _theS; }
59 
61  double gammahad() const { return _theGH; }
62 
63 
64 
66  const LorentzTransform& boostHCM() const {
67  return _hcm;
68  }
69 
71  const LorentzTransform& boostBreit() const {
72  return _breit;
73  }
74 
76  const Particle& beamHadron() const {
77  return _inHadron;
78  }
79 
81  const Particle& beamLepton() const {
82  return _inLepton;
83  }
84 
86  const Particle& scatteredLepton() const {
87  return _outLepton;
88  }
89 
93  int orientation() const {
94  return sign(_inHadron.pz());
95  }
96 
97 
98  private:
99 
101  double _theQ2;
102 
104  double _theW2;
105 
107  double _theX;
108 
110  double _theY;
111 
113  double _theS;
115  double _theGH;
116 
117 
119  Particle _inHadron, _inLepton, _outLepton;
120 
122  LorentzTransform _hcm;
123 
125  LorentzTransform _breit;
126 
127  };
128 
129 
130 }
131 
132 #endif
void setName(const std::string &name)
Used by derived classes to set their name.
Definition: Projection.hh:142
Definition: MC_Cent_pPb.hh:10
double s() const
The centre of mass energy .
Definition: DISKinematics.hh:58
double y() const
The inelasticity .
Definition: DISKinematics.hh:55
double gammahad() const
The angle $f{had}.
Definition: DISKinematics.hh:61
Get the DIS kinematic variables and relevant boosts for an event.
Definition: DISKinematics.hh:15
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
Particle representation, either from a HepMC::GenEvent or reconstructed.
Definition: Particle.hh:18
double Q2() const
The .
Definition: DISKinematics.hh:46
virtual CmpState compare(const Projection &p) const
Compare with other projections.
double x() const
The Bjorken .
Definition: DISKinematics.hh:52
double W2() const
The .
Definition: DISKinematics.hh:49
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
DEFAULT_RIVET_PROJ_CLONE(DISKinematics)
Clone on the heap.
const PROJ & declare(const PROJ &proj, const std::string &name)
Register a contained projection (user-facing version)
Definition: ProjectionApplier.hh:170
Project out the incoming beams.
Definition: Beam.hh:129
Get the incoming and outgoing leptons in a DIS event.
Definition: DISLepton.hh:17
virtual void project(const Event &e)
Perform the projection operation on the supplied event.
Base class for all Rivet projections.
Definition: Projection.hh:29
std::enable_if< std::is_arithmetic< NUM >::value, int >::type sign(NUM val)
Find the sign of a number.
Definition: MathUtils.hh:266