herwig
is hosted by
Hepforge
,
IPPP Durham
Herwig
7.2.1
Decay
ScalarMeson
DtoKPiPiE691.h
1
// -*- C++ -*-
2
//
3
// DtoKPiPiE691.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_DtoKPiPiE691_H
10
#define HERWIG_DtoKPiPiE691_H
11
//
12
// This is the declaration of the DtoKPiPiE691 class.
13
//
14
15
#include "Herwig/Decay/DecayIntegrator.h"
16
#include "Herwig/Decay/DecayPhaseSpaceMode.h"
17
18
namespace
Herwig
{
19
20
using namespace
ThePEG
;
21
28
class
DtoKPiPiE691
:
public
DecayIntegrator
{
29
30
public
:
34
DtoKPiPiE691
();
35
42
virtual
int
modeNumber(
bool
& cc,
tcPDPtr
parent,
43
const
tPDVector
& children)
const
;
44
53
double
me2(
const
int
ichan,
const
Particle
& part,
54
const
ParticleVector
& decay,
MEOption
meopt)
const
;
55
61
virtual
void
dataBaseOutput(ofstream & os,
bool
header)
const
;
62
63
public
:
64
71
void
persistentOutput(
PersistentOStream
& os)
const
;
72
78
void
persistentInput(
PersistentIStream
& is,
int
version);
80
87
static
void
Init();
88
89
protected
:
90
103
double
decayAngle
(
const
Lorentz5Momentum & pparent,
104
const
Lorentz5Momentum & pres,
105
const
Lorentz5Momentum & p1)
const
{
106
Energy2 dot = pparent*p1, mREp = pres*pparent;
107
Energy2 mRE1 = pres*p1, mp2 = pparent.mass2();
108
Energy2 mres2 = pres.mass2(), m12 = p1.mass2();
109
return
(dot*mres2-mREp*mRE1)/
110
sqrt
((mREp*mREp-mres2*mp2)*(mRE1*mRE1-mres2*m12));
111
}
112
121
Complex
amplitude
(
int
ispin,
double
costheta,Energy mAB,
122
Energy wres, Energy mres)
const
{
123
double
s = 0.;
124
switch
(ispin) {
125
case
0: s = 1.;
break
;
126
case
1: s = costheta;
break
;
127
case
2: s = 1.5*
sqr
(costheta)-0.5;
break
;
128
default
: assert(
false
);
129
}
130
Complex
bw =
sqrt
(0.5*wres/GeV/Constants::pi)*GeV/
131
(mAB-mres-complex<Energy>(
ZERO
,0.5*wres));
132
return
s*bw;
133
}
135
136
protected
:
137
144
virtual
IBPtr
clone
()
const
{
return
new_ptr(*
this
);}
145
150
virtual
IBPtr
fullclone
()
const
{
return
new_ptr(*
this
);}
152
153
protected
:
154
162
virtual
void
doinit();
163
168
virtual
void
doinitrun();
170
171
private
:
172
177
DtoKPiPiE691
& operator=(
const
DtoKPiPiE691
&) =
delete
;
178
179
private
:
180
188
double
_a1NR
;
189
193
double
_phi1NR
;
194
198
double
_a1K892
;
199
203
double
_phi1K892
;
204
208
double
_a1K1430
;
209
213
double
_phi1K1430
;
214
218
double
_a1K1680
;
219
223
double
_phi1K1680
;
224
228
double
_a2NR
;
229
233
double
_phi2NR
;
234
238
double
_a2K8920
;
239
243
double
_phi2K8920
;
244
248
double
_a2K892m
;
249
253
double
_phi2K892m
;
254
258
double
_a2rho
;
259
263
double
_phi2rho
;
264
268
double
_a3NR
;
269
273
double
_phi3NR
;
274
278
double
_a3K892
;
279
283
double
_phi3K892
;
284
288
double
_a3rho
;
289
293
double
_phi3rho
;
295
303
Complex
_c1NR
;
304
308
Complex
_c1K892
;
309
313
Complex
_c1K1430
;
314
318
Complex
_c1K1680
;
319
323
Complex
_c2NR
;
324
328
Complex
_c2K8920
;
329
333
Complex
_c2K892m
;
334
338
Complex
_c2rho
;
339
343
Complex
_c3NR
;
344
348
Complex
_c3K892
;
349
353
Complex
_c3rho
;
355
363
bool
_localparameters
;
364
368
Energy
_mK8920
;
369
373
Energy
_wK8920
;
374
378
Energy
_mK892m
;
379
383
Energy
_wK892m
;
384
388
Energy
_mK1680
;
389
393
Energy
_wK1680
;
394
398
Energy
_mK1430
;
399
403
Energy
_wK1430
;
404
408
Energy
_mrho0
;
409
413
Energy
_wrho0
;
414
418
Energy
_mrhop
;
419
423
Energy
_wrhop
;
425
433
vector<double>
_maxwgt
;
434
438
vector<double>
_weights
;
440
444
mutable
RhoDMatrix
_rho
;
445
};
446
447
}
448
449
#endif
/* HERWIG_DtoKPiPiE691_H */
Herwig::DtoKPiPiE691::_a1K1430
double _a1K1430
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:208
ThePEG::PersistentIStream
Herwig::DtoKPiPiE691::_phi3K892
double _phi3K892
Phase of the component for .
Definition:
DtoKPiPiE691.h:283
Herwig::DtoKPiPiE691::_mrho0
Energy _mrho0
Mass of the .
Definition:
DtoKPiPiE691.h:408
Herwig::DtoKPiPiE691::_phi2NR
double _phi2NR
Phase of the non-resonant component for .
Definition:
DtoKPiPiE691.h:233
sqrt
double sqrt(int x)
Herwig::DtoKPiPiE691::_c1K1680
Complex _c1K1680
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:318
Herwig::DtoKPiPiE691::_c2NR
Complex _c2NR
Amplitude of the non-resonant component for .
Definition:
DtoKPiPiE691.h:323
ThePEG::Complex
std::complex< double > Complex
ThePEG::Particle
Herwig::DtoKPiPiE691::_c2K8920
Complex _c2K8920
Amplitude of the for .
Definition:
DtoKPiPiE691.h:328
Herwig::DtoKPiPiE691::_wrhop
Energy _wrhop
Width of the .
Definition:
DtoKPiPiE691.h:423
Herwig::DtoKPiPiE691::_mK8920
Energy _mK8920
Mass of the .
Definition:
DtoKPiPiE691.h:368
Herwig::DtoKPiPiE691
Here is the documentation of the DtoKPiPiE691 class.
Definition:
DtoKPiPiE691.h:28
Herwig::DtoKPiPiE691::_a1NR
double _a1NR
Amplitudes and phases for the different components.
Definition:
DtoKPiPiE691.h:188
Herwig::DtoKPiPiE691::_phi1K892
double _phi1K892
Phase of the component for .
Definition:
DtoKPiPiE691.h:203
Herwig::DtoKPiPiE691::_phi3rho
double _phi3rho
Phase of the component for .
Definition:
DtoKPiPiE691.h:293
Herwig::DtoKPiPiE691::_c1NR
Complex _c1NR
Complex amplitudes for use in the matrix element.
Definition:
DtoKPiPiE691.h:303
Herwig::DtoKPiPiE691::_a3K892
double _a3K892
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:278
ThePEG::PersistentOStream
ThePEG::IBPtr
ThePEG::Ptr< InterfacedBase >::pointer IBPtr
Herwig::DtoKPiPiE691::_localparameters
bool _localparameters
Masses and widths of the various resonances.
Definition:
DtoKPiPiE691.h:363
Herwig::DtoKPiPiE691::_c3NR
Complex _c3NR
Amplitude of the non-resonant component for .
Definition:
DtoKPiPiE691.h:343
sqr
constexpr auto sqr(const T &x) -> decltype(x *x)
Herwig::DtoKPiPiE691::_mrhop
Energy _mrhop
Mass of the .
Definition:
DtoKPiPiE691.h:418
ThePEG
Herwig::DtoKPiPiE691::_a2rho
double _a2rho
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:258
Herwig::DtoKPiPiE691::_wK1680
Energy _wK1680
Width of the .
Definition:
DtoKPiPiE691.h:393
Herwig::DtoKPiPiE691::decayAngle
double decayAngle(const Lorentz5Momentum &pparent, const Lorentz5Momentum &pres, const Lorentz5Momentum &p1) const
Methods to calculate the amplitudes for a given channel.
Definition:
DtoKPiPiE691.h:103
Herwig::DtoKPiPiE691::_phi3NR
double _phi3NR
Phase of the non-resonant component for .
Definition:
DtoKPiPiE691.h:273
Herwig::DtoKPiPiE691::_mK1430
Energy _mK1430
Mass of the .
Definition:
DtoKPiPiE691.h:398
Herwig::DtoKPiPiE691::_wK1430
Energy _wK1430
Width of the .
Definition:
DtoKPiPiE691.h:403
Herwig::DtoKPiPiE691::_a3rho
double _a3rho
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:288
Herwig::DecayIntegrator::MEOption
MEOption
Enum for the matrix element option.
Definition:
DecayIntegrator.h:79
Herwig::DtoKPiPiE691::_phi2rho
double _phi2rho
Phase of the component for .
Definition:
DtoKPiPiE691.h:263
Herwig::DtoKPiPiE691::_c1K892
Complex _c1K892
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:308
Herwig::DtoKPiPiE691::amplitude
Complex amplitude(int ispin, double costheta, Energy mAB, Energy wres, Energy mres) const
Calculate the amplitude.
Definition:
DtoKPiPiE691.h:121
ThePEG::tPDVector
vector< tPDPtr > tPDVector
Herwig::DtoKPiPiE691::_weights
vector< double > _weights
Weights for the different integration channels.
Definition:
DtoKPiPiE691.h:438
Herwig::DecayIntegrator
Main class for Decayers implementing multi-channel phase space integration.
Definition:
DecayIntegrator.h:62
ThePEG::tcPDPtr
ThePEG::Ptr< ParticleData >::transient_const_pointer tcPDPtr
Herwig::DtoKPiPiE691::_mK892m
Energy _mK892m
Mass of the .
Definition:
DtoKPiPiE691.h:378
Herwig::DtoKPiPiE691::_wrho0
Energy _wrho0
Width of the .
Definition:
DtoKPiPiE691.h:413
Herwig::DtoKPiPiE691::_c1K1430
Complex _c1K1430
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:313
Herwig::DtoKPiPiE691::_phi2K892m
double _phi2K892m
Phase of the component for .
Definition:
DtoKPiPiE691.h:253
Herwig::DtoKPiPiE691::_a2NR
double _a2NR
Amplitude of the non-resonant component for .
Definition:
DtoKPiPiE691.h:228
Herwig::DtoKPiPiE691::_a2K892m
double _a2K892m
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:248
Herwig::DtoKPiPiE691::_phi1K1680
double _phi1K1680
Phase of the component for .
Definition:
DtoKPiPiE691.h:223
Herwig::DtoKPiPiE691::_mK1680
Energy _mK1680
Mass of the .
Definition:
DtoKPiPiE691.h:388
Herwig::DtoKPiPiE691::_c3K892
Complex _c3K892
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:348
Herwig::DtoKPiPiE691::_wK8920
Energy _wK8920
Width of the .
Definition:
DtoKPiPiE691.h:373
Herwig::DtoKPiPiE691::_c2K892m
Complex _c2K892m
Amplitude of the for .
Definition:
DtoKPiPiE691.h:333
Herwig::DtoKPiPiE691::clone
virtual IBPtr clone() const
Make a simple clone of this object.
Definition:
DtoKPiPiE691.h:144
Herwig::DtoKPiPiE691::_phi1NR
double _phi1NR
Phase of the non-resonant component for .
Definition:
DtoKPiPiE691.h:193
Herwig::DtoKPiPiE691::_c2rho
Complex _c2rho
Amplitude of the for .
Definition:
DtoKPiPiE691.h:338
Herwig::DtoKPiPiE691::_phi2K8920
double _phi2K8920
Phase of the component for .
Definition:
DtoKPiPiE691.h:243
Herwig
-*- C++ -*-
Definition:
BasicConsistency.h:17
Herwig::DtoKPiPiE691::_a1K892
double _a1K892
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:198
Herwig::DtoKPiPiE691::_c3rho
Complex _c3rho
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:353
Herwig::DtoKPiPiE691::_phi1K1430
double _phi1K1430
Phase of the component for .
Definition:
DtoKPiPiE691.h:213
Herwig::DtoKPiPiE691::_maxwgt
vector< double > _maxwgt
Parameters for the phase-space integration.
Definition:
DtoKPiPiE691.h:433
ThePEG::ParticleVector
vector< PPtr > ParticleVector
ThePEG::RhoDMatrix
Herwig::DtoKPiPiE691::_a2K8920
double _a2K8920
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:238
Herwig::DtoKPiPiE691::fullclone
virtual IBPtr fullclone() const
Make a clone of this object, possibly modifying the cloned object to make it sane.
Definition:
DtoKPiPiE691.h:150
Herwig::DtoKPiPiE691::_wK892m
Energy _wK892m
Width of the .
Definition:
DtoKPiPiE691.h:383
ZERO
constexpr ZeroUnit ZERO
Herwig::DtoKPiPiE691::_a1K1680
double _a1K1680
Amplitude of the component for .
Definition:
DtoKPiPiE691.h:218
Herwig::DtoKPiPiE691::_a3NR
double _a3NR
Amplitude of the non-resonant component for .
Definition:
DtoKPiPiE691.h:268
Herwig::DtoKPiPiE691::_rho
RhoDMatrix _rho
Spin density matrix.
Definition:
DtoKPiPiE691.h:444
Generated on Sat Apr 11 2020 14:50:29 for Herwig by
1.8.13