Aug 3rd 2013, 01:04 PM
|
#2 |
Forum Admin
Join Date: Apr 2008 Location: On the dance floor, baby!
Posts: 2,505
|
Originally Posted by SRS Hi all,
I try to calculate the Higgs production cross section by Mathematica using
for example the formula in [arXiv:hep-ph/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 |
I can't seem to get arXiv to reference that paper itself...it keeps sending me to the Cornell U. library and it won't let me view the article.
Question: The 50 pb...is this a measurement or is it a theoretical calculation? I would like to see how you did the calculation if so.
-Dan
__________________
Do not meddle in the affairs of dragons for you are crunchy and taste good with ketchup.
See the forum rules here.
|
| |