BOSS 7.0.7
BESIII Offline Software System
Loading...
Searching...
No Matches
EvtISGWFF.cc
Go to the documentation of this file.
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: EvtISGWFF.cc
12//
13// Description: Routine to implement semileptonic form factors
14// according to the model ISGW
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/EvtISGWFF.hh"
26#include "EvtGenBase/EvtPDL.hh"
27#include "EvtGenBase/EvtId.hh"
28#include <string>
29#include <math.h>
30using std::endl;
31
32void EvtISGWFF::getscalarff(EvtId parent,EvtId daught,
33 double t, double mass, double *fpf,
34 double *f0f ) {
35
36 //added by Lange Jan4,2000
37 static EvtId EM=EvtPDL::getId("e-");
38 static EvtId EP=EvtPDL::getId("e+");
39 static EvtId MUM=EvtPDL::getId("mu-");
40 static EvtId MUP=EvtPDL::getId("mu+");
41 static EvtId TAUM=EvtPDL::getId("tau-");
42 static EvtId TAUP=EvtPDL::getId("tau+");
43
44 static EvtId BP=EvtPDL::getId("B+");
45 static EvtId BM=EvtPDL::getId("B-");
46 static EvtId B0=EvtPDL::getId("B0");
47 static EvtId B0B=EvtPDL::getId("anti-B0");
48
49 static EvtId DST0=EvtPDL::getId("D*0");
50 static EvtId DSTB=EvtPDL::getId("anti-D*0");
51 static EvtId DSTP=EvtPDL::getId("D*+");
52 static EvtId DSTM=EvtPDL::getId("D*-");
53 static EvtId D0=EvtPDL::getId("D0");
54 static EvtId D0B=EvtPDL::getId("anti-D0");
55 static EvtId DP=EvtPDL::getId("D+");
56 static EvtId DM=EvtPDL::getId("D-");
57
58 static EvtId D1P1P=EvtPDL::getId("D_1+");
59 static EvtId D1P1N=EvtPDL::getId("D_1-");
60 static EvtId D1P10=EvtPDL::getId("D_10");
61 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
62
63 static EvtId D3P2P=EvtPDL::getId("D_2*+");
64 static EvtId D3P2N=EvtPDL::getId("D_2*-");
65 static EvtId D3P20=EvtPDL::getId("D_2*0");
66 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
67
68 static EvtId D3P1P=EvtPDL::getId("D'_1+");
69 static EvtId D3P1N=EvtPDL::getId("D'_1-");
70 static EvtId D3P10=EvtPDL::getId("D'_10");
71 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
72
73 static EvtId D3P0P=EvtPDL::getId("D_0*+");
74 static EvtId D3P0N=EvtPDL::getId("D_0*-");
75 static EvtId D3P00=EvtPDL::getId("D_0*0");
76 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
77
78 static EvtId D21S0P=EvtPDL::getId("hi");
79 static EvtId D21S0N=EvtPDL::getId("hi");
80 static EvtId D21S00=EvtPDL::getId("hi");
81 static EvtId D21S0B=EvtPDL::getId("hi");
82
83 static EvtId D23S1P=EvtPDL::getId("hi");
84 static EvtId D23S1N=EvtPDL::getId("hi");
85 static EvtId D23S10=EvtPDL::getId("hi");
86 static EvtId D23S1B=EvtPDL::getId("hi");
87
88 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
89 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
90 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
91 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
92 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
93
94 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
95 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
96 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
97
98 static EvtId PIP=EvtPDL::getId("pi+");
99 static EvtId PIM=EvtPDL::getId("pi-");
100 static EvtId PI0=EvtPDL::getId("pi0");
101
102 static EvtId RHOP=EvtPDL::getId("rho+");
103 static EvtId RHOM=EvtPDL::getId("rho-");
104 static EvtId RHO0=EvtPDL::getId("rho0");
105
106 static EvtId A2P=EvtPDL::getId("a_2+");
107 static EvtId A2M=EvtPDL::getId("a_2-");
108 static EvtId A20=EvtPDL::getId("a_20");
109
110 static EvtId A1P=EvtPDL::getId("a_1+");
111 static EvtId A1M=EvtPDL::getId("a_1-");
112 static EvtId A10=EvtPDL::getId("a_10");
113
114 static EvtId A0P=EvtPDL::getId("a_0+");
115 static EvtId A0M=EvtPDL::getId("a_0-");
116 static EvtId A00=EvtPDL::getId("a_00");
117
118 static EvtId B1P=EvtPDL::getId("b_1+");
119 static EvtId B1M=EvtPDL::getId("b_1-");
120 static EvtId B10=EvtPDL::getId("b_10");
121
122 static EvtId H1=EvtPDL::getId("h_1");
123 static EvtId H1PR=EvtPDL::getId("h'_1");
124
125 static EvtId F1=EvtPDL::getId("f_1");
126 static EvtId F1PR=EvtPDL::getId("f'_1");
127 static EvtId F0=EvtPDL::getId("f_0");
128 static EvtId F0PR=EvtPDL::getId("f'_0");
129 static EvtId F2=EvtPDL::getId("f_2");
130 static EvtId F2PR=EvtPDL::getId("f'_2");
131
132 static EvtId ETA=EvtPDL::getId("eta");
133 static EvtId ETAPR=EvtPDL::getId("eta'");
134 static EvtId OMEG=EvtPDL::getId("omega");
135
136 static EvtId KP=EvtPDL::getId("K+");
137 static EvtId KM=EvtPDL::getId("K-");
138 static EvtId K0=EvtPDL::getId("K0");
139 static EvtId KB=EvtPDL::getId("anti-K0");
140 static EvtId K0S=EvtPDL::getId("K_S0");
141 static EvtId K0L=EvtPDL::getId("K_L0");
142
143 static EvtId KSTP=EvtPDL::getId("K*+");
144 static EvtId KSTM=EvtPDL::getId("K*-");
145 static EvtId KST0=EvtPDL::getId("K*0");
146 static EvtId KSTB=EvtPDL::getId("anti-K*0");
147
148 static EvtId K1P=EvtPDL::getId("K_1+");
149 static EvtId K1M=EvtPDL::getId("K_1-");
150 static EvtId K10=EvtPDL::getId("K_10");
151 static EvtId K1B=EvtPDL::getId("anti-K_10");
152
153 static EvtId K1STP=EvtPDL::getId("K'_1+");
154 static EvtId K1STM=EvtPDL::getId("K'_1-");
155 static EvtId K1ST0=EvtPDL::getId("K'_10");
156 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
157
158 static EvtId K2STP=EvtPDL::getId("K_2*+");
159 static EvtId K2STM=EvtPDL::getId("K_2*-");
160 static EvtId K2ST0=EvtPDL::getId("K_2*0");
161 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
162
163 static EvtId K0STP=EvtPDL::getId("K_0*+");
164 static EvtId K0STM=EvtPDL::getId("K_0*-");
165 static EvtId K0ST0=EvtPDL::getId("K_0*0");
166 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
167
168 static EvtId PHI=EvtPDL::getId("phi");
169 static EvtId DSP=EvtPDL::getId("D_s+");
170 static EvtId DSM=EvtPDL::getId("D_s-");
171
172 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
173 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
174
175 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
176 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
177
178 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
179 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
180
181 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
182 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
183
184 static EvtId DSSTP=EvtPDL::getId("D_s*+");
185 static EvtId DSSTM=EvtPDL::getId("D_s*-");
186
187 static EvtId BSB=EvtPDL::getId("anti-B_s0");
188 static EvtId BS0=EvtPDL::getId("B_s0");
189 double fmf;
190 double mb=EvtPDL::getMeanMass(parent);
191
192 if (daught==PI0||daught==PIP||daught==PIM||daught==ETA||
193 daught==ETAPR||daught==D0||daught==D0B||daught==DP||
194 daught==DM||daught==KP||daught==KM||daught==K0||daught==K0L||
195 daught==KB||daught==DSP||daught==DSM||daught==K0S) {
196
197 EvtISGW1FF1S0(parent,daught,t,mass,fpf,&fmf);
198 }
199
200 if (daught==PI2S0||daught==PI2SP||daught==PI2SM||daught==ETA2S||
201 daught==D21S0P||daught==D21S0B||daught==D21S0N||daught==D21S00){
202 EvtISGW1FF21S0(parent,daught,t,mass,fpf,&fmf);
203 }
204
205 if (daught==A00||daught==A0P||daught==A0M||daught==F0||
206 daught==F0PR||daught==D3P0P||daught==D3P00||daught==D3P0B||
207 daught==D3P0N||daught==K0STM||daught==K0STB||daught==K0STP||
208 daught==D3P0SP||daught==D3P0SN||
209 daught==K0ST0) {
210 EvtISGW1FF3P0(parent,daught,t,mass,fpf,&fmf);
211 }
212
213 *f0f = (fmf/((mb*mb-mass*mass)/t))+(*fpf);
214
215 return ;
216}
217
218 void EvtISGWFF::gettensorff(EvtId parent,EvtId daught,
219 double t, double mass, double *hf,
220 double *kf, double *bpf, double *bmf ){
221
222 //added by Lange Jan4,2000
223 static EvtId EM=EvtPDL::getId("e-");
224 static EvtId EP=EvtPDL::getId("e+");
225 static EvtId MUM=EvtPDL::getId("mu-");
226 static EvtId MUP=EvtPDL::getId("mu+");
227 static EvtId TAUM=EvtPDL::getId("tau-");
228 static EvtId TAUP=EvtPDL::getId("tau+");
229
230 static EvtId BP=EvtPDL::getId("B+");
231 static EvtId BM=EvtPDL::getId("B-");
232 static EvtId B0=EvtPDL::getId("B0");
233 static EvtId B0B=EvtPDL::getId("anti-B0");
234
235 static EvtId DST0=EvtPDL::getId("D*0");
236 static EvtId DSTB=EvtPDL::getId("anti-D*0");
237 static EvtId DSTP=EvtPDL::getId("D*+");
238 static EvtId DSTM=EvtPDL::getId("D*-");
239 static EvtId D0=EvtPDL::getId("D0");
240 static EvtId D0B=EvtPDL::getId("anti-D0");
241 static EvtId DP=EvtPDL::getId("D+");
242 static EvtId DM=EvtPDL::getId("D-");
243
244 static EvtId D1P1P=EvtPDL::getId("D_1+");
245 static EvtId D1P1N=EvtPDL::getId("D_1-");
246 static EvtId D1P10=EvtPDL::getId("D_10");
247 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
248
249 static EvtId D3P2P=EvtPDL::getId("D_2*+");
250 static EvtId D3P2N=EvtPDL::getId("D_2*-");
251 static EvtId D3P20=EvtPDL::getId("D_2*0");
252 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
253
254 static EvtId D3P1P=EvtPDL::getId("D'_1+");
255 static EvtId D3P1N=EvtPDL::getId("D'_1-");
256 static EvtId D3P10=EvtPDL::getId("D'_10");
257 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
258
259 static EvtId D3P0P=EvtPDL::getId("D_0*+");
260 static EvtId D3P0N=EvtPDL::getId("D_0*-");
261 static EvtId D3P00=EvtPDL::getId("D_0*0");
262 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
263
264 static EvtId D21S0P=EvtPDL::getId("hi");
265 static EvtId D21S0N=EvtPDL::getId("hi");
266 static EvtId D21S00=EvtPDL::getId("hi");
267 static EvtId D21S0B=EvtPDL::getId("hi");
268
269 static EvtId D23S1P=EvtPDL::getId("hi");
270 static EvtId D23S1N=EvtPDL::getId("hi");
271 static EvtId D23S10=EvtPDL::getId("hi");
272 static EvtId D23S1B=EvtPDL::getId("hi");
273
274 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
275 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
276 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
277 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
278 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
279
280 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
281 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
282 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
283
284 static EvtId PIP=EvtPDL::getId("pi+");
285 static EvtId PIM=EvtPDL::getId("pi-");
286 static EvtId PI0=EvtPDL::getId("pi0");
287
288 static EvtId RHOP=EvtPDL::getId("rho+");
289 static EvtId RHOM=EvtPDL::getId("rho-");
290 static EvtId RHO0=EvtPDL::getId("rho0");
291
292 static EvtId A2P=EvtPDL::getId("a_2+");
293 static EvtId A2M=EvtPDL::getId("a_2-");
294 static EvtId A20=EvtPDL::getId("a_20");
295
296 static EvtId A1P=EvtPDL::getId("a_1+");
297 static EvtId A1M=EvtPDL::getId("a_1-");
298 static EvtId A10=EvtPDL::getId("a_10");
299
300 static EvtId A0P=EvtPDL::getId("a_0+");
301 static EvtId A0M=EvtPDL::getId("a_0-");
302 static EvtId A00=EvtPDL::getId("a_00");
303
304 static EvtId B1P=EvtPDL::getId("b_1+");
305 static EvtId B1M=EvtPDL::getId("b_1-");
306 static EvtId B10=EvtPDL::getId("b_10");
307
308 static EvtId H1=EvtPDL::getId("h_1");
309 static EvtId H1PR=EvtPDL::getId("h'_1");
310
311 static EvtId F1=EvtPDL::getId("f_1");
312 static EvtId F1PR=EvtPDL::getId("f'_1");
313 static EvtId F0=EvtPDL::getId("f_0");
314 static EvtId F0PR=EvtPDL::getId("f'_0");
315 static EvtId F2=EvtPDL::getId("f_2");
316 static EvtId F2PR=EvtPDL::getId("f'_2");
317
318 static EvtId ETA=EvtPDL::getId("eta");
319 static EvtId ETAPR=EvtPDL::getId("eta'");
320 static EvtId OMEG=EvtPDL::getId("omega");
321
322 static EvtId KP=EvtPDL::getId("K+");
323 static EvtId KM=EvtPDL::getId("K-");
324 static EvtId K0=EvtPDL::getId("K0");
325 static EvtId KB=EvtPDL::getId("anti-K0");
326 static EvtId K0S=EvtPDL::getId("K_S0");
327 static EvtId K0L=EvtPDL::getId("K_L0");
328
329 static EvtId KSTP=EvtPDL::getId("K*+");
330 static EvtId KSTM=EvtPDL::getId("K*-");
331 static EvtId KST0=EvtPDL::getId("K*0");
332 static EvtId KSTB=EvtPDL::getId("anti-K*0");
333
334 static EvtId K1P=EvtPDL::getId("K_1+");
335 static EvtId K1M=EvtPDL::getId("K_1-");
336 static EvtId K10=EvtPDL::getId("K_10");
337 static EvtId K1B=EvtPDL::getId("anti-K_10");
338
339 static EvtId K1STP=EvtPDL::getId("K'_1+");
340 static EvtId K1STM=EvtPDL::getId("K'_1-");
341 static EvtId K1ST0=EvtPDL::getId("K'_10");
342 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
343
344 static EvtId K2STP=EvtPDL::getId("K_2*+");
345 static EvtId K2STM=EvtPDL::getId("K_2*-");
346 static EvtId K2ST0=EvtPDL::getId("K_2*0");
347 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
348
349 static EvtId K0STP=EvtPDL::getId("K_0*+");
350 static EvtId K0STM=EvtPDL::getId("K_0*-");
351 static EvtId K0ST0=EvtPDL::getId("K_0*0");
352 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
353
354 static EvtId PHI=EvtPDL::getId("phi");
355 static EvtId DSP=EvtPDL::getId("D_s+");
356 static EvtId DSM=EvtPDL::getId("D_s-");
357
358 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
359 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
360
361 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
362 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
363
364 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
365 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
366
367 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
368 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
369
370 static EvtId DSSTP=EvtPDL::getId("D_s*+");
371 static EvtId DSSTM=EvtPDL::getId("D_s*-");
372
373 static EvtId BSB=EvtPDL::getId("anti-B_s0");
374 static EvtId BS0=EvtPDL::getId("B_s0");
375
376 EvtISGW1FF3P2(parent,daught,t,mass,hf,kf,bpf,bmf);
377
378 return;
379
380 }
381
382
383 void EvtISGWFF::getvectorff(EvtId parent,EvtId daught,
384 double t, double mass, double *a1f,
385 double *a2f, double *vf, double *a0f ){
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("hi");
430 static EvtId D21S0N=EvtPDL::getId("hi");
431 static EvtId D21S00=EvtPDL::getId("hi");
432 static EvtId D21S0B=EvtPDL::getId("hi");
433
434 static EvtId D23S1P=EvtPDL::getId("hi");
435 static EvtId D23S1N=EvtPDL::getId("hi");
436 static EvtId D23S10=EvtPDL::getId("hi");
437 static EvtId D23S1B=EvtPDL::getId("hi");
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 double ff,gf,apf,amf;
542
543 if (daught==DST0||daught==DSTP||daught==DSTM||daught==DSTB||
544 daught==OMEG||daught==RHO0||daught==RHOM||daught==RHOP||
545 daught==KSTP||daught==KSTM||daught==KST0||daught==KSTB||
546 daught==PHI||daught==DSSTP||daught==DSSTM) {
547 EvtISGW1FF3S1(parent,daught,t,mass,&ff,&gf,&apf,&amf);
548 }
549 if (daught==B10||daught==B1P||daught==B1M||daught==H1||
550 daught==H1PR||daught==D1P1P||daught==D1P10||daught==D1P1B||
551 daught==D1P1SP||daught==D1P1SN||
552 daught==D1P1N||daught==K10||daught==K1B||daught==K1P||
553 daught==K1M) {
554 EvtISGW1FF1P1(parent,daught,t,mass,&ff,&gf,&apf,&amf);
555 }
556 if (daught==RHO2S0||daught==RHO2SP||daught==RHO2SM||daught==OMEG2S||
557 daught==D23S1P||daught==D23S1B||daught==D23S1N||daught==D23S10){
558 EvtISGW1FF23S1(parent,daught,t,mass,&ff,&gf,&apf,&amf);
559 }
560 if (daught==A10||daught==A1P||daught==A1M||daught==F1||
561 daught==F1PR||daught==D3P1P||daught==D3P10||daught==D3P1B||
562 daught==D3P1N||daught==K1STM||daught==K1STB||daught==K1STP||
563 daught==D3P1SP||daught==D3P1SN||
564 daught==K1ST0) {
565 EvtISGW1FF3P1(parent,daught,t,mass,&ff,&gf,&apf,&amf);
566 }
567
568 // Need to stuff in some factors to make these the ffs that
569 // is used elsewhere...
570
571 double mb=EvtPDL::getMeanMass(parent);
572
573
574 *vf = (gf)*(mb+mass);
575 *a1f = (ff)/(mb+mass);
576 *a2f = -1.0*(apf)*(mb+mass);
577 double a3f = ((mb+mass)/(2.0*mass))*(*a1f) -
578 ((mb-mass)/(2.0*mass))*(*a2f);
579
580 *a0f = a3f - ( (t*amf)/(2.0*mass));
581
582
583
584 return;
585 }
586
587
588
589void EvtISGWFF::EvtISGW1FF3P2 (EvtId parent,EvtId daugt,
590 double t, double mass,
591 double *hf, double *kf, double *bpf, double *bmf ) {
592
593 //added by Lange Jan4,2000
594 static EvtId EM=EvtPDL::getId("e-");
595 static EvtId EP=EvtPDL::getId("e+");
596 static EvtId MUM=EvtPDL::getId("mu-");
597 static EvtId MUP=EvtPDL::getId("mu+");
598 static EvtId TAUM=EvtPDL::getId("tau-");
599 static EvtId TAUP=EvtPDL::getId("tau+");
600
601 static EvtId BP=EvtPDL::getId("B+");
602 static EvtId BM=EvtPDL::getId("B-");
603 static EvtId B0=EvtPDL::getId("B0");
604 static EvtId B0B=EvtPDL::getId("anti-B0");
605
606 static EvtId DST0=EvtPDL::getId("D*0");
607 static EvtId DSTB=EvtPDL::getId("anti-D*0");
608 static EvtId DSTP=EvtPDL::getId("D*+");
609 static EvtId DSTM=EvtPDL::getId("D*-");
610 static EvtId D0=EvtPDL::getId("D0");
611 static EvtId D0B=EvtPDL::getId("anti-D0");
612 static EvtId DP=EvtPDL::getId("D+");
613 static EvtId DM=EvtPDL::getId("D-");
614
615 static EvtId D1P1P=EvtPDL::getId("D_1+");
616 static EvtId D1P1N=EvtPDL::getId("D_1-");
617 static EvtId D1P10=EvtPDL::getId("D_10");
618 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
619
620 static EvtId D3P2P=EvtPDL::getId("D_2*+");
621 static EvtId D3P2N=EvtPDL::getId("D_2*-");
622 static EvtId D3P20=EvtPDL::getId("D_2*0");
623 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
624
625 static EvtId D3P1P=EvtPDL::getId("D'_1+");
626 static EvtId D3P1N=EvtPDL::getId("D'_1-");
627 static EvtId D3P10=EvtPDL::getId("D'_10");
628 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
629
630 static EvtId D3P0P=EvtPDL::getId("D_0*+");
631 static EvtId D3P0N=EvtPDL::getId("D_0*-");
632 static EvtId D3P00=EvtPDL::getId("D_0*0");
633 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
634
635 static EvtId D21S0P=EvtPDL::getId("hi");
636 static EvtId D21S0N=EvtPDL::getId("hi");
637 static EvtId D21S00=EvtPDL::getId("hi");
638 static EvtId D21S0B=EvtPDL::getId("hi");
639
640 static EvtId D23S1P=EvtPDL::getId("hi");
641 static EvtId D23S1N=EvtPDL::getId("hi");
642 static EvtId D23S10=EvtPDL::getId("hi");
643 static EvtId D23S1B=EvtPDL::getId("hi");
644
645 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
646 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
647 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
648 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
649 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
650
651 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
652 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
653 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
654
655 static EvtId PIP=EvtPDL::getId("pi+");
656 static EvtId PIM=EvtPDL::getId("pi-");
657 static EvtId PI0=EvtPDL::getId("pi0");
658
659 static EvtId RHOP=EvtPDL::getId("rho+");
660 static EvtId RHOM=EvtPDL::getId("rho-");
661 static EvtId RHO0=EvtPDL::getId("rho0");
662
663 static EvtId A2P=EvtPDL::getId("a_2+");
664 static EvtId A2M=EvtPDL::getId("a_2-");
665 static EvtId A20=EvtPDL::getId("a_20");
666
667 static EvtId A1P=EvtPDL::getId("a_1+");
668 static EvtId A1M=EvtPDL::getId("a_1-");
669 static EvtId A10=EvtPDL::getId("a_10");
670
671 static EvtId A0P=EvtPDL::getId("a_0+");
672 static EvtId A0M=EvtPDL::getId("a_0-");
673 static EvtId A00=EvtPDL::getId("a_00");
674
675 static EvtId B1P=EvtPDL::getId("b_1+");
676 static EvtId B1M=EvtPDL::getId("b_1-");
677 static EvtId B10=EvtPDL::getId("b_10");
678
679 static EvtId H1=EvtPDL::getId("h_1");
680 static EvtId H1PR=EvtPDL::getId("h'_1");
681
682 static EvtId F1=EvtPDL::getId("f_1");
683 static EvtId F1PR=EvtPDL::getId("f'_1");
684 static EvtId F0=EvtPDL::getId("f_0");
685 static EvtId F0PR=EvtPDL::getId("f'_0");
686 static EvtId F2=EvtPDL::getId("f_2");
687 static EvtId F2PR=EvtPDL::getId("f'_2");
688
689 static EvtId ETA=EvtPDL::getId("eta");
690 static EvtId ETAPR=EvtPDL::getId("eta'");
691 static EvtId OMEG=EvtPDL::getId("omega");
692
693 static EvtId KP=EvtPDL::getId("K+");
694 static EvtId KM=EvtPDL::getId("K-");
695 static EvtId K0=EvtPDL::getId("K0");
696 static EvtId KB=EvtPDL::getId("anti-K0");
697 static EvtId K0S=EvtPDL::getId("K_S0");
698 static EvtId K0L=EvtPDL::getId("K_L0");
699
700 static EvtId KSTP=EvtPDL::getId("K*+");
701 static EvtId KSTM=EvtPDL::getId("K*-");
702 static EvtId KST0=EvtPDL::getId("K*0");
703 static EvtId KSTB=EvtPDL::getId("anti-K*0");
704
705 static EvtId K1P=EvtPDL::getId("K_1+");
706 static EvtId K1M=EvtPDL::getId("K_1-");
707 static EvtId K10=EvtPDL::getId("K_10");
708 static EvtId K1B=EvtPDL::getId("anti-K_10");
709
710 static EvtId K1STP=EvtPDL::getId("K'_1+");
711 static EvtId K1STM=EvtPDL::getId("K'_1-");
712 static EvtId K1ST0=EvtPDL::getId("K'_10");
713 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
714
715 static EvtId K2STP=EvtPDL::getId("K_2*+");
716 static EvtId K2STM=EvtPDL::getId("K_2*-");
717 static EvtId K2ST0=EvtPDL::getId("K_2*0");
718 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
719
720 static EvtId K0STP=EvtPDL::getId("K_0*+");
721 static EvtId K0STM=EvtPDL::getId("K_0*-");
722 static EvtId K0ST0=EvtPDL::getId("K_0*0");
723 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
724
725 static EvtId PHI=EvtPDL::getId("phi");
726 static EvtId DSP=EvtPDL::getId("D_s+");
727 static EvtId DSM=EvtPDL::getId("D_s-");
728
729 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
730 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
731
732 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
733 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
734
735 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
736 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
737
738 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
739 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
740
741 static EvtId DSSTP=EvtPDL::getId("D_s*+");
742 static EvtId DSSTM=EvtPDL::getId("D_s*-");
743
744 static EvtId BSB=EvtPDL::getId("anti-B_s0");
745 static EvtId BS0=EvtPDL::getId("B_s0");
746
747 double mtb;
748 double msd(0.0), mx(0.0), mb(0.0);
749 double msq(0.0), bx2(0.0),mtx,f5;
750 double mum,mup,tm,bb2(0.0),bbx2;
751 double msb(0.0), kap;
752
753 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
754 msb=5.2;
755 msd=0.33;
756 bb2=0.41*0.41;
757 if (daugt==A20||daugt==A2P||daugt==A2M||daugt==F2||daugt==F2PR) {
758 msq=0.33;
759 bx2=0.27*0.27;
760 }
761 else{
762 if (daugt==D3P2P||daugt==D3P2N||daugt==D3P2B||daugt==D3P20) {
763 msq=1.82;
764 bx2=0.34*0.34;
765 }
766 else{
767 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P1.\n";
768 }
769 }
770 }
771 else{
772 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P1.\n";
773 }
774
775 mtb = msb + msd;
776 mtx = msq + msd;
777
778 mb = EvtPDL::getMeanMass( parent );
779 mx = mass;
780
781 mup=1.0/(1.0/msq+1.0/msb);
782 mum=1.0/(1.0/msq-1.0/msb);
783 bbx2=0.5*(bb2+bx2);
784
785 tm=(mb-mx)*(mb-mx);
786 if (t>tm) t = 0.99*tm;
787 kap = 0.7*0.7;
788
789 f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0)*
790 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
791
792 *hf = f5*(msd/(sqrt(8.0*bb2)*mtb))*((1.0/msq)-(msd*bb2/(2.0*mum*
793 mtx*bbx2)));
794
795 *kf = f5*msd*sqrt(2.0/bb2);
796
797 *bpf = (-1.0*f5*msd/(sqrt(8.0*bb2)*msb*mtx))*(1.0-(msd*msb*bx2/(
798 2.0*mup*mtb*bbx2))+(msd*msb*bx2*(1.0-(msd*bx2/(2.0*mtb*bbx2)))/
799 (4.0*mtb*mum*bbx2)));
800 *bmf = 0.0;
801 return;
802} //get_ff_isgw_1p1
803
804void EvtISGWFF::EvtISGW1FF1S0 ( EvtId parent, EvtId daugt,
805 double t, double mass, double *fpf, double *fmf ) {
806
807 //added by Lange Jan4,2000
808 static EvtId EM=EvtPDL::getId("e-");
809 static EvtId EP=EvtPDL::getId("e+");
810 static EvtId MUM=EvtPDL::getId("mu-");
811 static EvtId MUP=EvtPDL::getId("mu+");
812 static EvtId TAUM=EvtPDL::getId("tau-");
813 static EvtId TAUP=EvtPDL::getId("tau+");
814
815 static EvtId BP=EvtPDL::getId("B+");
816 static EvtId BM=EvtPDL::getId("B-");
817 static EvtId B0=EvtPDL::getId("B0");
818 static EvtId B0B=EvtPDL::getId("anti-B0");
819
820 static EvtId DST0=EvtPDL::getId("D*0");
821 static EvtId DSTB=EvtPDL::getId("anti-D*0");
822 static EvtId DSTP=EvtPDL::getId("D*+");
823 static EvtId DSTM=EvtPDL::getId("D*-");
824 static EvtId D0=EvtPDL::getId("D0");
825 static EvtId D0B=EvtPDL::getId("anti-D0");
826 static EvtId DP=EvtPDL::getId("D+");
827 static EvtId DM=EvtPDL::getId("D-");
828
829 static EvtId D1P1P=EvtPDL::getId("D_1+");
830 static EvtId D1P1N=EvtPDL::getId("D_1-");
831 static EvtId D1P10=EvtPDL::getId("D_10");
832 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
833
834 static EvtId D3P2P=EvtPDL::getId("D_2*+");
835 static EvtId D3P2N=EvtPDL::getId("D_2*-");
836 static EvtId D3P20=EvtPDL::getId("D_2*0");
837 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
838
839 static EvtId D3P1P=EvtPDL::getId("D'_1+");
840 static EvtId D3P1N=EvtPDL::getId("D'_1-");
841 static EvtId D3P10=EvtPDL::getId("D'_10");
842 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
843
844 static EvtId D3P0P=EvtPDL::getId("D_0*+");
845 static EvtId D3P0N=EvtPDL::getId("D_0*-");
846 static EvtId D3P00=EvtPDL::getId("D_0*0");
847 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
848
849 static EvtId D21S0P=EvtPDL::getId("hi");
850 static EvtId D21S0N=EvtPDL::getId("hi");
851 static EvtId D21S00=EvtPDL::getId("hi");
852 static EvtId D21S0B=EvtPDL::getId("hi");
853
854 static EvtId D23S1P=EvtPDL::getId("hi");
855 static EvtId D23S1N=EvtPDL::getId("hi");
856 static EvtId D23S10=EvtPDL::getId("hi");
857 static EvtId D23S1B=EvtPDL::getId("hi");
858
859 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
860 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
861 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
862 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
863 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
864
865 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
866 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
867 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
868
869 static EvtId PIP=EvtPDL::getId("pi+");
870 static EvtId PIM=EvtPDL::getId("pi-");
871 static EvtId PI0=EvtPDL::getId("pi0");
872
873 static EvtId RHOP=EvtPDL::getId("rho+");
874 static EvtId RHOM=EvtPDL::getId("rho-");
875 static EvtId RHO0=EvtPDL::getId("rho0");
876
877 static EvtId A2P=EvtPDL::getId("a_2+");
878 static EvtId A2M=EvtPDL::getId("a_2-");
879 static EvtId A20=EvtPDL::getId("a_20");
880
881 static EvtId A1P=EvtPDL::getId("a_1+");
882 static EvtId A1M=EvtPDL::getId("a_1-");
883 static EvtId A10=EvtPDL::getId("a_10");
884
885 static EvtId A0P=EvtPDL::getId("a_0+");
886 static EvtId A0M=EvtPDL::getId("a_0-");
887 static EvtId A00=EvtPDL::getId("a_00");
888
889 static EvtId B1P=EvtPDL::getId("b_1+");
890 static EvtId B1M=EvtPDL::getId("b_1-");
891 static EvtId B10=EvtPDL::getId("b_10");
892
893 static EvtId H1=EvtPDL::getId("h_1");
894 static EvtId H1PR=EvtPDL::getId("h'_1");
895
896 static EvtId F1=EvtPDL::getId("f_1");
897 static EvtId F1PR=EvtPDL::getId("f'_1");
898 static EvtId F0=EvtPDL::getId("f_0");
899 static EvtId F0PR=EvtPDL::getId("f'_0");
900 static EvtId F2=EvtPDL::getId("f_2");
901 static EvtId F2PR=EvtPDL::getId("f'_2");
902
903 static EvtId ETA=EvtPDL::getId("eta");
904 static EvtId ETAPR=EvtPDL::getId("eta'");
905 static EvtId OMEG=EvtPDL::getId("omega");
906
907 static EvtId KP=EvtPDL::getId("K+");
908 static EvtId KM=EvtPDL::getId("K-");
909 static EvtId K0=EvtPDL::getId("K0");
910 static EvtId KB=EvtPDL::getId("anti-K0");
911 static EvtId K0S=EvtPDL::getId("K_S0");
912 static EvtId K0L=EvtPDL::getId("K_L0");
913
914 static EvtId KSTP=EvtPDL::getId("K*+");
915 static EvtId KSTM=EvtPDL::getId("K*-");
916 static EvtId KST0=EvtPDL::getId("K*0");
917 static EvtId KSTB=EvtPDL::getId("anti-K*0");
918
919 static EvtId K1P=EvtPDL::getId("K_1+");
920 static EvtId K1M=EvtPDL::getId("K_1-");
921 static EvtId K10=EvtPDL::getId("K_10");
922 static EvtId K1B=EvtPDL::getId("anti-K_10");
923
924 static EvtId K1STP=EvtPDL::getId("K'_1+");
925 static EvtId K1STM=EvtPDL::getId("K'_1-");
926 static EvtId K1ST0=EvtPDL::getId("K'_10");
927 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
928
929 static EvtId K2STP=EvtPDL::getId("K_2*+");
930 static EvtId K2STM=EvtPDL::getId("K_2*-");
931 static EvtId K2ST0=EvtPDL::getId("K_2*0");
932 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
933
934 static EvtId K0STP=EvtPDL::getId("K_0*+");
935 static EvtId K0STM=EvtPDL::getId("K_0*-");
936 static EvtId K0ST0=EvtPDL::getId("K_0*0");
937 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
938
939 static EvtId PHI=EvtPDL::getId("phi");
940 static EvtId DSP=EvtPDL::getId("D_s+");
941 static EvtId DSM=EvtPDL::getId("D_s-");
942
943 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
944 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
945
946 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
947 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
948
949 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
950 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
951
952 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
953 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
954
955 static EvtId DSSTP=EvtPDL::getId("D_s*+");
956 static EvtId DSSTM=EvtPDL::getId("D_s*-");
957
958 static EvtId BSB=EvtPDL::getId("anti-B_s0");
959 static EvtId BS0=EvtPDL::getId("B_s0");
960
961 double mtb;
962 double msd(0.0), mx(0.0), mb(0.0);
963 double msq(0.0), bx2(0.0),mtx;
964 double f3,kap;
965 double msb(0.0),bb2(0.0),mup,mum,bbx2,tm;
966
967 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
968 msb=5.2;
969 msd=0.33;
970 bb2=0.41*0.41;
971 if (daugt==PIP||daugt==PIM||daugt==PI0||daugt==ETA||daugt==ETAPR) {
972 msq=0.33;
973 bx2=0.31*0.31;
974 }
975 else{
976 if (daugt==D0||daugt==DP||daugt==DM||daugt==D0B) {
977 msq=1.82;
978 bx2=0.39*0.39;
979 }
980 else{
981 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_1S0.\n";
982 }
983 }
984 }
985 else{
986 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_1S0.\n";
987 report(ERROR,"EvtGen") << "Parent:"<<parent.getId()<<endl;
988 }
989
990 mtb = msb + msd;
991 mtx = msq + msd;
992 mb = EvtPDL::getMeanMass( parent );
993 mx = mass;
994 mup=1.0/(1.0/msq+1.0/msb);
995 mum=1.0/(1.0/msq-1.0/msb);
996 bbx2=0.5*(bb2+bx2);
997 tm=(mb-mx)*(mb-mx);
998 if ( t>tm ) t=0.99*tm;
999
1000 kap = 0.7*0.7;
1001 f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0)*
1002 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1003
1004 *fpf = f3*(1+(msb/(2.0*mum))-(msb*msq*msd*bb2/(4.0*mup*mum*mtx*bbx2)));
1005 *fmf = f3*(1.0-(mtb+mtx)*(0.5/msq-(msd*bb2/(4.0*mup*mtx*bbx2))));
1006
1007 return;
1008} //get_ff_isgw_1s0
1009
1010
1011
1012void EvtISGWFF::EvtISGW1FF3S1(EvtId parent,EvtId daugt,double t,
1013 double mass, double *f,double *g,double *ap,double *am){
1014
1015 //added by Lange Jan4,2000
1016 static EvtId EM=EvtPDL::getId("e-");
1017 static EvtId EP=EvtPDL::getId("e+");
1018 static EvtId MUM=EvtPDL::getId("mu-");
1019 static EvtId MUP=EvtPDL::getId("mu+");
1020 static EvtId TAUM=EvtPDL::getId("tau-");
1021 static EvtId TAUP=EvtPDL::getId("tau+");
1022
1023 static EvtId BP=EvtPDL::getId("B+");
1024 static EvtId BM=EvtPDL::getId("B-");
1025 static EvtId B0=EvtPDL::getId("B0");
1026 static EvtId B0B=EvtPDL::getId("anti-B0");
1027
1028 static EvtId DST0=EvtPDL::getId("D*0");
1029 static EvtId DSTB=EvtPDL::getId("anti-D*0");
1030 static EvtId DSTP=EvtPDL::getId("D*+");
1031 static EvtId DSTM=EvtPDL::getId("D*-");
1032 static EvtId D0=EvtPDL::getId("D0");
1033 static EvtId D0B=EvtPDL::getId("anti-D0");
1034 static EvtId DP=EvtPDL::getId("D+");
1035 static EvtId DM=EvtPDL::getId("D-");
1036
1037 static EvtId D1P1P=EvtPDL::getId("D_1+");
1038 static EvtId D1P1N=EvtPDL::getId("D_1-");
1039 static EvtId D1P10=EvtPDL::getId("D_10");
1040 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
1041
1042 static EvtId D3P2P=EvtPDL::getId("D_2*+");
1043 static EvtId D3P2N=EvtPDL::getId("D_2*-");
1044 static EvtId D3P20=EvtPDL::getId("D_2*0");
1045 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
1046
1047 static EvtId D3P1P=EvtPDL::getId("D'_1+");
1048 static EvtId D3P1N=EvtPDL::getId("D'_1-");
1049 static EvtId D3P10=EvtPDL::getId("D'_10");
1050 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
1051
1052 static EvtId D3P0P=EvtPDL::getId("D_0*+");
1053 static EvtId D3P0N=EvtPDL::getId("D_0*-");
1054 static EvtId D3P00=EvtPDL::getId("D_0*0");
1055 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
1056
1057 static EvtId D21S0P=EvtPDL::getId("hi");
1058 static EvtId D21S0N=EvtPDL::getId("hi");
1059 static EvtId D21S00=EvtPDL::getId("hi");
1060 static EvtId D21S0B=EvtPDL::getId("hi");
1061
1062 static EvtId D23S1P=EvtPDL::getId("hi");
1063 static EvtId D23S1N=EvtPDL::getId("hi");
1064 static EvtId D23S10=EvtPDL::getId("hi");
1065 static EvtId D23S1B=EvtPDL::getId("hi");
1066
1067 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
1068 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
1069 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
1070 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
1071 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
1072
1073 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
1074 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
1075 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
1076
1077 static EvtId PIP=EvtPDL::getId("pi+");
1078 static EvtId PIM=EvtPDL::getId("pi-");
1079 static EvtId PI0=EvtPDL::getId("pi0");
1080
1081 static EvtId RHOP=EvtPDL::getId("rho+");
1082 static EvtId RHOM=EvtPDL::getId("rho-");
1083 static EvtId RHO0=EvtPDL::getId("rho0");
1084
1085 static EvtId A2P=EvtPDL::getId("a_2+");
1086 static EvtId A2M=EvtPDL::getId("a_2-");
1087 static EvtId A20=EvtPDL::getId("a_20");
1088
1089 static EvtId A1P=EvtPDL::getId("a_1+");
1090 static EvtId A1M=EvtPDL::getId("a_1-");
1091 static EvtId A10=EvtPDL::getId("a_10");
1092
1093 static EvtId A0P=EvtPDL::getId("a_0+");
1094 static EvtId A0M=EvtPDL::getId("a_0-");
1095 static EvtId A00=EvtPDL::getId("a_00");
1096
1097 static EvtId B1P=EvtPDL::getId("b_1+");
1098 static EvtId B1M=EvtPDL::getId("b_1-");
1099 static EvtId B10=EvtPDL::getId("b_10");
1100
1101 static EvtId H1=EvtPDL::getId("h_1");
1102 static EvtId H1PR=EvtPDL::getId("h'_1");
1103
1104 static EvtId F1=EvtPDL::getId("f_1");
1105 static EvtId F1PR=EvtPDL::getId("f'_1");
1106 static EvtId F0=EvtPDL::getId("f_0");
1107 static EvtId F0PR=EvtPDL::getId("f'_0");
1108 static EvtId F2=EvtPDL::getId("f_2");
1109 static EvtId F2PR=EvtPDL::getId("f'_2");
1110
1111 static EvtId ETA=EvtPDL::getId("eta");
1112 static EvtId ETAPR=EvtPDL::getId("eta'");
1113 static EvtId OMEG=EvtPDL::getId("omega");
1114
1115 static EvtId KP=EvtPDL::getId("K+");
1116 static EvtId KM=EvtPDL::getId("K-");
1117 static EvtId K0=EvtPDL::getId("K0");
1118 static EvtId KB=EvtPDL::getId("anti-K0");
1119 static EvtId K0S=EvtPDL::getId("K_S0");
1120 static EvtId K0L=EvtPDL::getId("K_L0");
1121
1122 static EvtId KSTP=EvtPDL::getId("K*+");
1123 static EvtId KSTM=EvtPDL::getId("K*-");
1124 static EvtId KST0=EvtPDL::getId("K*0");
1125 static EvtId KSTB=EvtPDL::getId("anti-K*0");
1126
1127 static EvtId K1P=EvtPDL::getId("K_1+");
1128 static EvtId K1M=EvtPDL::getId("K_1-");
1129 static EvtId K10=EvtPDL::getId("K_10");
1130 static EvtId K1B=EvtPDL::getId("anti-K_10");
1131
1132 static EvtId K1STP=EvtPDL::getId("K'_1+");
1133 static EvtId K1STM=EvtPDL::getId("K'_1-");
1134 static EvtId K1ST0=EvtPDL::getId("K'_10");
1135 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
1136
1137 static EvtId K2STP=EvtPDL::getId("K_2*+");
1138 static EvtId K2STM=EvtPDL::getId("K_2*-");
1139 static EvtId K2ST0=EvtPDL::getId("K_2*0");
1140 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
1141
1142 static EvtId K0STP=EvtPDL::getId("K_0*+");
1143 static EvtId K0STM=EvtPDL::getId("K_0*-");
1144 static EvtId K0ST0=EvtPDL::getId("K_0*0");
1145 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
1146
1147 static EvtId PHI=EvtPDL::getId("phi");
1148 static EvtId DSP=EvtPDL::getId("D_s+");
1149 static EvtId DSM=EvtPDL::getId("D_s-");
1150
1151 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
1152 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
1153
1154 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
1155 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
1156
1157 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
1158 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
1159
1160 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
1161 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
1162
1163 static EvtId DSSTP=EvtPDL::getId("D_s*+");
1164 static EvtId DSSTM=EvtPDL::getId("D_s*-");
1165
1166 static EvtId BSB=EvtPDL::getId("anti-B_s0");
1167 static EvtId BS0=EvtPDL::getId("B_s0");
1168
1169 double msd(0.0),mup,msq(0.0),bb2(0.0),mum,mtx,bbx2;
1170 double bx2(0.0),msb(0.0),tm;
1171 double mb,mx,f3, kap;
1172
1173 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
1174 msb=5.2;
1175 msd=0.33;
1176 bb2=0.41*0.41;
1177 if (daugt==DSTP||daugt==DSTM||daugt==DSTB||daugt==DST0) {
1178 msq=1.82;
1179 bx2=0.39*0.39;
1180 }
1181 else{
1182 if (daugt==RHOP||daugt==RHOM||daugt==RHO0||daugt==OMEG) {
1183 msq=0.33;
1184 bx2=0.31*0.31;
1185 }
1186 else{
1187 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3S1.\n";
1188 }
1189 }
1190 }
1191 else{
1192 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3S1.\n";
1193 }
1194
1195 double mtb;
1196
1197 mtb=msb+msd;
1198 mtx=msq+msd;
1199
1200 mup=1.0/(1.0/msq+1.0/msb);
1201 mum=1.0/(1.0/msq-1.0/msb);
1202 bbx2=0.5*(bb2+bx2);
1203 mb=EvtPDL::getMeanMass(parent);
1204 mx=mass;
1205 tm=(mb-mx)*(mb-mx);
1206 if ( t > tm ) t = 0.99*tm;
1207
1208 kap = 0.7*0.7;
1209 f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0)*
1210 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1211
1212 *f=2.0*mtb*f3;
1213 *g=0.5*f3*((1/msq)-(msd*bb2/(2.0*mum*mtx*bbx2)));
1214 *ap=(-1.0*f3/(2.0*mtx))*(1.0+(msd*(bb2-bx2)/(msb
1215 *(bb2+bx2)))-(msd*msd*bx2*bx2/(4.0*mum*mtb*bbx2*bbx2)));
1216 *am=0.0;
1217
1218}
1219
1220void EvtISGWFF::EvtISGW1FF23S1 (EvtId parent,EvtId daugt,
1221 double t, double mass, double *fpf, double *gpf,
1222 double *appf, double *apmf ) {
1223 //added by Lange Jan4,2000
1224 static EvtId EM=EvtPDL::getId("e-");
1225 static EvtId EP=EvtPDL::getId("e+");
1226 static EvtId MUM=EvtPDL::getId("mu-");
1227 static EvtId MUP=EvtPDL::getId("mu+");
1228 static EvtId TAUM=EvtPDL::getId("tau-");
1229 static EvtId TAUP=EvtPDL::getId("tau+");
1230
1231 static EvtId BP=EvtPDL::getId("B+");
1232 static EvtId BM=EvtPDL::getId("B-");
1233 static EvtId B0=EvtPDL::getId("B0");
1234 static EvtId B0B=EvtPDL::getId("anti-B0");
1235
1236 static EvtId DST0=EvtPDL::getId("D*0");
1237 static EvtId DSTB=EvtPDL::getId("anti-D*0");
1238 static EvtId DSTP=EvtPDL::getId("D*+");
1239 static EvtId DSTM=EvtPDL::getId("D*-");
1240 static EvtId D0=EvtPDL::getId("D0");
1241 static EvtId D0B=EvtPDL::getId("anti-D0");
1242 static EvtId DP=EvtPDL::getId("D+");
1243 static EvtId DM=EvtPDL::getId("D-");
1244
1245 static EvtId D1P1P=EvtPDL::getId("D_1+");
1246 static EvtId D1P1N=EvtPDL::getId("D_1-");
1247 static EvtId D1P10=EvtPDL::getId("D_10");
1248 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
1249
1250 static EvtId D3P2P=EvtPDL::getId("D_2*+");
1251 static EvtId D3P2N=EvtPDL::getId("D_2*-");
1252 static EvtId D3P20=EvtPDL::getId("D_2*0");
1253 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
1254
1255 static EvtId D3P1P=EvtPDL::getId("D'_1+");
1256 static EvtId D3P1N=EvtPDL::getId("D'_1-");
1257 static EvtId D3P10=EvtPDL::getId("D'_10");
1258 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
1259
1260 static EvtId D3P0P=EvtPDL::getId("D_0*+");
1261 static EvtId D3P0N=EvtPDL::getId("D_0*-");
1262 static EvtId D3P00=EvtPDL::getId("D_0*0");
1263 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
1264
1265 static EvtId D21S0P=EvtPDL::getId("hi");
1266 static EvtId D21S0N=EvtPDL::getId("hi");
1267 static EvtId D21S00=EvtPDL::getId("hi");
1268 static EvtId D21S0B=EvtPDL::getId("hi");
1269
1270 static EvtId D23S1P=EvtPDL::getId("hi");
1271 static EvtId D23S1N=EvtPDL::getId("hi");
1272 static EvtId D23S10=EvtPDL::getId("hi");
1273 static EvtId D23S1B=EvtPDL::getId("hi");
1274
1275 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
1276 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
1277 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
1278 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
1279 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
1280
1281 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
1282 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
1283 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
1284
1285 static EvtId PIP=EvtPDL::getId("pi+");
1286 static EvtId PIM=EvtPDL::getId("pi-");
1287 static EvtId PI0=EvtPDL::getId("pi0");
1288
1289 static EvtId RHOP=EvtPDL::getId("rho+");
1290 static EvtId RHOM=EvtPDL::getId("rho-");
1291 static EvtId RHO0=EvtPDL::getId("rho0");
1292
1293 static EvtId A2P=EvtPDL::getId("a_2+");
1294 static EvtId A2M=EvtPDL::getId("a_2-");
1295 static EvtId A20=EvtPDL::getId("a_20");
1296
1297 static EvtId A1P=EvtPDL::getId("a_1+");
1298 static EvtId A1M=EvtPDL::getId("a_1-");
1299 static EvtId A10=EvtPDL::getId("a_10");
1300
1301 static EvtId A0P=EvtPDL::getId("a_0+");
1302 static EvtId A0M=EvtPDL::getId("a_0-");
1303 static EvtId A00=EvtPDL::getId("a_00");
1304
1305 static EvtId B1P=EvtPDL::getId("b_1+");
1306 static EvtId B1M=EvtPDL::getId("b_1-");
1307 static EvtId B10=EvtPDL::getId("b_10");
1308
1309 static EvtId H1=EvtPDL::getId("h_1");
1310 static EvtId H1PR=EvtPDL::getId("h'_1");
1311
1312 static EvtId F1=EvtPDL::getId("f_1");
1313 static EvtId F1PR=EvtPDL::getId("f'_1");
1314 static EvtId F0=EvtPDL::getId("f_0");
1315 static EvtId F0PR=EvtPDL::getId("f'_0");
1316 static EvtId F2=EvtPDL::getId("f_2");
1317 static EvtId F2PR=EvtPDL::getId("f'_2");
1318
1319 static EvtId ETA=EvtPDL::getId("eta");
1320 static EvtId ETAPR=EvtPDL::getId("eta'");
1321 static EvtId OMEG=EvtPDL::getId("omega");
1322
1323 static EvtId KP=EvtPDL::getId("K+");
1324 static EvtId KM=EvtPDL::getId("K-");
1325 static EvtId K0=EvtPDL::getId("K0");
1326 static EvtId KB=EvtPDL::getId("anti-K0");
1327 static EvtId K0S=EvtPDL::getId("K_S0");
1328 static EvtId K0L=EvtPDL::getId("K_L0");
1329
1330 static EvtId KSTP=EvtPDL::getId("K*+");
1331 static EvtId KSTM=EvtPDL::getId("K*-");
1332 static EvtId KST0=EvtPDL::getId("K*0");
1333 static EvtId KSTB=EvtPDL::getId("anti-K*0");
1334
1335 static EvtId K1P=EvtPDL::getId("K_1+");
1336 static EvtId K1M=EvtPDL::getId("K_1-");
1337 static EvtId K10=EvtPDL::getId("K_10");
1338 static EvtId K1B=EvtPDL::getId("anti-K_10");
1339
1340 static EvtId K1STP=EvtPDL::getId("K'_1+");
1341 static EvtId K1STM=EvtPDL::getId("K'_1-");
1342 static EvtId K1ST0=EvtPDL::getId("K'_10");
1343 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
1344
1345 static EvtId K2STP=EvtPDL::getId("K_2*+");
1346 static EvtId K2STM=EvtPDL::getId("K_2*-");
1347 static EvtId K2ST0=EvtPDL::getId("K_2*0");
1348 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
1349
1350 static EvtId K0STP=EvtPDL::getId("K_0*+");
1351 static EvtId K0STM=EvtPDL::getId("K_0*-");
1352 static EvtId K0ST0=EvtPDL::getId("K_0*0");
1353 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
1354
1355 static EvtId PHI=EvtPDL::getId("phi");
1356 static EvtId DSP=EvtPDL::getId("D_s+");
1357 static EvtId DSM=EvtPDL::getId("D_s-");
1358
1359 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
1360 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
1361
1362 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
1363 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
1364
1365 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
1366 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
1367
1368 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
1369 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
1370
1371 static EvtId DSSTP=EvtPDL::getId("D_s*+");
1372 static EvtId DSSTM=EvtPDL::getId("D_s*-");
1373
1374 static EvtId BSB=EvtPDL::getId("anti-B_s0");
1375 static EvtId BS0=EvtPDL::getId("B_s0");
1376
1377 double mtb;
1378 double msd(0.0), mx(0.0), mb(0.0);
1379 double msq(0.0), bx2(0.0),mtx;
1380 double f3,f5,tt;
1381 double mum,mup,bb2(0.0),bbx2,tm,msb(0.0);
1382
1383 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
1384 msb=5.2;
1385 msd=0.33;
1386 bb2=0.41*0.41;
1387 if (daugt==RHO2SP||daugt==RHO2SM||daugt==RHO2S0||daugt==OMEG2S) {
1388 msq=0.33;
1389 bx2=0.31*0.31;
1390 }
1391 else{
1392 if (daugt==D23S1N||daugt==D23S10||daugt==D23S1P||daugt==D23S1B) {
1393 msq=1.82;
1394 bx2=0.39*0.39;
1395 }
1396 else{
1397 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_23P1.\n";
1398 }
1399 }
1400 }
1401 else{
1402 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_23P1.\n";
1403 }
1404
1405 mtb = msb + msd;
1406 mtx = msq + msd;
1407 mb = EvtPDL::getMeanMass( parent );
1408 mx = mass;
1409 mup=1.0/(1.0/msq+1.0/msb);
1410 mum=1.0/(1.0/msq-1.0/msb);
1411 bbx2=0.5*(bb2+bx2);
1412 tm=(mb-mx)*(mb-mx);
1413 if (t>tm) t = 0.99*tm;
1414
1415 double kap = 0.7*0.7;
1416 f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0)*
1417 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1418
1419 f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0)*
1420 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1421
1422 *fpf = sqrt(6.0)*f3*mtb*( ((bb2-bx2)/(bb2+bx2)) + ((msd*msd*bx2*
1423 (tm-t))/(6.0*mtx*mtb*bbx2*kap*bbx2)));
1424
1425 *gpf = sqrt(3.0/8.0)*f3*( ((((bb2-bx2)/(bb2+bx2)) + ((msd*msd*bx2*
1426 (tm-t))/(6.0*mtx*mtb*bbx2*kap*bbx2)))*
1427 ((1.0/msq)-((msd*bb2)/(2.0*mum*mtx*bbx2)))) +
1428 ((msd*bb2*bx2)/(3.0*mum*mtx*bbx2*bbx2)));
1429
1430 tt = (msd*msd*bx2*(tm-t))/(mtx*mtb*bb2*kap*bbx2);
1431
1432 *appf = (f5/(sqrt(6.0)*mtx))* ( ((3.0*mtb*bbx2/(2.0*msb*sqrt(bb2*bx2)))*
1433 (1.0 - ( (msd*msd*msb*bx2*bx2)/(4.0*mtb*mtb*mum*bbx2*bbx2)))) -
1434 ( (3.0*msd*sqrt(bx2/bb2))/(2.0*msb)) +
1435 ( (5.0*msd*sqrt(bx2*bb2)*(1.0 + 0.1*tt))/(2.0*msb*bbx2)) -
1436 ( (3.0*mtb*sqrt(bb2/bx2)*(1.0 + (tt/6.0)))/(2.0*msb)) +
1437 ( (7.0*msd*msd*sqrt(bb2/bx2)*bx2*bx2*(1.0 + (tt/14.0))) /
1438 (8.0*mtb*mum*bbx2*bbx2)));
1439
1440 *apmf = 0.0;
1441 return;
1442} //get_ff_isgw_23s1
1443
1444
1445void EvtISGWFF::EvtISGW1FF3P1 (EvtId parent,EvtId daugt,
1446 double t, double mass,
1447 double *lf, double *qf, double *cpf, double *cmf ) {
1448 //added by Lange Jan4,2000
1449 static EvtId EM=EvtPDL::getId("e-");
1450 static EvtId EP=EvtPDL::getId("e+");
1451 static EvtId MUM=EvtPDL::getId("mu-");
1452 static EvtId MUP=EvtPDL::getId("mu+");
1453 static EvtId TAUM=EvtPDL::getId("tau-");
1454 static EvtId TAUP=EvtPDL::getId("tau+");
1455
1456 static EvtId BP=EvtPDL::getId("B+");
1457 static EvtId BM=EvtPDL::getId("B-");
1458 static EvtId B0=EvtPDL::getId("B0");
1459 static EvtId B0B=EvtPDL::getId("anti-B0");
1460
1461 static EvtId DST0=EvtPDL::getId("D*0");
1462 static EvtId DSTB=EvtPDL::getId("anti-D*0");
1463 static EvtId DSTP=EvtPDL::getId("D*+");
1464 static EvtId DSTM=EvtPDL::getId("D*-");
1465 static EvtId D0=EvtPDL::getId("D0");
1466 static EvtId D0B=EvtPDL::getId("anti-D0");
1467 static EvtId DP=EvtPDL::getId("D+");
1468 static EvtId DM=EvtPDL::getId("D-");
1469
1470 static EvtId D1P1P=EvtPDL::getId("D_1+");
1471 static EvtId D1P1N=EvtPDL::getId("D_1-");
1472 static EvtId D1P10=EvtPDL::getId("D_10");
1473 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
1474
1475 static EvtId D3P2P=EvtPDL::getId("D_2*+");
1476 static EvtId D3P2N=EvtPDL::getId("D_2*-");
1477 static EvtId D3P20=EvtPDL::getId("D_2*0");
1478 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
1479
1480 static EvtId D3P1P=EvtPDL::getId("D'_1+");
1481 static EvtId D3P1N=EvtPDL::getId("D'_1-");
1482 static EvtId D3P10=EvtPDL::getId("D'_10");
1483 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
1484
1485 static EvtId D3P0P=EvtPDL::getId("D_0*+");
1486 static EvtId D3P0N=EvtPDL::getId("D_0*-");
1487 static EvtId D3P00=EvtPDL::getId("D_0*0");
1488 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
1489
1490 static EvtId D21S0P=EvtPDL::getId("hi");
1491 static EvtId D21S0N=EvtPDL::getId("hi");
1492 static EvtId D21S00=EvtPDL::getId("hi");
1493 static EvtId D21S0B=EvtPDL::getId("hi");
1494
1495 static EvtId D23S1P=EvtPDL::getId("hi");
1496 static EvtId D23S1N=EvtPDL::getId("hi");
1497 static EvtId D23S10=EvtPDL::getId("hi");
1498 static EvtId D23S1B=EvtPDL::getId("hi");
1499
1500 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
1501 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
1502 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
1503 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
1504 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
1505
1506 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
1507 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
1508 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
1509
1510 static EvtId PIP=EvtPDL::getId("pi+");
1511 static EvtId PIM=EvtPDL::getId("pi-");
1512 static EvtId PI0=EvtPDL::getId("pi0");
1513
1514 static EvtId RHOP=EvtPDL::getId("rho+");
1515 static EvtId RHOM=EvtPDL::getId("rho-");
1516 static EvtId RHO0=EvtPDL::getId("rho0");
1517
1518 static EvtId A2P=EvtPDL::getId("a_2+");
1519 static EvtId A2M=EvtPDL::getId("a_2-");
1520 static EvtId A20=EvtPDL::getId("a_20");
1521
1522 static EvtId A1P=EvtPDL::getId("a_1+");
1523 static EvtId A1M=EvtPDL::getId("a_1-");
1524 static EvtId A10=EvtPDL::getId("a_10");
1525
1526 static EvtId A0P=EvtPDL::getId("a_0+");
1527 static EvtId A0M=EvtPDL::getId("a_0-");
1528 static EvtId A00=EvtPDL::getId("a_00");
1529
1530 static EvtId B1P=EvtPDL::getId("b_1+");
1531 static EvtId B1M=EvtPDL::getId("b_1-");
1532 static EvtId B10=EvtPDL::getId("b_10");
1533
1534 static EvtId H1=EvtPDL::getId("h_1");
1535 static EvtId H1PR=EvtPDL::getId("h'_1");
1536
1537 static EvtId F1=EvtPDL::getId("f_1");
1538 static EvtId F1PR=EvtPDL::getId("f'_1");
1539 static EvtId F0=EvtPDL::getId("f_0");
1540 static EvtId F0PR=EvtPDL::getId("f'_0");
1541 static EvtId F2=EvtPDL::getId("f_2");
1542 static EvtId F2PR=EvtPDL::getId("f'_2");
1543
1544 static EvtId ETA=EvtPDL::getId("eta");
1545 static EvtId ETAPR=EvtPDL::getId("eta'");
1546 static EvtId OMEG=EvtPDL::getId("omega");
1547
1548 static EvtId KP=EvtPDL::getId("K+");
1549 static EvtId KM=EvtPDL::getId("K-");
1550 static EvtId K0=EvtPDL::getId("K0");
1551 static EvtId KB=EvtPDL::getId("anti-K0");
1552 static EvtId K0S=EvtPDL::getId("K_S0");
1553 static EvtId K0L=EvtPDL::getId("K_L0");
1554
1555 static EvtId KSTP=EvtPDL::getId("K*+");
1556 static EvtId KSTM=EvtPDL::getId("K*-");
1557 static EvtId KST0=EvtPDL::getId("K*0");
1558 static EvtId KSTB=EvtPDL::getId("anti-K*0");
1559
1560 static EvtId K1P=EvtPDL::getId("K_1+");
1561 static EvtId K1M=EvtPDL::getId("K_1-");
1562 static EvtId K10=EvtPDL::getId("K_10");
1563 static EvtId K1B=EvtPDL::getId("anti-K_10");
1564
1565 static EvtId K1STP=EvtPDL::getId("K'_1+");
1566 static EvtId K1STM=EvtPDL::getId("K'_1-");
1567 static EvtId K1ST0=EvtPDL::getId("K'_10");
1568 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
1569
1570 static EvtId K2STP=EvtPDL::getId("K_2*+");
1571 static EvtId K2STM=EvtPDL::getId("K_2*-");
1572 static EvtId K2ST0=EvtPDL::getId("K_2*0");
1573 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
1574
1575 static EvtId K0STP=EvtPDL::getId("K_0*+");
1576 static EvtId K0STM=EvtPDL::getId("K_0*-");
1577 static EvtId K0ST0=EvtPDL::getId("K_0*0");
1578 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
1579
1580 static EvtId PHI=EvtPDL::getId("phi");
1581 static EvtId DSP=EvtPDL::getId("D_s+");
1582 static EvtId DSM=EvtPDL::getId("D_s-");
1583
1584 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
1585 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
1586
1587 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
1588 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
1589
1590 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
1591 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
1592
1593 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
1594 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
1595
1596 static EvtId DSSTP=EvtPDL::getId("D_s*+");
1597 static EvtId DSSTM=EvtPDL::getId("D_s*-");
1598
1599 static EvtId BSB=EvtPDL::getId("anti-B_s0");
1600 static EvtId BS0=EvtPDL::getId("B_s0");
1601
1602 double mtb;
1603 double msd(0.0), mx(0.0), mb(0.0);
1604 double msq(0.0), bx2(0.0),mtx,f5;
1605 double msb(0.0),bb2(0.0),mup,mum,bbx2,tm;
1606 double kap;
1607
1608 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
1609 msb=5.2;
1610 msd=0.33;
1611 bb2=0.41*0.41;
1612 if (daugt==A10||daugt==A1P||daugt==A1M||daugt==F1||daugt==F1PR) {
1613 msq=0.33;
1614 bx2=0.27*0.27;
1615 }
1616 else{
1617 if (daugt==D3P1P||daugt==D3P1N||daugt==D3P1B||daugt==D3P10) {
1618 msq=1.82;
1619 bx2=0.34*0.34;
1620 }
1621 else{
1622 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P1.\n";
1623 }
1624 }
1625 }
1626 else{
1627 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P1.\n";
1628 }
1629
1630 mtb = msb + msd;
1631 mtx = msq + msd;
1632
1633 mb = EvtPDL::getMeanMass( parent );
1634 mx = mass;
1635
1636 mup=1.0/(1.0/msq+1.0/msb);
1637 mum=1.0/(1.0/msq-1.0/msb);
1638 bbx2=0.5*(bb2+bx2);
1639 tm=(mb-mx)*(mb-mx);
1640 if (t>tm) t = 0.99*tm;
1641
1642 kap = 0.7*0.7;
1643 f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0)*
1644 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1645
1646 *qf = (f5*msd)/(2.0*mtx*sqrt(bb2));
1647
1648 *lf = -1.0*mtb*sqrt(bb2)*f5*(1/mum+(msd*(tm-t)/(2.0*mtb*
1649 kap*bb2))*((1.0/msq)-(1.0*msd*bb2/(2.0*mum*mtx*bbx2))));
1650
1651 *cpf = (f5*msd*msb/(4.0*mtb*sqrt(bb2)*mum))*(1.0-(msd*msq*bb2/(
1652 2.0*mtx*mum*bbx2)));
1653 *cmf = 0.0;
1654 return;
1655} //get_ff_isgw_3p1
1656
1657
1658
1659void EvtISGWFF::EvtISGW1FF3P0 (EvtId parent,EvtId daugt,
1660 double t, double mass, double *upf, double *umf ) {
1661 //added by Lange Jan4,2000
1662 static EvtId EM=EvtPDL::getId("e-");
1663 static EvtId EP=EvtPDL::getId("e+");
1664 static EvtId MUM=EvtPDL::getId("mu-");
1665 static EvtId MUP=EvtPDL::getId("mu+");
1666 static EvtId TAUM=EvtPDL::getId("tau-");
1667 static EvtId TAUP=EvtPDL::getId("tau+");
1668
1669 static EvtId BP=EvtPDL::getId("B+");
1670 static EvtId BM=EvtPDL::getId("B-");
1671 static EvtId B0=EvtPDL::getId("B0");
1672 static EvtId B0B=EvtPDL::getId("anti-B0");
1673
1674 static EvtId DST0=EvtPDL::getId("D*0");
1675 static EvtId DSTB=EvtPDL::getId("anti-D*0");
1676 static EvtId DSTP=EvtPDL::getId("D*+");
1677 static EvtId DSTM=EvtPDL::getId("D*-");
1678 static EvtId D0=EvtPDL::getId("D0");
1679 static EvtId D0B=EvtPDL::getId("anti-D0");
1680 static EvtId DP=EvtPDL::getId("D+");
1681 static EvtId DM=EvtPDL::getId("D-");
1682
1683 static EvtId D1P1P=EvtPDL::getId("D_1+");
1684 static EvtId D1P1N=EvtPDL::getId("D_1-");
1685 static EvtId D1P10=EvtPDL::getId("D_10");
1686 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
1687
1688 static EvtId D3P2P=EvtPDL::getId("D_2*+");
1689 static EvtId D3P2N=EvtPDL::getId("D_2*-");
1690 static EvtId D3P20=EvtPDL::getId("D_2*0");
1691 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
1692
1693 static EvtId D3P1P=EvtPDL::getId("D'_1+");
1694 static EvtId D3P1N=EvtPDL::getId("D'_1-");
1695 static EvtId D3P10=EvtPDL::getId("D'_10");
1696 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
1697
1698 static EvtId D3P0P=EvtPDL::getId("D_0*+");
1699 static EvtId D3P0N=EvtPDL::getId("D_0*-");
1700 static EvtId D3P00=EvtPDL::getId("D_0*0");
1701 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
1702
1703 static EvtId D21S0P=EvtPDL::getId("hi");
1704 static EvtId D21S0N=EvtPDL::getId("hi");
1705 static EvtId D21S00=EvtPDL::getId("hi");
1706 static EvtId D21S0B=EvtPDL::getId("hi");
1707
1708 static EvtId D23S1P=EvtPDL::getId("hi");
1709 static EvtId D23S1N=EvtPDL::getId("hi");
1710 static EvtId D23S10=EvtPDL::getId("hi");
1711 static EvtId D23S1B=EvtPDL::getId("hi");
1712
1713 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
1714 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
1715 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
1716 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
1717 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
1718
1719 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
1720 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
1721 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
1722
1723 static EvtId PIP=EvtPDL::getId("pi+");
1724 static EvtId PIM=EvtPDL::getId("pi-");
1725 static EvtId PI0=EvtPDL::getId("pi0");
1726
1727 static EvtId RHOP=EvtPDL::getId("rho+");
1728 static EvtId RHOM=EvtPDL::getId("rho-");
1729 static EvtId RHO0=EvtPDL::getId("rho0");
1730
1731 static EvtId A2P=EvtPDL::getId("a_2+");
1732 static EvtId A2M=EvtPDL::getId("a_2-");
1733 static EvtId A20=EvtPDL::getId("a_20");
1734
1735 static EvtId A1P=EvtPDL::getId("a_1+");
1736 static EvtId A1M=EvtPDL::getId("a_1-");
1737 static EvtId A10=EvtPDL::getId("a_10");
1738
1739 static EvtId A0P=EvtPDL::getId("a_0+");
1740 static EvtId A0M=EvtPDL::getId("a_0-");
1741 static EvtId A00=EvtPDL::getId("a_00");
1742
1743 static EvtId B1P=EvtPDL::getId("b_1+");
1744 static EvtId B1M=EvtPDL::getId("b_1-");
1745 static EvtId B10=EvtPDL::getId("b_10");
1746
1747 static EvtId H1=EvtPDL::getId("h_1");
1748 static EvtId H1PR=EvtPDL::getId("h'_1");
1749
1750 static EvtId F1=EvtPDL::getId("f_1");
1751 static EvtId F1PR=EvtPDL::getId("f'_1");
1752 static EvtId F0=EvtPDL::getId("f_0");
1753 static EvtId F0PR=EvtPDL::getId("f'_0");
1754 static EvtId F2=EvtPDL::getId("f_2");
1755 static EvtId F2PR=EvtPDL::getId("f'_2");
1756
1757 static EvtId ETA=EvtPDL::getId("eta");
1758 static EvtId ETAPR=EvtPDL::getId("eta'");
1759 static EvtId OMEG=EvtPDL::getId("omega");
1760
1761 static EvtId KP=EvtPDL::getId("K+");
1762 static EvtId KM=EvtPDL::getId("K-");
1763 static EvtId K0=EvtPDL::getId("K0");
1764 static EvtId KB=EvtPDL::getId("anti-K0");
1765 static EvtId K0S=EvtPDL::getId("K_S0");
1766 static EvtId K0L=EvtPDL::getId("K_L0");
1767
1768 static EvtId KSTP=EvtPDL::getId("K*+");
1769 static EvtId KSTM=EvtPDL::getId("K*-");
1770 static EvtId KST0=EvtPDL::getId("K*0");
1771 static EvtId KSTB=EvtPDL::getId("anti-K*0");
1772
1773 static EvtId K1P=EvtPDL::getId("K_1+");
1774 static EvtId K1M=EvtPDL::getId("K_1-");
1775 static EvtId K10=EvtPDL::getId("K_10");
1776 static EvtId K1B=EvtPDL::getId("anti-K_10");
1777
1778 static EvtId K1STP=EvtPDL::getId("K'_1+");
1779 static EvtId K1STM=EvtPDL::getId("K'_1-");
1780 static EvtId K1ST0=EvtPDL::getId("K'_10");
1781 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
1782
1783 static EvtId K2STP=EvtPDL::getId("K_2*+");
1784 static EvtId K2STM=EvtPDL::getId("K_2*-");
1785 static EvtId K2ST0=EvtPDL::getId("K_2*0");
1786 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
1787
1788 static EvtId K0STP=EvtPDL::getId("K_0*+");
1789 static EvtId K0STM=EvtPDL::getId("K_0*-");
1790 static EvtId K0ST0=EvtPDL::getId("K_0*0");
1791 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
1792
1793 static EvtId PHI=EvtPDL::getId("phi");
1794 static EvtId DSP=EvtPDL::getId("D_s+");
1795 static EvtId DSM=EvtPDL::getId("D_s-");
1796
1797 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
1798 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
1799
1800 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
1801 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
1802
1803 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
1804 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
1805
1806 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
1807 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
1808
1809 static EvtId DSSTP=EvtPDL::getId("D_s*+");
1810 static EvtId DSSTM=EvtPDL::getId("D_s*-");
1811
1812 static EvtId BSB=EvtPDL::getId("anti-B_s0");
1813 static EvtId BS0=EvtPDL::getId("B_s0");
1814
1815 double mtb;
1816 double msd(0.0), mx(0.0), mb(0.0);
1817 double msq(0.0), bx2(0.0),mtx;
1818 double f5;
1819 double mum,mup,bb2(0.0),bbx2,msb(0.0),tm;
1820
1821 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
1822 msb=5.2;
1823 msd=0.33;
1824 bb2=0.41*0.41;
1825 if (daugt==A00||daugt==A0P||daugt==A0M||daugt==F0||daugt==F0PR) {
1826 msq=0.33;
1827 bx2=0.27*0.27;
1828 }
1829 else{
1830 if (daugt==D3P0P||daugt==D3P0N||daugt==D3P0B||daugt==D3P00) {
1831 msq=1.82;
1832 bx2=0.34*0.34;
1833 }
1834 else{
1835 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P0.\n";
1836 }
1837 }
1838 }
1839 else{
1840 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P0.\n";
1841 }
1842
1843 mtb = msb + msd;
1844 mtx = msq + msd;
1845
1846 mb = EvtPDL::getMeanMass( parent );
1847 mx = mass;
1848
1849 mup=1.0/(1.0/msq+1.0/msb);
1850 mum=1.0/(1.0/msq-1.0/msb);
1851 bbx2=0.5*(bb2+bx2);
1852 tm=(mb-mx)*(mb-mx);
1853 if (t>tm) t = 0.99*tm;
1854
1855 double kap = 0.7*0.7;
1856 f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0)*
1857 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
1858
1859 *upf = f5*msd*msq*msb/(sqrt(6.0*bb2)*mtx*mum);
1860 *umf = 0.0;
1861 return;
1862} //get_ff_isgw_3p0
1863
1864
1865
1866void EvtISGWFF::EvtISGW1FF1P1 (EvtId parent,EvtId daugt,
1867 double t, double mass,
1868 double *vf, double *rf, double *spf, double *smf ) {
1869 //added by Lange Jan4,2000
1870 static EvtId EM=EvtPDL::getId("e-");
1871 static EvtId EP=EvtPDL::getId("e+");
1872 static EvtId MUM=EvtPDL::getId("mu-");
1873 static EvtId MUP=EvtPDL::getId("mu+");
1874 static EvtId TAUM=EvtPDL::getId("tau-");
1875 static EvtId TAUP=EvtPDL::getId("tau+");
1876
1877 static EvtId BP=EvtPDL::getId("B+");
1878 static EvtId BM=EvtPDL::getId("B-");
1879 static EvtId B0=EvtPDL::getId("B0");
1880 static EvtId B0B=EvtPDL::getId("anti-B0");
1881
1882 static EvtId DST0=EvtPDL::getId("D*0");
1883 static EvtId DSTB=EvtPDL::getId("anti-D*0");
1884 static EvtId DSTP=EvtPDL::getId("D*+");
1885 static EvtId DSTM=EvtPDL::getId("D*-");
1886 static EvtId D0=EvtPDL::getId("D0");
1887 static EvtId D0B=EvtPDL::getId("anti-D0");
1888 static EvtId DP=EvtPDL::getId("D+");
1889 static EvtId DM=EvtPDL::getId("D-");
1890
1891 static EvtId D1P1P=EvtPDL::getId("D_1+");
1892 static EvtId D1P1N=EvtPDL::getId("D_1-");
1893 static EvtId D1P10=EvtPDL::getId("D_10");
1894 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
1895
1896 static EvtId D3P2P=EvtPDL::getId("D_2*+");
1897 static EvtId D3P2N=EvtPDL::getId("D_2*-");
1898 static EvtId D3P20=EvtPDL::getId("D_2*0");
1899 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
1900
1901 static EvtId D3P1P=EvtPDL::getId("D'_1+");
1902 static EvtId D3P1N=EvtPDL::getId("D'_1-");
1903 static EvtId D3P10=EvtPDL::getId("D'_10");
1904 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
1905
1906 static EvtId D3P0P=EvtPDL::getId("D_0*+");
1907 static EvtId D3P0N=EvtPDL::getId("D_0*-");
1908 static EvtId D3P00=EvtPDL::getId("D_0*0");
1909 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
1910
1911 static EvtId D21S0P=EvtPDL::getId("hi");
1912 static EvtId D21S0N=EvtPDL::getId("hi");
1913 static EvtId D21S00=EvtPDL::getId("hi");
1914 static EvtId D21S0B=EvtPDL::getId("hi");
1915
1916 static EvtId D23S1P=EvtPDL::getId("hi");
1917 static EvtId D23S1N=EvtPDL::getId("hi");
1918 static EvtId D23S10=EvtPDL::getId("hi");
1919 static EvtId D23S1B=EvtPDL::getId("hi");
1920
1921 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
1922 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
1923 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
1924 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
1925 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
1926
1927 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
1928 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
1929 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
1930
1931 static EvtId PIP=EvtPDL::getId("pi+");
1932 static EvtId PIM=EvtPDL::getId("pi-");
1933 static EvtId PI0=EvtPDL::getId("pi0");
1934
1935 static EvtId RHOP=EvtPDL::getId("rho+");
1936 static EvtId RHOM=EvtPDL::getId("rho-");
1937 static EvtId RHO0=EvtPDL::getId("rho0");
1938
1939 static EvtId A2P=EvtPDL::getId("a_2+");
1940 static EvtId A2M=EvtPDL::getId("a_2-");
1941 static EvtId A20=EvtPDL::getId("a_20");
1942
1943 static EvtId A1P=EvtPDL::getId("a_1+");
1944 static EvtId A1M=EvtPDL::getId("a_1-");
1945 static EvtId A10=EvtPDL::getId("a_10");
1946
1947 static EvtId A0P=EvtPDL::getId("a_0+");
1948 static EvtId A0M=EvtPDL::getId("a_0-");
1949 static EvtId A00=EvtPDL::getId("a_00");
1950
1951 static EvtId B1P=EvtPDL::getId("b_1+");
1952 static EvtId B1M=EvtPDL::getId("b_1-");
1953 static EvtId B10=EvtPDL::getId("b_10");
1954
1955 static EvtId H1=EvtPDL::getId("h_1");
1956 static EvtId H1PR=EvtPDL::getId("h'_1");
1957
1958 static EvtId F1=EvtPDL::getId("f_1");
1959 static EvtId F1PR=EvtPDL::getId("f'_1");
1960 static EvtId F0=EvtPDL::getId("f_0");
1961 static EvtId F0PR=EvtPDL::getId("f'_0");
1962 static EvtId F2=EvtPDL::getId("f_2");
1963 static EvtId F2PR=EvtPDL::getId("f'_2");
1964
1965 static EvtId ETA=EvtPDL::getId("eta");
1966 static EvtId ETAPR=EvtPDL::getId("eta'");
1967 static EvtId OMEG=EvtPDL::getId("omega");
1968
1969 static EvtId KP=EvtPDL::getId("K+");
1970 static EvtId KM=EvtPDL::getId("K-");
1971 static EvtId K0=EvtPDL::getId("K0");
1972 static EvtId KB=EvtPDL::getId("anti-K0");
1973 static EvtId K0S=EvtPDL::getId("K_S0");
1974 static EvtId K0L=EvtPDL::getId("K_L0");
1975
1976 static EvtId KSTP=EvtPDL::getId("K*+");
1977 static EvtId KSTM=EvtPDL::getId("K*-");
1978 static EvtId KST0=EvtPDL::getId("K*0");
1979 static EvtId KSTB=EvtPDL::getId("anti-K*0");
1980
1981 static EvtId K1P=EvtPDL::getId("K_1+");
1982 static EvtId K1M=EvtPDL::getId("K_1-");
1983 static EvtId K10=EvtPDL::getId("K_10");
1984 static EvtId K1B=EvtPDL::getId("anti-K_10");
1985
1986 static EvtId K1STP=EvtPDL::getId("K'_1+");
1987 static EvtId K1STM=EvtPDL::getId("K'_1-");
1988 static EvtId K1ST0=EvtPDL::getId("K'_10");
1989 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
1990
1991 static EvtId K2STP=EvtPDL::getId("K_2*+");
1992 static EvtId K2STM=EvtPDL::getId("K_2*-");
1993 static EvtId K2ST0=EvtPDL::getId("K_2*0");
1994 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
1995
1996 static EvtId K0STP=EvtPDL::getId("K_0*+");
1997 static EvtId K0STM=EvtPDL::getId("K_0*-");
1998 static EvtId K0ST0=EvtPDL::getId("K_0*0");
1999 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
2000
2001 static EvtId PHI=EvtPDL::getId("phi");
2002 static EvtId DSP=EvtPDL::getId("D_s+");
2003 static EvtId DSM=EvtPDL::getId("D_s-");
2004
2005 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
2006 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
2007
2008 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
2009 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
2010
2011 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
2012 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
2013
2014 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
2015 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
2016
2017 static EvtId DSSTP=EvtPDL::getId("D_s*+");
2018 static EvtId DSSTM=EvtPDL::getId("D_s*-");
2019
2020 static EvtId BSB=EvtPDL::getId("anti-B_s0");
2021 static EvtId BS0=EvtPDL::getId("B_s0");
2022
2023 double mtb;
2024 double msd(0.0), mx(0.0), mb(0.0);
2025 double msq(0.0), bx2(0.0),mtx,f5;
2026 double mup,mum,kap;
2027 double msb(0.0),bb2(0.0),bbx2,tm;
2028
2029 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
2030
2031 msb=5.2;
2032 msd=0.33;
2033 bb2=0.41*0.41;
2034 if (daugt==H1||daugt==H1PR||daugt==B10||daugt==B1P||daugt==B1M) {
2035 msq=0.33;
2036 bx2=0.27*0.27;
2037 }
2038 else{
2039 if (daugt==D1P1P||daugt==D1P1N||daugt==D1P10||daugt==D1P1B) {
2040 msq=1.82;
2041 bx2=0.34*0.34;
2042 }
2043 else{
2044 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw_ff_3P1.\n";
2045 }
2046 }
2047 }
2048 else{
2049 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw_ff_3P1.\n";
2050 }
2051
2052 mtb = msb + msd;
2053 mtx = msq + msd;
2054
2055 mb = EvtPDL::getMeanMass( parent );
2056 mx = mass;
2057
2058 mup=1.0/(1.0/msq+1.0/msb);
2059 mum=1.0/(1.0/msq-1.0/msb);
2060 bbx2=0.5*(bb2+bx2);
2061 tm=(mb-mx)*(mb-mx);
2062 if (t>tm) t = 0.99*tm;
2063
2064 kap = 0.7*0.7;
2065 f5 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,5.0/2.0)*
2066 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
2067
2068 *vf = f5*(((mtb*sqrt(bb2))/(4.0*sqrt(2.0)*msb*msq*mtx)));
2069 *rf = f5*mtb*sqrt(bb2/2)*((1.0/mup));
2070
2071 *spf = (f5*msd/(sqrt(2.0*bb2)*mtb))*(1.0+(msb/(2.0*mum))-
2072 (msb*msq*msd*bb2/(4.0*mup*mum*mtx*bbx2)));
2073 *smf = 0.0;
2074
2075 return;
2076//get_ff_isgw_1p1
2077
2078}
2079
2080void EvtISGWFF::EvtISGW1FF21S0 (EvtId parent,EvtId daugt,
2081 double t, double mass, double *fppf, double *fpmf ) {
2082 //added by Lange Jan4,2000
2083 static EvtId EM=EvtPDL::getId("e-");
2084 static EvtId EP=EvtPDL::getId("e+");
2085 static EvtId MUM=EvtPDL::getId("mu-");
2086 static EvtId MUP=EvtPDL::getId("mu+");
2087 static EvtId TAUM=EvtPDL::getId("tau-");
2088 static EvtId TAUP=EvtPDL::getId("tau+");
2089
2090 static EvtId BP=EvtPDL::getId("B+");
2091 static EvtId BM=EvtPDL::getId("B-");
2092 static EvtId B0=EvtPDL::getId("B0");
2093 static EvtId B0B=EvtPDL::getId("anti-B0");
2094
2095 static EvtId DST0=EvtPDL::getId("D*0");
2096 static EvtId DSTB=EvtPDL::getId("anti-D*0");
2097 static EvtId DSTP=EvtPDL::getId("D*+");
2098 static EvtId DSTM=EvtPDL::getId("D*-");
2099 static EvtId D0=EvtPDL::getId("D0");
2100 static EvtId D0B=EvtPDL::getId("anti-D0");
2101 static EvtId DP=EvtPDL::getId("D+");
2102 static EvtId DM=EvtPDL::getId("D-");
2103
2104 static EvtId D1P1P=EvtPDL::getId("D_1+");
2105 static EvtId D1P1N=EvtPDL::getId("D_1-");
2106 static EvtId D1P10=EvtPDL::getId("D_10");
2107 static EvtId D1P1B=EvtPDL::getId("anti-D_10");
2108
2109 static EvtId D3P2P=EvtPDL::getId("D_2*+");
2110 static EvtId D3P2N=EvtPDL::getId("D_2*-");
2111 static EvtId D3P20=EvtPDL::getId("D_2*0");
2112 static EvtId D3P2B=EvtPDL::getId("anti-D_2*0");
2113
2114 static EvtId D3P1P=EvtPDL::getId("D'_1+");
2115 static EvtId D3P1N=EvtPDL::getId("D'_1-");
2116 static EvtId D3P10=EvtPDL::getId("D'_10");
2117 static EvtId D3P1B=EvtPDL::getId("anti-D'_10");
2118
2119 static EvtId D3P0P=EvtPDL::getId("D_0*+");
2120 static EvtId D3P0N=EvtPDL::getId("D_0*-");
2121 static EvtId D3P00=EvtPDL::getId("D_0*0");
2122 static EvtId D3P0B=EvtPDL::getId("anti-D_0*0");
2123
2124 static EvtId D21S0P=EvtPDL::getId("hi");
2125 static EvtId D21S0N=EvtPDL::getId("hi");
2126 static EvtId D21S00=EvtPDL::getId("hi");
2127 static EvtId D21S0B=EvtPDL::getId("hi");
2128
2129 static EvtId D23S1P=EvtPDL::getId("hi");
2130 static EvtId D23S1N=EvtPDL::getId("hi");
2131 static EvtId D23S10=EvtPDL::getId("hi");
2132 static EvtId D23S1B=EvtPDL::getId("hi");
2133
2134 static EvtId RHO2S0=EvtPDL::getId("rho(2S)0");
2135 static EvtId RHO2SP=EvtPDL::getId("rho(2S)+");
2136 static EvtId RHO2SM=EvtPDL::getId("rho(2S)-");
2137 static EvtId OMEG2S=EvtPDL::getId("omega(2S)");
2138 static EvtId ETA2S=EvtPDL::getId("eta(2S)");
2139
2140 static EvtId PI2S0=EvtPDL::getId("pi(2S)0");
2141 static EvtId PI2SP=EvtPDL::getId("pi(2S)+");
2142 static EvtId PI2SM=EvtPDL::getId("pi(2S)-");
2143
2144 static EvtId PIP=EvtPDL::getId("pi+");
2145 static EvtId PIM=EvtPDL::getId("pi-");
2146 static EvtId PI0=EvtPDL::getId("pi0");
2147
2148 static EvtId RHOP=EvtPDL::getId("rho+");
2149 static EvtId RHOM=EvtPDL::getId("rho-");
2150 static EvtId RHO0=EvtPDL::getId("rho0");
2151
2152 static EvtId A2P=EvtPDL::getId("a_2+");
2153 static EvtId A2M=EvtPDL::getId("a_2-");
2154 static EvtId A20=EvtPDL::getId("a_20");
2155
2156 static EvtId A1P=EvtPDL::getId("a_1+");
2157 static EvtId A1M=EvtPDL::getId("a_1-");
2158 static EvtId A10=EvtPDL::getId("a_10");
2159
2160 static EvtId A0P=EvtPDL::getId("a_0+");
2161 static EvtId A0M=EvtPDL::getId("a_0-");
2162 static EvtId A00=EvtPDL::getId("a_00");
2163
2164 static EvtId B1P=EvtPDL::getId("b_1+");
2165 static EvtId B1M=EvtPDL::getId("b_1-");
2166 static EvtId B10=EvtPDL::getId("b_10");
2167
2168 static EvtId H1=EvtPDL::getId("h_1");
2169 static EvtId H1PR=EvtPDL::getId("h'_1");
2170
2171 static EvtId F1=EvtPDL::getId("f_1");
2172 static EvtId F1PR=EvtPDL::getId("f'_1");
2173 static EvtId F0=EvtPDL::getId("f_0");
2174 static EvtId F0PR=EvtPDL::getId("f'_0");
2175 static EvtId F2=EvtPDL::getId("f_2");
2176 static EvtId F2PR=EvtPDL::getId("f'_2");
2177
2178 static EvtId ETA=EvtPDL::getId("eta");
2179 static EvtId ETAPR=EvtPDL::getId("eta'");
2180 static EvtId OMEG=EvtPDL::getId("omega");
2181
2182 static EvtId KP=EvtPDL::getId("K+");
2183 static EvtId KM=EvtPDL::getId("K-");
2184 static EvtId K0=EvtPDL::getId("K0");
2185 static EvtId KB=EvtPDL::getId("anti-K0");
2186 static EvtId K0S=EvtPDL::getId("K_S0");
2187 static EvtId K0L=EvtPDL::getId("K_L0");
2188
2189 static EvtId KSTP=EvtPDL::getId("K*+");
2190 static EvtId KSTM=EvtPDL::getId("K*-");
2191 static EvtId KST0=EvtPDL::getId("K*0");
2192 static EvtId KSTB=EvtPDL::getId("anti-K*0");
2193
2194 static EvtId K1P=EvtPDL::getId("K_1+");
2195 static EvtId K1M=EvtPDL::getId("K_1-");
2196 static EvtId K10=EvtPDL::getId("K_10");
2197 static EvtId K1B=EvtPDL::getId("anti-K_10");
2198
2199 static EvtId K1STP=EvtPDL::getId("K'_1+");
2200 static EvtId K1STM=EvtPDL::getId("K'_1-");
2201 static EvtId K1ST0=EvtPDL::getId("K'_10");
2202 static EvtId K1STB=EvtPDL::getId("anti-K'_10");
2203
2204 static EvtId K2STP=EvtPDL::getId("K_2*+");
2205 static EvtId K2STM=EvtPDL::getId("K_2*-");
2206 static EvtId K2ST0=EvtPDL::getId("K_2*0");
2207 static EvtId K2STB=EvtPDL::getId("anti-K_2*0");
2208
2209 static EvtId K0STP=EvtPDL::getId("K_0*+");
2210 static EvtId K0STM=EvtPDL::getId("K_0*-");
2211 static EvtId K0ST0=EvtPDL::getId("K_0*0");
2212 static EvtId K0STB=EvtPDL::getId("anti-K_0*0");
2213
2214 static EvtId PHI=EvtPDL::getId("phi");
2215 static EvtId DSP=EvtPDL::getId("D_s+");
2216 static EvtId DSM=EvtPDL::getId("D_s-");
2217
2218 static EvtId D1P1SP=EvtPDL::getId("D_s1+");
2219 static EvtId D1P1SN=EvtPDL::getId("D_s1-");
2220
2221 static EvtId D3P0SP=EvtPDL::getId("D_s0*+");
2222 static EvtId D3P0SN=EvtPDL::getId("D_s0*-");
2223
2224 static EvtId D3P1SP=EvtPDL::getId("D'_s1*+");
2225 static EvtId D3P1SN=EvtPDL::getId("D'_s1*-");
2226
2227 static EvtId D3P2SP=EvtPDL::getId("D_s2*+");
2228 static EvtId D3P2SN=EvtPDL::getId("D_s2*-");
2229
2230 static EvtId DSSTP=EvtPDL::getId("D_s*+");
2231 static EvtId DSSTM=EvtPDL::getId("D_s*-");
2232
2233 static EvtId BSB=EvtPDL::getId("anti-B_s0");
2234 static EvtId BS0=EvtPDL::getId("B_s0");
2235
2236 double mtb;
2237 double msd(0.0), mx(0.0), mb(0.0);
2238 double msq(0.0), bx2(0.0),mtx;
2239 double f3;
2240 double msb(0.0);
2241 double mum,mup,tm,bb2(0.0),bbx2;
2242
2243 if (parent==BM||parent==BP||parent==B0||parent==B0B) {
2244 msb=5.2;
2245 msd=0.33;
2246 bb2=0.41*0.41;
2247 if (daugt==PI2S0||daugt==PI2SP||daugt==PI2SM||daugt==ETA2S){
2248 msq=0.33;
2249 bx2=0.31*0.31;
2250 }
2251 else{
2252 if (daugt==D21S00||daugt==D21S0P||daugt==D21S0N||daugt==D21S0B) {
2253 msq=1.82;
2254 bx2=0.39*0.39;
2255 }
2256 else{
2257 report(ERROR,"EvtGen") << "Not implemented daugt in get_isgw1_ff_21S0.\n";
2258 }
2259 }
2260 }
2261 else{
2262 report(ERROR,"EvtGen") << "Not implemented parent in get_isgw1_ff_21S0.\n";
2263 }
2264
2265 mtb = msb + msd;
2266 mtx = msq + msd;
2267
2268 mb = EvtPDL::getMeanMass( parent );
2269 mx = mass;
2270
2271 mup=1.0/(1.0/msq+1.0/msb);
2272 mum=1.0/(1.0/msq-1.0/msb);
2273 bbx2=0.5*(bb2+bx2);
2274 tm=(mb-mx)*(mb-mx);
2275 if (t>tm) t = 0.99*tm;
2276
2277 double kap = 0.7*0.7;
2278 f3 = sqrt(mtx/mtb)*pow(sqrt(bx2*bb2)/bbx2,3.0/2.0)*
2279 exp(-1.0*((msd*msd*(tm-t)/(4.0*mtb*mtx*kap*bbx2))));
2280
2281 *fppf = f3*sqrt(3.0/8.0)*(msb/mup)*( ((bb2-bx2)/(bb2+bx2)) +
2282 (((msq*msd*bb2)/(3.0*mum*mtx*bbx2))*((7.0*bx2-3.0*bb2)/
2283 (4.0*bbx2))) +
2284 (((msd*msd*bx2*(tm-t))/(6.0*mtx*mtb*bbx2*kap*bbx2))*
2285 (1.0 - ((msq*msd*bb2)/(2.0*mum*mtx*bbx2)))));
2286
2287 *fpmf = 0.0;
2288 return;
2289} //get_ff_isgw_21s0
2290
2291
mass
double mass
Definition: CosmicGenerator.cxx:138
exp
EvtComplex exp(const EvtComplex &c)
Definition: EvtComplex.hh:252
report
ostream & report(Severity severity, const char *facility)
Definition: EvtReport.cc:36
ERROR
@ ERROR
Definition: EvtReport.hh:49
t
TTree * t
Definition: binning.cxx:23
EvtId
Definition: EvtId.hh:27
EvtId::getId
int getId() const
Definition: EvtId.hh:41
EvtPDL::getMeanMass
static double getMeanMass(EvtId i)
Definition: EvtPDL.hh:45
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cc:287
f
TFile f("ana_bhabha660a_dqa_mcPat_zy_old.root")