herwig is hosted by Hepforge, IPPP Durham
Herwig 7.3.0
FIMassiveKinematics.h
1// -*- C++ -*-
2//
3// FILightKinematics.h is a part of Herwig - A multi-purpose Monte Carlo event generator
4// Copyright (C) 2002-2019 The Herwig Collaboration
5//
6// Herwig is licenced under version 3 of the GPL, see COPYING for details.
7// Please respect the MCnet academic guidelines, see GUIDELINES for details.
8//
9#ifndef HERWIG_FILightKinematics_H
10#define HERWIG_FILightKinematics_H
11//
12// This is the declaration of the FILightKinematics class.
13//
14
15#include "DipoleSplittingKinematics.h"
16
17namespace Herwig {
18
19using namespace ThePEG;
20
30
31public:
32
37
38public:
39
45 virtual pair<double,double> kappaSupport(const DipoleSplittingInfo& dIndex) const;
46
51 virtual pair<double,double> xiSupport(const DipoleSplittingInfo& dIndex) const;
52
56 virtual pair<double,double> zBoundaries(Energy,
58 const DipoleSplittingKernel&) const {
59 return {0.0,1.0};
60 }
61
70 virtual Energy dipoleScale(const Lorentz5Momentum& pEmitter,
71 const Lorentz5Momentum& pSpectator) const;
72
76 virtual Energy ptMax(Energy dScale,
77 double emX, double specX,
78 const DipoleSplittingInfo& dInfo,
79 const DipoleSplittingKernel& split) const;
80
84 virtual Energy ptMax(Energy dScale,
85 double, double,
86 const DipoleIndex& dIndex,
87 const DipoleSplittingKernel& split,
88 tPPtr emitter, tPPtr) const;
89
93 virtual Energy ptMax(Energy,
94 double, double,
95 const DipoleIndex&,
96 const DipoleSplittingKernel&) const {
97 // Only the DipoleSplittingInfo version should be used for massive
98 // dipoles, for now anyway.
99 assert(false);
100 return ZERO;
101 }
102
106 virtual Energy QMax(Energy dScale,
107 double emX, double specX,
108 const DipoleSplittingInfo& dInfo,
109 const DipoleSplittingKernel& split) const;
110
114 virtual Energy QMax(Energy,
115 double, double,
116 const DipoleIndex&,
117 const DipoleSplittingKernel&) const {
118 // Only the DipoleSplittingInfo version should be used for massive
119 // dipoles, for now anyway.
120 assert(false);
121 return ZERO;
122 }
123
127 virtual Energy PtFromQ(Energy scale, const DipoleSplittingInfo&) const;
128
132 virtual Energy QFromPt(Energy scale, const DipoleSplittingInfo&) const;
133
138 virtual double ptToRandom(Energy pt, Energy dScale,
139 double emX, double specX,
140 const DipoleIndex& dIndex,
141 const DipoleSplittingKernel&) const;
142
148 virtual bool generateSplitting(double kappa, double xi, double phi,
149 DipoleSplittingInfo& dIndex,
150 const DipoleSplittingKernel&);
151
157 virtual void generateKinematics(const Lorentz5Momentum& pEmitter,
158 const Lorentz5Momentum& pSpectator,
159 const DipoleSplittingInfo& dInfo);
160
161public:
162
166 template <class T>
167 inline double rootOfKallen (T a, T b, T c) const {
168 double sres=a*a + b*b + c*c - 2.*( a*b+a*c+b*c );
169 return sres>0.?sqrt( sres ):0.; }
170
171
172public:
173
181
187 void persistentInput(PersistentIStream & is, int version);
189
196 static void Init();
197
198protected:
199
206 virtual IBPtr clone() const;
207
212 virtual IBPtr fullclone() const;
214
215
216// If needed, insert declarations of virtual function defined in the
217// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
218
219
220private:
221
227
233
234};
235
236}
237
238#include "ThePEG/Utilities/ClassTraits.h"
239
240namespace ThePEG {
241
246template <>
247struct BaseClassTrait<Herwig::FIMassiveKinematics,1> {
250};
251
254template <>
255struct ClassTraits<Herwig::FIMassiveKinematics>
256 : public ClassTraitsBase<Herwig::FIMassiveKinematics> {
258 static string className() { return "Herwig::FIMassiveKinematics"; }
266 static string library() { return "HwDipoleShower.so"; }
267};
268
271}
272
273#endif /* HERWIG_FIMassiveKinematics_H */
DipoleIndex is used to index splitting generators for a particular dipole.
DipoleSplittingInfo contains all parameters to generate a full dipole splitting.
DipoleSplittingKernel is the base class for all kernels used within the dipole shower.
DipoleSplittingKinematics is the base class for dipole splittings as performed in the dipole shower.
FIMassiveKinematics implements massless splittings off a final-initial dipole.
virtual Energy dipoleScale(const Lorentz5Momentum &pEmitter, const Lorentz5Momentum &pSpectator) const
Return the dipole scale associated to the given pair of emitter and spectator.
virtual Energy QFromPt(Energy scale, const DipoleSplittingInfo &) const
Return the virtuality given a pt.
static ClassDescription< FIMassiveKinematics > initFIMassiveKinematics
The static object used to initialize the description of this class.
virtual Energy PtFromQ(Energy scale, const DipoleSplittingInfo &) const
Return the pt given a virtuality.
virtual double ptToRandom(Energy pt, Energy dScale, double emX, double specX, const DipoleIndex &dIndex, const DipoleSplittingKernel &) const
Return the random number associated to the given pt.
virtual pair< double, double > xiSupport(const DipoleSplittingInfo &dIndex) const
Return the boundaries in between the momentum fraction random number is to be sampled.
virtual Energy ptMax(Energy dScale, double emX, double specX, const DipoleSplittingInfo &dInfo, const DipoleSplittingKernel &split) const
Return the maximum pt for the given dipole scale.
virtual Energy QMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum virtuality for the given dipole scale.
virtual IBPtr clone() const
Make a simple clone of this object.
virtual pair< double, double > zBoundaries(Energy, const DipoleSplittingInfo &, const DipoleSplittingKernel &) const
Return the boundaries on the momentum fraction.
virtual Energy ptMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum pt for the given dipole scale.
FIMassiveKinematics & operator=(const FIMassiveKinematics &)=delete
The assignment operator is private and must never be called.
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
virtual void generateKinematics(const Lorentz5Momentum &pEmitter, const Lorentz5Momentum &pSpectator, const DipoleSplittingInfo &dInfo)
Generate the full kinematics given emitter and spectator momentum and a previously completeted Dipole...
double rootOfKallen(T a, T b, T c) const
Triangular / Kallen function.
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
virtual Energy ptMax(Energy dScale, double, double, const DipoleIndex &dIndex, const DipoleSplittingKernel &split, tPPtr emitter, tPPtr) const
Return the maximum pt for the given dipole scale.
static void Init()
The standard Init function used to initialize the interfaces.
FIMassiveKinematics()
The default constructor.
virtual pair< double, double > kappaSupport(const DipoleSplittingInfo &dIndex) const
Return the boundaries in between the evolution variable random number is to be sampled; the lower cuo...
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
virtual bool generateSplitting(double kappa, double xi, double phi, DipoleSplittingInfo &dIndex, const DipoleSplittingKernel &)
Generate splitting variables given three random numbers and the momentum fractions of the emitter and...
virtual Energy QMax(Energy dScale, double emX, double specX, const DipoleSplittingInfo &dInfo, const DipoleSplittingKernel &split) const
Return the maximum virtuality for the given dipole scale.
-*- C++ -*-
ThePEG::Ptr< Particle >::transient_pointer tPPtr
double sqrt(int x)
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
constexpr ZeroUnit ZERO
static string className()
static string library()