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

#include <D0ToKSpipi.h>

Public Member Functions

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

Detailed Description

Definition at line 9 of file D0ToKSpipi.h.

Constructor & Destructor Documentation

◆ D0ToKSpipi()

D0ToKSpipi::D0ToKSpipi ( )
inline

Definition at line 12 of file D0ToKSpipi.h.

12{}

◆ ~D0ToKSpipi()

D0ToKSpipi::~D0ToKSpipi ( )
virtual

Definition at line 23 of file D0ToKSpipi.cxx.

23{}

Member Function Documentation

◆ Amp_PFT()

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

Definition at line 183 of file D0ToKSpipi.cxx.

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

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0ToKSpipi::init ( )

Definition at line 25 of file D0ToKSpipi.cxx.

25 {
26 // std::cout << "D0ToKSpipi ==> 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
33 mass_R[0] = 0.77526;
34 width_R[0] = 0.14740;
35 spin_R[0] = 1;
36 ar[0] = 1;
37 phir[0] = 0;
38 mass_R[1] = 0.78266;
39 width_R[1] = 0.00868;
40 spin_R[1] = 1;
41 ar[1] = 0.037606;
42 phir[1] = 109.677;
43 mass_R[2] = 1.27550;
44 width_R[2] = 0.18670;
45 spin_R[2] = 2;
46 ar[2] = 1.54909;
47 phir[2] = -42.7425;
48 mass_R[3] = 1.46500;
49 width_R[3] = 0.40000;
50 spin_R[3] = 1;
51 ar[3] = 3.70735;
52 phir[3] = 103.644;
53 mass_R[4] = 0.89167;
54 width_R[4] = 0.0514;
55 spin_R[4] = 1;
56 ar[4] = 1.86093;
57 phir[4] = 136.529;
58 mass_R[5] = 1.42730;
59 width_R[5] = 0.10000;
60 spin_R[5] = 2;
61 ar[5] = 1.74288;
62 phir[5] = -48.0968;
63 mass_R[6] = 1.71800;
64 width_R[6] = 0.3220;
65 spin_R[6] = 1;
66 ar[6] = 3.31;
67 phir[6] = -118.2;
68 mass_R[7] = 1.41400;
69 width_R[7] = 0.2320;
70 spin_R[7] = 1;
71 ar[7] = 0.171672;
72 phir[7] = -68.41;
73 mass_R[8] = 0.89167;
74 width_R[8] = 0.0514;
75 spin_R[8] = 1;
76 ar[8] = 0.164;
77 phir[8] = -42.2;
78 mass_R[9] = 1.42730;
79 width_R[9] = 0.1000;
80 spin_R[9] = 2;
81 ar[9] = 0.1;
82 phir[9] = -89.6;
83 mass_R[10] = 1.41400;
84 width_R[10] = 0.2320;
85 spin_R[10] = 1;
86 ar[10] = 0.21;
87 phir[10] = 150.2;
88 mass_R[11] = 1.42500;
89 width_R[11] = 0.2700;
90 spin_R[11] = 1;
91 ar[11] = 2.78276;
92 phir[11] = 97.9608;
93 mass_R[12] = 1.42500;
94 width_R[12] = 0.2700;
95 spin_R[12] = 1;
96 ar[12] = 0.11;
97 phir[12] = 162.3;
98
99 beta[0] = complex<double>( 0.255303 * cos( 47.8861 * pi180inv ),
100 0.255303 * sin( 47.8861 * pi180inv ) );
101 beta[1] = complex<double>( 13.4446 * cos( -5.11127 * pi180inv ),
102 13.4446 * sin( -5.11127 * pi180inv ) );
103 beta[2] = complex<double>( 38.8496 * cos( -30.06 * pi180inv ),
104 38.8496 * sin( -30.06 * pi180inv ) );
105 beta[3] = complex<double>( 13.1086 * cos( -81.4148 * pi180inv ),
106 13.1086 * sin( -81.4148 * pi180inv ) );
107 beta[4] = complex<double>( 0., 0. );
108
109 fprod[0] = complex<double>( 5.08049 * cos( -182.312 * pi180inv ),
110 5.08049 * sin( -182.312 * pi180inv ) );
111 fprod[1] = complex<double>( 17.2388 * cos( -219.209 * pi180inv ),
112 17.2388 * sin( -219.209 * pi180inv ) );
113 fprod[2] = complex<double>( 19.0145 * cos( -76.9884 * pi180inv ),
114 19.0145 * sin( -76.9884 * pi180inv ) );
115 fprod[3] = complex<double>( 11.9875 * cos( -190.502 * pi180inv ),
116 11.9875 * sin( -190.502 * pi180inv ) );
117 fprod[4] = complex<double>( 0., 0. );
118 // beta.push_back( complex<double>( 0.255303*cos( 47.8861 *pi180inv), 0.255303*sin( 47.8861
119 // *pi180inv)) ); beta.push_back( complex<double>(13.4446 *cos( -5.11127*pi180inv), 13.4446
120 // *sin( -5.11127*pi180inv)) ); beta.push_back( complex<double>(38.8496 *cos(-30.06
121 // *pi180inv), 38.8496 *sin(-30.06 *pi180inv)) ); beta.push_back( complex<double>(13.1086
122 // *cos(-81.4148 *pi180inv), 13.1086 *sin(-81.4148 *pi180inv)) ); beta.push_back(
123 // complex<double>(0., 0.) );
124
125 // fprod.push_back(
126 // complex<double>(5.08049*cos(-182.312*pi180inv), 5.08049*sin(-182.312*pi180inv)));
127 // fprod.push_back(
128 // complex<double>(17.2388*cos(-219.209*pi180inv), 17.2388*sin(-219.209*pi180inv)));
129 // fprod.push_back(
130 // complex<double>(19.0145*cos(-76.9884*pi180inv), 19.0145*sin(-76.9884*pi180inv)));
131 // fprod.push_back(
132 // complex<double>(11.9875*cos(-190.502*pi180inv), 11.9875*sin(-190.502*pi180inv)));
133 // fprod.push_back( complex<double>(0., 0.));
134
135 ma[0] = 0.651;
136 g[0][0] = 0.22889;
137 g[0][1] = -0.55377;
138 g[0][2] = 0;
139 g[0][3] = -0.39899;
140 g[0][4] = -0.34639;
141 ma[1] = 1.20360;
142 g[1][0] = 0.94128;
143 g[1][1] = 0.55095;
144 g[1][2] = 0;
145 g[1][3] = 0.39065;
146 g[1][4] = 0.31503;
147 ma[2] = 1.55817;
148 g[2][0] = 0.36856;
149 g[2][1] = 0.23888;
150 g[2][2] = 0.55639;
151 g[2][3] = 0.18340;
152 g[2][4] = 0.18681;
153 ma[3] = 1.21000;
154 g[3][0] = 0.33650;
155 g[3][1] = 0.40907;
156 g[3][2] = 0.85679;
157 g[3][3] = 0.19906;
158 g[3][4] = -0.00984;
159 ma[4] = 1.82206;
160 g[4][0] = 0.18171;
161 g[4][1] = -0.17558;
162 g[4][2] = -0.79658;
163 g[4][3] = -0.00355;
164 g[4][4] = 0.22358;
165
166 // Hadronic parameters for tag modes: 0=no-specified, 1=Kpi, 2=Kpipi0, 3=K3pi
167 rd[0] = 0.0;
168 rd[1] = 0.0586;
169 rd[2] = 0.0440;
170 rd[3] = 0.0546;
171 deltad[0] = 0.0;
172 deltad[1] = 194.7 * pi180inv;
173 deltad[2] = 196.0 * pi180inv;
174 deltad[3] = 167.0 * pi180inv;
175 Rf[0] = 0.0;
176 Rf[1] = 1.0;
177 Rf[2] = 0.78;
178 Rf[3] = 0.52;
179
180 return;
181}
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

Referenced by QCMCFilter::findD0Decay().

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