BOSS 8.0.0
BESIII Offline Software System
Loading...
Searching...
No Matches
D0ToKSpipi2018 Class Reference

#include <D0ToKSpipi2018.h>

Public Member Functions

 D0ToKSpipi2018 ()
virtual ~D0ToKSpipi2018 ()
void init ()
complex< double > Amp_PFT (vector< double > k0l, vector< double > pip, vector< double > pim)

Detailed Description

Definition at line 9 of file D0ToKSpipi2018.h.

Constructor & Destructor Documentation

◆ D0ToKSpipi2018()

D0ToKSpipi2018::D0ToKSpipi2018 ( )
inline

Definition at line 12 of file D0ToKSpipi2018.h.

12{}

◆ ~D0ToKSpipi2018()

D0ToKSpipi2018::~D0ToKSpipi2018 ( )
virtual

Definition at line 23 of file D0ToKSpipi2018.cxx.

23{}

Member Function Documentation

◆ Amp_PFT()

complex< double > D0ToKSpipi2018::Amp_PFT ( vector< double > k0l,
vector< double > pip,
vector< double > pim )

Definition at line 177 of file D0ToKSpipi2018.cxx.

178 {
179 // Breit-Wigner lineshapes
180 vector<double> pD;
181 pD.clear();
182 if ( k0l.size() != 4 || pip.size() != 4 || pim.size() != 4 )
183 cout << "ERROR in KSPIPI daughter 4 momentum" << endl;
184 for ( int i = 0; i < k0l.size(); i++ ) { pD.push_back( k0l[i] + pip[i] + pim[i] ); }
185
186 complex<double> DK2piRes0 = Resonance2( pD, pip, pim, ar[0], phir[0], width_R[0], mass_R[0],
187 spin_R[0] ); // ar, phir, width, mass, spin Rho770
188 complex<double> DK2piRes1 = Resonance2( pD, pip, pim, ar[1], phir[1], width_R[1], mass_R[1],
189 spin_R[1] ); // ar, phir, width, mass, spin Omega782
190 complex<double> DK2piRes2 = Resonance2( pD, pip, pim, ar[2], phir[2], width_R[2], mass_R[2],
191 spin_R[2] ); // ar, phir, width, mass, spin ftwo1270
192 complex<double> DK2piRes3 = Resonance2( pD, pip, pim, ar[3], phir[3], width_R[3], mass_R[3],
193 spin_R[3] ); // ar, phir, width, mass, spin Rho1450
194 complex<double> DK2piRes4 = Resonance2( pD, k0l, pim, ar[4], phir[4], width_R[4], mass_R[4],
195 spin_R[4] ); // ar, phir, width, mass, spin Kstar892-
196 complex<double> DK2piRes5 =
197 Resonance2( pD, k0l, pim, ar[5], phir[5], width_R[5], mass_R[5],
198 spin_R[5] ); // ar, phir, width, mass, spin K2star1430-
199 complex<double> DK2piRes6 =
200 Resonance2( pD, k0l, pim, ar[6], phir[6], width_R[6], mass_R[6],
201 spin_R[6] ); // ar, phir, width, mass, spin Kstar1680-
202 complex<double> DK2piRes7 =
203 Resonance2( pD, k0l, pim, ar[7], phir[7], width_R[7], mass_R[7],
204 spin_R[7] ); // ar, phir, width, mass, spin Kstar1410-
205 complex<double> DK2piRes8 = Resonance2( pD, k0l, pip, ar[8], phir[8], width_R[8], mass_R[8],
206 spin_R[8] ); // ar, phir, width, mass, spin Kstar892+
207 complex<double> DK2piRes9 =
208 Resonance2( pD, k0l, pip, ar[9], phir[9], width_R[9], mass_R[9],
209 spin_R[9] ); // ar, phir, width, mass, spin K2star1430+
210 complex<double> DK2piRes10 =
211 Resonance2( pD, k0l, pip, ar[10], phir[10], width_R[10], mass_R[10],
212 spin_R[10] ); // ar, phir, width, mass, spin Kstar1410+
213 // K-matrix for pipi S-wave
214 complex<double> pipi_s_wave = K_matrix( pip, pim );
215 if ( pipi_s_wave == complex<double>( 9999., 9999. ) ) return 1e-20;
216 // LASS parametrization for Kpi S-wave
217 complex<double> kpi_s_wave =
218 amplitude_LASS( k0l, pip, pim, "k0spim", ar[11], phir[11] * pi180inv );
219 complex<double> dcs_kpi_s_wave =
220 amplitude_LASS( k0l, pip, pim, "k0spip", ar[12], phir[12] * pi180inv );
221
222 complex<double> _tmpAmp = DK2piRes0 + DK2piRes1 + DK2piRes2 + DK2piRes3 + pipi_s_wave;
223 // complex<double> TOT_PFT_AMP = DK2piRes0+ DK2piRes1+ DK2piRes2+ DK2piRes3+ DK2piRes4+
224 // DK2piRes5+ DK2piRes6+ DK2piRes7+ DK2piRes8+ DK2piRes9+ DK2piRes10+ pipi_s_wave +
225 // kpi_s_wave+ dcs_kpi_s_wave ;
226 complex<double> TOT_PFT_AMP = _tmpAmp + DK2piRes4 + DK2piRes5 + DK2piRes6 + DK2piRes7 +
227 DK2piRes8 + DK2piRes9 + DK2piRes10 + kpi_s_wave +
228 dcs_kpi_s_wave;
229 // Coherent sum for pure-flavor-tagged amplitudes (PFT)
230 return TOT_PFT_AMP;
231}

