herwig
is hosted by
Hepforge
,
IPPP Durham
Herwig
7.3.0
Shower
Dipole
Kinematics
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
17
namespace
Herwig
{
18
19
using namespace
ThePEG
;
20
29
class
FIMassiveKinematics
:
public
DipoleSplittingKinematics
{
30
31
public
:
32
36
FIMassiveKinematics
();
37
38
public
:
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,
57
const
DipoleSplittingInfo
&,
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
161
public
:
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
172
public
:
173
180
void
persistentOutput
(
PersistentOStream
& os)
const
;
181
187
void
persistentInput
(
PersistentIStream
& is,
int
version);
189
196
static
void
Init
();
197
198
protected
:
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
220
private
:
221
226
static
ClassDescription<FIMassiveKinematics>
initFIMassiveKinematics
;
227
232
FIMassiveKinematics
&
operator=
(
const
FIMassiveKinematics
&) =
delete
;
233
234
};
235
236
}
237
238
#include "ThePEG/Utilities/ClassTraits.h"
239
240
namespace
ThePEG
{
241
246
template
<>
247
struct
BaseClassTrait
<
Herwig
::FIMassiveKinematics,1> {
249
typedef
Herwig::DipoleSplittingKinematics
NthBase
;
250
};
251
254
template
<>
255
struct
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 */
Herwig::DipoleIndex
DipoleIndex is used to index splitting generators for a particular dipole.
Definition:
DipoleSplittingInfo.h:35
Herwig::DipoleSplittingInfo
DipoleSplittingInfo contains all parameters to generate a full dipole splitting.
Definition:
DipoleSplittingInfo.h:209
Herwig::DipoleSplittingKernel
DipoleSplittingKernel is the base class for all kernels used within the dipole shower.
Definition:
DipoleSplittingKernel.h:41
Herwig::DipoleSplittingKinematics
DipoleSplittingKinematics is the base class for dipole splittings as performed in the dipole shower.
Definition:
DipoleSplittingKinematics.h:41
Herwig::FIMassiveKinematics
FIMassiveKinematics implements massless splittings off a final-initial dipole.
Definition:
FIMassiveKinematics.h:29
Herwig::FIMassiveKinematics::dipoleScale
virtual Energy dipoleScale(const Lorentz5Momentum &pEmitter, const Lorentz5Momentum &pSpectator) const
Return the dipole scale associated to the given pair of emitter and spectator.
Herwig::FIMassiveKinematics::QFromPt
virtual Energy QFromPt(Energy scale, const DipoleSplittingInfo &) const
Return the virtuality given a pt.
Herwig::FIMassiveKinematics::initFIMassiveKinematics
static ClassDescription< FIMassiveKinematics > initFIMassiveKinematics
The static object used to initialize the description of this class.
Definition:
FIMassiveKinematics.h:226
Herwig::FIMassiveKinematics::PtFromQ
virtual Energy PtFromQ(Energy scale, const DipoleSplittingInfo &) const
Return the pt given a virtuality.
Herwig::FIMassiveKinematics::ptToRandom
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.
Herwig::FIMassiveKinematics::xiSupport
virtual pair< double, double > xiSupport(const DipoleSplittingInfo &dIndex) const
Return the boundaries in between the momentum fraction random number is to be sampled.
Herwig::FIMassiveKinematics::ptMax
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.
Herwig::FIMassiveKinematics::QMax
virtual Energy QMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum virtuality for the given dipole scale.
Definition:
FIMassiveKinematics.h:114
Herwig::FIMassiveKinematics::clone
virtual IBPtr clone() const
Make a simple clone of this object.
Herwig::FIMassiveKinematics::zBoundaries
virtual pair< double, double > zBoundaries(Energy, const DipoleSplittingInfo &, const DipoleSplittingKernel &) const
Return the boundaries on the momentum fraction.
Definition:
FIMassiveKinematics.h:56
Herwig::FIMassiveKinematics::ptMax
virtual Energy ptMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum pt for the given dipole scale.
Definition:
FIMassiveKinematics.h:93
Herwig::FIMassiveKinematics::operator=
FIMassiveKinematics & operator=(const FIMassiveKinematics &)=delete
The assignment operator is private and must never be called.
Herwig::FIMassiveKinematics::persistentInput
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
Herwig::FIMassiveKinematics::generateKinematics
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...
Herwig::FIMassiveKinematics::rootOfKallen
double rootOfKallen(T a, T b, T c) const
Triangular / Kallen function.
Definition:
FIMassiveKinematics.h:167
Herwig::FIMassiveKinematics::fullclone
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
Herwig::FIMassiveKinematics::ptMax
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.
Herwig::FIMassiveKinematics::Init
static void Init()
The standard Init function used to initialize the interfaces.
Herwig::FIMassiveKinematics::FIMassiveKinematics
FIMassiveKinematics()
The default constructor.
Herwig::FIMassiveKinematics::kappaSupport
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...
Herwig::FIMassiveKinematics::persistentOutput
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
Herwig::FIMassiveKinematics::generateSplitting
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...
Herwig::FIMassiveKinematics::QMax
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.
ThePEG::ClassDescription
ThePEG::PersistentIStream
ThePEG::PersistentOStream
Herwig
-*- C++ -*-
Definition:
BasicConsistency.h:17
ThePEG
ThePEG::tPPtr
ThePEG::Ptr< Particle >::transient_pointer tPPtr
ThePEG::sqrt
double sqrt(int x)
ThePEG::IBPtr
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
ThePEG::ZERO
constexpr ZeroUnit ZERO
ThePEG::BaseClassTrait
ThePEG::BaseClassTrait::NthBase
int NthBase
ThePEG::ClassTraitsBase
ThePEG::ClassTraitsBase::className
static string className()
ThePEG::ClassTraitsBase::library
static string library()
ThePEG::ClassTraits
Generated on Thu Jun 20 2024 17:50:53 for Herwig by
1.9.6