Page MenuHomeHEPForge
Differential D90 Diff 388 src/EvtGenModels/EvtLb2Lll.cpp

Changeset View

Standalone View

View Options

src/EvtGenModels/EvtLb2Lll.cpp

Show First 20 Lines • Show All 286 Lines • ▼ Show 20 Linesvoid EvtLb2Lll::initProbMax()
if ( m_maxProbability < 1e-10 ) { if ( m_maxProbability < 1e-10 ) {
EvtDiracParticle* parent = new EvtDiracParticle; EvtDiracParticle* parent = new EvtDiracParticle;
parent->noLifeTime(); parent->noLifeTime();
parent->init( getParentId(), parent->init( getParentId(),
EvtVector4R( EvtPDL::getMass( getParentId() ), 0, 0, 0 ) ); EvtVector4R( EvtPDL::getMass( getParentId() ), 0, 0, 0 ) );
parent->setDiagonalSpinDensity(); parent->setDiagonalSpinDensity();
EvtAmp amp; EvtAmp amp;
EvtId daughters[3] = {getDaug( 0 ), getDaug( 1 ), getDaug( 2 )}; EvtId daughters[3] = { getDaug( 0 ), getDaug( 1 ), getDaug( 2 ) };
amp.init( getParentId(), 3, daughters ); amp.init( getParentId(), 3, daughters );
parent->makeDaughters( 3, daughters ); parent->makeDaughters( 3, daughters );
EvtParticle* lambda = parent->getDaug( 0 ); EvtParticle* lambda = parent->getDaug( 0 );
EvtParticle* lep1 = parent->getDaug( 1 ); EvtParticle* lep1 = parent->getDaug( 1 );
EvtParticle* lep2 = parent->getDaug( 2 ); EvtParticle* lep2 = parent->getDaug( 2 );
lambda->noLifeTime(); lambda->noLifeTime();
lep1->noLifeTime(); lep1->noLifeTime();
lep2->noLifeTime(); lep2->noLifeTime();
▲ Show 20 Lines • Show All 59 Lines • ▼ Show 20 Lines if ( m_maxProbability < 1e-10 ) {
//::abort(); //::abort();
} }
//m_poleSize = 0.04*q2min; //m_poleSize = 0.04*q2min;
m_maxProbability *= 1.2; m_maxProbability *= 1.2;
delete parent; delete parent;
} }
if (m_maxProbability < 1e-10) { if ( m_maxProbability < 1e-10 ) {
EvtGenReport( EVTGEN_INFO, "EvtLb2Lll") EvtGenReport( EVTGEN_INFO, "EvtLb2Lll" )
<< "EvtLb2Lll found zero max prob, setting to 1.5e6" << std::endl; << "EvtLb2Lll found zero max prob, setting to 1.5e6" << std::endl;
m_maxProbability = 1.5e6; m_maxProbability = 1.5e6;
} }
setProbMax( m_maxProbability ); setProbMax( m_maxProbability );
EvtGenReport( EVTGEN_INFO, "EvtGen" ) EvtGenReport( EVTGEN_INFO, "EvtGen" )
<< " EvtLb2Lll set up maximum probability to " << m_maxProbability << " EvtLb2Lll set up maximum probability to " << m_maxProbability
<< std::endl; << std::endl;
} }
Show All 12 Linesvoid EvtLb2Lll::calcAmp( EvtAmp* amp, EvtParticle* parent )
EvtComplex Matrix[2][2][2][2]; EvtComplex Matrix[2][2][2][2];
EvtComplex i1( 0, 1 ); EvtComplex i1( 0, 1 );
int i, j, spins[4]; int i, j, spins[4];
char ch; char ch;
double r, M_L, M_Lb, M_s, M_c, M_b, q2, alpha, M_W, M_t; double r, M_L, M_Lb, M_s, M_c, M_b, q2, alpha, M_W, M_t;
double M_psi[2] = {0, 0}, Gamma_psi[2] = {0, 0}, k_psi[2] = {0, 0}; double M_psi[2] = { 0, 0 }, Gamma_psi[2] = { 0, 0 }, k_psi[2] = { 0, 0 };
double F0_1, F0_2, a_F1, a_F2, b_F1, b_F2, F1, F2; double F0_1, F0_2, a_F1, a_F2, b_F1, b_F2, F1, F2;
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, 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,
f_TV, f_TS, g_TV( 0.0 ), g_TS, f_T, g_T; f_TV, f_TS, g_TV( 0.0 ), g_TS, f_T, g_T;
EvtComplex A1, A2, A3, B1, B2, B3, D1, D2, D3, E1, E2, E3, N1, N2, H1, H2; EvtComplex A1, A2, A3, B1, B2, B3, D1, D2, D3, E1, E2, E3, N1, N2, H1, H2;
EvtComplex C_SL, C_BR, C_LLtot, C_LRtot, C_LL, C_LR, C_RL, C_RR, C_LRLR, EvtComplex C_SL, C_BR, C_LLtot, C_LRtot, C_LL, C_LR, C_RL, C_RR, C_LRLR,
C_RLLR, C_LRRL, C_RLRL, C_T, C_TE; C_RLLR, C_LRRL, C_RLRL, C_T, C_TE;
EvtComplex Yld, C_7eff, C_9eff; EvtComplex Yld, C_7eff, C_9eff;
EvtComplex V_ts, V_tb; EvtComplex V_ts, V_tb;
▲ Show 20 Lines • Show All 498 Lines • Show Last 20 Lines