BOSS 7.0.9
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtDalitzPoint Class Reference

#include <EvtDalitzPoint.hh>

Public Member Functions

 EvtDalitzPoint ()
 
 EvtDalitzPoint (double mA, double mB, double mC, double qAB, double qBC, double qCA)
 
 EvtDalitzPoint (double mA, double mB, double mC, EvtCyclic3::Pair i, double qres, double qhel, double qsum)
 
 EvtDalitzPoint (const EvtDalitzPlot &, const EvtDalitzCoord &)
 
 EvtDalitzPoint (const EvtDalitzPoint &other)
 
 ~EvtDalitzPoint ()
 
EvtDalitzCoord getDalitzPoint (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
EvtDalitzPlot getDalitzPlot () const
 
double q (EvtCyclic3::Pair) const
 
double bigM () const
 
double m (EvtCyclic3::Index) const
 
double qres (EvtCyclic3::Pair i) const
 
double qhel (EvtCyclic3::Pair i) const
 
double qsum () const
 
double qMin (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
double qMax (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
double pp (EvtCyclic3::Index i, EvtCyclic3::Index j) const
 
double e (EvtCyclic3::Index i, EvtCyclic3::Pair j) const
 
double p (EvtCyclic3::Index i, EvtCyclic3::Pair j) const
 
double cosTh (EvtCyclic3::Pair pairAng, EvtCyclic3::Pair pairRes) const
 
bool isValid () const
 
void print () const
 
 EvtDalitzPoint ()
 
 EvtDalitzPoint (double mA, double mB, double mC, double qAB, double qBC, double qCA)
 
 EvtDalitzPoint (double mA, double mB, double mC, EvtCyclic3::Pair i, double qres, double qhel, double qsum)
 
 EvtDalitzPoint (const EvtDalitzPlot &, const EvtDalitzCoord &)
 
 EvtDalitzPoint (const EvtDalitzPoint &other)
 
 ~EvtDalitzPoint ()
 
EvtDalitzCoord getDalitzPoint (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
EvtDalitzPlot getDalitzPlot () const
 
double q (EvtCyclic3::Pair) const
 
double bigM () const
 
double m (EvtCyclic3::Index) const
 
double qres (EvtCyclic3::Pair i) const
 
double qhel (EvtCyclic3::Pair i) const
 
double qsum () const
 
double qMin (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
double qMax (EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
 
double pp (EvtCyclic3::Index i, EvtCyclic3::Index j) const
 
double e (EvtCyclic3::Index i, EvtCyclic3::Pair j) const
 
double p (EvtCyclic3::Index i, EvtCyclic3::Pair j) const
 
double cosTh (EvtCyclic3::Pair pairAng, EvtCyclic3::Pair pairRes) const
 
bool isValid () const
 
void print () const
 

Detailed Description

Constructor & Destructor Documentation

◆ EvtDalitzPoint() [1/10]

EvtDalitzPoint::EvtDalitzPoint ( )

Definition at line 18 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

19 : _mA(-1.), _mB(-1.), _mC(-1.), _qAB(-1.), _qBC(-1.), _qCA(-1.)
20{}

◆ EvtDalitzPoint() [2/10]

EvtDalitzPoint::EvtDalitzPoint ( double  mA,
double  mB,
double  mC,
double  qAB,
double  qBC,
double  qCA 
)

Definition at line 22 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

23 : _mA(mA), _mB(mB), _mC(mC), _qAB(qAB), _qBC(qBC), _qCA(qCA)
24{}

◆ EvtDalitzPoint() [3/10]

EvtDalitzPoint::EvtDalitzPoint ( double  mA,
double  mB,
double  mC,
EvtCyclic3::Pair  i,
double  qres,
double  qhel,
double  qsum 
)

Definition at line 28 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

31 : _mA(mA), _mB(mB), _mC(mC)
32{
33 double qi = qres + qsum/3.;
34 double qj = -qres/2. + qhel + qsum/3.;
35 double qk = -qres/2. - qhel + qsum/3.;
36
37 if(i == AB) { _qAB = qi; _qBC = qj; _qCA = qk; }
38 else if(i == BC) { _qAB = qk; _qBC = qi; _qCA = qj; }
39 else if(i == CA) { _qAB = qj; _qBC = qk; _qCA = qi; }
40}

◆ EvtDalitzPoint() [4/10]

EvtDalitzPoint::EvtDalitzPoint ( const EvtDalitzPlot dp,
const EvtDalitzCoord x 
)

Definition at line 42 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

43 : _mA(dp.m(A)), _mB(dp.m(B)), _mC(dp.m(C))
44{
45 if(x.pair1() == AB) _qAB = x.q1();
46 else
47 if(x.pair2() == AB) _qAB = x.q2();
48 else _qAB = dp.sum() - x.q1() - x.q2();
49
50 if(x.pair1() == BC) _qBC = x.q1();
51 else
52 if(x.pair2() == BC) _qBC = x.q2();
53 else _qBC = dp.sum() - x.q1() - x.q2();
54
55 if(x.pair1() == CA) _qCA = x.q1();
56 else
57 if(x.pair2() == CA) _qCA = x.q2();
58 else _qCA = dp.sum() - x.q1() - x.q2();
59
60}
Double_t x[10]

◆ EvtDalitzPoint() [5/10]

EvtDalitzPoint::EvtDalitzPoint ( const EvtDalitzPoint other)

Definition at line 62 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

63 : _mA(other._mA), _mB(other._mB), _mC(other._mC),
64 _qAB(other._qAB), _qBC(other._qBC), _qCA(other._qCA)
65{}

◆ ~EvtDalitzPoint() [1/2]

EvtDalitzPoint::~EvtDalitzPoint ( )

◆ EvtDalitzPoint() [6/10]

EvtDalitzPoint::EvtDalitzPoint ( )

◆ EvtDalitzPoint() [7/10]

EvtDalitzPoint::EvtDalitzPoint ( double  mA,
double  mB,
double  mC,
double  qAB,
double  qBC,
double  qCA 
)

◆ EvtDalitzPoint() [8/10]

EvtDalitzPoint::EvtDalitzPoint ( double  mA,
double  mB,
double  mC,
EvtCyclic3::Pair  i,
double  qres,
double  qhel,
double  qsum 
)

◆ EvtDalitzPoint() [9/10]

EvtDalitzPoint::EvtDalitzPoint ( const EvtDalitzPlot ,
const EvtDalitzCoord  
)

◆ EvtDalitzPoint() [10/10]

EvtDalitzPoint::EvtDalitzPoint ( const EvtDalitzPoint other)

◆ ~EvtDalitzPoint() [2/2]

EvtDalitzPoint::~EvtDalitzPoint ( )

Member Function Documentation

◆ bigM() [1/2]

double EvtDalitzPoint::bigM ( ) const

Definition at line 176 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

177{
178 return sqrt(_qAB+_qBC+_qCA - _mA*_mA - _mB*_mB - _mC*_mC);
179}

Referenced by getDalitzPlot(), and isValid().

◆ bigM() [2/2]

double EvtDalitzPoint::bigM ( ) const

◆ cosTh() [1/2]

double EvtDalitzPoint::cosTh ( EvtCyclic3::Pair  pairAng,
EvtCyclic3::Pair  pairRes 
) const

Definition at line 138 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

139{
141 return dp.cosTh(pairAng,q(pairAng),pairRes,q(pairRes));
142}
****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
Definition: KKsem.h:33
double cosTh(EvtCyclic3::Pair i1, double q1, EvtCyclic3::Pair i2, double q2) const

Referenced by EvtPto3P::initDaughters().

◆ cosTh() [2/2]

double EvtDalitzPoint::cosTh ( EvtCyclic3::Pair  pairAng,
EvtCyclic3::Pair  pairRes 
) const

◆ e() [1/2]

double EvtDalitzPoint::e ( EvtCyclic3::Index  i,
EvtCyclic3::Pair  j 
) const

Definition at line 126 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

127{
129 return dp.e(i,j,q(j));
130}
double e(EvtCyclic3::Index i, EvtCyclic3::Pair j, double q) const

Referenced by EvtPto3P::initDaughters().

◆ e() [2/2]

double EvtDalitzPoint::e ( EvtCyclic3::Index  i,
EvtCyclic3::Pair  j 
) const

◆ getDalitzPlot() [1/2]

EvtDalitzPlot EvtDalitzPoint::getDalitzPlot ( ) const

◆ getDalitzPlot() [2/2]

EvtDalitzPlot EvtDalitzPoint::getDalitzPlot ( ) const

◆ getDalitzPoint() [1/2]

◆ getDalitzPoint() [2/2]

EvtDalitzCoord EvtDalitzPoint::getDalitzPoint ( EvtCyclic3::Pair  i,
EvtCyclic3::Pair  j 
) const

◆ isValid() [1/2]

bool EvtDalitzPoint::isValid ( ) const

Definition at line 156 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

157{
158 // Check masses
159
160 double M = bigM();
161 if(_mA < 0 || _mB < 0 || _mC < 0 || M <= 0) return false;
162 if(M < _mA + _mB + _mC) return false;
163
164 // Check that first coordinate is within absolute limits
165
166 bool inside = false;
168
169 if(dp.qAbsMin(AB) <= _qAB && _qAB <= dp.qAbsMax(AB))
170 if(qMin(BC,AB) <= _qBC && _qBC <= qMax(BC,AB))
171 inside = true;
172
173 return inside;
174}
double qMin(EvtCyclic3::Pair i, EvtCyclic3::Pair j) const
double qMax(EvtCyclic3::Pair i, EvtCyclic3::Pair j) const

Referenced by EvtPto3P::initDaughters(), and EvtDalitzResPdf::randomPoint().

◆ isValid() [2/2]

bool EvtDalitzPoint::isValid ( ) const

◆ m() [1/2]

double EvtDalitzPoint::m ( EvtCyclic3::Index  i) const

Definition at line 80 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

81{
82 double ret = _mA;
83 if(B == i) ret = _mB;
84 else
85 if(C == i) ret = _mC;
86
87 return ret;
88}

Referenced by pp().

◆ m() [2/2]

double EvtDalitzPoint::m ( EvtCyclic3::Index  ) const

◆ p() [1/2]

double EvtDalitzPoint::p ( EvtCyclic3::Index  i,
EvtCyclic3::Pair  j 
) const

Definition at line 132 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

133{
135 return dp.p(i,j,q(j));
136}
double p(EvtCyclic3::Index i, EvtCyclic3::Pair j, double q) const

Referenced by EvtPto3P::initDaughters().

◆ p() [2/2]

double EvtDalitzPoint::p ( EvtCyclic3::Index  i,
EvtCyclic3::Pair  j 
) const

◆ pp() [1/2]

double EvtDalitzPoint::pp ( EvtCyclic3::Index  i,
EvtCyclic3::Index  j 
) const

Definition at line 120 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

121{
122 if(i == j) return m(i)*m(i);
123 else return (q(combine(i,j)) - m(i)*m(i) - m(j)*m(j))/2.;
124}

◆ pp() [2/2]

double EvtDalitzPoint::pp ( EvtCyclic3::Index  i,
EvtCyclic3::Index  j 
) const

◆ print() [1/2]

void EvtDalitzPoint::print ( ) const

Definition at line 182 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

183{
185 printf("%f %f %f\n",_qAB,_qBC,_qCA);
186}

◆ print() [2/2]

void EvtDalitzPoint::print ( ) const

◆ q() [1/2]

double EvtDalitzPoint::q ( EvtCyclic3::Pair  i) const

Definition at line 70 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

71{
72 double ret = _qAB;
73 if(BC == i) ret = _qBC;
74 else
75 if(CA == i) ret = _qCA;
76
77 return ret;
78}

◆ q() [2/2]

double EvtDalitzPoint::q ( EvtCyclic3::Pair  ) const

◆ qhel() [1/2]

double EvtDalitzPoint::qhel ( EvtCyclic3::Pair  i) const

◆ qhel() [2/2]

double EvtDalitzPoint::qhel ( EvtCyclic3::Pair  i) const

◆ qMax() [1/2]

double EvtDalitzPoint::qMax ( EvtCyclic3::Pair  i,
EvtCyclic3::Pair  j 
) const

Definition at line 114 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

115{
117 return dp.qMax(i,j,q(j));
118}
double qMax(EvtCyclic3::Pair i, EvtCyclic3::Pair j, double q) const

Referenced by isValid().

◆ qMax() [2/2]

double EvtDalitzPoint::qMax ( EvtCyclic3::Pair  i,
EvtCyclic3::Pair  j 
) const

◆ qMin() [1/2]

double EvtDalitzPoint::qMin ( EvtCyclic3::Pair  i,
EvtCyclic3::Pair  j 
) const

Definition at line 108 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

109{
111 return dp.qMin(i,j,q(j));
112}
double qMin(EvtCyclic3::Pair i, EvtCyclic3::Pair j, double q) const

Referenced by isValid().

◆ qMin() [2/2]

double EvtDalitzPoint::qMin ( EvtCyclic3::Pair  i,
EvtCyclic3::Pair  j 
) const

◆ qres() [1/2]

double EvtDalitzPoint::qres ( EvtCyclic3::Pair  i) const

Definition at line 92 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

93{
94 return (2.*q(i) - q(EvtCyclic3::prev(i)) - q(EvtCyclic3::next(i)))/3.;
95}

Referenced by EvtDalitzPoint().

◆ qres() [2/2]

double EvtDalitzPoint::qres ( EvtCyclic3::Pair  i) const

◆ qsum() [1/2]

double EvtDalitzPoint::qsum ( ) const

Definition at line 102 of file bak-BesEvtGen-00-04-08/src/EvtGen/EvtGenBase/EvtDalitzPoint.cc.

103{
104 return _qAB + _qBC + _qCA;
105}

Referenced by EvtDalitzPoint().

◆ qsum() [2/2]

double EvtDalitzPoint::qsum ( ) const

The documentation for this class was generated from the following files: