Changeset View

Standalone View

Show First 20 Lines • Show All 150 Lines • ▼ Show 20 Lines	if ( fabs( dSqDenom ) < 1e-10 ) {
validAmp = false;		validAmp = false;
}		}

double factor( 1.0 );		double factor( 1.0 );
if ( validAmp ) {		if ( validAmp ) {
factor = dSq / ( dSqDenom * kSq );		factor = dSq / ( dSqDenom * kSq );
}		}

int iPols[4] = {0, 0, 0, 0};		int iPols[4] = { 0, 0, 0, 0 };

// Calculate the amplitude terms, looping over the chi, psi and lepton states		// Calculate the amplitude terms, looping over the chi, psi and lepton states
for ( int iChi = 0; iChi < 3; iChi++ ) {		for ( int iChi = 0; iChi < 3; iChi++ ) {
iPols[0] = iChi;		iPols[0] = iChi;
EvtVector4C epsChi = root->epsParent( iChi );		EvtVector4C epsChi = root->epsParent( iChi );

for ( int iPsi = 0; iPsi < 3; iPsi++ ) {		for ( int iPsi = 0; iPsi < 3; iPsi++ ) {
iPols[1] = iPsi;		iPols[1] = iPsi;
Show All 29 Lines