Differential D98 Diff 424 src/EvtGenModels/EvtbsToLLLLAmp.cpp

Changeset View

Standalone View

Show First 20 Lines • Show All 481 Lines • ▼ Show 20 Lines	void EvtbsToLLLLAmp::CalcAmp( EvtParticle* parent, EvtAmp& amp,
g_b2qsecond = c10a * Fa;		g_b2qsecond = c10a * Fa;
g_barb2barqsecond = g_b2qsecond;		g_barb2barqsecond = g_b2qsecond;

EvtTensor4C T1, T2; // Tensor structures for		EvtTensor4C T1, T2; // Tensor structures for
EvtTensor4C T1second,		EvtTensor4C T1second,
T2second; // the hadronic matrix element calculation		T2second; // the hadronic matrix element calculation

// B - and barB - mesons descriptors		// B - and barB - mesons descriptors
static EvtIdSet bmesons( "anti-B0", "anti-B_s0" );		static EvtIdSet bmesons{ "anti-B0", "anti-B_s0" };
static EvtIdSet bbarmesons( "B0", "B_s0" );		static EvtIdSet bbarmesons{ "B0", "B_s0" };

EvtId parentID = parent->getId();		EvtId parentID = parent->getId();

if ( bmesons.contains( parentID ) ) {		if ( bmesons.contains( parentID ) ) {
// The amplitude for the decay barB -> gamma ell^+ ell^- or		// The amplitude for the decay barB -> gamma ell^+ ell^- or
// b \bar q -> gamma ell^+ ell^-		// b \bar q -> gamma ell^+ ell^-

T1 = a_b2q * unit1 * dual( EvtGenFunctions::directProd( hatq, hatk ) ) -		T1 = a_b2q * unit1 * dual( EvtGenFunctions::directProd( hatq, hatk ) ) -
▲ Show 20 Lines • Show All 323 Lines • Show Last 20 Lines