Name: Francesco Date: 07/02/16-12:40:40 AM Z


Hi everyone.
I am trying to compute amplitudes of processes involving operators of dimension-6 (with four fermion fields) of the form QLbar Gamma[mu] T[a] QL QLbar Gamma[mu] T[a] QL.

I work with the SM.fr model of FeynRules (2.3.22) where I added a 4-fermion operator.

L81qq1:=C81qq /Lambda^2 Module[{ii,jj,kk,kk1,jj1,mu,aa,cc,sp1,sp2,sp3,sp4,cc1,cc2,cc3}, (QLbar[sp1,1,jj,cc].QL[sp2,1,jj,cc2] Ga[mu,sp1,sp2]) T[aa,cc,cc2]Taa,cc3,cc4];

The generation of the Model and GenericModel files goes well and I use them with the Dev version of FeynCalc 9.0.1 (and FeynArts 3.9).

You can download my model from here https://www.dropbox.com/sh/hm77oyedcb3p055/AABzBlEPZvcPjAstjoclCR3ya?dl=0

I read this <1042.html> and then I did

$LoadFeynArts = True;
«FeynCalc`;
$FAVerbose = 0;
FAPatch[PatchModelsOnly -> True]

top = CreateTopologies[0, 2 -> 2];
diags = InsertFields[
   top, {F[3, {1}], F[3, {1}]} -> {F[3, {1}], F[3, {1}]},
   InsertionLevel -> {Particles}, Model -> “Standard_Model_dijets_FA”,
    GenericModel -> “Standard_Model_dijets_FA”,
   ExcludeParticles -> {V[1], V[2], V[3]}];
Paint[%]

AmpList =
  FCFAConvert[CreateFeynAmp[diags], ChangeDimension -> 4,
    IncomingMomenta -> {p1, p2}, OutgoingMomenta -> {k1, k2},
    DropSumOver -> True, List -> False, UndoChiralSplittings -> True,
    FinalSubstitutions ->
     Flatten[Join[M$FACouplings, {SUNFDelta[__] :> 1}]]] // Contract;

Cmp = AmpList // ComplexConjugate // FCRenameDummyIndices;
sqAmp = AmpList Cmp // Simplify;
sqAmp // FermionSpinSum[#, ExtraFactor -> 1] &

and at this point I get the error

Something went wrong while computing the fermions spin sum! Returning \
unevaluated expression:

Specially, If I see the expression of the amplitude I get

AmpList
-I ((1/(3 Lambda^2))
   2 I C81qq IndexSum[
     DiracObject[DiracGamma[a], DiracGamma[7]][SequenceForm[“Dir”, 3],
        SequenceForm[
       ”Dir”, 2]] DiracObject[DiracGamma[a], DiracGamma[7]][
       SequenceForm[“Dir”, 4], SequenceForm[“Dir”, 1]], {a, 1, 4}] + (
   1/(3 Lambda^2))
   2 I C81qq IndexSum[
     DiracObject[DiracGamma[a], DiracGamma[7]][SequenceForm[“Dir”, 3],
        SequenceForm[
       ”Dir”, 1]] DiracObject[DiracGamma[a], DiracGamma[7]][
       SequenceForm[“Dir”, 4], SequenceForm[“Dir”, 2]], {a, 1,
      4}]) DiracObject[
   Spinor[Momentum[k1], FCGV[”]” “FCGV[” FCGV[”]” “FCGV[” MU]], 1]][
SequenceForm[
  ”Dir”, 3]] DiracObject[
   Spinor[Momentum[k2], FCGV[”]” “FCGV[” FCGV[”]” “FCGV[” MU]], 1]][
SequenceForm[
  ”Dir”, 4]] DiracObject[
   Spinor[-Momentum[p1], FCGV[”]” “FCGV[” FCGV[”]” “FCGV[” MU]], 1]][
SequenceForm[
  ”Dir”, 1]] DiracObject[
   Spinor[-Momentum[p2], FCGV[”]” “FCGV[” FCGV[”]” “FCGV[” MU]], 1]][
SequenceForm[“Dir”, 2]]

and there is a DiracObject function which is not defined.

Where is the problem?

Regards,
Francesco