Higgs production and pDF
Hi all,
I try to calculate the Higgs production cross section by Mathematica using
for example the formula in [arXiv:hepph/0503172].
The problem comes when I calculate " σ_0 " it's of order 10^10 pb at mh = 125 GeV, then when I convolute with the parton distribution function to get the pp > H production according to equ.( 3.58) it gives me of order 10^7 pb !! while in the plot (3.18) pp > H is around 50 pb.
Here is my calculation;
for σ0:
σ0[mh_] : = ( ( c^2)/ (256 *Pi)) *(1/2)^2* ( ((Yht)/(mT))* At [mh])^2
where
c = 0.117, Yht= mt/v, v= 174 and At is the fermion loop function.
For σ(pp > h) I used:
LogPlot[ {σ0[mh]*(q)* NIntegrate[(1/x)* pdf[iset, iparton, x, 14000] * pdf[iset, iparton, (q/x), 14000], {x, q, 1}], {mh, 100, 200}]
With
iset = 1; (* CTEQ5M*)
iparton = 0;
q := mh^2/14000^2;
I don't know what is the problem I have, any help please,
Thanks,
Safinaz
