herwig
is hosted by
Hepforge
,
IPPP Durham
Herwig
7.3.0
Shower
Dipole
Kinematics
FFMassiveKinematics.h
1
// -*- C++ -*-
2
//
3
// FFMassiveKinematics.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_FFMassiveKinematics_H
10
#define HERWIG_FFMassiveKinematics_H
11
//
12
// This is the declaration of the FFMassiveKinematics class.
13
//
14
15
#include "DipoleSplittingKinematics.h"
16
17
namespace
Herwig
{
18
19
using namespace
ThePEG
;
20
29
class
FFMassiveKinematics
:
public
DipoleSplittingKinematics
{
30
31
public
:
32
36
FFMassiveKinematics
();
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
spectator)
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
& split)
const
;
142
148
virtual
bool
generateSplitting
(
double
kappa,
double
xi,
double
phi,
149
DipoleSplittingInfo
& dIndex,
150
const
DipoleSplittingKernel
& split);
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
T
rootOfKallen
(T a, T b, T c)
const
{
168
if
( a*a + b*b + c*c - 2.*(a*b + a*c + b*c) >
ZERO
)
169
return
sqrt
(a*a + b*b + c*c - 2.*(a*b + a*c + b*c) ) ;
170
else
171
return
ZERO
; }
172
178
ThreeVector<double>
rotateToZ
(Lorentz5Momentum& pTarget, Lorentz5Momentum& p1){
179
ThreeVector<double>
oldAxis = pTarget.vect().
unit
();
180
double
ct = oldAxis.z();
double
st =
sqrt
( 1.-
sqr
(ct) );
// cos,sin(theta)
181
double
cp = oldAxis.x()/st;
double
sp = oldAxis.y()/st;
// cos,sin(phi)
182
pTarget.setZ( pTarget.vect().mag() ); pTarget.setX( 0.*GeV ); pTarget.setY( 0.*GeV );
183
Lorentz5Momentum p1old = p1;
184
p1.setX( sp*p1old.x() - cp*p1old.y() );
185
p1.setY( ct*cp*p1old.x() + ct*sp*p1old.y() - st*p1old.z() );
186
p1.setZ( st*cp*p1old.x() + st*sp*p1old.y() + ct*p1old.z() );
187
return
oldAxis;
188
}
189
190
public
:
191
198
void
persistentOutput
(
PersistentOStream
& os)
const
;
199
205
void
persistentInput
(
PersistentIStream
& is,
int
version);
207
214
static
void
Init
();
215
216
protected
:
217
224
virtual
IBPtr
clone
()
const
;
225
230
virtual
IBPtr
fullclone
()
const
;
232
233
234
// If needed, insert declarations of virtual function defined in the
235
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
236
237
238
private
:
239
244
static
ClassDescription<FFMassiveKinematics>
initFFMassiveKinematics
;
245
250
FFMassiveKinematics
&
operator=
(
const
FFMassiveKinematics
&) =
delete
;
251
252
};
253
254
}
255
256
#include "ThePEG/Utilities/ClassTraits.h"
257
258
namespace
ThePEG
{
259
264
template
<>
265
struct
BaseClassTrait
<
Herwig
::FFMassiveKinematics,1> {
267
typedef
Herwig::DipoleSplittingKinematics
NthBase
;
268
};
269
272
template
<>
273
struct
ClassTraits
<
Herwig
::FFMassiveKinematics>
274
:
public
ClassTraitsBase
<Herwig::FFMassiveKinematics> {
276
static
string
className
() {
return
"Herwig::FFMassiveKinematics"
; }
284
static
string
library
() {
return
"HwDipoleShower.so"
; }
285
};
286
289
}
290
291
#endif
/* HERWIG_FFMassiveKinematics_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::FFMassiveKinematics
FFMassiveKinematics implements massive splittings off a final-final dipole.
Definition:
FFMassiveKinematics.h:29
Herwig::FFMassiveKinematics::PtFromQ
virtual Energy PtFromQ(Energy scale, const DipoleSplittingInfo &) const
Return the pt given a virtuality.
Herwig::FFMassiveKinematics::persistentOutput
void persistentOutput(PersistentOStream &os) const
Function used to write out object persistently.
Herwig::FFMassiveKinematics::clone
virtual IBPtr clone() const
Make a simple clone of this object.
Herwig::FFMassiveKinematics::initFFMassiveKinematics
static ClassDescription< FFMassiveKinematics > initFFMassiveKinematics
The static object used to initialize the description of this class.
Definition:
FFMassiveKinematics.h:244
Herwig::FFMassiveKinematics::QFromPt
virtual Energy QFromPt(Energy scale, const DipoleSplittingInfo &) const
Return the virtuality given a pt.
Herwig::FFMassiveKinematics::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::FFMassiveKinematics::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::FFMassiveKinematics::generateSplitting
virtual bool generateSplitting(double kappa, double xi, double phi, DipoleSplittingInfo &dIndex, const DipoleSplittingKernel &split)
Generate splitting variables given three random numbers and the momentum fractions of the emitter and...
Herwig::FFMassiveKinematics::FFMassiveKinematics
FFMassiveKinematics()
The default constructor.
Herwig::FFMassiveKinematics::operator=
FFMassiveKinematics & operator=(const FFMassiveKinematics &)=delete
The assignment operator is private and must never be called.
Herwig::FFMassiveKinematics::zBoundaries
virtual pair< double, double > zBoundaries(Energy, const DipoleSplittingInfo &, const DipoleSplittingKernel &) const
Return the boundaries on the momentum fraction.
Definition:
FFMassiveKinematics.h:56
Herwig::FFMassiveKinematics::QMax
virtual Energy QMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum virtuality for the given dipole scale.
Definition:
FFMassiveKinematics.h:114
Herwig::FFMassiveKinematics::Init
static void Init()
The standard Init function used to initialize the interfaces.
Herwig::FFMassiveKinematics::ptToRandom
virtual double ptToRandom(Energy pt, Energy dScale, double emX, double specX, const DipoleIndex &dIndex, const DipoleSplittingKernel &split) const
Return the random number associated to the given pt.
Herwig::FFMassiveKinematics::persistentInput
void persistentInput(PersistentIStream &is, int version)
Function used to read in object persistently.
Herwig::FFMassiveKinematics::fullclone
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
Herwig::FFMassiveKinematics::rotateToZ
ThreeVector< double > rotateToZ(Lorentz5Momentum &pTarget, Lorentz5Momentum &p1)
Perform a rotation on both momenta such that the first one will point along the (positive) z axis.
Definition:
FFMassiveKinematics.h:178
Herwig::FFMassiveKinematics::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.
Herwig::FFMassiveKinematics::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::FFMassiveKinematics::rootOfKallen
T rootOfKallen(T a, T b, T c) const
Triangular / Kallen function.
Definition:
FFMassiveKinematics.h:167
Herwig::FFMassiveKinematics::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::FFMassiveKinematics::ptMax
virtual Energy ptMax(Energy dScale, double, double, const DipoleIndex &dIndex, const DipoleSplittingKernel &split, tPPtr emitter, tPPtr spectator) const
Return the maximum pt for the given dipole scale.
Herwig::FFMassiveKinematics::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::FFMassiveKinematics::ptMax
virtual Energy ptMax(Energy, double, double, const DipoleIndex &, const DipoleSplittingKernel &) const
Return the maximum pt for the given dipole scale.
Definition:
FFMassiveKinematics.h:93
ThePEG::ClassDescription
ThePEG::PersistentIStream
ThePEG::PersistentOStream
ThePEG::ThreeVector
ThePEG::ThreeVector::unit
ThreeVector< double > unit() const
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::sqr
constexpr auto sqr(const T &x) -> decltype(x *x)
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