Go to the source code of this file.
|
double | EvtDecayAngle (const EvtVector4R &, const EvtVector4R &, const EvtVector4R &) |
|
double | EvtDecayAngleChi (const EvtVector4R &, const EvtVector4R &, const EvtVector4R &, const EvtVector4R &, const EvtVector4R &) |
|
double | EvtDecayPlaneNormalAngle (const EvtVector4R &p, const EvtVector4R &q, const EvtVector4R &d1, const EvtVector4R &d2) |
|
double | EvtDecayAnglePhi (const EvtVector4R &g, const EvtVector4R &p, const EvtVector4R &q, const EvtVector4R &d) |
|
EvtComplex | wignerD (int j, int m1, int m2, double phi, double theta, double gamma) |
|
◆ EvtDecayAngle()
Definition at line 33 of file EvtKine.cc.
34 {
35
36 double pd=p*d;
41 double md2=d.mass2();
42
43 double cost=(pd*mq2-pq*qd)/sqrt((pq*pq-mq2*mp2)*(qd*qd-mq2*md2));
44
45 return cost;
46
47}
****INTEGER imax DOUBLE PRECISION m_pi *DOUBLE PRECISION m_amfin DOUBLE PRECISION m_Chfin DOUBLE PRECISION m_Xenph DOUBLE PRECISION m_sinw2 DOUBLE PRECISION m_GFermi DOUBLE PRECISION m_MfinMin DOUBLE PRECISION m_ta2 INTEGER m_out INTEGER m_KeyFSR INTEGER m_KeyQCD *COMMON c_Semalib $ !copy of input $ !CMS energy $ !beam mass $ !final mass $ !beam charge $ !final charge $ !smallest final mass $ !Z mass $ !Z width $ !EW mixing angle $ !Gmu Fermi $ alphaQED at q
Referenced by EvtMHelAmp::amplitude(), and EvtResonance::resAmpl().
◆ EvtDecayAngleChi()
Definition at line 49 of file EvtKine.cc.
51 {
52
54
55
56
57
62
63
66
67 D=p4_d1p+p4_d2p;
68
69 d1_perp=p4_d1p-(D.
dot(p4_d1p)/D.
dot(D))*D;
70 h1_perp=p4_h1p-(D.
dot(p4_h1p)/D.
dot(D))*D;
71
72
73
74 d1_prime=D.
cross(d1_perp);
75
76 d1_perp= d1_perp/d1_perp.
d3mag();
77 d1_prime= d1_prime/d1_prime.
d3mag();
78
80
81 x=d1_perp.
dot(h1_perp);
82 y=d1_prime.
dot(h1_perp);
83
84 double chi=atan2(
y,x);
85
87
88 return chi;
89
90}
EvtDiracSpinor boostTo(const EvtDiracSpinor &sp, const EvtVector4R p4)
static const double twoPi
double dot(const EvtVector4R &v2) const
EvtVector4R cross(const EvtVector4R &v2)
◆ EvtDecayAnglePhi()
Definition at line 111 of file EvtKine.cc.
113{
114 double eq = (p *
q) / p.
mass();
115 double ed = (p * d) / p.
mass();
116 double mq =
q.mass();
120 double zd = p.
dotr3(z,d);
121 double alpha = (eq - mq)/(q2 * mq) * qd - ed/mq;
122
124 double x = (zq * (qd +
alpha * q2) - q2 * (zd +
alpha * zq)) / sqrt(q2);
125
126 double phi = atan2(
y,x);
127
129}
double dotr3(const EvtVector4R &p1, const EvtVector4R &p2) const
double scalartripler3(const EvtVector4R &p1, const EvtVector4R &p2, const EvtVector4R &p3) const
double mag2r3(const EvtVector4R &p1) const
Referenced by EvtMHelAmp::amplitude().
◆ EvtDecayPlaneNormalAngle()
Definition at line 94 of file EvtKine.cc.
95 {
96
98
101
103
104 return q.
mass()*(p*l)/sqrt(-(pq*pq-p.
mass2()*
q.mass2())*l.mass2());
105
106
107}
Evt3Rank3C directProd(const EvtVector3C &c1, const EvtVector3C &c2, const EvtVector3C &c3)
EvtTensor4C dual(const EvtTensor4C &t2)
EvtVector4C cont2(const EvtVector4C &v4) const
const EvtComplex & get(int) const
◆ wignerD()
EvtComplex wignerD |
( |
int | j, |
|
|
int | m1, |
|
|
int | m2, |
|
|
double | phi, |
|
|
double | theta, |
|
|
double | gamma ) |
Definition at line 132 of file EvtKine.cc.
134{
137
138
140}
EvtComplex exp(const EvtComplex &c)
static double d(int j, int m1, int m2, double theta)
Referenced by EvtMHelAmp::amplitude().