]>
Commit | Line | Data |
---|---|---|
da0e9ce3 | 1 | //-------------------------------------------------------------------------- |
2 | // | |
3 | // Environment: | |
4 | // This software is part of the EvtGen package developed jointly | |
5 | // for the BaBar and CLEO collaborations. If you use all or part | |
6 | // of it, please give an appropriate acknowledgement. | |
7 | // | |
8 | // Copyright Information: See EvtGen/COPYRIGHT | |
9 | // Copyright (C) 1998 Caltech, UCSB | |
10 | // | |
11 | // Module: EvtISGW2FF.cc | |
12 | // | |
13 | // Description: Routine to implement semileptonic form factors | |
14 | // according to the model ISGW2 | |
15 | // | |
16 | // Modification history: | |
17 | // | |
18 | // DJL April 17, 1998 Module created | |
19 | // | |
20 | //------------------------------------------------------------------------ | |
21 | // | |
22 | #include "EvtGenBase/EvtPatches.hh" | |
23 | #include "EvtGenBase/EvtPatches.hh" | |
24 | #include "EvtGenBase/EvtReport.hh" | |
25 | #include "EvtGenModels/EvtISGW2FF.hh" | |
26 | #include "EvtGenBase/EvtPDL.hh" | |
27 | #include "EvtGenBase/EvtId.hh" | |
28 | #include "EvtGenBase/EvtIdSet.hh" | |
29 | #include "EvtGenBase/EvtConst.hh" | |
30 | #include <string> | |
31 | #include <math.h> | |
32 | #include <stdlib.h> | |
33 | using std::endl; | |
34 | ||
35 | void EvtISGW2FF::getscalarff(EvtId parent,EvtId daught, | |
36 | double t, double mass, double *fpf, | |
37 | double *f0f ) { | |
38 | ||
39 | //added by Lange Jan4,2000 | |
40 | ||
41 | ||
42 | static EvtId BP=EvtPDL::getId("B+"); | |
43 | static EvtId BM=EvtPDL::getId("B-"); | |
44 | static EvtId B0=EvtPDL::getId("B0"); | |
45 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
46 | ||
47 | static EvtId DST0=EvtPDL::getId("D*0"); | |
48 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
49 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
50 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
51 | static EvtId D0=EvtPDL::getId("D0"); | |
52 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
53 | static EvtId DP=EvtPDL::getId("D+"); | |
54 | static EvtId DM=EvtPDL::getId("D-"); | |
55 | ||
56 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
57 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
58 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
59 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
60 | ||
61 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
62 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
63 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
64 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
65 | ||
66 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
67 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
68 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
69 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
70 | ||
71 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
72 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
73 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
74 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
75 | ||
76 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
77 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
78 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
79 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
80 | ||
81 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
82 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
83 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
84 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
85 | ||
86 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
87 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
88 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
89 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
90 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
91 | ||
92 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
93 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
94 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
95 | ||
96 | static EvtId PIP=EvtPDL::getId("pi+"); | |
97 | static EvtId PIM=EvtPDL::getId("pi-"); | |
98 | static EvtId PI0=EvtPDL::getId("pi0"); | |
99 | ||
100 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
101 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
102 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
103 | ||
104 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
105 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
106 | static EvtId A20=EvtPDL::getId("a_20"); | |
107 | ||
108 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
109 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
110 | static EvtId A10=EvtPDL::getId("a_10"); | |
111 | ||
112 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
113 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
114 | static EvtId A00=EvtPDL::getId("a_00"); | |
115 | ||
116 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
117 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
118 | static EvtId B10=EvtPDL::getId("b_10"); | |
119 | ||
120 | static EvtId H1=EvtPDL::getId("h_1"); | |
121 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
122 | ||
123 | static EvtId F1=EvtPDL::getId("f_1"); | |
124 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
125 | static EvtId F0=EvtPDL::getId("f_0"); | |
126 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
127 | static EvtId F2=EvtPDL::getId("f_2"); | |
128 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
129 | ||
130 | static EvtId ETA=EvtPDL::getId("eta"); | |
131 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
132 | static EvtId OMEG=EvtPDL::getId("omega"); | |
133 | ||
134 | static EvtId KP=EvtPDL::getId("K+"); | |
135 | static EvtId KM=EvtPDL::getId("K-"); | |
136 | static EvtId K0=EvtPDL::getId("K0"); | |
137 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
138 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
139 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
140 | ||
141 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
142 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
143 | static EvtId KST0=EvtPDL::getId("K*0"); | |
144 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
145 | ||
146 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
147 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
148 | static EvtId K10=EvtPDL::getId("K_10"); | |
149 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
150 | ||
151 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
152 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
153 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
154 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
155 | ||
156 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
157 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
158 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
159 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
160 | ||
161 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
162 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
163 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
164 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
165 | ||
166 | static EvtId PHI=EvtPDL::getId("phi"); | |
167 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
168 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
169 | ||
170 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
171 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
172 | ||
173 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
174 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
175 | ||
176 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
177 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
178 | ||
179 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
180 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
181 | ||
182 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
183 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
184 | ||
185 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
186 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
187 | ||
188 | double fmf; | |
189 | double mb=EvtPDL::getMeanMass(parent); | |
190 | ||
191 | if (daught==PI0||daught==PIP||daught==PIM||daught==ETA|| | |
192 | daught==ETAPR||daught==D0||daught==D0B||daught==DP|| | |
193 | daught==DM||daught==KP||daught==KM||daught==K0||daught==K0S|| | |
194 | daught==K0L||daught==KB||daught==DSP||daught==DSM) { | |
195 | ||
196 | EvtISGW2FF1S0(parent,daught,t,mass,fpf,&fmf); | |
197 | } | |
198 | ||
199 | if (daught==PI2S0||daught==PI2SP||daught==PI2SM||daught==ETA2S|| | |
200 | daught==D21S0P||daught==D21S0B||daught==D21S0N||daught==D21S00){ | |
201 | EvtISGW2FF21S0(parent,daught,t,mass,fpf,&fmf); | |
202 | } | |
203 | ||
204 | if (daught==A00||daught==A0P||daught==A0M||daught==F0|| | |
205 | daught==F0PR||daught==D3P0P||daught==D3P00||daught==D3P0B|| | |
206 | daught==D3P0N||daught==K0STM||daught==K0STB||daught==K0STP|| | |
207 | daught==D3P0SP||daught==D3P0SN|| | |
208 | daught==K0ST0) { | |
209 | EvtISGW2FF3P0(parent,daught,t,mass,fpf,&fmf); | |
210 | } | |
211 | ||
212 | *f0f = (fmf/((mb*mb-mass*mass)/t))+(*fpf); | |
213 | ||
214 | return ; | |
215 | } | |
216 | ||
217 | void EvtISGW2FF::gettensorff(EvtId parent,EvtId daught, | |
218 | double t, double mass, double *hf, | |
219 | double *kf, double *bpf, double *bmf ){ | |
220 | ||
221 | //added by Lange Jan4,2000 | |
222 | static EvtId EM=EvtPDL::getId("e-"); | |
223 | static EvtId EP=EvtPDL::getId("e+"); | |
224 | static EvtId MUM=EvtPDL::getId("mu-"); | |
225 | static EvtId MUP=EvtPDL::getId("mu+"); | |
226 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
227 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
228 | ||
229 | static EvtId BP=EvtPDL::getId("B+"); | |
230 | static EvtId BM=EvtPDL::getId("B-"); | |
231 | static EvtId B0=EvtPDL::getId("B0"); | |
232 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
233 | ||
234 | static EvtId DST0=EvtPDL::getId("D*0"); | |
235 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
236 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
237 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
238 | static EvtId D0=EvtPDL::getId("D0"); | |
239 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
240 | static EvtId DP=EvtPDL::getId("D+"); | |
241 | static EvtId DM=EvtPDL::getId("D-"); | |
242 | ||
243 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
244 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
245 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
246 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
247 | ||
248 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
249 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
250 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
251 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
252 | ||
253 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
254 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
255 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
256 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
257 | ||
258 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
259 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
260 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
261 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
262 | ||
263 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
264 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
265 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
266 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
267 | ||
268 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
269 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
270 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
271 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
272 | ||
273 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
274 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
275 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
276 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
277 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
278 | ||
279 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
280 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
281 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
282 | ||
283 | static EvtId PIP=EvtPDL::getId("pi+"); | |
284 | static EvtId PIM=EvtPDL::getId("pi-"); | |
285 | static EvtId PI0=EvtPDL::getId("pi0"); | |
286 | ||
287 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
288 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
289 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
290 | ||
291 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
292 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
293 | static EvtId A20=EvtPDL::getId("a_20"); | |
294 | ||
295 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
296 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
297 | static EvtId A10=EvtPDL::getId("a_10"); | |
298 | ||
299 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
300 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
301 | static EvtId A00=EvtPDL::getId("a_00"); | |
302 | ||
303 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
304 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
305 | static EvtId B10=EvtPDL::getId("b_10"); | |
306 | ||
307 | static EvtId H1=EvtPDL::getId("h_1"); | |
308 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
309 | ||
310 | static EvtId F1=EvtPDL::getId("f_1"); | |
311 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
312 | static EvtId F0=EvtPDL::getId("f_0"); | |
313 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
314 | static EvtId F2=EvtPDL::getId("f_2"); | |
315 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
316 | ||
317 | static EvtId ETA=EvtPDL::getId("eta"); | |
318 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
319 | static EvtId OMEG=EvtPDL::getId("omega"); | |
320 | ||
321 | static EvtId KP=EvtPDL::getId("K+"); | |
322 | static EvtId KM=EvtPDL::getId("K-"); | |
323 | static EvtId K0=EvtPDL::getId("K0"); | |
324 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
325 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
326 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
327 | ||
328 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
329 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
330 | static EvtId KST0=EvtPDL::getId("K*0"); | |
331 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
332 | ||
333 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
334 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
335 | static EvtId K10=EvtPDL::getId("K_10"); | |
336 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
337 | ||
338 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
339 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
340 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
341 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
342 | ||
343 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
344 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
345 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
346 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
347 | ||
348 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
349 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
350 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
351 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
352 | ||
353 | static EvtId PHI=EvtPDL::getId("phi"); | |
354 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
355 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
356 | ||
357 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
358 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
359 | ||
360 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
361 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
362 | ||
363 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
364 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
365 | ||
366 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
367 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
368 | ||
369 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
370 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
371 | ||
372 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
373 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
374 | ||
375 | EvtISGW2FF3P2(parent,daught,t,mass,hf,kf,bpf,bmf); | |
376 | ||
377 | return; | |
378 | ||
379 | } | |
380 | ||
381 | ||
382 | void EvtISGW2FF::getvectorff(EvtId parent,EvtId daught, | |
383 | double t, double mass, double *a1f, | |
384 | double *a2f, double *vf, double *a0f ){ | |
385 | double ff,gf,apf,amf; | |
386 | ||
387 | //added by Lange Jan4,2000 | |
388 | static EvtId EM=EvtPDL::getId("e-"); | |
389 | static EvtId EP=EvtPDL::getId("e+"); | |
390 | static EvtId MUM=EvtPDL::getId("mu-"); | |
391 | static EvtId MUP=EvtPDL::getId("mu+"); | |
392 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
393 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
394 | ||
395 | static EvtId BP=EvtPDL::getId("B+"); | |
396 | static EvtId BM=EvtPDL::getId("B-"); | |
397 | static EvtId B0=EvtPDL::getId("B0"); | |
398 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
399 | ||
400 | static EvtId DST0=EvtPDL::getId("D*0"); | |
401 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
402 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
403 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
404 | static EvtId D0=EvtPDL::getId("D0"); | |
405 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
406 | static EvtId DP=EvtPDL::getId("D+"); | |
407 | static EvtId DM=EvtPDL::getId("D-"); | |
408 | ||
409 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
410 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
411 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
412 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
413 | ||
414 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
415 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
416 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
417 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
418 | ||
419 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
420 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
421 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
422 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
423 | ||
424 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
425 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
426 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
427 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
428 | ||
429 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
430 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
431 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
432 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
433 | ||
434 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
435 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
436 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
437 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
438 | ||
439 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
440 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
441 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
442 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
443 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
444 | ||
445 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
446 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
447 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
448 | ||
449 | static EvtId PIP=EvtPDL::getId("pi+"); | |
450 | static EvtId PIM=EvtPDL::getId("pi-"); | |
451 | static EvtId PI0=EvtPDL::getId("pi0"); | |
452 | ||
453 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
454 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
455 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
456 | ||
457 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
458 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
459 | static EvtId A20=EvtPDL::getId("a_20"); | |
460 | ||
461 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
462 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
463 | static EvtId A10=EvtPDL::getId("a_10"); | |
464 | ||
465 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
466 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
467 | static EvtId A00=EvtPDL::getId("a_00"); | |
468 | ||
469 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
470 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
471 | static EvtId B10=EvtPDL::getId("b_10"); | |
472 | ||
473 | static EvtId H1=EvtPDL::getId("h_1"); | |
474 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
475 | ||
476 | static EvtId F1=EvtPDL::getId("f_1"); | |
477 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
478 | static EvtId F0=EvtPDL::getId("f_0"); | |
479 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
480 | static EvtId F2=EvtPDL::getId("f_2"); | |
481 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
482 | ||
483 | static EvtId ETA=EvtPDL::getId("eta"); | |
484 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
485 | static EvtId OMEG=EvtPDL::getId("omega"); | |
486 | ||
487 | static EvtId KP=EvtPDL::getId("K+"); | |
488 | static EvtId KM=EvtPDL::getId("K-"); | |
489 | static EvtId K0=EvtPDL::getId("K0"); | |
490 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
491 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
492 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
493 | ||
494 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
495 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
496 | static EvtId KST0=EvtPDL::getId("K*0"); | |
497 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
498 | ||
499 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
500 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
501 | static EvtId K10=EvtPDL::getId("K_10"); | |
502 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
503 | ||
504 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
505 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
506 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
507 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
508 | ||
509 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
510 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
511 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
512 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
513 | ||
514 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
515 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
516 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
517 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
518 | ||
519 | static EvtId PHI=EvtPDL::getId("phi"); | |
520 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
521 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
522 | ||
523 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
524 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
525 | ||
526 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
527 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
528 | ||
529 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
530 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
531 | ||
532 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
533 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
534 | ||
535 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
536 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
537 | ||
538 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
539 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
540 | ||
541 | if (daught==DST0||daught==DSTP||daught==DSTM||daught==DSTB|| | |
542 | daught==OMEG||daught==RHO0||daught==RHOM||daught==RHOP|| | |
543 | daught==KSTP||daught==KSTM||daught==KST0||daught==KSTB|| | |
544 | daught==PHI||daught==DSSTP||daught==DSSTM) { | |
545 | EvtISGW2FF3S1(parent,daught,t,mass,&ff,&gf,&apf,&amf); | |
546 | } | |
547 | if (daught==B10||daught==B1P||daught==B1M||daught==H1|| | |
548 | daught==H1PR||daught==D1P1P||daught==D1P10||daught==D1P1B|| | |
549 | daught==D1P1SP||daught==D1P1SN|| | |
550 | daught==D1P1N||daught==K10||daught==K1B||daught==K1P|| | |
551 | daught==K1M) { | |
552 | EvtISGW2FF1P1(parent,daught,t,mass,&ff,&gf,&apf,&amf); | |
553 | } | |
554 | if (daught==RHO2S0||daught==RHO2SP||daught==RHO2SM||daught==OMEG2S|| | |
555 | daught==D23S1P||daught==D23S1B||daught==D23S1N||daught==D23S10){ | |
556 | EvtISGW2FF23S1(parent,daught,t,mass,&ff,&gf,&apf,&amf); | |
557 | } | |
558 | if (daught==A10||daught==A1P||daught==A1M||daught==F1|| | |
559 | daught==F1PR||daught==D3P1P||daught==D3P10||daught==D3P1B|| | |
560 | daught==D3P1N||daught==K1STM||daught==K1STB||daught==K1STP|| | |
561 | daught==D3P1SP||daught==D3P1SN|| | |
562 | daught==K1ST0) { | |
563 | EvtISGW2FF3P1(parent,daught,t,mass,&ff,&gf,&apf,&amf); | |
564 | } | |
565 | ||
566 | // Need to stuff in some factors to make these the ffs that | |
567 | // is used elsewhere... | |
568 | ||
569 | double mb=EvtPDL::getMeanMass(parent); | |
570 | ||
571 | ||
572 | *vf = (gf)*(mb+mass); | |
573 | *a1f = (ff)/(mb+mass); | |
574 | *a2f = -1.0*(apf)*(mb+mass); | |
575 | ||
576 | double a3f = ((mb+mass)/(2.0*mass))*(*a1f) - | |
577 | ((mb-mass)/(2.0*mass))*(*a2f); | |
578 | ||
579 | *a0f = a3f + ( (t*amf)/(2.0*mass)); | |
580 | ||
581 | return; | |
582 | } | |
583 | ||
584 | ||
585 | ||
586 | void EvtISGW2FF::EvtISGW2FF1S0 (EvtId parent,EvtId daugt, | |
587 | double t, double mass, double *fpf, double *fmf ) { | |
588 | ||
589 | double mtb, mbb(0.0); | |
590 | double msd(0.0), mx,mb,nf(0.0),nfp(0.0); | |
591 | double msq(0.0),bx2(0.0),mbx(0.0),mtx; | |
592 | double zji,cji,w,gammaji,chiji,betaji_fppfm; | |
593 | double rfppfm,rfpmfm,f3fppfm,f3fpmfm,fppfm,fpmfm,al,ai,rcji,f3; | |
594 | double mqm,msb(0.0),bb2(0.0),mup,mum,bbx2,tm,wt,r2,betaji_fpmfm; | |
595 | ||
596 | EvtId prnt=parent; | |
597 | EvtId dgt=daugt; | |
598 | ||
599 | //added by Lange Jan4,2000 | |
600 | static EvtId EM=EvtPDL::getId("e-"); | |
601 | static EvtId EP=EvtPDL::getId("e+"); | |
602 | static EvtId MUM=EvtPDL::getId("mu-"); | |
603 | static EvtId MUP=EvtPDL::getId("mu+"); | |
604 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
605 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
606 | ||
607 | static EvtIdSet theB("B+","B-","B0","anti-B0"); | |
608 | ||
609 | static EvtId DST0=EvtPDL::getId("D*0"); | |
610 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
611 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
612 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
613 | static EvtId D0=EvtPDL::getId("D0"); | |
614 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
615 | static EvtId DP=EvtPDL::getId("D+"); | |
616 | static EvtId DM=EvtPDL::getId("D-"); | |
617 | ||
618 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
619 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
620 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
621 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
622 | ||
623 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
624 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
625 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
626 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
627 | ||
628 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
629 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
630 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
631 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
632 | ||
633 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
634 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
635 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
636 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
637 | ||
638 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
639 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
640 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
641 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
642 | ||
643 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
644 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
645 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
646 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
647 | ||
648 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
649 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
650 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
651 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
652 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
653 | ||
654 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
655 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
656 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
657 | ||
658 | static EvtId PIP=EvtPDL::getId("pi+"); | |
659 | static EvtId PIM=EvtPDL::getId("pi-"); | |
660 | static EvtId PI0=EvtPDL::getId("pi0"); | |
661 | ||
662 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
663 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
664 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
665 | ||
666 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
667 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
668 | static EvtId A20=EvtPDL::getId("a_20"); | |
669 | ||
670 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
671 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
672 | static EvtId A10=EvtPDL::getId("a_10"); | |
673 | ||
674 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
675 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
676 | static EvtId A00=EvtPDL::getId("a_00"); | |
677 | ||
678 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
679 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
680 | static EvtId B10=EvtPDL::getId("b_10"); | |
681 | ||
682 | static EvtId H1=EvtPDL::getId("h_1"); | |
683 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
684 | ||
685 | static EvtId F1=EvtPDL::getId("f_1"); | |
686 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
687 | static EvtId F0=EvtPDL::getId("f_0"); | |
688 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
689 | static EvtId F2=EvtPDL::getId("f_2"); | |
690 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
691 | ||
692 | static EvtId ETA=EvtPDL::getId("eta"); | |
693 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
694 | static EvtId OMEG=EvtPDL::getId("omega"); | |
695 | ||
696 | static EvtId KP=EvtPDL::getId("K+"); | |
697 | static EvtId KM=EvtPDL::getId("K-"); | |
698 | static EvtId K0=EvtPDL::getId("K0"); | |
699 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
700 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
701 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
702 | ||
703 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
704 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
705 | static EvtId KST0=EvtPDL::getId("K*0"); | |
706 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
707 | ||
708 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
709 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
710 | static EvtId K10=EvtPDL::getId("K_10"); | |
711 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
712 | ||
713 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
714 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
715 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
716 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
717 | ||
718 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
719 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
720 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
721 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
722 | ||
723 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
724 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
725 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
726 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
727 | ||
728 | static EvtId PHI=EvtPDL::getId("phi"); | |
729 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
730 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
731 | ||
732 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
733 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
734 | ||
735 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
736 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
737 | ||
738 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
739 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
740 | ||
741 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
742 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
743 | ||
744 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
745 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
746 | ||
747 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
748 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
749 | ||
750 | if (theB.contains(prnt)) { | |
751 | ||
752 | msb=5.2; | |
753 | msd=0.33; | |
754 | bb2=0.431*0.431; | |
755 | mbb=5.31; | |
756 | nf = 4.0; | |
757 | ||
758 | if (dgt==PI0||dgt==PIP||dgt==PIM||dgt==ETA||dgt==ETAPR) { | |
759 | ||
760 | msq=0.33; | |
761 | bx2=0.406*0.406; | |
762 | mbx=0.75*0.770+0.25*0.14; | |
763 | nfp = 0.0; | |
764 | } | |
765 | else{ | |
766 | if (dgt==D0||dgt==D0B||dgt==DP||dgt==DM) { | |
767 | msq=1.82; | |
768 | bx2=0.45*0.45; | |
769 | mbx=0.75*2.01+0.25*1.87; | |
770 | nfp = 3.0; | |
771 | } | |
772 | else{ | |
773 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
774 | } | |
775 | } | |
776 | } | |
777 | else{ | |
778 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
779 | msb=1.82; | |
780 | msd=0.33; | |
781 | bb2=0.45*0.45; | |
782 | mbb=1.963; | |
783 | nf = 3.0; | |
784 | ||
785 | if (dgt==PIP||dgt==PIM||dgt==PI0||dgt==ETA||dgt==ETAPR) { | |
786 | msq=0.33; | |
787 | bx2=0.406*0.406; | |
788 | mbx=0.75*0.770+0.25*0.14; | |
789 | nfp = 0.0; | |
790 | } | |
791 | else{ | |
792 | if (dgt==K0||dgt==K0S||dgt==K0L||dgt==KB||dgt==KP||dgt==KM) { | |
793 | msq=0.55; | |
794 | bx2=0.44*0.44; | |
795 | mbx=0.75*0.892+0.25*0.49767; | |
796 | nfp = 2.0; | |
797 | } | |
798 | else{ | |
799 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
800 | } | |
801 | } | |
802 | } | |
803 | else{ | |
804 | if (prnt==DSP||prnt==DSM){ | |
805 | msb=1.82; | |
806 | msd=0.55; | |
807 | bb2=0.56*0.56; | |
808 | mbb=1.968; | |
809 | nf = 3.0; | |
810 | ||
811 | if (dgt==K0||dgt==K0S||dgt==K0L||dgt==KB) { | |
812 | ||
813 | msq=0.33; | |
814 | bx2=0.44*0.44; | |
815 | mbx=0.75*0.770+0.25*0.14; | |
816 | nfp = 0.0; | |
817 | } | |
818 | else{ | |
819 | if (dgt==PI0||dgt==ETA||dgt==ETAPR) { | |
820 | msq=0.33; | |
821 | bx2=0.53*0.53; | |
822 | mbx=0.75*0.892+0.25*0.49767; | |
823 | nfp = 0.0; | |
824 | } | |
825 | else{ | |
826 | ||
827 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
828 | } | |
829 | } | |
830 | } | |
831 | else{ | |
832 | //BS -> cs constants added by djl on Jan. 21,1998 | |
833 | if (prnt==BS0||prnt==BSB){ | |
834 | ||
835 | msb=5.2; | |
836 | msd=0.55; | |
837 | bb2=0.54*0.54; | |
838 | mbb=5.38; | |
839 | nf = 4.0; | |
840 | ||
841 | if (dgt==DSP||dgt==DSM) { | |
842 | ||
843 | msq=1.82; | |
844 | bx2=0.56*0.56; | |
845 | mbx=0.75*2.11+0.25*1.97; | |
846 | nfp = 3.0; | |
847 | } | |
848 | else if (dgt==KP||dgt==KM) { | |
849 | ||
850 | msq=0.55; | |
851 | bx2=0.44*0.44; | |
852 | mbx=0.75*0.892+0.25*0.49767; | |
853 | nfp = 2.0; | |
854 | } | |
855 | else{ | |
856 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
857 | } | |
858 | } | |
859 | ||
860 | else{ | |
861 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_1S0.\n"; | |
862 | report(ERROR,"EvtGen") << "Parent:"<<parent.getId()<<endl; | |
863 | } | |
864 | } | |
865 | } | |
866 | } | |
867 | ||
868 | mtb = msb + msd; | |
869 | mtx = msq + msd; | |
870 | mb=EvtPDL::getMeanMass(parent); | |
871 | mx=mass; | |
872 | ||
873 | mup=1.0/(1.0/msq+1.0/msb); | |
874 | mum=1.0/(1.0/msq-1.0/msb); | |
875 | bbx2=0.5*(bb2+bx2); | |
876 | tm=(mb-mx)*(mb-mx); | |
877 | if ( t>tm ) t=0.99*tm; | |
878 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
879 | ||
880 | mqm = 0.1; | |
881 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
882 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
883 | log(EvtGetas(mqm,mqm)/EvtGetas(msq,msq)); | |
884 | ||
885 | f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,1.5) / | |
886 | ((1.0+r2*(tm-t)/12.0)*(1.0+r2*(tm-t)/12.0)); | |
887 | ||
888 | // for w use wt def with physical masses. | |
889 | // report(ERROR,"EvtGen") << "before w\n"; | |
890 | ||
891 | w = 1.0 + (( tm - t ) / ( 2.0* mb * mx )); | |
892 | rcji = ( 1/sqrt(w*w -1 ))*log( w + sqrt( w*w -1 )); | |
893 | al = (8.0 / ( 33.0 - 2.0*nfp ))*(w*rcji -1.0 ); | |
894 | ai = -1.0* ( 6.0/( 33.0 - 2.0*nf)); | |
895 | cji = pow(( EvtGetas( msb,msb ) / EvtGetas( msq,msq ) ),ai); | |
896 | ||
897 | zji = msq / msb; | |
898 | ||
899 | gammaji = EvtGetGammaji( zji ); | |
900 | chiji = -1.0 - ( gammaji / ( 1- zji )); | |
901 | betaji_fppfm = gammaji - (2.0/3.0)*chiji; | |
902 | betaji_fpmfm = gammaji + (2.0/3.0)*chiji; | |
903 | rfppfm = cji *(1.0 + betaji_fppfm*EvtGetas( msq,sqrt(msb*msq) )/EvtConst::pi); | |
904 | rfpmfm = cji *(1.0 + betaji_fpmfm*EvtGetas( msq,sqrt(msb*msq) )/EvtConst::pi); | |
905 | f3fppfm = f3*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),0.5); | |
906 | f3fpmfm = f3*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),-0.5); | |
907 | fppfm = f3fppfm* rfppfm * ( 2.0 - ( ( mtx/msq)*(1- ( (msd*msq*bb2) | |
908 | /(2.0*mup*mtx*bbx2))))); | |
909 | fpmfm = f3fpmfm* rfpmfm * ( mtb/msq) * ( 1 - ( ( msd*msq*bb2)/ | |
910 | ( 2.0*mup*mtx*bbx2))); | |
911 | ||
912 | *fpf = (fppfm + fpmfm)/2.0; | |
913 | *fmf = (fppfm - fpmfm)/2.0; | |
914 | ||
915 | return; | |
916 | } //get_ff_isgw_1s0 | |
917 | ||
918 | ||
919 | ||
920 | ||
921 | ||
922 | void EvtISGW2FF::EvtISGW2FF3S1(EvtId parent,EvtId daugt,double t,double mass, | |
923 | double *f,double *g,double *ap,double *am){ | |
924 | ||
925 | //added by Lange Jan4,2000 | |
926 | static EvtId EM=EvtPDL::getId("e-"); | |
927 | static EvtId EP=EvtPDL::getId("e+"); | |
928 | static EvtId MUM=EvtPDL::getId("mu-"); | |
929 | static EvtId MUP=EvtPDL::getId("mu+"); | |
930 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
931 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
932 | ||
933 | static EvtId BP=EvtPDL::getId("B+"); | |
934 | static EvtId BM=EvtPDL::getId("B-"); | |
935 | static EvtId B0=EvtPDL::getId("B0"); | |
936 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
937 | ||
938 | static EvtId DST0=EvtPDL::getId("D*0"); | |
939 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
940 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
941 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
942 | static EvtId D0=EvtPDL::getId("D0"); | |
943 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
944 | static EvtId DP=EvtPDL::getId("D+"); | |
945 | static EvtId DM=EvtPDL::getId("D-"); | |
946 | ||
947 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
948 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
949 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
950 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
951 | ||
952 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
953 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
954 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
955 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
956 | ||
957 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
958 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
959 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
960 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
961 | ||
962 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
963 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
964 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
965 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
966 | ||
967 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
968 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
969 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
970 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
971 | ||
972 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
973 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
974 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
975 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
976 | ||
977 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
978 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
979 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
980 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
981 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
982 | ||
983 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
984 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
985 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
986 | ||
987 | static EvtId PIP=EvtPDL::getId("pi+"); | |
988 | static EvtId PIM=EvtPDL::getId("pi-"); | |
989 | static EvtId PI0=EvtPDL::getId("pi0"); | |
990 | ||
991 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
992 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
993 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
994 | ||
995 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
996 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
997 | static EvtId A20=EvtPDL::getId("a_20"); | |
998 | ||
999 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
1000 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
1001 | static EvtId A10=EvtPDL::getId("a_10"); | |
1002 | ||
1003 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
1004 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
1005 | static EvtId A00=EvtPDL::getId("a_00"); | |
1006 | ||
1007 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
1008 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
1009 | static EvtId B10=EvtPDL::getId("b_10"); | |
1010 | ||
1011 | static EvtId H1=EvtPDL::getId("h_1"); | |
1012 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
1013 | ||
1014 | static EvtId F1=EvtPDL::getId("f_1"); | |
1015 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
1016 | static EvtId F0=EvtPDL::getId("f_0"); | |
1017 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
1018 | static EvtId F2=EvtPDL::getId("f_2"); | |
1019 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
1020 | ||
1021 | static EvtId ETA=EvtPDL::getId("eta"); | |
1022 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
1023 | static EvtId OMEG=EvtPDL::getId("omega"); | |
1024 | ||
1025 | static EvtId KP=EvtPDL::getId("K+"); | |
1026 | static EvtId KM=EvtPDL::getId("K-"); | |
1027 | static EvtId K0=EvtPDL::getId("K0"); | |
1028 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
1029 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
1030 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
1031 | ||
1032 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
1033 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
1034 | static EvtId KST0=EvtPDL::getId("K*0"); | |
1035 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
1036 | ||
1037 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
1038 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
1039 | static EvtId K10=EvtPDL::getId("K_10"); | |
1040 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
1041 | ||
1042 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
1043 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
1044 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
1045 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
1046 | ||
1047 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
1048 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
1049 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
1050 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
1051 | ||
1052 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
1053 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
1054 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
1055 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
1056 | ||
1057 | static EvtId PHI=EvtPDL::getId("phi"); | |
1058 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
1059 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
1060 | ||
1061 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
1062 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
1063 | ||
1064 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
1065 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
1066 | ||
1067 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
1068 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
1069 | ||
1070 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
1071 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
1072 | ||
1073 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
1074 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
1075 | ||
1076 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
1077 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
1078 | ||
1079 | double cf(0.0),mtb,wt,msd(0.0),mup,f3f,msq(0.0),bb2(0.0),mum,mtx,bbx2,f3g; | |
1080 | double cji,bx2(0.0),f3appam,msb(0.0),tm,mbb(0.0),mbx(0.0); | |
1081 | double f3apmam,appam,apmam,mb,mx,f3; | |
1082 | double r_f,r_g,r_appam,r_apmam, betaji_f,betaji_g; | |
1083 | double betaji_appam, betaji_apmam; | |
1084 | double w,mqm,r2,chiji,zji,ai,al,rcji,nf(0.0),nfp(0.0),gammaji; | |
1085 | ||
1086 | EvtId prnt=parent; | |
1087 | EvtId dgt=daugt; | |
1088 | ||
1089 | if (parent==B0||parent==B0B||parent==BP||parent==BM) { | |
1090 | ||
1091 | msb=5.2; | |
1092 | msd=0.33; | |
1093 | bb2=0.431*0.431; | |
1094 | mbb=5.31; | |
1095 | nf = 4.0; | |
1096 | ||
1097 | if (dgt==DST0||dgt==DSTP||dgt==DSTM||dgt==DSTB) { | |
1098 | ||
1099 | cf=0.989; | |
1100 | msq=1.82; | |
1101 | bx2=0.38*0.38; | |
1102 | mbx=0.75*2.01+0.25*1.87; | |
1103 | nfp = 3.0; | |
1104 | } | |
1105 | else{ | |
1106 | if (dgt==OMEG||dgt==RHO0||dgt==RHOM||dgt==RHOP) { | |
1107 | ||
1108 | cf=0.905; | |
1109 | msq=0.33; | |
1110 | bx2=0.299*0.299; | |
1111 | mbx=0.75*0.770+0.25*0.14; | |
1112 | nfp = 0.0; | |
1113 | } | |
1114 | else{ | |
1115 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_3S1.\n"; | |
1116 | } | |
1117 | } | |
1118 | } | |
1119 | else{ | |
1120 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
1121 | ||
1122 | msb=1.82; | |
1123 | msd=0.33; | |
1124 | bb2=0.45*0.45; | |
1125 | mbb=1.963; | |
1126 | nf = 3.0; | |
1127 | ||
1128 | if (dgt==KSTP||dgt==KSTM||dgt==KST0||dgt==KSTB) { | |
1129 | ||
1130 | cf=0.928; | |
1131 | msq=0.55; | |
1132 | bx2=0.33*0.33; | |
1133 | mbx=0.75*0.892+0.25*0.494; | |
1134 | nfp = 2.0; | |
1135 | } | |
1136 | else{ | |
1137 | if (dgt==RHO0||dgt==OMEG||dgt==RHOM||dgt==RHOP) { | |
1138 | cf=0.889; | |
1139 | msq=0.33; | |
1140 | bx2=0.299*0.299; | |
1141 | mbx=0.75*0.770+0.25*0.14; | |
1142 | nfp = 0.0; | |
1143 | } | |
1144 | else{ | |
1145 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_3S1.\n"; | |
1146 | } | |
1147 | } | |
1148 | } | |
1149 | else{ | |
1150 | if (prnt==DSP||prnt==DSM){ | |
1151 | ||
1152 | msb=1.82; | |
1153 | msd=0.55; | |
1154 | bb2=0.56*0.56; | |
1155 | mbb=1.968; | |
1156 | nf = 3.0; | |
1157 | ||
1158 | if (dgt==KSTB||dgt==KST0) { | |
1159 | ||
1160 | cf=0.8731; | |
1161 | msq=0.55; | |
1162 | bx2=0.33*0.33; | |
1163 | mbx=0.87; | |
1164 | nfp = 2.0; | |
1165 | } | |
1166 | else{ | |
1167 | if(dgt==PHI){ | |
1168 | cf=0.911; | |
1169 | msq=0.55; | |
1170 | bx2=0.37*0.37; | |
1171 | mbx=0.97; | |
1172 | nfp = 2.0; | |
1173 | } | |
1174 | else{ | |
1175 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_3S1.\n"; | |
1176 | } | |
1177 | } | |
1178 | } | |
1179 | else{ | |
1180 | //BS -> cs constants added by djl on Jan. 21,1998 | |
1181 | if (prnt==BS0||prnt==BSB){ | |
1182 | ||
1183 | msb=5.2; | |
1184 | msd=0.55; | |
1185 | bb2=0.54*0.54; | |
1186 | mbb=5.38; | |
1187 | nf = 4.0; | |
1188 | ||
1189 | if (dgt==DSSTP||dgt==DSSTM) { | |
1190 | ||
1191 | cf=0.984; | |
1192 | msq=1.82; | |
1193 | bx2=0.49*0.49; | |
1194 | mbx=0.75*2.11+0.25*1.97; | |
1195 | nfp = 3.0; | |
1196 | } | |
1197 | else if (dgt==KSTP||dgt==KSTM||dgt==KST0||dgt==KSTB) { | |
1198 | ||
1199 | cf=0.928; | |
1200 | msq=0.55; | |
1201 | bx2=0.33*0.33; | |
1202 | mbx=0.75*0.892+0.25*0.494; | |
1203 | nfp = 2.0; | |
1204 | } | |
1205 | else{ | |
1206 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
1207 | } | |
1208 | } | |
1209 | ||
1210 | else{ | |
1211 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw2_ff_3S1.\n"; | |
1212 | } | |
1213 | } | |
1214 | } | |
1215 | } | |
1216 | ||
1217 | ||
1218 | mtb=msb+msd; | |
1219 | mtx=msq+msd; | |
1220 | ||
1221 | mup=1.0/(1.0/msq+1.0/msb); | |
1222 | mum=1.0/(1.0/msq-1.0/msb); | |
1223 | bbx2=0.5*(bb2+bx2); | |
1224 | mb=EvtPDL::getMeanMass(parent); | |
1225 | mx=mass; | |
1226 | tm=(mb-mx)*(mb-mx); | |
1227 | if ( t > tm ) t = 0.99*tm; | |
1228 | ||
1229 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
1230 | mqm = 0.1; | |
1231 | ||
1232 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
1233 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
1234 | log(EvtGetas(mqm,mqm)/EvtGetas(msq,msq)); | |
1235 | ||
1236 | w = 1.0 + (( tm - t ) / ( 2.0* mb * mx )); | |
1237 | ||
1238 | rcji = ( 1/sqrt(w*w -1 ))*log( w + sqrt( w*w -1 )); | |
1239 | ||
1240 | al = (8.0 / ( 33.0 - 2.0*nfp ))*(w*rcji -1.0 ); | |
1241 | ||
1242 | ai = -1.0* ( 6.0/( 33.0 - 2.0*nf)); | |
1243 | ||
1244 | cji = pow(( EvtGetas( msb,msb ) / EvtGetas( msq,msq ) ),ai); | |
1245 | zji = msq / msb; | |
1246 | ||
1247 | gammaji = EvtGetGammaji( zji ); | |
1248 | ||
1249 | chiji = -1.0 - ( gammaji / ( 1- zji )); | |
1250 | ||
1251 | betaji_g = (2.0/3.0)+gammaji; | |
1252 | betaji_f = (-2.0/3.0)+gammaji; | |
1253 | betaji_appam = -1.0-chiji+(4.0/(3.0*(1.0-zji)))+ | |
1254 | (2.0*(1+zji)*gammaji/(3.0*(1.0-zji)*(1.0-zji))); | |
1255 | ||
1256 | betaji_apmam = (1.0/3.0)-chiji-(4.0/(3.0*(1.0-zji)))- | |
1257 | (2.0*(1+zji)*gammaji/(3.0*(1.0-zji)*(1.0-zji)))+ | |
1258 | gammaji; | |
1259 | ||
1260 | r_g = cji*(1+(betaji_g*EvtGetas( msq,sqrt(mb*msq) )/(EvtConst::pi))); | |
1261 | r_f = cji*(1+(betaji_f*EvtGetas( msq,sqrt(mb*msq) )/(EvtConst::pi))); | |
1262 | r_appam = cji*(1+(betaji_appam*EvtGetas( msq,sqrt(mb*msq) )/(EvtConst::pi))); | |
1263 | r_apmam = cji*(1+(betaji_apmam*EvtGetas( msq,sqrt(mb*msq) )/(EvtConst::pi))); | |
1264 | ||
1265 | ||
1266 | f3=sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,1.5)/ | |
1267 | ((1.0+r2*(tm-t)/12.0)*(1.0+r2*(tm-t)/12.0)); | |
1268 | ||
1269 | f3f=sqrt(mbx*mbb/(mtx*mtb))*f3; | |
1270 | f3g=sqrt(mtx*mtb/(mbx*mbb))*f3; | |
1271 | f3appam=sqrt(mtb*mtb*mtb*mbx/(mbb*mbb*mbb*mtx))*f3; | |
1272 | f3apmam=sqrt(mtx*mtb/(mbx*mbb))*f3; | |
1273 | *f=cf*mtb*(1+wt+msd*(wt-1)/(2*mup))*f3f*r_f; | |
1274 | *g=0.5*(1/msq-msd*bb2/(2*mum*mtx*bbx2))*f3g*r_g; | |
1275 | ||
1276 | appam=cji*(msd*bx2*(1-msd*bx2/(2*mtb*bbx2))/ | |
1277 | ((1+wt)*msq*msb*bbx2)- | |
1278 | betaji_appam*EvtGetas( msq,sqrt(msq*mb) )/ | |
1279 | (mtb*EvtConst::pi))*f3appam; | |
1280 | ||
1281 | apmam=-1.0*(mtb/msb-msd*bx2/(2*mup*bbx2)+wt*msd*mtb*bx2* | |
1282 | (1-msd*bx2/(2*mtb*bbx2))/((wt+1)*msq*msb*bbx2))* | |
1283 | f3apmam*r_apmam/mtx; | |
1284 | ||
1285 | *ap=0.5*(appam+apmam); | |
1286 | *am=0.5*(appam-apmam); | |
1287 | return; | |
1288 | } | |
1289 | ||
1290 | ||
1291 | void EvtISGW2FF::EvtISGW2FF21S0 (EvtId parent,EvtId daugt, | |
1292 | double t, double mass, double *fppf, double *fpmf ) { | |
1293 | ||
1294 | //added by Lange Jan4,2000 | |
1295 | static EvtId EM=EvtPDL::getId("e-"); | |
1296 | static EvtId EP=EvtPDL::getId("e+"); | |
1297 | static EvtId MUM=EvtPDL::getId("mu-"); | |
1298 | static EvtId MUP=EvtPDL::getId("mu+"); | |
1299 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
1300 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
1301 | ||
1302 | static EvtId BP=EvtPDL::getId("B+"); | |
1303 | static EvtId BM=EvtPDL::getId("B-"); | |
1304 | static EvtId B0=EvtPDL::getId("B0"); | |
1305 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
1306 | ||
1307 | static EvtId DST0=EvtPDL::getId("D*0"); | |
1308 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
1309 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
1310 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
1311 | static EvtId D0=EvtPDL::getId("D0"); | |
1312 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
1313 | static EvtId DP=EvtPDL::getId("D+"); | |
1314 | static EvtId DM=EvtPDL::getId("D-"); | |
1315 | ||
1316 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
1317 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
1318 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
1319 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
1320 | ||
1321 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
1322 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
1323 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
1324 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
1325 | ||
1326 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
1327 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
1328 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
1329 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
1330 | ||
1331 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
1332 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
1333 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
1334 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
1335 | ||
1336 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
1337 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
1338 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
1339 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
1340 | ||
1341 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
1342 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
1343 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
1344 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
1345 | ||
1346 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
1347 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
1348 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
1349 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
1350 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
1351 | ||
1352 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
1353 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
1354 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
1355 | ||
1356 | static EvtId PIP=EvtPDL::getId("pi+"); | |
1357 | static EvtId PIM=EvtPDL::getId("pi-"); | |
1358 | static EvtId PI0=EvtPDL::getId("pi0"); | |
1359 | ||
1360 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
1361 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
1362 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
1363 | ||
1364 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
1365 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
1366 | static EvtId A20=EvtPDL::getId("a_20"); | |
1367 | ||
1368 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
1369 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
1370 | static EvtId A10=EvtPDL::getId("a_10"); | |
1371 | ||
1372 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
1373 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
1374 | static EvtId A00=EvtPDL::getId("a_00"); | |
1375 | ||
1376 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
1377 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
1378 | static EvtId B10=EvtPDL::getId("b_10"); | |
1379 | ||
1380 | static EvtId H1=EvtPDL::getId("h_1"); | |
1381 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
1382 | ||
1383 | static EvtId F1=EvtPDL::getId("f_1"); | |
1384 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
1385 | static EvtId F0=EvtPDL::getId("f_0"); | |
1386 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
1387 | static EvtId F2=EvtPDL::getId("f_2"); | |
1388 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
1389 | ||
1390 | static EvtId ETA=EvtPDL::getId("eta"); | |
1391 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
1392 | static EvtId OMEG=EvtPDL::getId("omega"); | |
1393 | ||
1394 | static EvtId KP=EvtPDL::getId("K+"); | |
1395 | static EvtId KM=EvtPDL::getId("K-"); | |
1396 | static EvtId K0=EvtPDL::getId("K0"); | |
1397 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
1398 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
1399 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
1400 | ||
1401 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
1402 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
1403 | static EvtId KST0=EvtPDL::getId("K*0"); | |
1404 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
1405 | ||
1406 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
1407 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
1408 | static EvtId K10=EvtPDL::getId("K_10"); | |
1409 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
1410 | ||
1411 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
1412 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
1413 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
1414 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
1415 | ||
1416 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
1417 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
1418 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
1419 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
1420 | ||
1421 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
1422 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
1423 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
1424 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
1425 | ||
1426 | static EvtId PHI=EvtPDL::getId("phi"); | |
1427 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
1428 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
1429 | ||
1430 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
1431 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
1432 | ||
1433 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
1434 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
1435 | ||
1436 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
1437 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
1438 | ||
1439 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
1440 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
1441 | ||
1442 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
1443 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
1444 | ||
1445 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
1446 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
1447 | ||
1448 | double mtb, mbb(0.0); | |
1449 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
1450 | double msq(0.0),bx2(0.0),mbx(0.0),mtx; | |
1451 | double f3fppfm,f3fpmfm,fppfm,fpmfm,f3; | |
1452 | double mqm,msb(0.0); | |
1453 | double mum,mup,r2,wt,tm,bb2(0.0),bbx2; | |
1454 | double tau,udef,vdef; | |
1455 | ||
1456 | EvtId prnt=parent; | |
1457 | EvtId dgt=daugt; | |
1458 | ||
1459 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
1460 | ||
1461 | msb=5.2; | |
1462 | msd=0.33; | |
1463 | bb2=0.431*0.431; | |
1464 | mbb=0.75*5.325+0.25*5.279; | |
1465 | nf = 4.0; | |
1466 | ||
1467 | if (dgt==PI2S0||dgt==PI2SP||dgt==PI2SM||dgt==ETA2S) { | |
1468 | ||
1469 | msq=0.33; | |
1470 | bx2=0.406*0.406; | |
1471 | mbx=0.75*1.45+0.25*1.300; | |
1472 | nfp = 0.0; | |
1473 | } | |
1474 | else{ | |
1475 | if (dgt==D21S0P||dgt==D21S0B||dgt==D21S0N||dgt==D21S00) { | |
1476 | msq=1.82; | |
1477 | bx2=0.45*0.45; | |
1478 | mbx=0.75*2.64+0.25*2.58; | |
1479 | nfp=3.0; | |
1480 | } | |
1481 | else{ | |
1482 | ||
1483 | report(ERROR,"EvtGen") << "Not implemented daugt in get_EvtISGW2_ff_21S0.\n"; | |
1484 | } | |
1485 | } | |
1486 | } | |
1487 | else{ | |
1488 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
1489 | msb=1.82; | |
1490 | msd=0.33; | |
1491 | bb2=0.45*0.45; | |
1492 | mbb=1.963; | |
1493 | nf = 3.0; | |
1494 | if (dgt==PI2SP||dgt==PI2SM||dgt==PI2S0||dgt==ETA2S) { | |
1495 | msq=0.33; | |
1496 | bx2=0.406*0.406; | |
1497 | mbx=0.75*1.45+0.25*1.300; | |
1498 | nfp = 0.0; | |
1499 | } | |
1500 | else{ | |
1501 | report(ERROR,"EvtGen") << "Not implemented daugt in get_EvtISGW2_ff_21S0.\n"; | |
1502 | } | |
1503 | } | |
1504 | else{ | |
1505 | report(ERROR,"EvtGen") << "Not implemented parent in get_EvtISGW2_ff_21S0.\n"; | |
1506 | } | |
1507 | } | |
1508 | ||
1509 | mtb = msb + msd; | |
1510 | mtx = msq + msd; | |
1511 | ||
1512 | mb = EvtPDL::getMeanMass( parent ); | |
1513 | mx = mass; | |
1514 | ||
1515 | mup=1.0/(1.0/msq+1.0/msb); | |
1516 | mum=1.0/(1.0/msq-1.0/msb); | |
1517 | bbx2=0.5*(bb2+bx2); | |
1518 | tm=(mb-mx)*(mb-mx); | |
1519 | if (t>tm) t = 0.99*tm; | |
1520 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
1521 | ||
1522 | mqm = 0.1; | |
1523 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
1524 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
1525 | log(EvtGetas(mqm)/EvtGetas(msq)); | |
1526 | ||
1527 | f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0) / | |
1528 | (pow((1.0+r2*(tm-t)/24.0),4.0)); | |
1529 | ||
1530 | f3fppfm = f3*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),0.5); | |
1531 | f3fpmfm = f3*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),-0.5); | |
1532 | ||
1533 | tau = msd*msd*bx2*(wt-1)/(bb2*bbx2); | |
1534 | udef = (( bb2-bx2)/(2.0*bbx2)) + ((bb2*tau)/(3.0*bbx2)); | |
1535 | vdef = (bb2*(1.0+(msq/msb))/(6.0*bbx2))*(7.0 - ((bb2/bbx2)*(5+tau))); | |
1536 | ||
1537 | fppfm = f3fppfm*sqrt(1.5)*((1.0-(msd/msq))*udef-(msd*vdef/msq)); | |
1538 | fpmfm = f3fpmfm*sqrt(1.5)*(mtb/msq)*(udef+(msd*vdef/mtx)); | |
1539 | ||
1540 | *fppf = (fppfm + fpmfm) /2.0; | |
1541 | *fpmf = (fppfm - fpmfm) /2.0; | |
1542 | return; | |
1543 | ||
1544 | } //get_ff_isgw_21s0 | |
1545 | ||
1546 | ||
1547 | void EvtISGW2FF::EvtISGW2FF23S1 (EvtId parent,EvtId daugt, | |
1548 | double t, double mass, double *fpf, double *gpf, | |
1549 | double *appf, double *apmf ) { | |
1550 | ||
1551 | //added by Lange Jan4,2000 | |
1552 | static EvtId EM=EvtPDL::getId("e-"); | |
1553 | static EvtId EP=EvtPDL::getId("e+"); | |
1554 | static EvtId MUM=EvtPDL::getId("mu-"); | |
1555 | static EvtId MUP=EvtPDL::getId("mu+"); | |
1556 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
1557 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
1558 | ||
1559 | static EvtId BP=EvtPDL::getId("B+"); | |
1560 | static EvtId BM=EvtPDL::getId("B-"); | |
1561 | static EvtId B0=EvtPDL::getId("B0"); | |
1562 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
1563 | ||
1564 | static EvtId DST0=EvtPDL::getId("D*0"); | |
1565 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
1566 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
1567 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
1568 | static EvtId D0=EvtPDL::getId("D0"); | |
1569 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
1570 | static EvtId DP=EvtPDL::getId("D+"); | |
1571 | static EvtId DM=EvtPDL::getId("D-"); | |
1572 | ||
1573 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
1574 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
1575 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
1576 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
1577 | ||
1578 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
1579 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
1580 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
1581 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
1582 | ||
1583 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
1584 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
1585 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
1586 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
1587 | ||
1588 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
1589 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
1590 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
1591 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
1592 | ||
1593 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
1594 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
1595 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
1596 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
1597 | ||
1598 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
1599 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
1600 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
1601 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
1602 | ||
1603 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
1604 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
1605 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
1606 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
1607 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
1608 | ||
1609 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
1610 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
1611 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
1612 | ||
1613 | static EvtId PIP=EvtPDL::getId("pi+"); | |
1614 | static EvtId PIM=EvtPDL::getId("pi-"); | |
1615 | static EvtId PI0=EvtPDL::getId("pi0"); | |
1616 | ||
1617 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
1618 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
1619 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
1620 | ||
1621 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
1622 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
1623 | static EvtId A20=EvtPDL::getId("a_20"); | |
1624 | ||
1625 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
1626 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
1627 | static EvtId A10=EvtPDL::getId("a_10"); | |
1628 | ||
1629 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
1630 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
1631 | static EvtId A00=EvtPDL::getId("a_00"); | |
1632 | ||
1633 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
1634 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
1635 | static EvtId B10=EvtPDL::getId("b_10"); | |
1636 | ||
1637 | static EvtId H1=EvtPDL::getId("h_1"); | |
1638 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
1639 | ||
1640 | static EvtId F1=EvtPDL::getId("f_1"); | |
1641 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
1642 | static EvtId F0=EvtPDL::getId("f_0"); | |
1643 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
1644 | static EvtId F2=EvtPDL::getId("f_2"); | |
1645 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
1646 | ||
1647 | static EvtId ETA=EvtPDL::getId("eta"); | |
1648 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
1649 | static EvtId OMEG=EvtPDL::getId("omega"); | |
1650 | ||
1651 | static EvtId KP=EvtPDL::getId("K+"); | |
1652 | static EvtId KM=EvtPDL::getId("K-"); | |
1653 | static EvtId K0=EvtPDL::getId("K0"); | |
1654 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
1655 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
1656 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
1657 | ||
1658 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
1659 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
1660 | static EvtId KST0=EvtPDL::getId("K*0"); | |
1661 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
1662 | ||
1663 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
1664 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
1665 | static EvtId K10=EvtPDL::getId("K_10"); | |
1666 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
1667 | ||
1668 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
1669 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
1670 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
1671 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
1672 | ||
1673 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
1674 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
1675 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
1676 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
1677 | ||
1678 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
1679 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
1680 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
1681 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
1682 | ||
1683 | static EvtId PHI=EvtPDL::getId("phi"); | |
1684 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
1685 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
1686 | ||
1687 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
1688 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
1689 | ||
1690 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
1691 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
1692 | ||
1693 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
1694 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
1695 | ||
1696 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
1697 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
1698 | ||
1699 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
1700 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
1701 | ||
1702 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
1703 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
1704 | ||
1705 | double mtb, mbb(0.0); | |
1706 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
1707 | double msq(0.0),bx2(0.0),mbx(0.0),mtx; | |
1708 | double f3appam,f3apmam,f3,appam,apmam,f3fp,f3gp; | |
1709 | double udef,tau,mum,mup,bb2(0.0),bbx2,tm,wt,mqm,r2,msb(0.0); | |
1710 | double cfp(0.0); | |
1711 | ||
1712 | EvtId prnt=parent; | |
1713 | EvtId dgt=daugt; | |
1714 | ||
1715 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
1716 | ||
1717 | msb=5.2; | |
1718 | msd=0.33; | |
1719 | bb2=0.431*0.431; | |
1720 | mbb=0.75*5.325+0.25*5.279; | |
1721 | nf = 4.0; | |
1722 | ||
1723 | if (dgt==RHO2S0||dgt==RHO2SP||dgt==RHO2SM||dgt==OMEG2S) { | |
1724 | ||
1725 | cfp=0.776; | |
1726 | msq=0.33; | |
1727 | bx2=0.299*0.299; | |
1728 | mbx=0.75*1.45+0.25*1.300; | |
1729 | nfp = 0.0; | |
1730 | ||
1731 | } | |
1732 | else{ | |
1733 | if (dgt==D23S1N||dgt==D23S1P||dgt==D23S1B||dgt==D23S10) { | |
1734 | cfp=0.929; | |
1735 | msq=1.82; | |
1736 | bx2=0.38*0.38; | |
1737 | mbx=0.75*2.64+0.25*2.58; | |
1738 | nfp=3.0; | |
1739 | } | |
1740 | else{ | |
1741 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_23P1.\n"; | |
1742 | } | |
1743 | } | |
1744 | } | |
1745 | else{ | |
1746 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
1747 | msb=1.82; | |
1748 | msd=0.33; | |
1749 | bb2=0.45*0.45; | |
1750 | mbb=1.963; | |
1751 | nf = 3.0; | |
1752 | ||
1753 | if (dgt==RHO2S0||dgt==RHO2SP||dgt==RHO2SM||dgt==OMEG2S) { | |
1754 | cfp=0.74; | |
1755 | msq=0.33; | |
1756 | bx2=0.299*0.299; | |
1757 | mbx=0.75*1.45+0.25*1.300; | |
1758 | nfp = 0.0; | |
1759 | } | |
1760 | else{ | |
1761 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_23P1.\n"; | |
1762 | } | |
1763 | } | |
1764 | else{ | |
1765 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_23P1.\n"; | |
1766 | } | |
1767 | } | |
1768 | ||
1769 | mtb = msb + msd; | |
1770 | mtx = msq + msd; | |
1771 | mb = EvtPDL::getMeanMass( parent ); | |
1772 | mx = mass; | |
1773 | ||
1774 | mup=1.0/(1.0/msq+1.0/msb); | |
1775 | mum=1.0/(1.0/msq-1.0/msb); | |
1776 | bbx2=0.5*(bb2+bx2); | |
1777 | tm=(mb-mx)*(mb-mx); | |
1778 | ||
1779 | if (t>tm) t = 0.99*tm; | |
1780 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
1781 | ||
1782 | mqm = 0.1; | |
1783 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
1784 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
1785 | log(EvtGetas(mqm)/EvtGetas(msq)); | |
1786 | ||
1787 | f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0) / | |
1788 | (pow((1.0+r2*(tm-t)/24.0),4.0)); | |
1789 | ||
1790 | f3fp = f3*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),0.5); | |
1791 | f3gp = f3*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
1792 | f3appam = f3*pow(( mbb / mtb ),-1.5)*pow((mbx/mtx),0.5); | |
1793 | f3apmam = f3*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
1794 | ||
1795 | tau = msd*msd*bx2*(wt-1.0)/(bb2*bbx2); | |
1796 | udef = (( bb2-bx2)/(2.0*bbx2)); | |
1797 | udef = udef + ((bb2*tau)/(3.0*bbx2)); | |
1798 | ||
1799 | *fpf = cfp*sqrt(1.5)*mtb*(1.0+wt)*udef*f3fp; | |
1800 | ||
1801 | *gpf = sqrt(3.0/8.0)*f3gp*(((1.0/msq)-((msd*bb2)/(2.0*mum*mtx*bbx2)))* | |
1802 | udef + ( (msd*bb2*bx2)/(3.0*mum*mtx*bbx2*bbx2))); | |
1803 | ||
1804 | appam = f3appam*sqrt(2.0/3.0)*(bb2/(msq*msb*bbx2))*((-7.0*msd*msd*bx2* | |
1805 | bx2*(1.0+(tau/7.0))/(8.0*mtb*bbx2*bbx2))+(5.0*msd*bx2*(1.0+ | |
1806 | (tau/5.0))/(4.0*bbx2))+(3.0*msd*msd*bx2*bx2/(8.0*mtb*bb2*bbx2))- | |
1807 | (3.0*msd*bx2/(4.0*bb2))); | |
1808 | ||
1809 | apmam = f3apmam*sqrt(3.0/2.0)*(mtb/(msb*mtx))*(1.0-(bb2*(1.0+(tau/7.0))/ | |
1810 | bbx2)-(msd*bx2*(1.0-(5.0*bb2*(1.0+(tau/5.0))/(3.0*bbx2))) | |
1811 | /(2.0*mtb*bbx2))-(7.0*msd*msd*bb2*bx2/(12.0*msq*mtb*bbx2*bbx2))* | |
1812 | (1.0-(bx2/bbx2)+(bb2*tau/(7.0*bbx2)))); | |
1813 | ||
1814 | *appf = (appam + apmam) /2.0; | |
1815 | *apmf = (appam - apmam) /2.0; | |
1816 | return; | |
1817 | } //get_ff_isgw_23s1 | |
1818 | ||
1819 | void EvtISGW2FF::EvtISGW2FF1P1 (EvtId parent,EvtId daugt, | |
1820 | double t, double mass, double *rf, double *vf, | |
1821 | double *spf, double *smf ) { | |
1822 | ||
1823 | //added by Lange Jan4,2000 | |
1824 | static EvtId EM=EvtPDL::getId("e-"); | |
1825 | static EvtId EP=EvtPDL::getId("e+"); | |
1826 | static EvtId MUM=EvtPDL::getId("mu-"); | |
1827 | static EvtId MUP=EvtPDL::getId("mu+"); | |
1828 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
1829 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
1830 | ||
1831 | static EvtId BP=EvtPDL::getId("B+"); | |
1832 | static EvtId BM=EvtPDL::getId("B-"); | |
1833 | static EvtId B0=EvtPDL::getId("B0"); | |
1834 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
1835 | ||
1836 | static EvtId DST0=EvtPDL::getId("D*0"); | |
1837 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
1838 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
1839 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
1840 | static EvtId D0=EvtPDL::getId("D0"); | |
1841 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
1842 | static EvtId DP=EvtPDL::getId("D+"); | |
1843 | static EvtId DM=EvtPDL::getId("D-"); | |
1844 | ||
1845 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
1846 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
1847 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
1848 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
1849 | ||
1850 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
1851 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
1852 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
1853 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
1854 | ||
1855 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
1856 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
1857 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
1858 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
1859 | ||
1860 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
1861 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
1862 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
1863 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
1864 | ||
1865 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
1866 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
1867 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
1868 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
1869 | ||
1870 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
1871 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
1872 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
1873 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
1874 | ||
1875 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
1876 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
1877 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
1878 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
1879 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
1880 | ||
1881 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
1882 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
1883 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
1884 | ||
1885 | static EvtId PIP=EvtPDL::getId("pi+"); | |
1886 | static EvtId PIM=EvtPDL::getId("pi-"); | |
1887 | static EvtId PI0=EvtPDL::getId("pi0"); | |
1888 | ||
1889 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
1890 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
1891 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
1892 | ||
1893 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
1894 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
1895 | static EvtId A20=EvtPDL::getId("a_20"); | |
1896 | ||
1897 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
1898 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
1899 | static EvtId A10=EvtPDL::getId("a_10"); | |
1900 | ||
1901 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
1902 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
1903 | static EvtId A00=EvtPDL::getId("a_00"); | |
1904 | ||
1905 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
1906 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
1907 | static EvtId B10=EvtPDL::getId("b_10"); | |
1908 | ||
1909 | static EvtId H1=EvtPDL::getId("h_1"); | |
1910 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
1911 | ||
1912 | static EvtId F1=EvtPDL::getId("f_1"); | |
1913 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
1914 | static EvtId F0=EvtPDL::getId("f_0"); | |
1915 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
1916 | static EvtId F2=EvtPDL::getId("f_2"); | |
1917 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
1918 | ||
1919 | static EvtId ETA=EvtPDL::getId("eta"); | |
1920 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
1921 | static EvtId OMEG=EvtPDL::getId("omega"); | |
1922 | ||
1923 | static EvtId KP=EvtPDL::getId("K+"); | |
1924 | static EvtId KM=EvtPDL::getId("K-"); | |
1925 | static EvtId K0=EvtPDL::getId("K0"); | |
1926 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
1927 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
1928 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
1929 | ||
1930 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
1931 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
1932 | static EvtId KST0=EvtPDL::getId("K*0"); | |
1933 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
1934 | ||
1935 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
1936 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
1937 | static EvtId K10=EvtPDL::getId("K_10"); | |
1938 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
1939 | ||
1940 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
1941 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
1942 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
1943 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
1944 | ||
1945 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
1946 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
1947 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
1948 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
1949 | ||
1950 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
1951 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
1952 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
1953 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
1954 | ||
1955 | static EvtId PHI=EvtPDL::getId("phi"); | |
1956 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
1957 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
1958 | ||
1959 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
1960 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
1961 | ||
1962 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
1963 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
1964 | ||
1965 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
1966 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
1967 | ||
1968 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
1969 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
1970 | ||
1971 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
1972 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
1973 | ||
1974 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
1975 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
1976 | ||
1977 | double mtb, mbb(0.0); | |
1978 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
1979 | double msq(0.0),bx2(0.0),mbx(0.0),mtx,f5; | |
1980 | double f5sppsm,f5spmsm; | |
1981 | double f5v,f5r,mup,mum,vv,rr,spmsm,sppsm; | |
1982 | double mqm,msb(0.0),bb2(0.0),bbx2,tm,wt,r2; | |
1983 | EvtId prnt=parent; | |
1984 | EvtId dgt=daugt; | |
1985 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
1986 | msb=5.2; | |
1987 | msd=0.33; | |
1988 | bb2=0.431*0.431; | |
1989 | mbb=5.31; | |
1990 | nf = 4.0; | |
1991 | if (dgt==B10||dgt==B1P||dgt==B1M||dgt==H1||dgt==H1PR) { | |
1992 | msq=0.33; | |
1993 | bx2=0.275*0.275; | |
1994 | mbx=(3.0*1.123+0.98+5.0*1.32+3.0*1.26)/12.0; | |
1995 | nfp = 0.0; | |
1996 | } | |
1997 | else{ | |
1998 | if (dgt==D1P1P||dgt==D1P10||dgt==D1P1B||dgt==D1P1N) { | |
1999 | msq=1.82; | |
2000 | bx2=0.33*0.33; | |
2001 | mbx=(5.0*2.46+3.0*2.42)/8.0; | |
2002 | nfp = 3.0; | |
2003 | } | |
2004 | else{ | |
2005 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_1P1.\n"; | |
2006 | } | |
2007 | } | |
2008 | } | |
2009 | else{ | |
2010 | if (prnt==DM||prnt==DP||prnt==D0B||prnt==D0) { | |
2011 | msb=1.82; | |
2012 | msd=0.33; | |
2013 | bb2=0.45*0.45; | |
2014 | mbb=1.963; | |
2015 | nf = 3.0; | |
2016 | if (dgt==B10||dgt==B1P||dgt==B1M||dgt==H1||dgt==H1PR) { | |
2017 | msq=0.33; | |
2018 | bx2=0.275*0.275; | |
2019 | mbx=(3.0*1.123+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2020 | nfp = 0.0; | |
2021 | } | |
2022 | else{ | |
2023 | if (dgt==K10||dgt==K1B||dgt==K1P||dgt==K1M) { | |
2024 | msq=0.55; | |
2025 | bx2=0.30*0.30; | |
2026 | mbx=(3.0*1.27+1.43+5.0*1.43+3.0*1.4)/12.0; | |
2027 | nfp = 2.0; | |
2028 | } | |
2029 | else{ | |
2030 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_1P1.\n"; | |
2031 | } | |
2032 | } | |
2033 | } | |
2034 | else{ | |
2035 | //BS -> cs constants added by djl on Jan. 21,1998 | |
2036 | if (prnt==BS0||prnt==BSB){ | |
2037 | ||
2038 | msb=5.2; | |
2039 | msd=0.55; | |
2040 | bb2=0.54*0.54; | |
2041 | mbb=5.38; | |
2042 | nf = 4.0; | |
2043 | ||
2044 | if (dgt==D1P1SP||dgt==D1P1SN) { | |
2045 | ||
2046 | msq=1.82; | |
2047 | bx2=0.41*0.41; | |
2048 | mbx=(5.0*2.61+3.0*2.54)/8.0; | |
2049 | nfp = 3.0; | |
2050 | } | |
2051 | else if (dgt==K10||dgt==K1B||dgt==K1P||dgt==K1M) { | |
2052 | msq=0.55; | |
2053 | bx2=0.30*0.30; | |
2054 | mbx=(3.0*1.27+1.43+5.0*1.43+3.0*1.4)/12.0; | |
2055 | nfp = 2.0; | |
2056 | } | |
2057 | else{ | |
2058 | report(ERROR,"EvtGen") << "Not implemented daugt:" | |
2059 | <<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
2060 | } | |
2061 | } | |
2062 | ||
2063 | else{ | |
2064 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_1P1.\n"; | |
2065 | } | |
2066 | } | |
2067 | } | |
2068 | ||
2069 | ||
2070 | mtb = msb + msd; | |
2071 | mtx = msq + msd; | |
2072 | ||
2073 | mb = EvtPDL::getMeanMass( parent ); | |
2074 | mx = mass; | |
2075 | ||
2076 | mup=1.0/(1.0/msq+1.0/msb); | |
2077 | mum=1.0/(1.0/msq-1.0/msb); | |
2078 | bbx2=0.5*(bb2+bx2); | |
2079 | tm=(mb-mx)*(mb-mx); | |
2080 | if (t>tm) t = 0.99*tm; | |
2081 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
2082 | ||
2083 | mqm = 0.1; | |
2084 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2)+ | |
2085 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
2086 | log(EvtGetas(mqm,mqm)/EvtGetas(msq,msq)); | |
2087 | ||
2088 | f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0) / | |
2089 | (pow((1.0+r2*(tm-t)/18.0),3.0)); | |
2090 | ||
2091 | f5v = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
2092 | f5r = f5*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),0.5); | |
2093 | f5sppsm = f5*pow(( mbb / mtb ),-1.5)*pow((mbx/mtx),0.5); | |
2094 | f5spmsm = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
2095 | ||
2096 | if (msq == msd) { | |
2097 | vv = f5v*(((mtb*sqrt(bb2))/(4.0*sqrt(2.0)*msb*msq*mtx)) + | |
2098 | (((wt-1)*msd)/(6.0*sqrt(2.0*bb2)*mtx))); | |
2099 | ||
2100 | rr = f5r*mtb*sqrt(bb2/2)*((1.0/mup)+((msd*mtx*(wt-1)*(wt-1))/ | |
2101 | (3.0*msq*bb2))); | |
2102 | ||
2103 | sppsm = msd*f5sppsm/(sqrt(2.0*bb2)*mtb)*(1.0-(msd/msq)+((msd*bb2)/ | |
2104 | (2.0*mup*bbx2))); | |
2105 | ||
2106 | spmsm = msd*f5spmsm/(sqrt(2.0*bb2)*msq)*(((4-wt)/3.0)- ( (msd*msq*bb2)/ | |
2107 | (2.0*mtx*mup*bbx2))); | |
2108 | ||
2109 | } else { | |
2110 | vv = -1.0*msd*f5v/(2.0*sqrt(3.0*bb2)*mtx)* | |
2111 | ((wt+1)/2.0+bb2*mtb/(2.0*msd*msq*msb)); | |
2112 | ||
2113 | rr = -2.0*mtb*sqrt(bb2/3.0)*f5r*(1.0/msq + mtx*msd*(wt-1)/(2.0*bb2)* | |
2114 | ((wt+1)/(2.0*msq)-msd*bb2/(2.0*mum*mtx*bbx2))); | |
2115 | ||
2116 | sppsm = -1.0*sqrt(3.0)*msd*f5sppsm/(2.0*sqrt(bb2)*mtb)*(1 - msd/(3.0*msq) - | |
2117 | msd*bb2/(3.0*bbx2)*(1.0/(2.0*mum)-1.0/mup)); | |
2118 | ||
2119 | spmsm = -1.0*msd*f5spmsm/(2.0*sqrt(3.0*bb2)*mtx)*((2-wt)*mtx/msq + | |
2120 | msd*bb2/bbx2*(1.0/(2.0*mum)-1.0/mup)); | |
2121 | ||
2122 | } | |
2123 | ||
2124 | //smooth out the mass(meson) dependence a little | |
2125 | double parMass=EvtPDL::getMeanMass(prnt); | |
2126 | double q2max = parMass*parMass + mass*mass - 2.0*parMass*mass; | |
2127 | double massNom= EvtPDL::getMeanMass(dgt); | |
2128 | double q2maxNom = parMass*parMass + massNom*massNom - 2.0*parMass*massNom; | |
2129 | double q2maxin=sqrt(q2maxNom/q2max); | |
2130 | if ( q2maxin > 1000. ) q2maxin=1000.; | |
2131 | ||
2132 | vv*=q2maxin; | |
2133 | rr*=q2maxin; | |
2134 | sppsm*=q2maxin; | |
2135 | spmsm*=q2maxin; | |
2136 | ||
2137 | *vf = vv; | |
2138 | *rf = rr; | |
2139 | *spf = (sppsm + spmsm)/2.0; | |
2140 | *smf = (sppsm - spmsm)/2.0; | |
2141 | return; | |
2142 | } //get_ff_isgw_1p1 | |
2143 | ||
2144 | ||
2145 | void EvtISGW2FF::EvtISGW2FF3P1 (EvtId parent,EvtId daugt, | |
2146 | double t, double mass, double *lf, double *qf, | |
2147 | double *cpf, double *cmf ) { | |
2148 | ||
2149 | //added by Lange Jan4,2000 | |
2150 | static EvtId EM=EvtPDL::getId("e-"); | |
2151 | static EvtId EP=EvtPDL::getId("e+"); | |
2152 | static EvtId MUM=EvtPDL::getId("mu-"); | |
2153 | static EvtId MUP=EvtPDL::getId("mu+"); | |
2154 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
2155 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
2156 | ||
2157 | static EvtId BP=EvtPDL::getId("B+"); | |
2158 | static EvtId BM=EvtPDL::getId("B-"); | |
2159 | static EvtId B0=EvtPDL::getId("B0"); | |
2160 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
2161 | ||
2162 | static EvtId DST0=EvtPDL::getId("D*0"); | |
2163 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
2164 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
2165 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
2166 | static EvtId D0=EvtPDL::getId("D0"); | |
2167 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
2168 | static EvtId DP=EvtPDL::getId("D+"); | |
2169 | static EvtId DM=EvtPDL::getId("D-"); | |
2170 | ||
2171 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
2172 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
2173 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
2174 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
2175 | ||
2176 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
2177 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
2178 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
2179 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
2180 | ||
2181 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
2182 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
2183 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
2184 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
2185 | ||
2186 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
2187 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
2188 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
2189 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
2190 | ||
2191 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
2192 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
2193 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
2194 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
2195 | ||
2196 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
2197 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
2198 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
2199 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
2200 | ||
2201 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
2202 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
2203 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
2204 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
2205 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
2206 | ||
2207 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
2208 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
2209 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
2210 | ||
2211 | static EvtId PIP=EvtPDL::getId("pi+"); | |
2212 | static EvtId PIM=EvtPDL::getId("pi-"); | |
2213 | static EvtId PI0=EvtPDL::getId("pi0"); | |
2214 | ||
2215 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
2216 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
2217 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
2218 | ||
2219 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
2220 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
2221 | static EvtId A20=EvtPDL::getId("a_20"); | |
2222 | ||
2223 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
2224 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
2225 | static EvtId A10=EvtPDL::getId("a_10"); | |
2226 | ||
2227 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
2228 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
2229 | static EvtId A00=EvtPDL::getId("a_00"); | |
2230 | ||
2231 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
2232 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
2233 | static EvtId B10=EvtPDL::getId("b_10"); | |
2234 | ||
2235 | static EvtId H1=EvtPDL::getId("h_1"); | |
2236 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
2237 | ||
2238 | static EvtId F1=EvtPDL::getId("f_1"); | |
2239 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
2240 | static EvtId F0=EvtPDL::getId("f_0"); | |
2241 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
2242 | static EvtId F2=EvtPDL::getId("f_2"); | |
2243 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
2244 | ||
2245 | static EvtId ETA=EvtPDL::getId("eta"); | |
2246 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
2247 | static EvtId OMEG=EvtPDL::getId("omega"); | |
2248 | ||
2249 | static EvtId KP=EvtPDL::getId("K+"); | |
2250 | static EvtId KM=EvtPDL::getId("K-"); | |
2251 | static EvtId K0=EvtPDL::getId("K0"); | |
2252 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
2253 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
2254 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
2255 | ||
2256 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
2257 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
2258 | static EvtId KST0=EvtPDL::getId("K*0"); | |
2259 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
2260 | ||
2261 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
2262 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
2263 | static EvtId K10=EvtPDL::getId("K_10"); | |
2264 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
2265 | ||
2266 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
2267 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
2268 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
2269 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
2270 | ||
2271 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
2272 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
2273 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
2274 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
2275 | ||
2276 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
2277 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
2278 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
2279 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
2280 | ||
2281 | static EvtId PHI=EvtPDL::getId("phi"); | |
2282 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
2283 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
2284 | ||
2285 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
2286 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
2287 | ||
2288 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
2289 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
2290 | ||
2291 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
2292 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
2293 | ||
2294 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
2295 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
2296 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
2297 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
2298 | ||
2299 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
2300 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
2301 | ||
2302 | double mtb, mbb(0.0); | |
2303 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
2304 | double msq(0.0),bx2(0.0),mbx(0.0),mtx; | |
2305 | double f5cppcm,f5cpmcm,f5,ql,ll,cppcm,cpmcm,f5q,f5l; | |
2306 | double mqm,msb(0.0),bb2(0.0),mup,mum,bbx2,tm,wt,r2; | |
2307 | EvtId prnt=parent; | |
2308 | EvtId dgt=daugt; | |
2309 | ||
2310 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
2311 | ||
2312 | msb=5.2; | |
2313 | msd=0.33; | |
2314 | bb2=0.431*0.431; | |
2315 | mbb=5.31; | |
2316 | nf = 4.0; | |
2317 | ||
2318 | if (dgt==A10||dgt==A1P||dgt==A1M||dgt==F1||dgt==F1PR) { | |
2319 | ||
2320 | msq=0.33; | |
2321 | bx2=0.275*0.275; | |
2322 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2323 | nfp = 0.0; | |
2324 | } | |
2325 | else{ | |
2326 | if (dgt==D3P1P||dgt==D3P1N||dgt==D3P10||dgt==D3P1B) { | |
2327 | msq=1.82; | |
2328 | bx2=0.33*0.33; | |
2329 | mbx=(3.0*2.49+2.40)/4.0; | |
2330 | nfp = 3.0; | |
2331 | } | |
2332 | else{ | |
2333 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_3P1.\n"; | |
2334 | } | |
2335 | } | |
2336 | } | |
2337 | else{ | |
2338 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
2339 | ||
2340 | msb=1.82; | |
2341 | msd=0.33; | |
2342 | bb2=0.45*0.45; | |
2343 | mbb=1.963; | |
2344 | nf = 3.0; | |
2345 | ||
2346 | if (dgt==F1||dgt==F1PR||dgt==A10||dgt==A1P||dgt==A1M) { | |
2347 | ||
2348 | msq=0.33; | |
2349 | bx2=0.275*0.275; | |
2350 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2351 | nfp = 0.0; | |
2352 | } | |
2353 | else{ | |
2354 | if (dgt==K1STM||dgt==K1STB||dgt==K1STP||dgt==K1ST0) { | |
2355 | msq=0.55; | |
2356 | bx2=0.30*0.30; | |
2357 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
2358 | nfp = 2.0; | |
2359 | } | |
2360 | else{ | |
2361 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_3P1.\n"; | |
2362 | } | |
2363 | } | |
2364 | } | |
2365 | else{ | |
2366 | //BS -> cs constants added by djl on Jan. 21,1998 | |
2367 | if (prnt==BS0||prnt==BSB){ | |
2368 | ||
2369 | msb=5.2; | |
2370 | msd=0.55; | |
2371 | bb2=0.54*0.54; | |
2372 | mbb=5.38; | |
2373 | nf = 4.0; | |
2374 | ||
2375 | if (dgt==D3P1SP||dgt==D3P1SN) { | |
2376 | ||
2377 | msq=1.82; | |
2378 | bx2=0.41*0.41; | |
2379 | mbx=(3.0*2.54+2.46)/4.0; | |
2380 | nfp = 3.0; | |
2381 | } | |
2382 | else if (dgt==K1STM||dgt==K1STB||dgt==K1STP||dgt==K1ST0) { | |
2383 | msq=0.55; | |
2384 | bx2=0.30*0.30; | |
2385 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
2386 | nfp = 2.0; | |
2387 | } | |
2388 | else{ | |
2389 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
2390 | } | |
2391 | } | |
2392 | ||
2393 | else{ | |
2394 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P1.\n"; | |
2395 | } | |
2396 | } | |
2397 | } | |
2398 | ||
2399 | ||
2400 | ||
2401 | mtb = msb + msd; | |
2402 | mtx = msq + msd; | |
2403 | ||
2404 | mb = EvtPDL::getMeanMass( parent ); | |
2405 | mx = mass; | |
2406 | ||
2407 | mup=1.0/(1.0/msq+1.0/msb); | |
2408 | mum=1.0/(1.0/msq-1.0/msb); | |
2409 | bbx2=0.5*(bb2+bx2); | |
2410 | tm=(mb-mx)*(mb-mx); | |
2411 | if (t>tm) t = 0.99*tm; | |
2412 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
2413 | ||
2414 | mqm = 0.1; | |
2415 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
2416 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
2417 | log(EvtGetas(mqm)/EvtGetas(msq)); | |
2418 | ||
2419 | f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0) / | |
2420 | (pow((1.0+r2*(tm-t)/18.0),3.0)); | |
2421 | ||
2422 | f5q = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
2423 | f5l = f5*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),0.5); | |
2424 | f5cppcm = f5*pow(( mbb / mtb ),-1.5)*pow((mbx/mtx),0.5); | |
2425 | f5cpmcm = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),-0.5); | |
2426 | ||
2427 | if (msq == msd) { | |
2428 | ||
2429 | ql = -1.0*(msd*(5.0+wt)*f5q/(2.0*mtx*sqrt(bb2)*6.0)); | |
2430 | ||
2431 | ll = -1.0*mtb*sqrt(bb2)*f5l*(1/mum+ ( (msd*mtx*(wt-1)/bb2)* | |
2432 | ( (5.0+wt)/(6.0*msq)-(msd*bb2)/(2.0*mum*mtx*bbx2)))); | |
2433 | ||
2434 | cppcm = (-1.0*(msd*mtx*f5cppcm/(2.0*msq*mtb*sqrt(bb2)))* | |
2435 | (1-(msd*msq*bb2)/(2.0*mtx*mum*bbx2))); | |
2436 | ||
2437 | cpmcm = 1.0*(msd*mtx*f5cpmcm/(2.0*msq*mtb*sqrt(bb2)))* | |
2438 | (((wt+2.0)/3.0)-(msd*msq*bb2)/(2.0*mtx*mum*bbx2)) | |
2439 | *(mtb/mtx); | |
2440 | } else { | |
2441 | ||
2442 | ql = f5q*sqrt(1.0/6.0)*msd/(sqrt(bb2)*mtx)* | |
2443 | (1.0-bb2*mtb/(4.0*msd*msq*msb)); | |
2444 | ll = f5l*sqrt(2.0/3.0)*mtb*sqrt(bb2)*(1.0/(2.0*msq) - 3.0/(2.0*msb) + | |
2445 | msd*mtx*(wt-1)/bb2*(1.0/msq-msd*bb2/(2.0*mum*mtx*bbx2))); | |
2446 | cppcm = msd*msd*bx2*f5cppcm/(sqrt(6.0)*mtb*msq*sqrt(bb2)*bbx2); | |
2447 | cpmcm = -sqrt(2.0/3.0)*msd*f5cpmcm/(sqrt(bb2)*mtx)* | |
2448 | (1+msd*bx2/(2.0*msq*bbx2)); | |
2449 | } | |
2450 | ||
2451 | //smooth out the mass(meson) dependence a little | |
2452 | double parMass=EvtPDL::getMeanMass(prnt); | |
2453 | double q2max = parMass*parMass + mass*mass - 2.0*parMass*mass; | |
2454 | double massNom= EvtPDL::getMeanMass(dgt); | |
2455 | double q2maxNom = parMass*parMass + massNom*massNom - 2.0*parMass*massNom; | |
2456 | double q2maxin=sqrt(q2maxNom/q2max); | |
2457 | if ( q2maxin > 1000. ) q2maxin=1000.; | |
2458 | ql*=q2maxin; | |
2459 | ll*=q2maxin; | |
2460 | cppcm*=q2maxin; | |
2461 | cpmcm*=q2maxin; | |
2462 | ||
2463 | *qf = ql; | |
2464 | *lf = ll; | |
2465 | *cpf = (cppcm + cpmcm)/2.0; | |
2466 | *cmf = (cppcm - cpmcm)/2.0; | |
2467 | return; | |
2468 | } //get_ff_isgw_3p1 | |
2469 | ||
2470 | ||
2471 | void EvtISGW2FF::EvtISGW2FF3P0 (EvtId parent,EvtId daugt, | |
2472 | double t, double mass, double *upf, double *umf ) { | |
2473 | ||
2474 | //added by Lange Jan4,2000 | |
2475 | static EvtId EM=EvtPDL::getId("e-"); | |
2476 | static EvtId EP=EvtPDL::getId("e+"); | |
2477 | static EvtId MUM=EvtPDL::getId("mu-"); | |
2478 | static EvtId MUP=EvtPDL::getId("mu+"); | |
2479 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
2480 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
2481 | ||
2482 | static EvtId BP=EvtPDL::getId("B+"); | |
2483 | static EvtId BM=EvtPDL::getId("B-"); | |
2484 | static EvtId B0=EvtPDL::getId("B0"); | |
2485 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
2486 | ||
2487 | static EvtId DST0=EvtPDL::getId("D*0"); | |
2488 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
2489 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
2490 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
2491 | static EvtId D0=EvtPDL::getId("D0"); | |
2492 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
2493 | static EvtId DP=EvtPDL::getId("D+"); | |
2494 | static EvtId DM=EvtPDL::getId("D-"); | |
2495 | ||
2496 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
2497 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
2498 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
2499 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
2500 | ||
2501 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
2502 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
2503 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
2504 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
2505 | ||
2506 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
2507 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
2508 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
2509 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
2510 | ||
2511 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
2512 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
2513 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
2514 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
2515 | ||
2516 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
2517 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
2518 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
2519 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
2520 | ||
2521 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
2522 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
2523 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
2524 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
2525 | ||
2526 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
2527 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
2528 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
2529 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
2530 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
2531 | ||
2532 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
2533 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
2534 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
2535 | ||
2536 | static EvtId PIP=EvtPDL::getId("pi+"); | |
2537 | static EvtId PIM=EvtPDL::getId("pi-"); | |
2538 | static EvtId PI0=EvtPDL::getId("pi0"); | |
2539 | ||
2540 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
2541 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
2542 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
2543 | ||
2544 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
2545 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
2546 | static EvtId A20=EvtPDL::getId("a_20"); | |
2547 | ||
2548 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
2549 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
2550 | static EvtId A10=EvtPDL::getId("a_10"); | |
2551 | ||
2552 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
2553 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
2554 | static EvtId A00=EvtPDL::getId("a_00"); | |
2555 | ||
2556 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
2557 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
2558 | static EvtId B10=EvtPDL::getId("b_10"); | |
2559 | ||
2560 | static EvtId H1=EvtPDL::getId("h_1"); | |
2561 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
2562 | ||
2563 | static EvtId F1=EvtPDL::getId("f_1"); | |
2564 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
2565 | static EvtId F0=EvtPDL::getId("f_0"); | |
2566 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
2567 | static EvtId F2=EvtPDL::getId("f_2"); | |
2568 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
2569 | ||
2570 | static EvtId ETA=EvtPDL::getId("eta"); | |
2571 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
2572 | static EvtId OMEG=EvtPDL::getId("omega"); | |
2573 | ||
2574 | static EvtId KP=EvtPDL::getId("K+"); | |
2575 | static EvtId KM=EvtPDL::getId("K-"); | |
2576 | static EvtId K0=EvtPDL::getId("K0"); | |
2577 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
2578 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
2579 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
2580 | ||
2581 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
2582 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
2583 | static EvtId KST0=EvtPDL::getId("K*0"); | |
2584 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
2585 | ||
2586 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
2587 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
2588 | static EvtId K10=EvtPDL::getId("K_10"); | |
2589 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
2590 | ||
2591 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
2592 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
2593 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
2594 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
2595 | ||
2596 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
2597 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
2598 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
2599 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
2600 | ||
2601 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
2602 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
2603 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
2604 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
2605 | ||
2606 | static EvtId PHI=EvtPDL::getId("phi"); | |
2607 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
2608 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
2609 | ||
2610 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
2611 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
2612 | ||
2613 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
2614 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
2615 | ||
2616 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
2617 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
2618 | ||
2619 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
2620 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
2621 | ||
2622 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
2623 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
2624 | ||
2625 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
2626 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
2627 | ||
2628 | double mtb, mbb(0.0); | |
2629 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
2630 | double msq(0.0),bx2(0.0),mbx(0.0),mtx; | |
2631 | double f5uppum,f5upmum,uppum,upmum,f5; | |
2632 | double mqm,mum,mup,wt,r2,bb2(0.0),bbx2,msb(0.0),tm; | |
2633 | ||
2634 | EvtId prnt=parent; | |
2635 | EvtId dgt=daugt; | |
2636 | ||
2637 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
2638 | ||
2639 | msb=5.2; | |
2640 | msd=0.33; | |
2641 | bb2=0.431*0.431; | |
2642 | mbb=5.31; | |
2643 | nf = 4.0; | |
2644 | if (dgt==A00||dgt==A0P||dgt==A0M||dgt==F0||dgt==F0PR) { | |
2645 | ||
2646 | msq=0.33; | |
2647 | bx2=0.275*0.275; | |
2648 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2649 | nfp = 0.0; | |
2650 | } | |
2651 | else{ | |
2652 | if (dgt==D3P0P||dgt==D3P0N||dgt==D3P00||dgt==D3P0B) { | |
2653 | msq=1.82; | |
2654 | bx2=0.33*0.33; | |
2655 | mbx=(3.0*2.49+2.40)/4.0; | |
2656 | nfp = 3.0; | |
2657 | } | |
2658 | else{ | |
2659 | report(ERROR,"EvtGen") << "Not implemented daugt in get_EvtISGW2_ff_3P0.\n"; | |
2660 | } | |
2661 | } | |
2662 | } | |
2663 | else{ | |
2664 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
2665 | ||
2666 | msb=1.82; | |
2667 | msd=0.33; | |
2668 | bb2=0.45*0.45; | |
2669 | mbb=1.963; | |
2670 | nf = 3.0; | |
2671 | if (dgt==F0||dgt==F0PR||dgt==A00||dgt==A0P||dgt==A0M) { | |
2672 | msq=0.33; | |
2673 | bx2=0.275*0.275; | |
2674 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2675 | nfp = 0.0; | |
2676 | } | |
2677 | else{ | |
2678 | if (dgt==K0STM||dgt==K0STB||dgt==K0STP||dgt==K0ST0) { | |
2679 | msq=0.55; | |
2680 | bx2=0.30*0.30; | |
2681 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
2682 | nfp = 2.0; | |
2683 | } | |
2684 | else{ | |
2685 | report(ERROR,"EvtGen") << "Not implemented daugt in get_EvtISGW2_ff_3P0.\n"; | |
2686 | } | |
2687 | } | |
2688 | } | |
2689 | else{ | |
2690 | if (prnt==DSP||prnt==DSM){ | |
2691 | msb=1.82; | |
2692 | msd=0.55; | |
2693 | bb2=0.56*0.56; | |
2694 | mbb=1.968; | |
2695 | nf = 3.0; | |
2696 | ||
2697 | if (dgt==F0||dgt==F0PR||dgt==A00||dgt==A0P||dgt==A0M) { | |
2698 | msq=0.55; | |
2699 | bx2=0.33*0.33; | |
2700 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
2701 | nfp = 2.0; | |
2702 | } | |
2703 | else{ | |
2704 | if (dgt==K0STM||dgt==K0STB||dgt==K0STP||dgt==K0ST0) { | |
2705 | msq=0.33; | |
2706 | bx2=0.30*0.30; | |
2707 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2708 | nfp = 0.0; | |
2709 | } | |
2710 | else{ | |
2711 | report(ERROR,"EvtGen") << "Not implemented daugt in get_EvtISGW2_ff_3P0.\n"; | |
2712 | } | |
2713 | } | |
2714 | } | |
2715 | else{ | |
2716 | //BS -> cs constants added by djl on Jan. 21,1998 | |
2717 | if (prnt==BS0||prnt==BSB){ | |
2718 | ||
2719 | msb=5.2; | |
2720 | msd=0.55; | |
2721 | bb2=0.54*0.54; | |
2722 | mbb=5.38; | |
2723 | nf = 4.0; | |
2724 | ||
2725 | if (dgt==D3P0SP||dgt==D3P0SN) { | |
2726 | ||
2727 | msq=1.82; | |
2728 | bx2=0.41*0.41; | |
2729 | mbx=(3.0*2.54+2.46)/4.0; | |
2730 | nfp = 3.0; | |
2731 | } | |
2732 | else if (dgt==K0STM||dgt==K0STB||dgt==K0STP||dgt==K0ST0) { | |
2733 | msq=0.55; | |
2734 | bx2=0.30*0.30; | |
2735 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
2736 | nfp = 2.0; | |
2737 | } | |
2738 | else{ | |
2739 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
2740 | } | |
2741 | } | |
2742 | else{ | |
2743 | report(ERROR,"EvtGen") << "Not implemented parent in get_EvtISGW2_ff_3P0.\n"; | |
2744 | } | |
2745 | } | |
2746 | } | |
2747 | } | |
2748 | ||
2749 | ||
2750 | mtb = msb + msd; | |
2751 | mtx = msq + msd; | |
2752 | ||
2753 | mb = EvtPDL::getMeanMass( parent ); | |
2754 | mx = mass; | |
2755 | ||
2756 | mup=1.0/(1.0/msq+1.0/msb); | |
2757 | mum=1.0/(1.0/msq-1.0/msb); | |
2758 | bbx2=0.5*(bb2+bx2); | |
2759 | tm=(mb-mx)*(mb-mx); | |
2760 | if (t>tm) t = 0.99*tm; | |
2761 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
2762 | ||
2763 | mqm = 0.1; | |
2764 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2) + | |
2765 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
2766 | log(EvtGetas(mqm)/EvtGetas(msq)); | |
2767 | ||
2768 | f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0) / | |
2769 | (pow((1.0+r2*(tm-t)/18.0),3.0)); | |
2770 | ||
2771 | f5uppum = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),0.5); | |
2772 | f5upmum = f5*pow(( mbb / mtb ),0.5)*pow((mbx/mtx),-0.5); | |
2773 | ||
2774 | uppum = -1.0*f5uppum*sqrt(2.0/(3.0*bb2))*msd; | |
2775 | upmum = 1.0*f5upmum*sqrt(2.0/(3.0*bb2))*msd*mtb/mtx; | |
2776 | ||
2777 | *upf = (uppum + upmum)/2.0; | |
2778 | *umf = (uppum - upmum)/2.0; | |
2779 | ||
2780 | return; | |
2781 | ||
2782 | } | |
2783 | ||
2784 | ||
2785 | void EvtISGW2FF::EvtISGW2FF3P2 (EvtId parent,EvtId daugt, | |
2786 | double t, double mass, double *hf, double *kf, | |
2787 | double *bpf, double *bmf ) { | |
2788 | ||
2789 | //added by Lange Jan4,2000 | |
2790 | static EvtId EM=EvtPDL::getId("e-"); | |
2791 | static EvtId EP=EvtPDL::getId("e+"); | |
2792 | static EvtId MUM=EvtPDL::getId("mu-"); | |
2793 | static EvtId MUP=EvtPDL::getId("mu+"); | |
2794 | static EvtId TAUM=EvtPDL::getId("tau-"); | |
2795 | static EvtId TAUP=EvtPDL::getId("tau+"); | |
2796 | ||
2797 | static EvtId BP=EvtPDL::getId("B+"); | |
2798 | static EvtId BM=EvtPDL::getId("B-"); | |
2799 | static EvtId B0=EvtPDL::getId("B0"); | |
2800 | static EvtId B0B=EvtPDL::getId("anti-B0"); | |
2801 | ||
2802 | static EvtId DST0=EvtPDL::getId("D*0"); | |
2803 | static EvtId DSTB=EvtPDL::getId("anti-D*0"); | |
2804 | static EvtId DSTP=EvtPDL::getId("D*+"); | |
2805 | static EvtId DSTM=EvtPDL::getId("D*-"); | |
2806 | static EvtId D0=EvtPDL::getId("D0"); | |
2807 | static EvtId D0B=EvtPDL::getId("anti-D0"); | |
2808 | static EvtId DP=EvtPDL::getId("D+"); | |
2809 | static EvtId DM=EvtPDL::getId("D-"); | |
2810 | ||
2811 | static EvtId D1P1P=EvtPDL::getId("D_1+"); | |
2812 | static EvtId D1P1N=EvtPDL::getId("D_1-"); | |
2813 | static EvtId D1P10=EvtPDL::getId("D_10"); | |
2814 | static EvtId D1P1B=EvtPDL::getId("anti-D_10"); | |
2815 | ||
2816 | static EvtId D3P2P=EvtPDL::getId("D_2*+"); | |
2817 | static EvtId D3P2N=EvtPDL::getId("D_2*-"); | |
2818 | static EvtId D3P20=EvtPDL::getId("D_2*0"); | |
2819 | static EvtId D3P2B=EvtPDL::getId("anti-D_2*0"); | |
2820 | ||
2821 | static EvtId D3P1P=EvtPDL::getId("D'_1+"); | |
2822 | static EvtId D3P1N=EvtPDL::getId("D'_1-"); | |
2823 | static EvtId D3P10=EvtPDL::getId("D'_10"); | |
2824 | static EvtId D3P1B=EvtPDL::getId("anti-D'_10"); | |
2825 | ||
2826 | static EvtId D3P0P=EvtPDL::getId("D_0*+"); | |
2827 | static EvtId D3P0N=EvtPDL::getId("D_0*-"); | |
2828 | static EvtId D3P00=EvtPDL::getId("D_0*0"); | |
2829 | static EvtId D3P0B=EvtPDL::getId("anti-D_0*0"); | |
2830 | ||
2831 | static EvtId D21S0P=EvtPDL::getId("D(2S)+"); | |
2832 | static EvtId D21S0N=EvtPDL::getId("D(2S)-"); | |
2833 | static EvtId D21S00=EvtPDL::getId("D(2S)0"); | |
2834 | static EvtId D21S0B=EvtPDL::getId("anti-D(2S)0"); | |
2835 | ||
2836 | static EvtId D23S1P=EvtPDL::getId("D*(2S)+"); | |
2837 | static EvtId D23S1N=EvtPDL::getId("D*(2S)-"); | |
2838 | static EvtId D23S10=EvtPDL::getId("D*(2S)0"); | |
2839 | static EvtId D23S1B=EvtPDL::getId("anti-D*(2S)0"); | |
2840 | ||
2841 | ||
2842 | static EvtId RHO2S0=EvtPDL::getId("rho(2S)0"); | |
2843 | static EvtId RHO2SP=EvtPDL::getId("rho(2S)+"); | |
2844 | static EvtId RHO2SM=EvtPDL::getId("rho(2S)-"); | |
2845 | static EvtId OMEG2S=EvtPDL::getId("omega(2S)"); | |
2846 | static EvtId ETA2S=EvtPDL::getId("eta(2S)"); | |
2847 | ||
2848 | static EvtId PI2S0=EvtPDL::getId("pi(2S)0"); | |
2849 | static EvtId PI2SP=EvtPDL::getId("pi(2S)+"); | |
2850 | static EvtId PI2SM=EvtPDL::getId("pi(2S)-"); | |
2851 | ||
2852 | static EvtId PIP=EvtPDL::getId("pi+"); | |
2853 | static EvtId PIM=EvtPDL::getId("pi-"); | |
2854 | static EvtId PI0=EvtPDL::getId("pi0"); | |
2855 | ||
2856 | static EvtId RHOP=EvtPDL::getId("rho+"); | |
2857 | static EvtId RHOM=EvtPDL::getId("rho-"); | |
2858 | static EvtId RHO0=EvtPDL::getId("rho0"); | |
2859 | ||
2860 | static EvtId A2P=EvtPDL::getId("a_2+"); | |
2861 | static EvtId A2M=EvtPDL::getId("a_2-"); | |
2862 | static EvtId A20=EvtPDL::getId("a_20"); | |
2863 | ||
2864 | static EvtId A1P=EvtPDL::getId("a_1+"); | |
2865 | static EvtId A1M=EvtPDL::getId("a_1-"); | |
2866 | static EvtId A10=EvtPDL::getId("a_10"); | |
2867 | ||
2868 | static EvtId A0P=EvtPDL::getId("a_0+"); | |
2869 | static EvtId A0M=EvtPDL::getId("a_0-"); | |
2870 | static EvtId A00=EvtPDL::getId("a_00"); | |
2871 | ||
2872 | static EvtId B1P=EvtPDL::getId("b_1+"); | |
2873 | static EvtId B1M=EvtPDL::getId("b_1-"); | |
2874 | static EvtId B10=EvtPDL::getId("b_10"); | |
2875 | ||
2876 | static EvtId H1=EvtPDL::getId("h_1"); | |
2877 | static EvtId H1PR=EvtPDL::getId("h'_1"); | |
2878 | ||
2879 | static EvtId F1=EvtPDL::getId("f_1"); | |
2880 | static EvtId F1PR=EvtPDL::getId("f'_1"); | |
2881 | static EvtId F0=EvtPDL::getId("f_0"); | |
2882 | static EvtId F0PR=EvtPDL::getId("f'_0"); | |
2883 | static EvtId F2=EvtPDL::getId("f_2"); | |
2884 | static EvtId F2PR=EvtPDL::getId("f'_2"); | |
2885 | ||
2886 | static EvtId ETA=EvtPDL::getId("eta"); | |
2887 | static EvtId ETAPR=EvtPDL::getId("eta'"); | |
2888 | static EvtId OMEG=EvtPDL::getId("omega"); | |
2889 | ||
2890 | static EvtId KP=EvtPDL::getId("K+"); | |
2891 | static EvtId KM=EvtPDL::getId("K-"); | |
2892 | static EvtId K0=EvtPDL::getId("K0"); | |
2893 | static EvtId KB=EvtPDL::getId("anti-K0"); | |
2894 | static EvtId K0S=EvtPDL::getId("K_S0"); | |
2895 | static EvtId K0L=EvtPDL::getId("K_L0"); | |
2896 | ||
2897 | static EvtId KSTP=EvtPDL::getId("K*+"); | |
2898 | static EvtId KSTM=EvtPDL::getId("K*-"); | |
2899 | static EvtId KST0=EvtPDL::getId("K*0"); | |
2900 | static EvtId KSTB=EvtPDL::getId("anti-K*0"); | |
2901 | ||
2902 | static EvtId K1P=EvtPDL::getId("K_1+"); | |
2903 | static EvtId K1M=EvtPDL::getId("K_1-"); | |
2904 | static EvtId K10=EvtPDL::getId("K_10"); | |
2905 | static EvtId K1B=EvtPDL::getId("anti-K_10"); | |
2906 | ||
2907 | static EvtId K1STP=EvtPDL::getId("K'_1+"); | |
2908 | static EvtId K1STM=EvtPDL::getId("K'_1-"); | |
2909 | static EvtId K1ST0=EvtPDL::getId("K'_10"); | |
2910 | static EvtId K1STB=EvtPDL::getId("anti-K'_10"); | |
2911 | ||
2912 | static EvtId K2STP=EvtPDL::getId("K_2*+"); | |
2913 | static EvtId K2STM=EvtPDL::getId("K_2*-"); | |
2914 | static EvtId K2ST0=EvtPDL::getId("K_2*0"); | |
2915 | static EvtId K2STB=EvtPDL::getId("anti-K_2*0"); | |
2916 | ||
2917 | static EvtId K0STP=EvtPDL::getId("K_0*+"); | |
2918 | static EvtId K0STM=EvtPDL::getId("K_0*-"); | |
2919 | static EvtId K0ST0=EvtPDL::getId("K_0*0"); | |
2920 | static EvtId K0STB=EvtPDL::getId("anti-K_0*0"); | |
2921 | ||
2922 | static EvtId PHI=EvtPDL::getId("phi"); | |
2923 | static EvtId DSP=EvtPDL::getId("D_s+"); | |
2924 | static EvtId DSM=EvtPDL::getId("D_s-"); | |
2925 | ||
2926 | static EvtId D1P1SP=EvtPDL::getId("D_s1+"); | |
2927 | static EvtId D1P1SN=EvtPDL::getId("D_s1-"); | |
2928 | ||
2929 | static EvtId D3P0SP=EvtPDL::getId("D_s0*+"); | |
2930 | static EvtId D3P0SN=EvtPDL::getId("D_s0*-"); | |
2931 | ||
2932 | static EvtId D3P1SP=EvtPDL::getId("D'_s1+"); | |
2933 | static EvtId D3P1SN=EvtPDL::getId("D'_s1-"); | |
2934 | ||
2935 | static EvtId D3P2SP=EvtPDL::getId("D_s2*+"); | |
2936 | static EvtId D3P2SN=EvtPDL::getId("D_s2*-"); | |
2937 | ||
2938 | static EvtId DSSTP=EvtPDL::getId("D_s*+"); | |
2939 | static EvtId DSSTM=EvtPDL::getId("D_s*-"); | |
2940 | ||
2941 | static EvtId BSB=EvtPDL::getId("anti-B_s0"); | |
2942 | static EvtId BS0=EvtPDL::getId("B_s0"); | |
2943 | ||
2944 | ||
2945 | double mtb, mbb(0.0); | |
2946 | double msd(0.0), mx,mb,nf,nfp(0.0); | |
2947 | double msq(0.0),bx2(0.0),mbx(0.0),mtx,f5; | |
2948 | double f5h,f5k,f5bppbm,f5bpmbm,bppbm,bpmbm; | |
2949 | double mqm,mum,mup,tm,wt,r2,bb2(0.0),bbx2; | |
2950 | double msb(0.0); | |
2951 | EvtId prnt=parent; | |
2952 | EvtId dgt=daugt; | |
2953 | ||
2954 | if (prnt==B0||prnt==B0B||prnt==BP||prnt==BM) { | |
2955 | ||
2956 | msb=5.2; | |
2957 | msd=0.33; | |
2958 | bb2=0.431*0.431; | |
2959 | mbb=5.31; | |
2960 | nf = 4.0; | |
2961 | ||
2962 | if (dgt==A20||dgt==A2P||dgt==A2M||dgt==F2||dgt==F2PR) { | |
2963 | ||
2964 | msq=0.33; | |
2965 | bx2=0.275*0.275; | |
2966 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2967 | nfp = 0.0; | |
2968 | ||
2969 | } | |
2970 | ||
2971 | else{ | |
2972 | if (dgt==D3P2P||dgt==D3P2N||dgt==D3P20||dgt==D3P2B) { | |
2973 | ||
2974 | msq=1.82; | |
2975 | bx2=0.33*0.33; | |
2976 | mbx=(5.0*2.46+3.0*2.42)/8.0; | |
2977 | nfp = 3.0; | |
2978 | } | |
2979 | else{ | |
2980 | ||
2981 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P2.\n"; | |
2982 | } | |
2983 | } | |
2984 | } | |
2985 | else{ | |
2986 | if (prnt==D0||prnt==D0B||prnt==DP||prnt==DM) { | |
2987 | ||
2988 | msb=1.82; | |
2989 | msd=0.33; | |
2990 | bb2=0.45*0.45; | |
2991 | mbb=1.963; | |
2992 | nf = 3.0; | |
2993 | if (dgt==F2||dgt==F2PR||dgt==A20||dgt==A2P||dgt==A2M) { | |
2994 | msq=0.33; | |
2995 | bx2=0.275*0.275; | |
2996 | mbx=(3.0*1.23+0.98+5.0*1.32+3.0*1.26)/12.0; | |
2997 | nfp = 0.0; | |
2998 | } | |
2999 | else{ | |
3000 | if (dgt==K2STM||dgt==K2STB||dgt==K2STP||dgt==K2ST0) { | |
3001 | msq=0.55; | |
3002 | bx2=0.30*0.30; | |
3003 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
3004 | nfp = 2.0; | |
3005 | } | |
3006 | else{ | |
3007 | report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P2.\n"; | |
3008 | } | |
3009 | } | |
3010 | } | |
3011 | else{ | |
3012 | //BS -> cs constants added by djl on Jan. 21,1998 | |
3013 | if (prnt==BS0||prnt==BSB){ | |
3014 | ||
3015 | msb=5.2; | |
3016 | msd=0.55; | |
3017 | bb2=0.54*0.54; | |
3018 | mbb=5.38; | |
3019 | nf = 4.0; | |
3020 | ||
3021 | if (dgt==D3P2SP||dgt==D3P2SN) { | |
3022 | ||
3023 | msq=1.82; | |
3024 | bx2=0.41*0.41; | |
3025 | mbx=(5.0*2.61+3.0*2.54)/8.0; | |
3026 | nfp = 3.0; | |
3027 | } | |
3028 | else if (dgt==K2STM||dgt==K2STB||dgt==K2STP||dgt==K2ST0) { | |
3029 | msq=0.55; | |
3030 | bx2=0.30*0.30; | |
3031 | mbx=(3.0*1.40+1.43+5.0*1.43+3.0*1.27)/12.0; | |
3032 | nfp = 2.0; | |
3033 | } | |
3034 | else{ | |
3035 | report(ERROR,"EvtGen") << "Not implemented daugt:"<<daugt.getId()<<" in get_isgw_ff_1S0.\n"; | |
3036 | } | |
3037 | } | |
3038 | ||
3039 | else{ | |
3040 | report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P2.\n"; | |
3041 | } | |
3042 | } | |
3043 | } | |
3044 | mtb = msb + msd; | |
3045 | mtx = msq + msd; | |
3046 | ||
3047 | mb = EvtPDL::getMeanMass( parent ); | |
3048 | mx = mass; | |
3049 | mup=1.0/(1.0/msq+1.0/msb); | |
3050 | mum=1.0/(1.0/msq-1.0/msb); | |
3051 | bbx2=0.5*(bb2+bx2); | |
3052 | tm=(mb-mx)*(mb-mx); | |
3053 | if (t>tm) t = 0.99*tm; | |
3054 | wt=1.0+(tm-t)/(2.0*mbb*mbx); | |
3055 | ||
3056 | mqm = 0.1; | |
3057 | r2=3.0/(4.0*msb*msq)+3*msd*msd/(2*mbb*mbx*bbx2)+ | |
3058 | (16.0/(mbb*mbx*(33.0-2.0*nfp)))* | |
3059 | log(EvtGetas(mqm)/EvtGetas(msq)); | |
3060 | ||
3061 | f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0) / | |
3062 | (pow((1.0+r2*(tm-t)/18.0),3.0)); | |
3063 | ||
3064 | f5h = f5*pow(( mbb / mtb ),-1.5)*pow((mbx/mtx),-0.5); | |
3065 | f5k = f5*pow(( mbb / mtb ),-0.5)*pow((mbx/mtx),0.5); | |
3066 | f5bppbm = f5*pow(( mbb / mtb ),-2.5)*pow((mbx/mtx),0.5); | |
3067 | f5bpmbm = f5*pow(( mbb / mtb ),-1.5)*pow((mbx/mtx),-0.5); | |
3068 | ||
3069 | *hf = f5h*(msd/(sqrt(8.0*bb2)*mtb))*((1.0/msq)-(msd*bb2/(2.0*mum* | |
3070 | mtx*bbx2))); | |
3071 | ||
3072 | *kf = f5k*(msd/(sqrt(2.0*bb2)))*(1.0+wt); | |
3073 | ||
3074 | bppbm = ((msd*msd*f5bppbm*bx2)/(sqrt(32.0*bb2)*msq*msb*mtb*bbx2))* | |
3075 | (1.0-(msd*bx2/(2.0*mtb*bbx2))); | |
3076 | ||
3077 | bpmbm = -1.0*(msd*f5bpmbm/(sqrt(2.0*bb2)*msb*mtx))*(1.0- | |
3078 | ((msd*msb*bx2)/(2.0*mup*mtb*bbx2))+((msd*bx2*(1.0- | |
3079 | ((msd*bx2)/(2.0*mtb*bbx2))))/(4.0*msq*bbx2))); | |
3080 | ||
3081 | *bpf = (bppbm + bpmbm)/2.0; | |
3082 | *bmf = (bppbm - bpmbm)/2.0; | |
3083 | return; | |
3084 | } //get_ff_isgw_1p1 | |
3085 | ||
3086 | ||
3087 | double EvtISGW2FF::EvtGetGammaji ( double z ) | |
3088 | ||
3089 | { | |
3090 | double temp; | |
3091 | ||
3092 | temp = 2+((2.0*z)/(1-z))*log(z); | |
3093 | temp = -1.0*temp; | |
3094 | ||
3095 | return temp; | |
3096 | ||
3097 | } //EvtGetGammaji | |
3098 | ||
3099 | ||
3100 | ||
3101 | double EvtISGW2FF::EvtGetas ( double massq, double massx ) | |
3102 | { | |
3103 | double lqcd2 = 0.04; | |
3104 | double nflav = 4; | |
3105 | double temp = 0.6; | |
3106 | ||
3107 | if ( massx > 0.6 ) { | |
3108 | if ( massq < 1.85 ) { | |
3109 | nflav = 3.0;} | |
3110 | ||
3111 | temp = 12.0*EvtConst::pi / ( 33.0 - 2.0*nflav) / | |
3112 | log( massx*massx/lqcd2); | |
3113 | } | |
3114 | return temp; | |
3115 | ||
3116 | } //EvtGetas | |
3117 | ||
3118 | double EvtISGW2FF::EvtGetas ( double mass ) | |
3119 | ||
3120 | { | |
3121 | double lqcd2 = 0.04; | |
3122 | double nflav = 4; | |
3123 | double temp = 0.6; | |
3124 | ||
3125 | if ( mass > 0.6 ) { | |
3126 | if ( mass < 1.85 ) { | |
3127 | nflav = 3.0;} | |
3128 | ||
3129 | temp = 12.0*EvtConst::pi / ( 33.0 - 2.0*nflav) / | |
3130 | log( mass*mass/lqcd2); | |
3131 | } | |
3132 | return temp; | |
3133 | ||
3134 | } //EvtGetas | |
3135 | ||
3136 | ||
3137 | void EvtISGW2FF::getbaryonff(EvtId, EvtId, double, double, double*, | |
3138 | double*, double*, double*){ | |
3139 | ||
3140 | report(ERROR,"EvtGen") << "Not implemented :getbaryonff in EvtISGW2FF.\n"; | |
3141 | ||
3142 | ::abort(); | |
3143 | ||
3144 | } | |
3145 |