MadGraph5_aMC@NLO

gg into 4 leptons

Asked by Giacomo Magni

Hi,
I'm trying to simulate that process:

     import model SMEFTatNLO_U2_2_U3_3_cG_4F_LO_UFO-NLO
     generate g g > mu- vm~ e+ ve [noborn=QCD] QED=4 NP=0 @1
     add process g g > mu+ vm e- ve~ [noborn=QCD] QED=4 NP=0 @2
     output myfolder
     launch myfolder

Everything goes fine, except for this warning:

 ^[[1;34mWARNING: Channel G1//data/theorie/gmagni/ATLAS_WW/gg_WW_SMlimit_massless_NLO_ATLAScut/SubProcesses/P2_gg_lvllvl has a very low efficiency of unweighting. Might not be possible to reach target ^[[0m

The simulation goes on but then I got this error:

INFO: Refine results to 100000
INFO: Generating 100000.0 unweighted events.
INFO: Effective Luminosity 1491346.36418 pb^-1
INFO: need to improve 2 channels
Current estimate of cross-section: 0.077110859531 +- 0.000209939353442
INFO: Idle: 0, Running: 2, Completed: 0 [ current time: 21h21 ]
INFO: Idle: 1, Running: 1, Completed: 1 [ 3s ]
^[[1;34mWARNING: float division by zero ^[[0m
^[[1;34mWARNING: float division by zero ^[[0m
INFO: Idle: 0, Running: 0, Completed: 4 [ 3.1s ]
float division by zero
Cluster Error detected. Do you want to clean the queue? ("c"=continue the run anyway) [^[[4my^[[0m, n, c]
^[[1mThe answer to the previous question is not set in your input file^[[0m
^[[1mUse y value^[[0m

The log file reports:

Traceback (most recent call last):
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/extended_cmd.py", line 1515, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/common_run_interface.py", line 6963, in new_fct
    original_fct(obj, *args, **opts)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/madevent_interface.py", line 2529, in run_generate_events
    self.exec_cmd('refine %s' % nb_event, postcmd=False)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/extended_cmd.py", line 1544, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/extended_cmd.py", line 1464, in onecmd_orig
    return func(arg, **opt)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/madevent_interface.py", line 3508, in do_refine
    html=True)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/interface/madevent_interface.py", line 5552, in monitor
    self.cluster.wait(self.me_dir, update_status, update_first=update_first)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/various/cluster.py", line 2807, in worker
    returncode = exe(*arg, **opt)
  File "/data/theorie/gmagni/ATLAS_WW/tmp/madgraph/madevent/gen_ximprove.py", line 1465, in combine_iteration
    n_target_combined = (needed_event-nunwgt) / efficiency
ZeroDivisionError: float division by zero

I'm running on multicore mode on a desktop.
What can I do to fix this?

Thank you for your help,
Giacomo

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

You should have more precise information in
/data/theorie/gmagni/ATLAS_WW/gg_WW_SMlimit_massless_NLO_ATLAScut/SubProcesses/P2_gg_lvllvl/G1
directory.

Without looking at such log file, it seems to be either a numerical issue or a pure model issue

Cheers,

Olivier

Revision history for this message
Giacomo Magni (gmagni) said : 		#2

Hi Olivier,
yes something was wrong in /SubProcesses/P2_gg_lvllvl/G1 because there is no .lhe file.

I attached you the link to the log file.

https://www.dropbox.com/s/heqdor36dcpjp5c/log.txt?dl=0

The result.dat file reports:

    0.23007E-02 0.19861E-03 0.19861E-03 2 113861683 9 2 0.86930E+03 0.39260E+05 0.23007E-02 0.0 0.0 0

And the input_app.txt:

        1000 1 1 !Number of events and max and min iterations
   0.02 !Accuracy
  -2 !Grid Adjustment 0=none, 2=adjust
  1 !Suppress Amplitude 1=yes
  1 !Helicity Sum/event 0=exact
  1

Thank you very much for your help.

Best,
Giacomo

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

G1 seems fine actually.

By the way which technical cut did you use to kill the 0 pt for the W in this channel?
If you do not put cut on the PT of the W this is likely divergent.

This might be related to your issue actually.

Cheers,

Olivier

Revision history for this message
Giacomo Magni (gmagni) said : 		#4

Hi,
so that may be the issue, I set only:

ptl 10
etal 2.5

Best,
Giacomo

Revision history for this message
Bryan (bryankortman) said : 		#5

Hello,

I have encountered a similar issue with the
generate g g > mu- vm~ e+ ve [noborn=QCD] QED=4 NP=0 @1
add process g g > mu+ vm e- ve~ [noborn=QCD] QED=4 NP=0 @2
generation.

Did the PTmin cut of the W solve the issue?
I am wondering since I am successfully generating the SM contribution with these settings.

However the ctW interference with the SM is running into the same problem. Even with the pTmin cut on the W employed like:
'pt_min_pdg':{24:1.0}

Best,
Bryan

1:39:15 INFO: Survey finished for P0_gg_llvlvl/G3.26 at 4.38e-07 +- 1.82e-08 pb
01:39:52 INFO: Idle: 0, Running: 0, Completed: 220 [ 38m 50s ]
01:39:52 INFO: Idle: 0, Running: 0, Completed: 220 [ 38m 50s ]
01:39:52 INFO: End survey
01:39:52 refine 15000
01:39:52 Creating Jobs
01:39:52 INFO: Refine results to 15000
01:39:52 INFO: Generating 15000.0 unweigthed events.
01:39:52 INFO: Effective Luminosity 51091714.4269 pb^-1
01:39:52 INFO: need to improve 9 channels
01:39:52 WARNING: Channel G1//data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/P0_gg_llvlvl has a very low efficiency of unweighting. Might not be possible to reach target 
01:39:52 Current estimate of cross-section: 0.000124792965929 +- 4.70401221541e-05

.........

03:35:33 INFO: Idle: 8, Running: 16, Completed: 60 [ 1h 55m ]
03:44:17 INFO: Idle: 7, Running: 16, Completed: 61 [ 2h 4m ]
03:46:36 INFO: Idle: 6, Running: 16, Completed: 62 [ 2h 6m ]
03:47:43 INFO: Idle: 5, Running: 16, Completed: 64 [ 2h 7m ]
03:48:43 INFO: Idle: 4, Running: 16, Completed: 65 [ 2h 8m ]
03:48:59 INFO: Idle: 4, Running: 16, Completed: 66 [ 2h 9m ]
04:10:13 INFO: Idle: 3, Running: 16, Completed: 67 [ 2h 30m ]
04:10:14 INFO: P0_gg_llvlvl/G3.22 is at 513/1003 ('51%') event. Resubmit 9 job at iteration 3.
04:10:14 WARNING: float division by zero 
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 97799 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 97800 | tee -a $k
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 98663 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 98664 | tee -a $k
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 101698 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 101699 | tee -a $k
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 103021 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 103022 | tee -a $k
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 103051 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 103052 | tee -a $k
04:10:14 /data/atlas/users/bkortman/MGEvntGen_output/ggWW_smefit/emununu_ci_WWcomb/mc.SMEFTatNLO_full_ggWW_emu_NPsq2_ctW_smefit.100012.000/test/PROC_ggWW_NLO_ctW/SubProcesses/refine_splitted.sh: line 61: 103179 Terminated ../madevent 2>&1 >> $k < input_app.txt
04:10:14 103180 | tee -a $k

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#6

Hi, "pt_min_pdg" only impacts final state particles.
With your syntax you do not have any W in the final state and therefore
'pt_min_pdg':{24:1.0}
does not cut anything

Cheers,

Olivier

Revision history for this message
Bryan (bryankortman) said : 		#7

Hello Olivier,

Thanks for for the comment. I did not know that.

Is there a way I can keep the 4 body final state as described above and put a cut on the V boson pT to avoid divergences?

Best,
Bryan

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#8

The only solution is to define such cuts at the fortran level.
(via the file dummy_fcts.f)

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Giacomo Magni for more information if necessary.

To post a message you must log in.