diff --git a/FeynCalc/Feynman/FCFADiracChainJoin.m b/FeynCalc/Feynman/FCFADiracChainJoin.m index 0dd548f64..4541f58f3 100644 --- a/FeynCalc/Feynman/FCFADiracChainJoin.m +++ b/FeynCalc/Feynman/FCFADiracChainJoin.m @@ -197,7 +197,7 @@ (* u_i v_j A_ij -> ubar.A.v or vbar.A.u *) diracChainEvalS[rest_. DiracChain[S: Spinor[_. m1_Momentum, ___], a_DiracIndex] DiracChain[z_, a_DiracIndex, b_DiracIndex] DiracChain[Spinor[s_. m2_Momentum, r___], b_DiracIndex]]:= - diracChainEvalS[rest ordering[First[m2],First[m1]]] FCCCT[DOT[S, z, Spinor[-s m2, r]], Explicit->True, FCDiracIsolate->True, FCI->True]/; MemberQ[optFirst, Spinor[s m2, r]]; + diracChainEvalS[rest ordering[First[m2],First[m1]]] SpinorChainTranspose[DOT[S, z, Spinor[-s m2, r]], FCI->True]/; MemberQ[optFirst, Spinor[s m2, r]]; diracChainEvalS[rest_. ordering[a__] ordering[b__]]:= diracChainEvalS[rest ordering[a,b]]; diff --git a/Tests/Feynman/FCFADiracChainJoin.test b/Tests/Feynman/FCFADiracChainJoin.test index 6bc7f9993..a5c1954b4 100644 --- a/Tests/Feynman/FCFADiracChainJoin.test +++ b/Tests/Feynman/FCFADiracChainJoin.test @@ -191,7 +191,25 @@ Col6] SUNTF[{Glu5}, colInB, Col7]"}, DiracIndex[Dir6]]}, FCE -> True]", "{Spinor[-Momentum[InMom2], mqu, 1] . GA[Lor1] . GA[7] . (mqu - GS[InMom2 - OutMom3]) . GS[-OutMom1 - OutMom2] . -GA[Lor2] . GA[7] . Spinor[Momentum[InMom1], mqu, 1]}"} +GA[Lor2] . GA[7] . Spinor[Momentum[InMom1], mqu, 1]}"}, +{"fcstFCFADiracChainJoin-ID7","FCFADiracChainJoin[ +DiracChain[Spinor[-Momentum[p1, D], mqu, 1], + DiracIndex[Dir1]] DiracChain[Spinor[-Momentum[p2, D], mqu, 1], + DiracIndex[Dir2]] DiracChain[DiracGamma[LorentzIndex[Lor2, D], D], + DiracIndex[Dir5], DiracIndex[Dir1]] DiracChain[ + DiracGamma[LorentzIndex[Lor3, D], D], DiracIndex[Dir2], + DiracIndex[Dir7]] DiracChain[ + mqu + DiracGamma[Momentum[-l - p2, D], D], DiracIndex[Dir7], + DiracIndex[Dir8]] DiracChain[ + mqu + DiracGamma[Momentum[-k - k1 - k2 + l + p2, D], D], + DiracIndex[Dir5], DiracIndex[Dir6]] DiracChain[ + DiracGamma[LorentzIndex[LCdummy2, D], D] . DiracGamma[6], + DiracIndex[Dir8], DiracIndex[Dir6]] // FCE, +First -> {Spinor[-Momentum[p2, D], mqu, 1]}, FCE -> True]", +"-Spinor[-Momentum[p2, D], mqu, 1] . +GAD[Lor3] . (mqu + GSD[-l - p2]) . GAD[LCdummy2] . +GA[6] . (mqu - GSD[-k - k1 - k2 + l + p2]) . GAD[Lor2] . +Spinor[Momentum[p1, D], mqu, 1]"} };