◆ init()

void D0ToKSpipi2018::init ( )

Definition at line 25 of file D0ToKSpipi2018.cxx.

25 {
26 // std::cout << "D0ToKSpipi2018 ==> Initialization !" << std::endl;
27
28 _nd = 3;
29 tan2thetaC = ( 0.22650 * 0.22650 ) /
30 ( 1. - ( 0.22650 * 0.22650 ) ); // sin(theta_C) = 0.22650 +/- 0.00048
31 pi180inv = 1.0 * 3.1415926 / 180;
32 double Pi = 3.1415926;
33
34 mass_R[0] = 0.77155;
35 width_R[0] = 0.13469;
36 spin_R[0] = 1;
37 ar[0] = 1;
38 phir[0] = 0;
39 mass_R[1] = 0.78265;
40 width_R[1] = 0.00849;
41 spin_R[1] = 1;
42 ar[1] = 0.038791;
43 phir[1] = ( 180. / Pi ) * 2.1073;
44 mass_R[2] = 1.27510;
45 width_R[2] = 0.18420;
46 spin_R[2] = 2;
47 ar[2] = 1.42887;
48 phir[2] = ( 180. / Pi ) * -0.633296;
49 mass_R[3] = 1.46500;
50 width_R[3] = 0.40000;
51 spin_R[3] = 1;
52 ar[3] = 2.85131;
53 phir[3] = ( 180. / Pi ) * 1.7820801;
54 mass_R[4] = 0.89371;
55 width_R[4] = 0.04719;
56 spin_R[4] = 1;
57 ar[4] = 1.72044;
58 phir[4] = ( 180. / Pi ) * 2.38835877;
59 mass_R[5] = 1.42560;
60 width_R[5] = 0.09850;
61 spin_R[5] = 2;
62 ar[5] = 1.27268;
63 phir[5] = ( 180. / Pi ) * -0.769095;
64 mass_R[6] = 1.71700;
65 width_R[6] = 0.3220;
66 spin_R[6] = 1;
67 ar[6] = 3.307642;
68 phir[6] = ( 180. / Pi ) * -2.062227;
69 mass_R[7] = 1.41400;
70 width_R[7] = 0.2320;
71 spin_R[7] = 1;
72 ar[7] = 0.286927;
73 phir[7] = ( 180. / Pi ) * 1.7346186;
74 mass_R[8] = 0.89371;
75 width_R[8] = 0.04719;
76 spin_R[8] = 1;
77 ar[8] = 0.1641792;
78 phir[8] = ( 180. / Pi ) * -0.735903;
79 mass_R[9] = 1.42560;
80 width_R[9] = 0.0985;
81 spin_R[9] = 2;
82 ar[9] = 0.1025736;
83 phir[9] = ( 180. / Pi ) * -1.56397;
84 mass_R[10] = 1.41400;
85 width_R[10] = 0.2320;
86 spin_R[10] = 1;
87 ar[10] = 0.2090326;
88 phir[10] = ( 180. / Pi ) * 2.6208986;
89 ////////////////////////////////////////
90 mass_R[11] = 1.42500;
91 width_R[11] = 0.2700;
92 spin_R[11] = 1;
93 ar[11] = 2.36;
94 phir[11] = 99.4; // not found
95 mass_R[12] = 1.42500;
96 width_R[12] = 0.2700;
97 spin_R[12] = 1;
98 ar[12] = 0.11267;
99 phir[12] = -162.3; // rad
100
101 beta[0] = complex<double>( 8.521486 * cos( 1.195641 ), 8.521486 * sin( 1.195641 ) ); //
102 beta[1] = complex<double>( 12.1895 * cos( 0.41802 ), 12.1895 * sin( 0.41802 ) );
103 beta[2] = complex<double>( 29.14616 * cos( -0.0018386 ), 29.14616 * sin( -0.0018386 ) );
104 beta[3] = complex<double>( 10.745735 * cos( -0.9057014 ), 10.745735 * sin( -0.9057014 ) );
105 beta[4] = complex<double>( 0., 0. );
106
107 fprod[0] = complex<double>( 8.04427 * cos( -2.19847 ), 8.04427 * sin( -2.19847 ) );
108 fprod[1] = complex<double>( 26.2986 * cos( -2.65853 ), 26.2986 * sin( -2.65853 ) );
109 fprod[2] = complex<double>( 33.0349 * cos( -1.62714 ), 33.0349 * sin( -1.62714 ) );
110 fprod[3] = complex<double>( 26.1741 * cos( -2.11891 ), 26.1741 * sin( -2.11891 ) );
111 fprod[4] = complex<double>( 0., 0. );
112 // beta.push_back( complex<double>( 0.255303*cos( 47.8861 *pi180inv), 0.255303*sin( 47.8861
113 // *pi180inv)) ); beta.push_back( complex<double>(13.4446 *cos( -5.11127*pi180inv), 13.4446
114 // *sin( -5.11127*pi180inv)) ); beta.push_back( complex<double>(38.8496 *cos(-30.06
115 // *pi180inv), 38.8496 *sin(-30.06 *pi180inv)) ); beta.push_back( complex<double>(13.1086
116 // *cos(-81.4148 *pi180inv), 13.1086 *sin(-81.4148 *pi180inv)) ); beta.push_back(
117 // complex<double>(0., 0.) );
118
119 // fprod.push_back(
120 // complex<double>(5.08049*cos(-182.312*pi180inv), 5.08049*sin(-182.312*pi180inv)));
121 // fprod.push_back(
122 // complex<double>(17.2388*cos(-219.209*pi180inv), 17.2388*sin(-219.209*pi180inv)));
123 // fprod.push_back(
124 // complex<double>(19.0145*cos(-76.9884*pi180inv), 19.0145*sin(-76.9884*pi180inv)));
125 // fprod.push_back(
126 // complex<double>(11.9875*cos(-190.502*pi180inv), 11.9875*sin(-190.502*pi180inv)));
127 // fprod.push_back( complex<double>(0., 0.));
128
129 ma[0] = 0.651;
130 g[0][0] = 0.22889;
131 g[0][1] = -0.55377;
132 g[0][2] = 0;
133 g[0][3] = -0.39899;
134 g[0][4] = -0.34639;
135 ma[1] = 1.20360;
136 g[1][0] = 0.94128;
137 g[1][1] = 0.55095;
138 g[1][2] = 0;
139 g[1][3] = 0.39065;
140 g[1][4] = 0.31503;
141 ma[2] = 1.55817;
142 g[2][0] = 0.36856;
143 g[2][1] = 0.23888;
144 g[2][2] = 0.55639;
145 g[2][3] = 0.18340;
146 g[2][4] = 0.18681;
147 ma[3] = 1.21000;
148 g[3][0] = 0.33650;
149 g[3][1] = 0.40907;
150 g[3][2] = 0.85679;
151 g[3][3] = 0.19906;
152 g[3][4] = -0.00984;
153 ma[4] = 1.82206;
154 g[4][0] = 0.18171;
155 g[4][1] = -0.17558;
156 g[4][2] = -0.79658;
157 g[4][3] = -0.00355;
158 g[4][4] = 0.22358;
159
160 // Hadronic parameters for tag modes: 0=no-specified, 1=Kpi, 2=Kpipi0, 3=K3pi
161 rd[0] = 0.0;
162 rd[1] = 0.0586;
163 rd[2] = 0.0440;
164 rd[3] = 0.0546;
165 deltad[0] = 0.0;
166 deltad[1] = 194.7 * pi180inv;
167 deltad[2] = 196.0 * pi180inv;
168 deltad[3] = 167.0 * pi180inv;
169 Rf[0] = 0.0;
170 Rf[1] = 1.0;
171 Rf[2] = 0.78;
172 Rf[3] = 0.52;
173
174 return;
175}
sin
double sin(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:185
cos
double cos(const BesAngle a)
Definition InstallArea/x86_64-el9-gcc13-dbg/include/MdcGeom/BesAngle.h:187
The documentation for this class was generated from the following files: