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

#include <D0ToKLpipi.h>

Public Member Functions

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

Detailed Description

Definition at line 9 of file D0ToKLpipi.h.

Constructor & Destructor Documentation

◆ D0ToKLpipi()

D0ToKLpipi::D0ToKLpipi ( )
inline

Definition at line 12 of file D0ToKLpipi.h.

12{}

◆ ~D0ToKLpipi()

D0ToKLpipi::~D0ToKLpipi ( )
virtual

Definition at line 24 of file D0ToKLpipi.cxx.

24{}

Member Function Documentation

◆ Amp_PFT()

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

Definition at line 191 of file D0ToKLpipi.cxx.

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

Referenced by QCMCFilter::findD0Decay().

◆ init()

void D0ToKLpipi::init ( )

Definition at line 26 of file D0ToKLpipi.cxx.

26 {
27 // std::cout << "D0ToKLpipi ==> 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 // spin[0] = 1 ; //Rho770
35 // spin[1] = 1 ; //Omega782
36 // spin[2] = 2 ; //ftwo1270
37 // spin[3] = 1 ; //Rho1450
38 // spin[4] = 1 ; //Kstar892minus
39 // spin[5] = 2 ; //K2star1430minus
40 // spin[6] = 1 ; //Kstar1680minus
41 // spin[7] = 1 ; //Kstar1410minus
42 // spin[8] = 1 ; //Kstar892plus
43 // spin[9] = 2 ; //K2star1430plus
44 // spin[10]= 1 ; //Kstar1410plus
45 mass_R[0] = 0.77526;
46 width_R[0] = 0.14740;
47 spin_R[0] = 1;
48 ar[0] = 1;
49 phir[0] = 0;
50 mass_R[1] = 0.78266;
51 width_R[1] = 0.00868;
52 spin_R[1] = 1;
53 ar[1] = 0.03928;
54 phir[1] = 106.1;
55 mass_R[2] = 1.27550;
56 width_R[2] = 0.18670;
57 spin_R[2] = 2;
58 ar[2] = 1.302;
59 phir[2] = -39.97;
60 mass_R[3] = 1.46500;
61 width_R[3] = 0.40000;
62 spin_R[3] = 1;
63 ar[3] = 1.656;
64 phir[3] = 113;
65 mass_R[4] = 0.89167;
66 width_R[4] = 0.0514;
67 spin_R[4] = 1;
68 ar[4] = 1.843;
69 phir[4] = 138;
70 mass_R[5] = 1.42730;
71 width_R[5] = 0.10000;
72 spin_R[5] = 2;
73 ar[5] = 1.537;
74 phir[5] = -49.45;
75 mass_R[6] = 1.71800;
76 width_R[6] = 0.3220;
77 spin_R[6] = 1;
78 ar[6] = 2.251;
79 phir[6] = -196.3;
80 mass_R[7] = 1.41400;
81 width_R[7] = 0.2320;
82 spin_R[7] = 1;
83 ar[7] = 0.4733;
84 phir[7] = 1659;
85 mass_R[8] = 0.89167;
86 width_R[8] = 0.0514;
87 spin_R[8] = 1;
88 ar[8] = 0.1617;
89 phir[8] = -36.25;
90 mass_R[9] = 1.42730;
91 width_R[9] = 0.1000;
92 spin_R[9] = 2;
93 ar[9] = 0.2267;
94 phir[9] = -84.64;
95 mass_R[10] = 1.41400;
96 width_R[10] = 0.2320;
97 spin_R[10] = 1;
98 ar[10] = 0.2185;
99 phir[10] = 73.12;
100 mass_R[11] = 1.42500;
101 width_R[11] = 0.2700;
102 spin_R[11] = 1;
103 ar[11] = 2.396;
104 phir[11] = 98.14;
105 // beta[kb] = EvtComplex(mag*cos(phase*pi180inv), mag*sin(phase*pi180inv)) ;
106 // fprod[kf] = EvtComplex(mag*cos(phase*pi180inv), mag*sin(phase*pi180inv)) ;
107 beta[0] = complex<double>( 8.5 * cos( 68.5 * pi180inv ), 8.5 * sin( 68.5 * pi180inv ) );
108 beta[1] = complex<double>( 12.2 * cos( 24 * pi180inv ), 12.2 * sin( 24 * pi180inv ) );
109 beta[2] = complex<double>( 29.2 * cos( -0.1 * pi180inv ), 29.2 * sin( -0.1 * pi180inv ) );
110 beta[3] = complex<double>( 10.8 * cos( -51.9 * pi180inv ), 10.8 * sin( -51.9 * pi180inv ) );
111 beta[4] = complex<double>( 0., 0. );
112
113 fprod[0] = complex<double>( 8 * cos( -126 * pi180inv ), 8 * sin( -126 * pi180inv ) );
114 fprod[1] =
115 complex<double>( 26.3 * cos( -152.3 * pi180inv ), 26.3 * sin( -152.3 * pi180inv ) );
116 fprod[2] = complex<double>( 33 * cos( -93.2 * pi180inv ), 33 * sin( -93.2 * pi180inv ) );
117 fprod[3] =
118 complex<double>( 26.2 * cos( -121.4 * pi180inv ), 26.2 * sin( -121.4 * pi180inv ) );
119 fprod[4] = complex<double>( 0., 0. );
120
121 // CP_mult[w] = (EvtComplex(1., 0.) - 2.*tan2thetaC*EvtComplex(r*cos(delta), r*sin(delta))) ;
122 CP_mult[0] =
123 complex<double>( 1., 0. ) - 2. * tan2thetaC *
124 complex<double>( 1.851 * cos( -94.07 * pi180inv ),
125 1.851 * sin( -94.07 * pi180inv ) );
126 CP_mult[1] =
127 complex<double>( 1., 0. ) -
128 2. * tan2thetaC *
129 complex<double>( 6.332 * cos( 2.103 * pi180inv ), 6.332 * sin( 2.103 * pi180inv ) );
130 CP_mult[2] =
131 complex<double>( 1., 0. ) - 2. * tan2thetaC *
132 complex<double>( 3.229 * cos( -60.05 * pi180inv ),
133 3.229 * sin( -60.05 * pi180inv ) );
134 CP_mult[3] =
135 complex<double>( 1., 0. ) -
136 2. * tan2thetaC *
137 complex<double>( 12.75 * cos( 73.62 * pi180inv ), 12.75 * sin( 73.62 * pi180inv ) );
138 CP_mult[4] =
139 complex<double>( 1., 0. ) - 2. * tan2thetaC *
140 complex<double>( 0.7116 * cos( -177.149 * pi180inv ),
141 0.7116 * sin( -177.149 * pi180inv ) );
142
143 ma[0] = 0.651;
144 g[0][0] = 0.22889;
145 g[0][1] = -0.55377;
146 g[0][2] = 0;
147 g[0][3] = -0.39899;
148 g[0][4] = -0.34639;
149 ma[1] = 1.20360;
150 g[1][0] = 0.94128;
151 g[1][1] = 0.55095;
152 g[1][2] = 0;
153 g[1][3] = 0.39065;
154 g[1][4] = 0.31503;
155 ma[2] = 1.55817;
156 g[2][0] = 0.36856;
157 g[2][1] = 0.23888;
158 g[2][2] = 0.55639;
159 g[2][3] = 0.18340;
160 g[2][4] = 0.18681;
161 ma[3] = 1.21000;
162 g[3][0] = 0.33650;
163 g[3][1] = 0.40907;
164 g[3][2] = 0.85679;
165 g[3][3] = 0.19906;
166 g[3][4] = -0.00984;
167 ma[4] = 1.82206;
168 g[4][0] = 0.18171;
169 g[4][1] = -0.17558;
170 g[4][2] = -0.79658;
171 g[4][3] = -0.00355;
172 g[4][4] = 0.22358;
173
174 // Hadronic parameters for tag modes: 0=no-specified, 1=Kpi, 2=Kpipi0, 3=K3pi
175 rd[0] = 0.0;
176 rd[1] = 0.0586;
177 rd[2] = 0.0440;
178 rd[3] = 0.0546;
179 deltad[0] = 0.0;
180 deltad[1] = 194.7 * pi180inv;
181 deltad[2] = 196.0 * pi180inv;
182 deltad[3] = 167.0 * pi180inv;
183 Rf[0] = 0.0;
184 Rf[1] = 1.0;
185 Rf[2] = 0.78;
186 Rf[3] = 0.52;
187
188 return;
189}
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: