**Next message:**D. Azevedo: "Re: Problem contracting Lorentz indexes"**Previous message:**D. Azevedo: "Problem contracting Lorentz indexes"**In reply to:**D. Azevedo: "Problem contracting Lorentz indexes"**Next in thread:**D. Azevedo: "Re: Problem contracting Lorentz indexes"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

Hi,

well if you have an expression like

SpinorUBar[p1].GA[mu].SpinorV[p2] SpinorVBar[p2].GA[mu].SpinorU[p1]

then there is nothing more you can contract at this stage, so the

Lorentz index will of course remain there.

If you are calculating an amplitude squared, the next step is to apply

the usual summation formulas for Dirac spinors

u ubar = pslash +m /2m,

v vbar = pslash -m /2m,

take the Dirac traces and contract the remaining indices (now that the

expression is free of Dirac spinors and matrices)

res=((M12 // FermionSpinSum) /. DiracTrace -> Tr) // Contract

FreeQ[res, LorentzIndex]

(*True*)

Cheers,

Vladyslav

Am 28.05.2018 um 22:14 schrieb D. Azevedo:

*> Hello FeynCalc team,
*

*>
*

*> first of all, thank you for such a powerful tool.
*

*>
*

*> I am beginning learning FeynCalc and before venturing into loops, I would like to compute the squared matrix element of e+ e- -> t bar{t} h at tree level. I am having some problems related to contracting lorentz indexes. My nb file is:
*

*>
*

*> SetOptions[InitializeModel,
*

*> ModelEdit :> (M$ClassesDescription =
*

*> M$ClassesDescription /. MZ -> MZc)]
*

*> InitializeModel[{SM, UnitarySM},
*

*> GenericModel -> {Lorentz, UnitaryLorentz}];
*

*>
*

*> NoElectronHCoupling =
*

*> ExcludeFieldPoints -> {FieldPoint[0][-F[2, {1}], F[1, {1}], S[3]],
*

*> FieldPoint[0][-F[2, {1}], F[2, {1}], S[1]],
*

*> FieldPoint[0][-F[2, {1}], F[2, {1}], S[2]]}
*

*>
*

*> part = InsertFields[
*

*> CreateTopologies[0,
*

*> 2 -> 3], {F[2, {1}], -F[2, {1}]} -> {F[3, {3}], -F[3, {3}], S[1]},
*

*> Restrictions -> NoElectronHCoupling, InsertionLevel -> {Classes},
*

*> Model -> {SM, UnitarySM},
*

*> GenericModel -> {Lorentz, UnitaryLorentz}];
*

*> Paint[part, PaintLevel -> {Classes}]
*

*>
*

*> listdiag =
*

*> FCFAConvert[CreateFeynAmp[part], SMP -> True, ChangeDimension -> 4,
*

*> IncomingMomenta -> {p1, p2}, OutgoingMomenta -> {k1, k2, k3},
*

*> DropSumOver -> True, List -> False, UndoChiralSplittings -> True];
*

*> diag1 = listdiag[[5]]
*

*> diag2 = listdiag[[4]]
*

*> diag3 = listdiag[[3]]
*

*> diag4 = listdiag[[2]]
*

*> diag5 = listdiag[[1]]
*

*>
*

*> I am splitting the diagrams just to follow the same idea as some tutorials I have been reading. diag[[5]] is the Higgsstrahung from the Z boson. In this diagram (as with all of the others) I haven't been able to contract the lorentz indexes, there is one that remains for some reason.
*

*>
*

*> (I set up my kinematics)
*

*> ClearScalarProducts;
*

*>
*

*> SP[p1, p1] = me^2;
*

*> SP[p2, p2] = me^2;
*

*> SP[p1, p2] = (s - 2 me^2)/2 ;
*

*>
*

*> SP[k1, k1] = mt^2;
*

*> SP[k2, k2] = mt^2;
*

*> SP[k3, k3] = mh^2;
*

*> SP[k1, k2] = (s3 - 2 mt^2)/2;
*

*> SP[k1, k3] = (s2 - mt^2 - mh^2)/2;
*

*> SP[k2, k3] = (s1 - mt^2 - mh^2)/2;
*

*>
*

*> diag1 = diag1 /. {Lor1 -> \[Mu], Lor2 -> \[Nu], Lor3 -> \[Alpha],
*

*> Lor4 -> \[Beta]};
*

*> diag1C = ComplexConjugate[
*

*> diag1] /. {\[Mu] -> \[Mu]lin, \[Nu] -> \[Nu]lin, \[Alpha] -> \
*

*> \[Alpha]lin, \[Beta] -> \[Beta]lin};
*

*> M12 = diag1 diag1C // Contract
*

*>
*

*> M12 still have mu and mulin indexes left to be contracted. What am I doing wrong?
*

*>
*

*> Thank you,
*

*> Duarte
*

*>
*

**Next message:**D. Azevedo: "Re: Problem contracting Lorentz indexes"**Previous message:**D. Azevedo: "Problem contracting Lorentz indexes"**In reply to:**D. Azevedo: "Problem contracting Lorentz indexes"**Next in thread:**D. Azevedo: "Re: Problem contracting Lorentz indexes"**Messages sorted by:**[ date ] [ thread ] [ subject ] [ author ]**Mail actions:**[ respond to this message ] [ mail a new topic ]

*
This archive was generated by hypermail 2b29
: 11/20/18-10:40:01 PM Z CET
*