Problem About Calculation of double Higgs production Cross Section (for gluon fusion) with FeynArts+FormCalc

+ 2 like - 0 dislike
144 views

I want to mention about my problem about FeynArts/FormCalc.I have tried to calculate LHC cross section for "gluon gluon -> H H" process at Standard Model. However, my result XS is quite different from the ones in the literature. I have no idea about this difference or where is the problem. I'm using FeynArts 3.9, FormCalc 9.3 and LoopTools 2.13. I attached the my FeynArts mathematica file.

In literature

XS(gluon gluon -> H H at 14 TeV at LHC)=~0,033 pb

but it is ~0,95 pb in my result.

Could you give me some suggestions? Please, check my mathematica file. Is there any mistakes or missing in my codes.

Thank you very much

<< FeynArts.m

<< FormCalc.m

InitialDirectory;name="glglHHSM"process=V[5],V[5]−>S[1],>S[1];SetOptions[InsertFields,Model−>"SMQCD",InsertionLevel−>>Particles,ExcludeFieldPoints−>FieldPoint[0][V[5],V[5],>V[5]]];