MadGraph prediction for neutrino trident production
Hi! I am trying to calculate the cross-section of a process
nu_mu + gamma -> mu+ + mu- + nu_mu
in the limit s_nu gamma >> m_mu^2. The cross-section of this process under the limit has been derived analytically in https:/
CV = 1/2 + sin^2(theta_W)*2, CA = 1/2
I have also derived it independently.
I assume that e.g. at s > 10 GeV^2 the approximation s >> 4*m_mu^2 works well (as the ratio is >200), and it makes sense to compare the predictions of MG with the analytic formula. To generate the cross-section in MG, I used (the loop_sm model is necessary since m_mu != 0.)
import model loop_sm
generate vm a > mu+ mu- vm
Then, I changed the run_card.dat accordingly (please see the comment below), where I assumed equal energies of nu and gamma, such that s_nu gamma = (Enu+Egamma)^2.
The values of the cross-section for various energies predicted by MG is shown in this screenshot (https:/
However, the problem is that there is a difference at any s from 10 GeV^2 to 300 GeV^2, please see the plot (https:/
Could you please tell me whether I made something wrong with the MG part? If not, then what may be a reason for the discrepancy?
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Maksym Ovchynnikov for more information if necessary.