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

#include <D0ToKLpipi2018.h>

Public Member Functions

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

Detailed Description

Definition at line 9 of file D0ToKLpipi2018.h.

Constructor & Destructor Documentation

◆ D0ToKLpipi2018()

D0ToKLpipi2018::D0ToKLpipi2018 ( )
inline

Definition at line 12 of file D0ToKLpipi2018.h.

12{}

◆ ~D0ToKLpipi2018()

D0ToKLpipi2018::~D0ToKLpipi2018 ( )
virtual

Definition at line 24 of file D0ToKLpipi2018.cxx.

24{}

Member Function Documentation

◆ Amp_PFT()

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

Definition at line 181 of file D0ToKLpipi2018.cxx.

182 {
183
184 vector<double> pD;
185 pD.clear();
186 if ( k0l.size() != 4 || pip.size() != 4 || pim.size() != 4 )
187 cout << "ERROR in KSPIPI daughter 4 momentum" << endl;
188 for ( int i = 0; i < k0l.size(); i++ ) { pD.push_back( k0l[i] + pip[i] + pim[i] ); }
189
190 // EvtResonance2 DK2piRes0(pD, pip, pim, ar[0], phir[0], width_R[0], mass_R[0], spin[0]);
191 // //ar, phir, width, mass, spin //Rho770 EvtResonance2 DK2piRes1(pD, pip, pim, ar[1],
192 // phir[1], width_R[1], mass_R[1], spin[1]); //ar, phir, width, mass, spin //Omega782
193 // EvtResonance2 DK2piRes2(pD, pip, pim, ar[2], phir[2], width_R[2], mass_R[2], spin[2]);
194 // //ar, phir, width, mass, spin //ftwo1270 EvtResonance2 DK2piRes3(pD, pip, pim, ar[3],
195 // phir[3], width_R[3], mass_R[3], spin[3]); //ar, phir, width, mass, spin //Rho1450
196 // EvtResonance2 DK2piRes4(pD, k0l, pim, ar[4], phir[4], width_R[4], mass_R[4], spin[4]);
197 // //ar, phir, width, mass, spin //Kstar892minus EvtResonance2 DK2piRes5(pD, k0l, pim,
198 // ar[5], phir[5], width_R[5], mass_R[5], spin[5]); //ar, phir, width, mass, spin
199 // //K2star1430minus EvtResonance2 DK2piRes6(pD, k0l, pim, ar[6], phir[6], width_R[6],
200 // mass_R[6], spin[6]); //ar, phir, width, mass, spin //Kstar1680minus EvtResonance2
201 // DK2piRes7(pD, k0l, pim, ar[7], phir[7], width_R[7], mass_R[7], spin[7]); //ar, phir,
202 // width, mass, spin //Kstar1410minus EvtResonance2 DK2piRes8(pD, k0l, pip, ar[8], phir[8],
203 // width_R[8], mass_R[8], spin[8]); //ar, phir, width, mass, spin //Kstar892plus
204 // EvtResonance2 DK2piRes9(pD, k0l, pip, ar[9], phir[9], width_R[9], mass_R[9], spin[9]);
205 // //ar, phir, width, mass, spin //K2star1430plus EvtResonance2 DK2piRes10(pD, k0l, pip,
206 // ar[10], phir[10], width_R[10], mass_R[10], spin[10]); //ar, phir, width, mass, spin
207 // //Kstar1410plus
208 complex<double> DK2piRes0 = Resonance2( pD, pip, pim, ar[0], phir[0], width_R[0], mass_R[0],
209 spin_R[0] ); // ar, phir, width, mass, spin Rho770
210 complex<double> DK2piRes1 = Resonance2( pD, pip, pim, ar[1], phir[1], width_R[1], mass_R[1],
211 spin_R[1] ); // ar, phir, width, mass, spin Omega782
212 complex<double> DK2piRes2 = Resonance2( pD, pip, pim, ar[2], phir[2], width_R[2], mass_R[2],
213 spin_R[2] ); // ar, phir, width, mass, spin ftwo1270
214 complex<double> DK2piRes3 = Resonance2( pD, pip, pim, ar[3], phir[3], width_R[3], mass_R[3],
215 spin_R[3] ); // ar, phir, width, mass, spin Rho1450
216 complex<double> DK2piRes4 = Resonance2( pD, k0l, pim, ar[4], phir[4], width_R[4], mass_R[4],
217 spin_R[4] ); // ar, phir, width, mass, spin Kstar892-
218 complex<double> DK2piRes5 =
219 Resonance2( pD, k0l, pim, ar[5], phir[5], width_R[5], mass_R[5],
220 spin_R[5] ); // ar, phir, width, mass, spin K2star1430-
221 complex<double> DK2piRes6 =
222 Resonance2( pD, k0l, pim, ar[6], phir[6], width_R[6], mass_R[6],
223 spin_R[6] ); // ar, phir, width, mass, spin Kstar1680-
224 complex<double> DK2piRes7 =
225 Resonance2( pD, k0l, pim, ar[7], phir[7], width_R[7], mass_R[7],
226 spin_R[7] ); // ar, phir, width, mass, spin Kstar1410-
227 complex<double> DK2piRes8 = Resonance2( pD, k0l, pip, ar[8], phir[8], width_R[8], mass_R[8],
228 spin_R[8] ); // ar, phir, width, mass, spin Kstar892+
229 complex<double> DK2piRes9 =
230 Resonance2( pD, k0l, pip, ar[9], phir[9], width_R[9], mass_R[9],
231 spin_R[9] ); // ar, phir, width, mass, spin K2star1430+
232 complex<double> DK2piRes10 =
233 Resonance2( pD, k0l, pip, ar[10], phir[10], width_R[10], mass_R[10],
234 spin_R[10] ); // ar, phir, width, mass, spin Kstar1410+
235
236 complex<double> pipi_s_wave = K_matrix( pip, pim );
237 if ( pipi_s_wave == complex<double>( 9999., 9999. ) ) return 1e-20;
238
239 complex<double> kpi_s_wave =
240 amplitude_LASS( k0l, pip, pim, "k0lpim", ar[11], phir[11] * pi180inv );
241 // complex<double> kpi_s_wave_dcs = amplitude_LASS(k0l, pip, pim, "k0spip", ar[12],
242 // phir[12]*pi180inv) ; should be there but not observed yet GUESS
243
244 complex<double> TOT_PFT_AMP = DK2piRes0 * CP_mult[0] + DK2piRes1 * CP_mult[1] +
245 DK2piRes2 * CP_mult[2] + DK2piRes3 * CP_mult[3] + DK2piRes4 +
246 DK2piRes5 + DK2piRes6 + DK2piRes7 + DK2piRes8 * ( -1. ) +
247 DK2piRes9 * ( -1. ) + DK2piRes10 * ( -1. ) +
248 pipi_s_wave * CP_mult[4] + kpi_s_wave;
249
250 return TOT_PFT_AMP;
251}

◆ init()

void D0ToKLpipi2018::init ( )

Definition at line 26 of file D0ToKLpipi2018.cxx.

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