50 <<
" ERROR: EvtLb2Lll generator expected max. 8 arguments but found: " 53 <<
" 1. Lambda_b0 polarization - zero is default" << std::endl;
55 <<
" 2. Model type - \"SM\" for Standard Model is default" 58 <<
" 3. Form-Factors - \"HQET\" is used by default" << std::endl;
60 <<
" 4. How to set polarization - \"ModifiedSpinors\" is default" 63 <<
" 5. Include long distance (LD) effects - \"SD\" (no) is default" 66 <<
" 6. NonFactorizable contribution (omega) to b->sg decay at q2=0 " 69 <<
" 7. Note on every x-th decay" << std::endl;
71 <<
" 8. Maximum probability - automatic by default" << std::endl;
77 <<
" ERROR: EvtLb2Lll generator expected 3 daughters but found: " 94 <<
" EvtLb2Lll generator found Lambda_b0" << std::endl;
102 <<
" EvtLb2Lll generator found decay: Lambda_b0 -> Lambda0 e- e+" 107 <<
" EvtLb2Lll generator found decay: Lambda_b0 -> Lambda0 mu- mu+" 112 <<
" EvtLb2Lll generator found decay: Lambda_b0 -> Lambda0 tau- tau+" 116 <<
" ERROR: EvtLb2Lll generator expected lepton pair daughters but found: " 124 <<
" EvtLb2Lll generator found anti-Lambda_b0" << std::endl;
130 m_decayName =
"anti-Lambda_b0 -> anti-Lambda0 e+ e-";
132 <<
" EvtLb2Lll generator found decay: anti-Lambda_b0 -> anti-Lambda0 e+ e-" 135 m_decayName =
"anti-Lambda_b0 -> anti-Lambda0 mu+ mu-";
137 <<
" EvtLb2Lll generator found decay: anti-Lambda_b0 -> anti-Lambda0 mu+ mu-" 140 m_decayName =
"anti-Lambda_b0 -> anti-Lambda0 tau+ tau-";
142 <<
" EvtLb2Lll generator found decay: anti-Lambda_b0 -> anti-Lambda0 tau+ tau-" 146 <<
" ERROR: EvtLb2Lll generator expected lepton pair daughters but found: " 153 <<
" ERROR: EvtLb2Lll generator expected (anti-)Lambda_b0 parent but found: " 162 <<
" ERROR: EvtLb2Lll expects polarization to be in interval <-1,1>, not " 163 <<
getArg( 0 ) << std::endl;
180 <<
" ERROR: EvtLb2Lll doesn't know this physics model: " 189 <<
" EvtLb2Lll will use this physics model: " <<
m_HEPmodel << std::endl;
194 getArgStr( 2 ).substr( 1, 11 ) !=
"HQET-delta=" ) {
196 <<
" ERROR: EvtLb2Lll doesn't know this Form-Factors model: " 205 <<
" EvtLb2Lll will use this Form-Factors model: " <<
m_FFtype 212 <<
" ERROR: EvtLb2Lll doesn't know kind of introducing polarization: " 222 <<
" EvtLb2Lll will use this kind of introducing polarization: " 229 <<
" ERROR: EvtLb2Lll didn't find SD or LD parameter: " 239 <<
" EvtLb2Lll will include contribution from these effects: " 243 if ( fabs(
getArg( 5 ) ) > 0.15 ) {
245 <<
" WARNING: EvtLb2Lll found very high contribution to b->sg decay at q2=0: " 246 <<
getArg( 5 ) << std::endl;
253 <<
" EvtLb2Lll will use this contribution to b->sg decay at q2=0: " 264 <<
" ERROR: EvtLb2Lll expects positive maximum probability not : " 265 <<
getArg( 7 ) << std::endl;
285 <<
" EvtLb2Lll is finding maximum probability ... " << std::endl;
313 double q2, pstar, elambda, theta;
314 double q2min = ( m1 + m2 ) * ( m1 + m2 );
315 double q2max = ( M0 - mL ) * ( M0 - mL );
320 <<
" EvtLb2Lll is probing whole phase space ..." << std::endl;
324 for ( i = 0; i <= 100; i++ ) {
325 q2 = q2min + i * ( q2max - q2min ) / 100.;
326 elambda = ( M0 * M0 + mL * mL - q2 ) / 2 / M0;
330 pstar = sqrt( q2 - ( m1 + m2 ) * ( m1 + m2 ) ) *
331 sqrt( q2 - ( m1 - m2 ) * ( m1 - m2 ) ) / 2 / sqrt( q2 );
332 boost.
set( M0 - elambda, 0, 0, +sqrt( elambda * elambda - mL * mL ) );
334 p4lambda.
set( elambda, 0, 0,
335 -sqrt( elambda * elambda - mL * mL ) );
337 p4lambda.
set( mL, 0, 0, 0 );
339 for ( j = 0; j <= 45; j++ ) {
341 p4lep1.
set( sqrt( pstar * pstar + m1 * m1 ), 0,
342 +pstar * sin( theta ), +pstar * cos( theta ) );
343 p4lep2.
set( sqrt( pstar * pstar + m2 * m2 ), 0,
344 -pstar * sin( theta ), -pstar * cos( theta ) );
347 p4lep1 =
boostTo( p4lep1, boost );
348 p4lep2 =
boostTo( p4lep2, boost );
356 <<
" - probability " << prob <<
" found at q2 = " << q2
357 <<
" (" << 100 * ( q2 - q2min ) / ( q2max - q2min )
386 static long noTries = 0;
387 static double delta = 0;
396 double r, M_L, M_Lb, M_s, M_c, M_b, q2, alpha, M_W, M_t;
397 double M_psi[2] = {0, 0}, Gamma_psi[2] = {0, 0}, k_psi[2] = {0, 0};
398 double F0_1, F0_2, a_F1, a_F2, b_F1, b_F2, F1, F2;
399 double f_1, f_2, f_3, g_1, g_2, g_3, f_1T, f_2T, f_3T, g_1T, g_2T, g_3T,
400 f_TV, f_TS, g_TV( 0.0 ), g_TS, f_T, g_T;
401 EvtComplex A1, A2, A3, B1, B2, B3, D1, D2, D3, E1, E2, E3, N1, N2, H1, H2;
402 EvtComplex C_SL, C_BR, C_LLtot, C_LRtot, C_LL, C_LR, C_RL, C_RR, C_LRLR,
403 C_RLLR, C_LRRL, C_RLRL, C_T, C_TE;
407 EvtVector4C lbar_Gmu_l[2][2], lbar_GmuG5_l[2][2], hbar_GmuPlusG5_h[2][2],
408 hbar_GmuMinusG5_h[2][2], hbar_Gmu_h[2][2];
409 EvtComplex lbar_l[2][2], lbar_G5_l[2][2], hbar_1PlusG5_h[2][2],
410 hbar_1MinusG5_h[2][2], hbar_G5_h[2][2], hbar_h[2][2];
411 EvtTensor4C lbar_Smunu_l[2][2], lbar_ESmunu_l[2][2],
412 hbar_SmunuPlusG5_h[2][2], hbar_SmunuMinusG5_h[2][2], hbar_Smunu_h[2][2];
417 M_Lb = parent->
mass();
422 alpha = 1. / 137.036;
427 if (
m_decayName ==
"Lambda_b0 -> Lambda0 e- e+" ||
428 m_decayName ==
"anti-Lambda_b0 -> anti-Lambda0 e+ e-" ) {
432 if (
m_decayName ==
"Lambda_b0 -> Lambda0 mu- mu+" ||
433 m_decayName ==
"anti-Lambda_b0 -> anti-Lambda0 mu+ mu-" ) {
437 if (
m_decayName ==
"Lambda_b0 -> Lambda0 tau- tau+" ||
438 m_decayName ==
"anti-Lambda_b0 -> anti-Lambda0 tau+ tau-" ) {
457 <<
" EvtLb2Lll already finished " << noTries
458 <<
" matrix element calculations" << std::endl;
461 r = M_L * M_L / M_Lb / M_Lb;
468 F1 = F0_1 / ( 1.0 - ( q2 / M_Lb / M_Lb ) *
469 ( a_F1 - b_F1 * ( q2 / M_Lb / M_Lb ) ) );
470 F2 = F0_2 / ( 1.0 - ( q2 / M_Lb / M_Lb ) *
471 ( a_F2 - b_F2 * ( q2 / M_Lb / M_Lb ) ) );
472 g_1 = f_1 = f_2T = g_2T = F1 + sqrt( r ) * F2;
475 g_2 = f_2 = g_3 = f_3 = g_TV = f_TV = F2 / M_Lb;
477 g_1T = f_1T = F2 / M_Lb * q2;
478 g_3T = +F2 / M_Lb * ( M_Lb + M_L );
479 f_3T = -F2 / M_Lb * ( M_Lb - M_L );
480 f_T = f_2T - f_TS * q2;
481 g_T = g_2T - g_TS * q2;
482 }
else if ( strstr(
m_FFtype.c_str(),
"HQET-delta=" ) ==
m_FFtype.c_str() ) {
485 sscanf(
m_FFtype.c_str(),
"%c%c%c%c%c%c%c%c%c%c%c%lf", &ch, &ch,
486 &ch, &ch, &ch, &ch, &ch, &ch, &ch, &ch, &ch, &delta );
487 r = M_L * M_L / M_Lb / M_Lb;
494 F1 = F0_1 / ( 1.0 - ( q2 / M_Lb / M_Lb ) *
495 ( a_F1 - b_F1 * ( q2 / M_Lb / M_Lb ) ) );
496 F2 = F0_2 / ( 1.0 - ( q2 / M_Lb / M_Lb ) *
497 ( a_F2 - b_F2 * ( q2 / M_Lb / M_Lb ) ) );
498 g_1 = f_1 = f_2T = g_2T = F1 + sqrt( r ) * F2;
501 g_2 = f_2 = g_3 = f_3 = g_TV = f_TV = F2 / M_Lb;
503 g_1T = f_1T = F2 / M_Lb * q2;
504 g_3T = +F2 / M_Lb * ( M_Lb + M_L );
505 f_3T = -F2 / M_Lb * ( M_Lb - M_L );
506 f_T = f_2T - f_TS * q2;
507 g_T = g_2T - g_TS * q2;
508 }
else if (
m_FFtype ==
"HQET-noF2" ) {
510 r = M_L * M_L / M_Lb / M_Lb;
514 F1 = F0_1 / ( 1.0 - ( q2 / M_Lb / M_Lb ) *
515 ( a_F1 - b_F1 * ( q2 / M_Lb / M_Lb ) ) );
516 g_1 = f_1 = f_2T = g_2T = F1;
517 g_2 = f_2 = g_3 = f_3 = g_TV = f_TV = 0;
522 f_T = f_2T - f_TS * q2;
523 g_T = g_2T - g_TS * q2;
525 f_1 = f_2 = f_3 = g_1 = g_2 = g_3 = f_3T = g_3T = f_TS = g_TS = f_T =
527 f_2T = f_T + f_TS * q2;
528 f_1T = ( f_TV + f_TS * ( M_L + M_Lb ) ) * q2;
529 f_1T = -q2 / ( M_Lb - M_L ) * f_3T;
530 g_2T = g_T + g_TS * q2;
531 g_1T = ( g_TV - g_TS * ( M_L - M_Lb ) ) * q2;
532 g_1T = +q2 / ( M_Lb + M_L ) * g_3T;
534 <<
" ERROR: EvtLb2Lll - unknown Form-Factors model: " <<
m_FFtype 535 <<
" - this should never happen !" << std::endl;
540 C_LL = C_LR = C_RL = C_RR = C_LRLR = C_RLLR = C_LRRL = C_RLRL = C_T =
549 (
m_WC.
hzs( M_c / M_b, q2 / M_b / M_b, M_b, M_b ) + Yld );
556 C_SL = -2 * M_s * C_7eff;
557 C_BR = -2 * M_b * C_7eff;
562 C_LL = C_LR = C_RL = C_RR = C_LRLR = C_RLLR = C_LRRL = C_RLRL = C_T =
571 (
m_WC.
hzs( M_c / M_b, q2 / M_b / M_b, M_b, M_b ) + Yld );
578 C_SL = +2 * M_s * C_7eff;
579 C_BR = +2 * M_b * C_7eff;
585 C_LL = C_LR = C_RL = C_RR = C_LRLR = C_RLLR = C_LRRL = C_RLRL = C_T =
591 EvtComplex C_7susy = -1.75 * d_u23LL - 0.25 * d_u23LR - 10.3 * d_d23LR;
593 EvtComplex C_10susy = -9.37 * d_u23LR + 1.4 * d_u23LR * d_u33RL +
602 (
m_WC.
hzs( M_c / M_b, q2 / M_b / M_b, M_b, M_b ) + Yld );
609 C_SL = -2 * M_s * C_7eff;
610 C_BR = -2 * M_b * C_7eff;
616 <<
" ERROR: EvtLb2Lll - unknown physics model: " <<
m_HEPmodel 617 <<
" - this should never happen !" << std::endl;
621 A1 = ( f_1T - g_1T ) * C_SL / q2 + ( f_1T + g_1T ) * C_BR / q2 +
622 0.5 * ( f_1 - g_1 ) * ( C_LLtot + C_LRtot ) +
623 0.5 * ( f_1 + g_1 ) * ( C_RL + C_RR );
625 A2 = ( f_2T - g_2T ) * C_SL / q2 + ( f_2T + g_2T ) * C_BR / q2 +
626 0.5 * ( f_2 - g_2 ) * ( C_LLtot + C_LRtot ) +
627 0.5 * ( f_2 + g_2 ) * ( C_RL + C_RR );
628 A3 = ( f_3T - g_3T ) * C_SL / q2 + ( f_3T + g_3T ) * C_BR / q2 +
629 0.5 * ( f_3 - g_3 ) * ( C_LLtot + C_LRtot ) +
630 0.5 * ( f_3 + g_3 ) * ( C_RL + C_RR );
632 B1 = ( f_1T + g_1T ) * C_SL / q2 + ( f_1T - g_1T ) * C_BR / q2 +
633 0.5 * ( f_1 + g_1 ) * ( C_LLtot + C_LRtot ) +
634 0.5 * ( f_1 - g_1 ) * ( C_RL + C_RR );
635 B2 = ( f_2T + g_2T ) * C_SL / q2 + ( f_2T - g_2T ) * C_BR / q2 +
636 0.5 * ( f_2 + g_2 ) * ( C_LLtot + C_LRtot ) +
637 0.5 * ( f_2 - g_2 ) * ( C_RL + C_RR );
638 B3 = ( f_3T + g_3T ) * C_SL / q2 + ( f_3T - g_3T ) * C_BR / q2 +
639 0.5 * ( f_3 + g_3 ) * ( C_LLtot + C_LRtot ) +
640 0.5 * ( f_3 - g_3 ) * ( C_RL + C_RR );
642 D1 = 0.5 * ( C_RR - C_RL ) * ( f_1 + g_1 ) +
643 0.5 * ( C_LRtot - C_LLtot ) * ( f_1 - g_1 );
644 D2 = 0.5 * ( C_RR - C_RL ) * ( f_2 + g_2 ) +
645 0.5 * ( C_LRtot - C_LLtot ) * ( f_2 - g_2 );
646 D3 = 0.5 * ( C_RR - C_RL ) * ( f_3 + g_3 ) +
647 0.5 * ( C_LRtot - C_LLtot ) * ( f_3 - g_3 );
649 E1 = 0.5 * ( C_RR - C_RL ) * ( f_1 - g_1 ) +
650 0.5 * ( C_LRtot - C_LLtot ) * ( f_1 + g_1 );
651 E2 = 0.5 * ( C_RR - C_RL ) * ( f_2 - g_2 ) +
652 0.5 * ( C_LRtot - C_LLtot ) * ( f_2 + g_2 );
653 E3 = 0.5 * ( C_RR - C_RL ) * ( f_3 - g_3 ) +
654 0.5 * ( C_LRtot - C_LLtot ) * ( f_3 + g_3 );
656 N1 = ( f_1 * ( M_Lb - M_L ) + f_3 * q2 ) / M_b *
657 ( C_LRLR + C_RLLR + C_LRRL + C_RLRL );
658 N2 = ( f_1 * ( M_Lb - M_L ) + f_3 * q2 ) / M_b *
659 ( C_LRLR + C_RLLR - C_LRRL - C_RLRL );
661 H1 = ( g_1 * ( M_Lb + M_L ) - g_3 * q2 ) / M_b *
662 ( C_LRLR - C_RLLR + C_LRRL - C_RLRL );
663 H2 = ( g_1 * ( M_Lb + M_L ) - g_3 * q2 ) / M_b *
664 ( C_LRLR - C_RLLR - C_LRRL + C_RLRL );
666 for ( i = 0; i < 4; i++ ) {
667 lbar_Gmu_l[i / 2][i % 2] =
670 lbar_GmuG5_l[i / 2][i % 2] =
673 lbar_l[i / 2][i % 2] =
676 lbar_G5_l[i / 2][i % 2] =
679 lbar_Smunu_l[i / 2][i % 2] =
682 lbar_ESmunu_l[i / 2][i % 2] =
dual(
690 parent__spParent[0] = parent->
sp( 0 );
691 parent__spParent[1] = parent->
sp( 1 );
693 parent__spParent[0] = parent->
sp( 0 );
694 parent__spParent[1] = parent->
sp( 1 );
696 parent__spParent[0] = parent->
sp( 0 );
697 parent__spParent[1] = parent->
sp( 1 );
700 <<
" ERROR: EvtLb2Lll - unknown polarization: " 706 for ( i = 0; i < 4; i++ ) {
707 hbar_GmuPlusG5_h[i / 2][i % 2] =
709 parent__spParent[i % 2] ) +
711 parent__spParent[i % 2] );
714 hbar_SmunuPlusG5_h[i / 2][i % 2] =
716 parent__spParent[i % 2] ) +
718 parent__spParent[i % 2] );
719 hbar_SmunuMinusG5_h[i / 2][i % 2] =
721 parent__spParent[i % 2] ) -
723 parent__spParent[i % 2] );
724 hbar_1PlusG5_h[i / 2][i % 2] =
726 parent__spParent[i % 2] ) +
728 parent__spParent[i % 2] );
729 hbar_1MinusG5_h[i / 2][i % 2] =
731 parent__spParent[i % 2] ) -
733 parent__spParent[i % 2] );
744 for ( i = 0; i < 4; i++ )
745 for ( j = 0; j < 4; j++ ) {
748 Matrix[j / 2][j % 2][i / 2][i % 2] +=
749 lbar_Gmu_l[i / 2][i % 2] *
750 ( A1 * hbar_GmuPlusG5_h[j / 2][j % 2] +
751 B1 * hbar_GmuMinusG5_h[j / 2][j % 2] );
756 Matrix[j / 2][j % 2][i / 2][i % 2] +=
757 lbar_Gmu_l[i / 2][i % 2] *
758 ( i1 * A2 * ( hbar_SmunuPlusG5_h[j / 2][j % 2].
cont2( q_mu ) ) +
759 B2 * ( hbar_SmunuMinusG5_h[j / 2][j % 2].cont2( q_mu ) ) );
761 Matrix[j / 2][j % 2][i / 2][i % 2] +=
762 lbar_Gmu_l[i / 2][i % 2] *
763 ( ( A3 * hbar_1PlusG5_h[j / 2][j % 2] +
764 B3 * hbar_1MinusG5_h[j / 2][j % 2] ) *
767 Matrix[j / 2][j % 2][i / 2][i % 2] +=
768 lbar_GmuG5_l[i / 2][i % 2] *
769 ( D1 * hbar_GmuPlusG5_h[j / 2][j % 2] +
770 E1 * hbar_GmuMinusG5_h[j / 2][j % 2] );
772 Matrix[j / 2][j % 2][i / 2][i % 2] +=
773 lbar_GmuG5_l[i / 2][i % 2] *
774 ( i1 * D2 * ( hbar_SmunuPlusG5_h[j / 2][j % 2].
cont2( q_mu ) ) +
775 E2 * ( hbar_SmunuMinusG5_h[j / 2][j % 2].cont2( q_mu ) ) );
777 Matrix[j / 2][j % 2][i / 2][i % 2] +=
778 lbar_GmuG5_l[i / 2][i % 2] *
779 ( ( D3 * hbar_1PlusG5_h[j / 2][j % 2] +
780 E3 * hbar_1MinusG5_h[j / 2][j % 2] ) *
783 Matrix[j / 2][j % 2][i / 2][i % 2] +=
784 lbar_l[i / 2][i % 2] *
785 ( N1 * hbar_h[j / 2][j % 2] + H1 * hbar_G5_h[j / 2][j % 2] );
787 Matrix[j / 2][j % 2][i / 2][i % 2] +=
788 lbar_G5_l[i / 2][i % 2] *
789 ( N2 * hbar_h[j / 2][j % 2] + H2 * hbar_G5_h[j / 2][j % 2] );
791 Matrix[j / 2][j % 2][i / 2][i % 2] +=
792 cont( lbar_Smunu_l[i / 2][i % 2],
793 4 * C_T * f_T * hbar_Smunu_h[j / 2][j % 2] );
795 Matrix[j / 2][j % 2][i / 2][i % 2] +=
cont(
796 lbar_Smunu_l[i / 2][i % 2],
797 -4 * C_T * f_TV * i1 *
802 Matrix[j / 2][j % 2][i / 2][i % 2] +=
803 cont( lbar_Smunu_l[i / 2][i % 2],
804 -4 * C_T * f_TS * i1 *
807 hbar_h[j / 2][j % 2] );
809 Matrix[j / 2][j % 2][i / 2][i % 2] +=
810 cont( lbar_ESmunu_l[i / 2][i % 2],
811 4 * C_TE * f_T * i1 * hbar_Smunu_h[j / 2][j % 2] );
813 Matrix[j / 2][j % 2][i / 2][i % 2] +=
cont(
814 lbar_ESmunu_l[i / 2][i % 2],
820 Matrix[j / 2][j % 2][i / 2][i % 2] +=
821 cont( lbar_ESmunu_l[i / 2][i % 2],
825 hbar_h[j / 2][j % 2] );
834 amp->
vertex( spins, Matrix[j / 2][j % 2][i / 2][i % 2] );
884 temp.
set( 0, 1, i2 * ( d * ( mat01 * dp ) ) );
885 temp.
set( 1, 0, -temp.
get( 0, 1 ) );
887 temp.
set( 0, 2, i2 * ( d * ( mat02 * dp ) ) );
888 temp.
set( 2, 0, -temp.
get( 0, 2 ) );
890 temp.
set( 0, 3, i2 * ( d * ( mat03 * dp ) ) );
891 temp.
set( 3, 0, -temp.
get( 0, 3 ) );
893 temp.
set( 1, 2, i2 * ( d * ( mat12 * dp ) ) );
894 temp.
set( 2, 1, -temp.
get( 1, 2 ) );
896 temp.
set( 1, 3, i2 * ( d * ( mat13 * dp ) ) );
897 temp.
set( 3, 1, -temp.
get( 1, 3 ) );
899 temp.
set( 2, 3, i2 * ( d * ( mat23 * dp ) ) );
900 temp.
set( 3, 2, -temp.
get( 2, 3 ) );
EvtTensor3C directProd(const EvtVector3C &c1, const EvtVector3C &c2)
EvtComplex Yld(double q2, double ki[], double Gi[], double Mi[], int ni, EvtComplex c1, EvtComplex c2, EvtComplex c3, EvtComplex c4, EvtComplex c5, EvtComplex c6, double ialpha)
void set(int i, int j, const EvtComplex &c)
static std::string name(EvtId i)
virtual EvtDiracSpinor sp(int) const
EvtTensor4C EvtLeptonTCurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtComplex hzs(double z, double shat, double mu, double M_b)
double getArg(unsigned int j)
std::string getArgStr(int j) const
static const EvtGammaMatrix & g0()
EvtComplex EvtLeptonPCurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtTensor4C dual(const EvtTensor4C &t2)
std::ostream & EvtGenReport(EvtGenSeverity severity, const char *facility=0)
double GetStrongCouplingConst()
EvtDecayBase * clone() override
EvtComplex cont(const EvtTensor4C &t1, const EvtTensor4C &t2)
std::string m_polarizationIntroduction
EvtVector4C cont2(const EvtVector4C &v4) const
EvtWilsonCoefficients m_WC
void decay(EvtParticle *p) override
int getSpinStates() const
void makeDaughters(unsigned int ndaug, EvtId *id)
std::string getName() override
void set(int i, double d)
EvtRaritaSchwinger boostTo(const EvtRaritaSchwinger &rs, const EvtVector4R p4)
void init(EvtId part_n, const EvtVector4R &p4) override
EvtComplex C9efftilda(double z, double shat, double alpha_S, EvtComplex c1, EvtComplex c2, EvtComplex c3, EvtComplex c4, EvtComplex c5, EvtComplex c6, EvtComplex c9tilda, int ksi)
EvtSpinDensity getSpinDensity()
void CalculateAllCoefficients()
virtual EvtDiracSpinor spParent(int) const
void setProbMax(double prbmx)
EvtVector4C EvtLeptonACurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
void vertex(const EvtComplex &)
void calcAmp(EvtAmp *amp, EvtParticle *parent)
EvtId getParentId() const
static const EvtGammaMatrix & g1()
double initializePhaseSpace(unsigned int numdaughter, EvtId *daughters, bool forceResetMasses=false, double poleSize=-1., int whichTwo1=0, int whichTwo2=1)
static EvtId getId(const std::string &name)
const EvtVector4R & getP4() const
static const EvtGammaMatrix & g2()
double lambda(double q, double m1, double m2)
static const EvtGammaMatrix & g3()
EvtComplex exp(const EvtComplex &c)
EvtTensor4C EvtLeptonTG5Current(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
std::string m_effectContribution
double m_polarizationLambdab0
EvtVector4C EvtLeptonVACurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
void setDiagonalSpinDensity()
EvtParticle * getDaug(int i)
EvtVector4C EvtLeptonVCurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
EvtComplex C7b2sg(double alpha_S, double et, EvtComplex c2, double M_t, double M_W)
static double getMass(EvtId i)
EvtComplex EvtLeptonSCurrent(const EvtDiracSpinor &d, const EvtDiracSpinor &dp)
void initProbMax() override
const EvtComplex & get(int i, int j) const
EvtVector4R getP4Restframe() const
void init(EvtId p, int ndaug, EvtId *daug)
static const EvtGammaMatrix & g5()
double normalizedProb(const EvtSpinDensity &d)
EvtId getDaug(int i) const