evtgen is hosted by Hepforge, IPPP Durham
EvtGen  2.0.0
Monte Carlo generator of particle decays, in particular the weak decays of heavy flavour particles such as B mesons.
EvtISGWFF.cpp
Go to the documentation of this file.
1 
2 /***********************************************************************
3 * Copyright 1998-2020 CERN for the benefit of the EvtGen authors *
4 * *
5 * This file is part of EvtGen. *
6 * *
7 * EvtGen is free software: you can redistribute it and/or modify *
8 * it under the terms of the GNU General Public License as published by *
9 * the Free Software Foundation, either version 3 of the License, or *
10 * (at your option) any later version. *
11 * *
12 * EvtGen is distributed in the hope that it will be useful, *
13 * but WITHOUT ANY WARRANTY; without even the implied warranty of *
14 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
15 * GNU General Public License for more details. *
16 * *
17 * You should have received a copy of the GNU General Public License *
18 * along with EvtGen. If not, see <https://www.gnu.org/licenses/>. *
19 ***********************************************************************/
20 
21 #include "EvtGenModels/EvtISGWFF.hh"
22 
23 #include "EvtGenBase/EvtId.hh"
24 #include "EvtGenBase/EvtPDL.hh"
25 #include "EvtGenBase/EvtPatches.hh"
26 #include "EvtGenBase/EvtReport.hh"
27 
28 #include <math.h>
29 #include <stdlib.h>
30 #include <string>
31 using std::endl;
32 
33 void EvtISGWFF::getscalarff( EvtId parent, EvtId daught, double t, double mass,
34  double* fpf, double* f0f )
35 {
36  //added by Lange Jan4,2000
37  static EvtId D0 = EvtPDL::getId( "D0" );
38  static EvtId D0B = EvtPDL::getId( "anti-D0" );
39  static EvtId DP = EvtPDL::getId( "D+" );
40  static EvtId DM = EvtPDL::getId( "D-" );
41 
42  static EvtId D3P0P = EvtPDL::getId( "D_0*+" );
43  static EvtId D3P0N = EvtPDL::getId( "D_0*-" );
44  static EvtId D3P00 = EvtPDL::getId( "D_0*0" );
45  static EvtId D3P0B = EvtPDL::getId( "anti-D_0*0" );
46 
47  static EvtId D21S0P = EvtPDL::getId( "hi" );
48  static EvtId D21S0N = EvtPDL::getId( "hi" );
49  static EvtId D21S00 = EvtPDL::getId( "hi" );
50  static EvtId D21S0B = EvtPDL::getId( "hi" );
51 
52  static EvtId ETA2S = EvtPDL::getId( "eta(2S)" );
53 
54  static EvtId PI2S0 = EvtPDL::getId( "pi(2S)0" );
55  static EvtId PI2SP = EvtPDL::getId( "pi(2S)+" );
56  static EvtId PI2SM = EvtPDL::getId( "pi(2S)-" );
57 
58  static EvtId PIP = EvtPDL::getId( "pi+" );
59  static EvtId PIM = EvtPDL::getId( "pi-" );
60  static EvtId PI0 = EvtPDL::getId( "pi0" );
61 
62  static EvtId A0P = EvtPDL::getId( "a_0+" );
63  static EvtId A0M = EvtPDL::getId( "a_0-" );
64  static EvtId A00 = EvtPDL::getId( "a_00" );
65 
66  static EvtId F0 = EvtPDL::getId( "f_0" );
67  static EvtId F0PR = EvtPDL::getId( "f'_0" );
68 
69  static EvtId ETA = EvtPDL::getId( "eta" );
70  static EvtId ETAPR = EvtPDL::getId( "eta'" );
71 
72  static EvtId KP = EvtPDL::getId( "K+" );
73  static EvtId KM = EvtPDL::getId( "K-" );
74  static EvtId K0 = EvtPDL::getId( "K0" );
75  static EvtId KB = EvtPDL::getId( "anti-K0" );
76  static EvtId K0S = EvtPDL::getId( "K_S0" );
77  static EvtId K0L = EvtPDL::getId( "K_L0" );
78 
79  static EvtId K0STP = EvtPDL::getId( "K_0*+" );
80  static EvtId K0STM = EvtPDL::getId( "K_0*-" );
81  static EvtId K0ST0 = EvtPDL::getId( "K_0*0" );
82  static EvtId K0STB = EvtPDL::getId( "anti-K_0*0" );
83 
84  static EvtId DSP = EvtPDL::getId( "D_s+" );
85  static EvtId DSM = EvtPDL::getId( "D_s-" );
86 
87  static EvtId D3P0SP = EvtPDL::getId( "D_s0*+" );
88  static EvtId D3P0SN = EvtPDL::getId( "D_s0*-" );
89 
90  double fmf;
91  double mb = EvtPDL::getMeanMass( parent );
92 
93  if ( daught == PI0 || daught == PIP || daught == PIM || daught == ETA ||
94  daught == ETAPR || daught == D0 || daught == D0B || daught == DP ||
95  daught == DM || daught == KP || daught == KM || daught == K0 ||
96  daught == K0L || daught == KB || daught == DSP || daught == DSM ||
97  daught == K0S ) {
98  EvtISGW1FF1S0( parent, daught, t, mass, fpf, &fmf );
99  }
100 
101  if ( daught == PI2S0 || daught == PI2SP || daught == PI2SM ||
102  daught == ETA2S || daught == D21S0P || daught == D21S0B ||
103  daught == D21S0N || daught == D21S00 ) {
104  EvtISGW1FF21S0( parent, daught, t, mass, fpf, &fmf );
105  }
106 
107  if ( daught == A00 || daught == A0P || daught == A0M || daught == F0 ||
108  daught == F0PR || daught == D3P0P || daught == D3P00 ||
109  daught == D3P0B || daught == D3P0N || daught == K0STM ||
110  daught == K0STB || daught == K0STP || daught == D3P0SP ||
111  daught == D3P0SN || daught == K0ST0 ) {
112  EvtISGW1FF3P0( parent, daught, t, mass, fpf, &fmf );
113  }
114 
115  *f0f = ( fmf / ( ( mb * mb - mass * mass ) / t ) ) + ( *fpf );
116 
117  return;
118 }
119 
120 void EvtISGWFF::gettensorff( EvtId parent, EvtId daught, double t, double mass,
121  double* hf, double* kf, double* bpf, double* bmf )
122 {
123  //added by Lange Jan4,2000
124  EvtISGW1FF3P2( parent, daught, t, mass, hf, kf, bpf, bmf );
125 
126  return;
127 }
128 
129 void EvtISGWFF::getvectorff( EvtId parent, EvtId daught, double t, double mass,
130  double* a1f, double* a2f, double* vf, double* a0f )
131 {
132  //added by Lange Jan4,2000
133  static EvtId DST0 = EvtPDL::getId( "D*0" );
134  static EvtId DSTB = EvtPDL::getId( "anti-D*0" );
135  static EvtId DSTP = EvtPDL::getId( "D*+" );
136  static EvtId DSTM = EvtPDL::getId( "D*-" );
137 
138  static EvtId D1P1P = EvtPDL::getId( "D_1+" );
139  static EvtId D1P1N = EvtPDL::getId( "D_1-" );
140  static EvtId D1P10 = EvtPDL::getId( "D_10" );
141  static EvtId D1P1B = EvtPDL::getId( "anti-D_10" );
142 
143  static EvtId D3P1P = EvtPDL::getId( "D'_1+" );
144  static EvtId D3P1N = EvtPDL::getId( "D'_1-" );
145  static EvtId D3P10 = EvtPDL::getId( "D'_10" );
146  static EvtId D3P1B = EvtPDL::getId( "anti-D'_10" );
147 
148  static EvtId D23S1P = EvtPDL::getId( "hi" );
149  static EvtId D23S1N = EvtPDL::getId( "hi" );
150  static EvtId D23S10 = EvtPDL::getId( "hi" );
151  static EvtId D23S1B = EvtPDL::getId( "hi" );
152 
153  static EvtId RHO2S0 = EvtPDL::getId( "rho(2S)0" );
154  static EvtId RHO2SP = EvtPDL::getId( "rho(2S)+" );
155  static EvtId RHO2SM = EvtPDL::getId( "rho(2S)-" );
156  static EvtId OMEG2S = EvtPDL::getId( "omega(2S)" );
157 
158  static EvtId RHOP = EvtPDL::getId( "rho+" );
159  static EvtId RHOM = EvtPDL::getId( "rho-" );
160  static EvtId RHO0 = EvtPDL::getId( "rho0" );
161 
162  static EvtId A1P = EvtPDL::getId( "a_1+" );
163  static EvtId A1M = EvtPDL::getId( "a_1-" );
164  static EvtId A10 = EvtPDL::getId( "a_10" );
165 
166  static EvtId B1P = EvtPDL::getId( "b_1+" );
167  static EvtId B1M = EvtPDL::getId( "b_1-" );
168  static EvtId B10 = EvtPDL::getId( "b_10" );
169 
170  static EvtId H1 = EvtPDL::getId( "h_1" );
171  static EvtId H1PR = EvtPDL::getId( "h'_1" );
172 
173  static EvtId F1 = EvtPDL::getId( "f_1" );
174  static EvtId F1PR = EvtPDL::getId( "f'_1" );
175 
176  static EvtId OMEG = EvtPDL::getId( "omega" );
177 
178  static EvtId KSTP = EvtPDL::getId( "K*+" );
179  static EvtId KSTM = EvtPDL::getId( "K*-" );
180  static EvtId KST0 = EvtPDL::getId( "K*0" );
181  static EvtId KSTB = EvtPDL::getId( "anti-K*0" );
182 
183  static EvtId K1P = EvtPDL::getId( "K_1+" );
184  static EvtId K1M = EvtPDL::getId( "K_1-" );
185  static EvtId K10 = EvtPDL::getId( "K_10" );
186  static EvtId K1B = EvtPDL::getId( "anti-K_10" );
187 
188  static EvtId K1STP = EvtPDL::getId( "K'_1+" );
189  static EvtId K1STM = EvtPDL::getId( "K'_1-" );
190  static EvtId K1ST0 = EvtPDL::getId( "K'_10" );
191  static EvtId K1STB = EvtPDL::getId( "anti-K'_10" );
192 
193  static EvtId PHI = EvtPDL::getId( "phi" );
194 
195  static EvtId D1P1SP = EvtPDL::getId( "D_s1+" );
196  static EvtId D1P1SN = EvtPDL::getId( "D_s1-" );
197 
198  static EvtId D3P1SP = EvtPDL::getId( "D'_s1*+" );
199  static EvtId D3P1SN = EvtPDL::getId( "D'_s1*-" );
200 
201  static EvtId DSSTP = EvtPDL::getId( "D_s*+" );
202  static EvtId DSSTM = EvtPDL::getId( "D_s*-" );
203 
204  double ff, gf, apf, amf;
205 
206  if ( daught == DST0 || daught == DSTP || daught == DSTM || daught == DSTB ||
207  daught == OMEG || daught == RHO0 || daught == RHOM || daught == RHOP ||
208  daught == KSTP || daught == KSTM || daught == KST0 || daught == KSTB ||
209  daught == PHI || daught == DSSTP || daught == DSSTM ) {
210  EvtISGW1FF3S1( parent, daught, t, mass, &ff, &gf, &apf, &amf );
211  }
212  if ( daught == B10 || daught == B1P || daught == B1M || daught == H1 ||
213  daught == H1PR || daught == D1P1P || daught == D1P10 || daught == D1P1B ||
214  daught == D1P1SP || daught == D1P1SN || daught == D1P1N ||
215  daught == K10 || daught == K1B || daught == K1P || daught == K1M ) {
216  EvtISGW1FF1P1( parent, daught, t, mass, &ff, &gf, &apf, &amf );
217  }
218  if ( daught == RHO2S0 || daught == RHO2SP || daught == RHO2SM ||
219  daught == OMEG2S || daught == D23S1P || daught == D23S1B ||
220  daught == D23S1N || daught == D23S10 ) {
221  EvtISGW1FF23S1( parent, daught, t, mass, &ff, &gf, &apf, &amf );
222  }
223  if ( daught == A10 || daught == A1P || daught == A1M || daught == F1 ||
224  daught == F1PR || daught == D3P1P || daught == D3P10 ||
225  daught == D3P1B || daught == D3P1N || daught == K1STM ||
226  daught == K1STB || daught == K1STP || daught == D3P1SP ||
227  daught == D3P1SN || daught == K1ST0 ) {
228  EvtISGW1FF3P1( parent, daught, t, mass, &ff, &gf, &apf, &amf );
229  }
230 
231  // Need to stuff in some factors to make these the ffs that
232  // is used elsewhere...
233 
234  double mb = EvtPDL::getMeanMass( parent );
235 
236  *vf = ( gf ) * ( mb + mass );
237  *a1f = ( ff ) / ( mb + mass );
238  *a2f = -1.0 * ( apf ) * ( mb + mass );
239  double a3f = ( ( mb + mass ) / ( 2.0 * mass ) ) * ( *a1f ) -
240  ( ( mb - mass ) / ( 2.0 * mass ) ) * ( *a2f );
241 
242  *a0f = a3f - ( ( t * amf ) / ( 2.0 * mass ) );
243 
244  return;
245 }
246 
247 void EvtISGWFF::EvtISGW1FF3P2( EvtId parent, EvtId daugt, double t, double mass,
248  double* hf, double* kf, double* bpf, double* bmf )
249 {
250  //added by Lange Jan4,2000
251  static EvtId BP = EvtPDL::getId( "B+" );
252  static EvtId BM = EvtPDL::getId( "B-" );
253  static EvtId B0 = EvtPDL::getId( "B0" );
254  static EvtId B0B = EvtPDL::getId( "anti-B0" );
255 
256  static EvtId D3P2P = EvtPDL::getId( "D_2*+" );
257  static EvtId D3P2N = EvtPDL::getId( "D_2*-" );
258  static EvtId D3P20 = EvtPDL::getId( "D_2*0" );
259  static EvtId D3P2B = EvtPDL::getId( "anti-D_2*0" );
260 
261  static EvtId A2P = EvtPDL::getId( "a_2+" );
262  static EvtId A2M = EvtPDL::getId( "a_2-" );
263  static EvtId A20 = EvtPDL::getId( "a_20" );
264 
265  static EvtId F2 = EvtPDL::getId( "f_2" );
266  static EvtId F2PR = EvtPDL::getId( "f'_2" );
267 
268  double mtb;
269  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
270  double msq( 0.0 ), bx2( 0.0 ), mtx, f5;
271  double mum, mup, tm, bb2( 0.0 ), bbx2;
272  double msb( 0.0 ), kap;
273 
274  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
275  msb = 5.2;
276  msd = 0.33;
277  bb2 = 0.41 * 0.41;
278  if ( daugt == A20 || daugt == A2P || daugt == A2M || daugt == F2 ||
279  daugt == F2PR ) {
280  msq = 0.33;
281  bx2 = 0.27 * 0.27;
282  } else {
283  if ( daugt == D3P2P || daugt == D3P2N || daugt == D3P2B ||
284  daugt == D3P20 ) {
285  msq = 1.82;
286  bx2 = 0.34 * 0.34;
287  } else {
288  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
289  << "Not implemented daugt in get_isgw_ff_3P1.\n";
290  }
291  }
292  } else {
293  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
294  << "Not implemented parent in get_isgw_ff_3P1.\n";
295  }
296 
297  mtb = msb + msd;
298  mtx = msq + msd;
299 
300  mb = EvtPDL::getMeanMass( parent );
301  mx = mass;
302 
303  mup = 1.0 / ( 1.0 / msq + 1.0 / msb );
304  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
305  bbx2 = 0.5 * ( bb2 + bx2 );
306 
307  tm = ( mb - mx ) * ( mb - mx );
308  if ( t > tm )
309  t = 0.99 * tm;
310  kap = 0.7 * 0.7;
311 
312  f5 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 5.0 / 2.0 ) *
313  exp( -1.0 *
314  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
315 
316  *hf = f5 * ( msd / ( sqrt( 8.0 * bb2 ) * mtb ) ) *
317  ( ( 1.0 / msq ) - ( msd * bb2 / ( 2.0 * mum * mtx * bbx2 ) ) );
318 
319  *kf = f5 * msd * sqrt( 2.0 / bb2 );
320 
321  *bpf = ( -1.0 * f5 * msd / ( sqrt( 8.0 * bb2 ) * msb * mtx ) ) *
322  ( 1.0 - ( msd * msb * bx2 / ( 2.0 * mup * mtb * bbx2 ) ) +
323  ( msd * msb * bx2 * ( 1.0 - ( msd * bx2 / ( 2.0 * mtb * bbx2 ) ) ) /
324  ( 4.0 * mtb * mum * bbx2 ) ) );
325  *bmf = 0.0;
326  return;
327 } //get_ff_isgw_1p1
328 
329 void EvtISGWFF::EvtISGW1FF1S0( EvtId parent, EvtId daugt, double t, double mass,
330  double* fpf, double* fmf )
331 {
332  //added by Lange Jan4,2000
333  static EvtId BP = EvtPDL::getId( "B+" );
334  static EvtId BM = EvtPDL::getId( "B-" );
335  static EvtId B0 = EvtPDL::getId( "B0" );
336  static EvtId B0B = EvtPDL::getId( "anti-B0" );
337 
338  static EvtId D0 = EvtPDL::getId( "D0" );
339  static EvtId D0B = EvtPDL::getId( "anti-D0" );
340  static EvtId DP = EvtPDL::getId( "D+" );
341  static EvtId DM = EvtPDL::getId( "D-" );
342 
343  static EvtId PIP = EvtPDL::getId( "pi+" );
344  static EvtId PIM = EvtPDL::getId( "pi-" );
345  static EvtId PI0 = EvtPDL::getId( "pi0" );
346 
347  static EvtId ETA = EvtPDL::getId( "eta" );
348  static EvtId ETAPR = EvtPDL::getId( "eta'" );
349 
350  double mtb;
351  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
352  double msq( 0.0 ), bx2( 0.0 ), mtx;
353  double f3, kap;
354  double msb( 0.0 ), bb2( 0.0 ), mup, mum, bbx2, tm;
355 
356  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
357  msb = 5.2;
358  msd = 0.33;
359  bb2 = 0.41 * 0.41;
360  if ( daugt == PIP || daugt == PIM || daugt == PI0 || daugt == ETA ||
361  daugt == ETAPR ) {
362  msq = 0.33;
363  bx2 = 0.31 * 0.31;
364  } else {
365  if ( daugt == D0 || daugt == DP || daugt == DM || daugt == D0B ) {
366  msq = 1.82;
367  bx2 = 0.39 * 0.39;
368  } else {
369  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
370  << "Not implemented daugt in get_isgw_ff_1S0.\n";
371  }
372  }
373  } else {
374  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
375  << "Not implemented parent in get_isgw_ff_1S0.\n";
376  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
377  << "Parent:" << parent.getId() << endl;
378  }
379 
380  mtb = msb + msd;
381  mtx = msq + msd;
382  mb = EvtPDL::getMeanMass( parent );
383  mx = mass;
384  mup = 1.0 / ( 1.0 / msq + 1.0 / msb );
385  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
386  bbx2 = 0.5 * ( bb2 + bx2 );
387  tm = ( mb - mx ) * ( mb - mx );
388  if ( t > tm )
389  t = 0.99 * tm;
390 
391  kap = 0.7 * 0.7;
392  f3 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 3.0 / 2.0 ) *
393  exp( -1.0 *
394  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
395 
396  *fpf = f3 * ( 1 + ( msb / ( 2.0 * mum ) ) -
397  ( msb * msq * msd * bb2 / ( 4.0 * mup * mum * mtx * bbx2 ) ) );
398  *fmf = f3 * ( 1.0 - ( mtb + mtx ) *
399  ( 0.5 / msq -
400  ( msd * bb2 / ( 4.0 * mup * mtx * bbx2 ) ) ) );
401 
402  return;
403 } //get_ff_isgw_1s0
404 
405 void EvtISGWFF::EvtISGW1FF3S1( EvtId parent, EvtId daugt, double t, double mass,
406  double* f, double* g, double* ap, double* am )
407 {
408  //added by Lange Jan4,2000
409  static EvtId BP = EvtPDL::getId( "B+" );
410  static EvtId BM = EvtPDL::getId( "B-" );
411  static EvtId B0 = EvtPDL::getId( "B0" );
412  static EvtId B0B = EvtPDL::getId( "anti-B0" );
413 
414  static EvtId DST0 = EvtPDL::getId( "D*0" );
415  static EvtId DSTB = EvtPDL::getId( "anti-D*0" );
416  static EvtId DSTP = EvtPDL::getId( "D*+" );
417  static EvtId DSTM = EvtPDL::getId( "D*-" );
418 
419  static EvtId RHOP = EvtPDL::getId( "rho+" );
420  static EvtId RHOM = EvtPDL::getId( "rho-" );
421  static EvtId RHO0 = EvtPDL::getId( "rho0" );
422 
423  static EvtId OMEG = EvtPDL::getId( "omega" );
424 
425  double msd( 0.0 ), msq( 0.0 ), bb2( 0.0 ), mum, mtx, bbx2;
426  double bx2( 0.0 ), msb( 0.0 ), tm;
427  double mb, mx, f3, kap;
428 
429  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
430  msb = 5.2;
431  msd = 0.33;
432  bb2 = 0.41 * 0.41;
433  if ( daugt == DSTP || daugt == DSTM || daugt == DSTB || daugt == DST0 ) {
434  msq = 1.82;
435  bx2 = 0.39 * 0.39;
436  } else {
437  if ( daugt == RHOP || daugt == RHOM || daugt == RHO0 ||
438  daugt == OMEG ) {
439  msq = 0.33;
440  bx2 = 0.31 * 0.31;
441  } else {
442  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
443  << "Not implemented daugt in get_isgw_ff_3S1.\n";
444  }
445  }
446  } else {
447  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
448  << "Not implemented parent in get_isgw_ff_3S1.\n";
449  }
450 
451  double mtb;
452 
453  mtb = msb + msd;
454  mtx = msq + msd;
455 
456  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
457  bbx2 = 0.5 * ( bb2 + bx2 );
458  mb = EvtPDL::getMeanMass( parent );
459  mx = mass;
460  tm = ( mb - mx ) * ( mb - mx );
461  if ( t > tm )
462  t = 0.99 * tm;
463 
464  kap = 0.7 * 0.7;
465  f3 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 3.0 / 2.0 ) *
466  exp( -1.0 *
467  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
468 
469  *f = 2.0 * mtb * f3;
470  *g = 0.5 * f3 * ( ( 1 / msq ) - ( msd * bb2 / ( 2.0 * mum * mtx * bbx2 ) ) );
471  *ap = ( -1.0 * f3 / ( 2.0 * mtx ) ) *
472  ( 1.0 + ( msd * ( bb2 - bx2 ) / ( msb * ( bb2 + bx2 ) ) ) -
473  ( msd * msd * bx2 * bx2 / ( 4.0 * mum * mtb * bbx2 * bbx2 ) ) );
474  *am = 0.0;
475 }
476 
477 void EvtISGWFF::EvtISGW1FF23S1( EvtId parent, EvtId daugt, double t,
478  double mass, double* fpf, double* gpf,
479  double* appf, double* apmf )
480 {
481  //added by Lange Jan4,2000
482  static EvtId BP = EvtPDL::getId( "B+" );
483  static EvtId BM = EvtPDL::getId( "B-" );
484  static EvtId B0 = EvtPDL::getId( "B0" );
485  static EvtId B0B = EvtPDL::getId( "anti-B0" );
486 
487  static EvtId D23S1P = EvtPDL::getId( "hi" );
488  static EvtId D23S1N = EvtPDL::getId( "hi" );
489  static EvtId D23S10 = EvtPDL::getId( "hi" );
490  static EvtId D23S1B = EvtPDL::getId( "hi" );
491 
492  static EvtId RHO2S0 = EvtPDL::getId( "rho(2S)0" );
493  static EvtId RHO2SP = EvtPDL::getId( "rho(2S)+" );
494  static EvtId RHO2SM = EvtPDL::getId( "rho(2S)-" );
495  static EvtId OMEG2S = EvtPDL::getId( "omega(2S)" );
496 
497  double mtb;
498  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
499  double msq( 0.0 ), bx2( 0.0 ), mtx;
500  double f3, f5, tt;
501  double mum, bb2( 0.0 ), bbx2, tm, msb( 0.0 );
502 
503  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
504  msb = 5.2;
505  msd = 0.33;
506  bb2 = 0.41 * 0.41;
507  if ( daugt == RHO2SP || daugt == RHO2SM || daugt == RHO2S0 ||
508  daugt == OMEG2S ) {
509  msq = 0.33;
510  bx2 = 0.31 * 0.31;
511  } else {
512  if ( daugt == D23S1N || daugt == D23S10 || daugt == D23S1P ||
513  daugt == D23S1B ) {
514  msq = 1.82;
515  bx2 = 0.39 * 0.39;
516  } else {
517  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
518  << "Not implemented daugt in get_isgw_ff_23P1.\n";
519  }
520  }
521  } else {
522  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
523  << "Not implemented parent in get_isgw_ff_23P1.\n";
524  }
525 
526  mtb = msb + msd;
527  mtx = msq + msd;
528  mb = EvtPDL::getMeanMass( parent );
529  mx = mass;
530  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
531  bbx2 = 0.5 * ( bb2 + bx2 );
532  tm = ( mb - mx ) * ( mb - mx );
533  if ( t > tm )
534  t = 0.99 * tm;
535 
536  double kap = 0.7 * 0.7;
537  f3 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 3.0 / 2.0 ) *
538  exp( -1.0 *
539  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
540 
541  f5 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 5.0 / 2.0 ) *
542  exp( -1.0 *
543  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
544 
545  *fpf = sqrt( 6.0 ) * f3 * mtb *
546  ( ( ( bb2 - bx2 ) / ( bb2 + bx2 ) ) +
547  ( ( msd * msd * bx2 * ( tm - t ) ) /
548  ( 6.0 * mtx * mtb * bbx2 * kap * bbx2 ) ) );
549 
550  *gpf = sqrt( 3.0 / 8.0 ) * f3 *
551  ( ( ( ( ( bb2 - bx2 ) / ( bb2 + bx2 ) ) +
552  ( ( msd * msd * bx2 * ( tm - t ) ) /
553  ( 6.0 * mtx * mtb * bbx2 * kap * bbx2 ) ) ) *
554  ( ( 1.0 / msq ) - ( ( msd * bb2 ) / ( 2.0 * mum * mtx * bbx2 ) ) ) ) +
555  ( ( msd * bb2 * bx2 ) / ( 3.0 * mum * mtx * bbx2 * bbx2 ) ) );
556 
557  tt = ( msd * msd * bx2 * ( tm - t ) ) / ( mtx * mtb * bb2 * kap * bbx2 );
558 
559  *appf = ( f5 / ( sqrt( 6.0 ) * mtx ) ) *
560  ( ( ( 3.0 * mtb * bbx2 / ( 2.0 * msb * sqrt( bb2 * bx2 ) ) ) *
561  ( 1.0 - ( ( msd * msd * msb * bx2 * bx2 ) /
562  ( 4.0 * mtb * mtb * mum * bbx2 * bbx2 ) ) ) ) -
563  ( ( 3.0 * msd * sqrt( bx2 / bb2 ) ) / ( 2.0 * msb ) ) +
564  ( ( 5.0 * msd * sqrt( bx2 * bb2 ) * ( 1.0 + 0.1 * tt ) ) /
565  ( 2.0 * msb * bbx2 ) ) -
566  ( ( 3.0 * mtb * sqrt( bb2 / bx2 ) * ( 1.0 + ( tt / 6.0 ) ) ) /
567  ( 2.0 * msb ) ) +
568  ( ( 7.0 * msd * msd * sqrt( bb2 / bx2 ) * bx2 * bx2 *
569  ( 1.0 + ( tt / 14.0 ) ) ) /
570  ( 8.0 * mtb * mum * bbx2 * bbx2 ) ) );
571 
572  *apmf = 0.0;
573  return;
574 } //get_ff_isgw_23s1
575 
576 void EvtISGWFF::EvtISGW1FF3P1( EvtId parent, EvtId daugt, double t, double mass,
577  double* lf, double* qf, double* cpf, double* cmf )
578 {
579  //added by Lange Jan4,2000
580  static EvtId BP = EvtPDL::getId( "B+" );
581  static EvtId BM = EvtPDL::getId( "B-" );
582  static EvtId B0 = EvtPDL::getId( "B0" );
583  static EvtId B0B = EvtPDL::getId( "anti-B0" );
584 
585  static EvtId D3P1P = EvtPDL::getId( "D'_1+" );
586  static EvtId D3P1N = EvtPDL::getId( "D'_1-" );
587  static EvtId D3P10 = EvtPDL::getId( "D'_10" );
588  static EvtId D3P1B = EvtPDL::getId( "anti-D'_10" );
589 
590  static EvtId A1P = EvtPDL::getId( "a_1+" );
591  static EvtId A1M = EvtPDL::getId( "a_1-" );
592  static EvtId A10 = EvtPDL::getId( "a_10" );
593 
594  static EvtId F1 = EvtPDL::getId( "f_1" );
595  static EvtId F1PR = EvtPDL::getId( "f'_1" );
596 
597  double mtb;
598  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
599  double msq( 0.0 ), bx2( 0.0 ), mtx, f5;
600  double msb( 0.0 ), bb2( 0.0 ), mum, bbx2, tm;
601  double kap;
602 
603  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
604  msb = 5.2;
605  msd = 0.33;
606  bb2 = 0.41 * 0.41;
607  if ( daugt == A10 || daugt == A1P || daugt == A1M || daugt == F1 ||
608  daugt == F1PR ) {
609  msq = 0.33;
610  bx2 = 0.27 * 0.27;
611  } else {
612  if ( daugt == D3P1P || daugt == D3P1N || daugt == D3P1B ||
613  daugt == D3P10 ) {
614  msq = 1.82;
615  bx2 = 0.34 * 0.34;
616  } else {
617  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
618  << "Not implemented daugt in get_isgw_ff_3P1.\n";
619  }
620  }
621  } else {
622  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
623  << "Not implemented parent in get_isgw_ff_3P1.\n";
624  }
625 
626  mtb = msb + msd;
627  mtx = msq + msd;
628 
629  mb = EvtPDL::getMeanMass( parent );
630  mx = mass;
631 
632  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
633  bbx2 = 0.5 * ( bb2 + bx2 );
634  tm = ( mb - mx ) * ( mb - mx );
635  if ( t > tm )
636  t = 0.99 * tm;
637 
638  kap = 0.7 * 0.7;
639  f5 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 5.0 / 2.0 ) *
640  exp( -1.0 *
641  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
642 
643  *qf = ( f5 * msd ) / ( 2.0 * mtx * sqrt( bb2 ) );
644 
645  *lf = -1.0 * mtb * sqrt( bb2 ) * f5 *
646  ( 1 / mum + ( msd * ( tm - t ) / ( 2.0 * mtb * kap * bb2 ) ) *
647  ( ( 1.0 / msq ) -
648  ( 1.0 * msd * bb2 / ( 2.0 * mum * mtx * bbx2 ) ) ) );
649 
650  *cpf = ( f5 * msd * msb / ( 4.0 * mtb * sqrt( bb2 ) * mum ) ) *
651  ( 1.0 - ( msd * msq * bb2 / ( 2.0 * mtx * mum * bbx2 ) ) );
652  *cmf = 0.0;
653  return;
654 } //get_ff_isgw_3p1
655 
656 void EvtISGWFF::EvtISGW1FF3P0( EvtId parent, EvtId daugt, double t, double mass,
657  double* upf, double* umf )
658 {
659  //added by Lange Jan4,2000
660  static EvtId BP = EvtPDL::getId( "B+" );
661  static EvtId BM = EvtPDL::getId( "B-" );
662  static EvtId B0 = EvtPDL::getId( "B0" );
663  static EvtId B0B = EvtPDL::getId( "anti-B0" );
664 
665  static EvtId D3P0P = EvtPDL::getId( "D_0*+" );
666  static EvtId D3P0N = EvtPDL::getId( "D_0*-" );
667  static EvtId D3P00 = EvtPDL::getId( "D_0*0" );
668  static EvtId D3P0B = EvtPDL::getId( "anti-D_0*0" );
669 
670  static EvtId A0P = EvtPDL::getId( "a_0+" );
671  static EvtId A0M = EvtPDL::getId( "a_0-" );
672  static EvtId A00 = EvtPDL::getId( "a_00" );
673 
674  static EvtId F0 = EvtPDL::getId( "f_0" );
675  static EvtId F0PR = EvtPDL::getId( "f'_0" );
676 
677  double mtb;
678  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
679  double msq( 0.0 ), bx2( 0.0 ), mtx;
680  double f5;
681  double mum, bb2( 0.0 ), bbx2, msb( 0.0 ), tm;
682 
683  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
684  msb = 5.2;
685  msd = 0.33;
686  bb2 = 0.41 * 0.41;
687  if ( daugt == A00 || daugt == A0P || daugt == A0M || daugt == F0 ||
688  daugt == F0PR ) {
689  msq = 0.33;
690  bx2 = 0.27 * 0.27;
691  } else {
692  if ( daugt == D3P0P || daugt == D3P0N || daugt == D3P0B ||
693  daugt == D3P00 ) {
694  msq = 1.82;
695  bx2 = 0.34 * 0.34;
696  } else {
697  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
698  << "Not implemented daugt in get_isgw_ff_3P0.\n";
699  }
700  }
701  } else {
702  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
703  << "Not implemented parent in get_isgw_ff_3P0.\n";
704  }
705 
706  mtb = msb + msd;
707  mtx = msq + msd;
708 
709  mb = EvtPDL::getMeanMass( parent );
710  mx = mass;
711 
712  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
713  bbx2 = 0.5 * ( bb2 + bx2 );
714  tm = ( mb - mx ) * ( mb - mx );
715  if ( t > tm )
716  t = 0.99 * tm;
717 
718  double kap = 0.7 * 0.7;
719  f5 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 5.0 / 2.0 ) *
720  exp( -1.0 *
721  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
722 
723  *upf = f5 * msd * msq * msb / ( sqrt( 6.0 * bb2 ) * mtx * mum );
724  *umf = 0.0;
725  return;
726 } //get_ff_isgw_3p0
727 
728 void EvtISGWFF::EvtISGW1FF1P1( EvtId parent, EvtId daugt, double t, double mass,
729  double* vf, double* rf, double* spf, double* smf )
730 {
731  //added by Lange Jan4,2000
732  static EvtId BP = EvtPDL::getId( "B+" );
733  static EvtId BM = EvtPDL::getId( "B-" );
734  static EvtId B0 = EvtPDL::getId( "B0" );
735  static EvtId B0B = EvtPDL::getId( "anti-B0" );
736 
737  static EvtId D1P1P = EvtPDL::getId( "D_1+" );
738  static EvtId D1P1N = EvtPDL::getId( "D_1-" );
739  static EvtId D1P10 = EvtPDL::getId( "D_10" );
740  static EvtId D1P1B = EvtPDL::getId( "anti-D_10" );
741 
742  static EvtId B1P = EvtPDL::getId( "b_1+" );
743  static EvtId B1M = EvtPDL::getId( "b_1-" );
744  static EvtId B10 = EvtPDL::getId( "b_10" );
745 
746  static EvtId H1 = EvtPDL::getId( "h_1" );
747  static EvtId H1PR = EvtPDL::getId( "h'_1" );
748 
749  double mtb;
750  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
751  double msq( 0.0 ), bx2( 0.0 ), mtx, f5;
752  double mup, mum, kap;
753  double msb( 0.0 ), bb2( 0.0 ), bbx2, tm;
754 
755  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
756  msb = 5.2;
757  msd = 0.33;
758  bb2 = 0.41 * 0.41;
759  if ( daugt == H1 || daugt == H1PR || daugt == B10 || daugt == B1P ||
760  daugt == B1M ) {
761  msq = 0.33;
762  bx2 = 0.27 * 0.27;
763  } else {
764  if ( daugt == D1P1P || daugt == D1P1N || daugt == D1P10 ||
765  daugt == D1P1B ) {
766  msq = 1.82;
767  bx2 = 0.34 * 0.34;
768  } else {
769  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
770  << "Not implemented daugt in get_isgw_ff_3P1.\n";
771  }
772  }
773  } else {
774  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
775  << "Not implemented parent in get_isgw_ff_3P1.\n";
776  }
777 
778  mtb = msb + msd;
779  mtx = msq + msd;
780 
781  mb = EvtPDL::getMeanMass( parent );
782  mx = mass;
783 
784  mup = 1.0 / ( 1.0 / msq + 1.0 / msb );
785  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
786  bbx2 = 0.5 * ( bb2 + bx2 );
787  tm = ( mb - mx ) * ( mb - mx );
788  if ( t > tm )
789  t = 0.99 * tm;
790 
791  kap = 0.7 * 0.7;
792  f5 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 5.0 / 2.0 ) *
793  exp( -1.0 *
794  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
795 
796  *vf = f5 *
797  ( ( ( mtb * sqrt( bb2 ) ) / ( 4.0 * sqrt( 2.0 ) * msb * msq * mtx ) ) );
798  *rf = f5 * mtb * sqrt( bb2 / 2 ) * ( ( 1.0 / mup ) );
799 
800  *spf = ( f5 * msd / ( sqrt( 2.0 * bb2 ) * mtb ) ) *
801  ( 1.0 + ( msb / ( 2.0 * mum ) ) -
802  ( msb * msq * msd * bb2 / ( 4.0 * mup * mum * mtx * bbx2 ) ) );
803  *smf = 0.0;
804 
805  return;
806  //get_ff_isgw_1p1
807 }
808 
809 void EvtISGWFF::EvtISGW1FF21S0( EvtId parent, EvtId daugt, double t,
810  double mass, double* fppf, double* fpmf )
811 {
812  //added by Lange Jan4,2000
813  static EvtId BP = EvtPDL::getId( "B+" );
814  static EvtId BM = EvtPDL::getId( "B-" );
815  static EvtId B0 = EvtPDL::getId( "B0" );
816  static EvtId B0B = EvtPDL::getId( "anti-B0" );
817 
818  static EvtId D21S0P = EvtPDL::getId( "hi" );
819  static EvtId D21S0N = EvtPDL::getId( "hi" );
820  static EvtId D21S00 = EvtPDL::getId( "hi" );
821  static EvtId D21S0B = EvtPDL::getId( "hi" );
822 
823  static EvtId ETA2S = EvtPDL::getId( "eta(2S)" );
824 
825  static EvtId PI2S0 = EvtPDL::getId( "pi(2S)0" );
826  static EvtId PI2SP = EvtPDL::getId( "pi(2S)+" );
827  static EvtId PI2SM = EvtPDL::getId( "pi(2S)-" );
828 
829  double mtb;
830  double msd( 0.0 ), mx( 0.0 ), mb( 0.0 );
831  double msq( 0.0 ), bx2( 0.0 ), mtx;
832  double f3;
833  double msb( 0.0 );
834  double mum, mup, tm, bb2( 0.0 ), bbx2;
835 
836  if ( parent == BM || parent == BP || parent == B0 || parent == B0B ) {
837  msb = 5.2;
838  msd = 0.33;
839  bb2 = 0.41 * 0.41;
840  if ( daugt == PI2S0 || daugt == PI2SP || daugt == PI2SM ||
841  daugt == ETA2S ) {
842  msq = 0.33;
843  bx2 = 0.31 * 0.31;
844  } else {
845  if ( daugt == D21S00 || daugt == D21S0P || daugt == D21S0N ||
846  daugt == D21S0B ) {
847  msq = 1.82;
848  bx2 = 0.39 * 0.39;
849  } else {
850  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
851  << "Not implemented daugt in get_isgw1_ff_21S0.\n";
852  }
853  }
854  } else {
855  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
856  << "Not implemented parent in get_isgw1_ff_21S0.\n";
857  }
858 
859  mtb = msb + msd;
860  mtx = msq + msd;
861 
862  mb = EvtPDL::getMeanMass( parent );
863  mx = mass;
864 
865  mup = 1.0 / ( 1.0 / msq + 1.0 / msb );
866  mum = 1.0 / ( 1.0 / msq - 1.0 / msb );
867  bbx2 = 0.5 * ( bb2 + bx2 );
868  tm = ( mb - mx ) * ( mb - mx );
869  if ( t > tm )
870  t = 0.99 * tm;
871 
872  double kap = 0.7 * 0.7;
873  f3 = sqrt( mtx / mtb ) * pow( sqrt( bx2 * bb2 ) / bbx2, 3.0 / 2.0 ) *
874  exp( -1.0 *
875  ( ( msd * msd * ( tm - t ) / ( 4.0 * mtb * mtx * kap * bbx2 ) ) ) );
876 
877  *fppf = f3 * sqrt( 3.0 / 8.0 ) * ( msb / mup ) *
878  ( ( ( bb2 - bx2 ) / ( bb2 + bx2 ) ) +
879  ( ( ( msq * msd * bb2 ) / ( 3.0 * mum * mtx * bbx2 ) ) *
880  ( ( 7.0 * bx2 - 3.0 * bb2 ) / ( 4.0 * bbx2 ) ) ) +
881  ( ( ( msd * msd * bx2 * ( tm - t ) ) /
882  ( 6.0 * mtx * mtb * bbx2 * kap * bbx2 ) ) *
883  ( 1.0 - ( ( msq * msd * bb2 ) / ( 2.0 * mum * mtx * bbx2 ) ) ) ) );
884 
885  *fpmf = 0.0;
886  return;
887 } //get_ff_isgw_21s0
888 
889 void EvtISGWFF::getbaryonff( EvtId, EvtId, double, double, double*, double*,
890  double*, double* )
891 {
892  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
893  << "Not implemented :getbaryonff in EvtISGWFF.\n";
894  ::abort();
895 }
896 
897 void EvtISGWFF::getdiracff( EvtId, EvtId, double, double, double*, double*,
898  double*, double*, double*, double* )
899 {
900  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
901  << "Not implemented :getdiracff in EvtISGWFF.\n";
902  ::abort();
903 }
904 
905 void EvtISGWFF::getraritaff( EvtId, EvtId, double, double, double*, double*,
906  double*, double*, double*, double*, double*, double* )
907 {
908  EvtGenReport( EVTGEN_ERROR, "EvtGen" )
909  << "Not implemented :getraritaff in EvtISGWFF.\n";
910  ::abort();
911 }
EvtGenReport
std::ostream & EvtGenReport(EvtGenSeverity severity, const char *facility=0)
Definition: EvtReport.cpp:33
EvtISGWFF::EvtISGW1FF21S0
void EvtISGW1FF21S0(EvtId parent, EvtId daught, double t, double mass, double *fppf, double *fpmf)
Definition: EvtISGWFF.cpp:809
EvtISGWFF::gettensorff
void gettensorff(EvtId parent, EvtId daught, double t, double mass, double *hf, double *kf, double *bpf, double *bmf) override
Definition: EvtISGWFF.cpp:120
EvtISGWFF::EvtISGW1FF1S0
void EvtISGW1FF1S0(EvtId parent, EvtId daught, double t, double mass, double *fpf, double *fmf)
Definition: EvtISGWFF.cpp:329
EvtPDL::getMeanMass
static double getMeanMass(EvtId i)
Definition: EvtPDL.cpp:314
EvtISGWFF::EvtISGW1FF3P0
void EvtISGW1FF3P0(EvtId parent, EvtId daught, double t, double mass, double *upf, double *umf)
Definition: EvtISGWFF.cpp:656
EvtId
Definition: EvtId.hh:27
EvtISGWFF::getdiracff
void getdiracff(EvtId, EvtId, double, double, double *, double *, double *, double *, double *, double *) override
Definition: EvtISGWFF.cpp:897
EvtPDL::getId
static EvtId getId(const std::string &name)
Definition: EvtPDL.cpp:287
EvtISGWFF::EvtISGW1FF3P1
void EvtISGW1FF3P1(EvtId parent, EvtId daught, double t, double mass, double *lf, double *qf, double *cpf, double *cmf)
Definition: EvtISGWFF.cpp:576
EvtISGWFF::getvectorff
void getvectorff(EvtId parent, EvtId daught, double t, double mass, double *a1f, double *a2f, double *vf, double *a0f) override
Definition: EvtISGWFF.cpp:129
EvtISGWFF::getscalarff
void getscalarff(EvtId parent, EvtId daught, double t, double mass, double *fpf, double *f0f) override
Definition: EvtISGWFF.cpp:33
EvtId::getId
int getId() const
Definition: EvtId.hh:42
EvtISGWFF::getraritaff
void getraritaff(EvtId, EvtId, double, double, double *, double *, double *, double *, double *, double *, double *, double *) override
Definition: EvtISGWFF.cpp:905
exp
EvtComplex exp(const EvtComplex &c)
Definition: EvtComplex.hh:240
EvtISGWFF::EvtISGW1FF3S1
void EvtISGW1FF3S1(EvtId parent, EvtId daught, double t, double mass, double *ff, double *gf, double *apf, double *amf)
Definition: EvtISGWFF.cpp:405
EvtISGWFF::EvtISGW1FF1P1
void EvtISGW1FF1P1(EvtId parent, EvtId daught, double t, double mass, double *vf, double *rf, double *spf, double *smf)
Definition: EvtISGWFF.cpp:728
EvtISGWFF::EvtISGW1FF3P2
void EvtISGW1FF3P2(EvtId parent, EvtId daught, double t, double mass, double *h, double *k, double *bp, double *bm)
Definition: EvtISGWFF.cpp:247
EvtISGWFF::EvtISGW1FF23S1
void EvtISGW1FF23S1(EvtId parent, EvtId daught, double t, double mass, double *fpf, double *gpf, double *app, double *apm)
Definition: EvtISGWFF.cpp:477
EvtISGWFF::getbaryonff
void getbaryonff(EvtId, EvtId, double, double, double *, double *, double *, double *) override
Definition: EvtISGWFF.cpp:889
Generated by   doxygen 1.8.15