Problem of convergence and cross section value for Bhabha scattering

Asked by Yehia Abdelaziz

I am trying to generate the process ee > ee but I get a very large cross section of O(10^20) fb. I don't think this is normal since other generators give much less values. This is my sindarin file:

sqrts= 500 GeV

?vis_diags=true
process DY = e1, E1 => e1, E1

beams = "e1" , "E1" => circe2 => isr

$circe2_file = "ilc500.circe"
$circe2_design = "ILC"
?circe2_polarized = false

beams_pol_density= @(-1), @(+1)
beams_pol_fraction = 80%, 30%

?ps_fsr_active = true
?hadronization_active = true
$shower_method = "PYTHIA6"
$hadronization_method= "PYTHIA6"

integrate(DY) {iterations = 10:40000, 6:80000}

$sample = "Drell_Yann_ee_1000000_events_ilc_sqrts_500GeV"
sample_format = lhef

simulate(DY) {n_events = 1000000}
show(results)

Question information

Language:
English Edit question
Status:
Solved
For:
WHIZARD Edit question
Assignee:
Juergen Reuter Edit question
Solved by:
Yehia Abdelaziz
Solved:
Last query:
Last reply:
Revision history for this message
Juergen Reuter (j.r.reuter) said :
#1

Juergen Reuter suggests this article as an answer to your question:
FAQ #3404: “Why does Whizard give a much larger cross section than other generators?”.

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#2

Dear Yehia,
this is actually a recurring question for which we created a FAQ now (cf. the link). Whizard doesn't set any default cuts, so you need to set a minimal angular cut here or a cut on the momentum transfer for the t-channel contribution. Other generators normally set a default minimal pT for final-state leptons which would affect Bhabha scattering as in your example as well.
Cheers,
    JRR (Juergen)

Revision history for this message
Yehia Abdelaziz (yehia95) said :
#3

Thank you Juergen.

I think the cuts are not the main reason for this since the difference orders of magnitude higher. In madgrap I get a cross section of O(10^5) fb while in whizard I get O(10^16) fb with an error of O(10^16)!.Even after I applied cuts on pt and eta in the sindarin file.

Revision history for this message
Juergen Reuter (j.r.reuter) said :
#4

Using this cut
cuts = all M < -10 GeV [incoming e1:E1,e1:E1]
on the momentum transfer between incoming and outgoing electrons/positrons, and then I get (with beam spectrum and ISR)
a perfectly converging integration of the right order of magnitude:
| VAMP: parameter mismatch, discarding grid file 'DY.m1.vg'
   1 1000 2.7341512E+06 3.62E+05 13.25 4.19* 4.11
   2 998 2.6873194E+06 1.23E+05 4.56 1.44* 11.44
   3 996 2.8615803E+06 5.99E+04 2.09 0.66* 25.26
|-----------------------------------------------------------------------------|
   3 2994 2.8259557E+06 5.33E+04 1.88 1.03 25.26 0.85 3
|-----------------------------------------------------------------------------|
   4 9997 2.7911859E+06 1.86E+04 0.67 0.67 20.97
   5 9997 2.7432835E+06 1.84E+04 0.67 0.67 20.37
   6 9997 2.7925499E+06 1.91E+04 0.68 0.68 13.08
|-----------------------------------------------------------------------------|
   6 29991 2.7750988E+06 1.08E+04 0.39 0.67 13.08 2.29 3
|=============================================================================

Revision history for this message
Yehia Abdelaziz (yehia95) said :
#5

Thank you this solved my problem.