MadGraph5_aMC@NLO

Error while generating events at NLO

Asked by Antara Sarkar

Hi,

I was trying to study pp to z gamma at SM NLO using the following commands:

import model loop_qcd_qed_sm-full
define p = p b b~
define j = j b b~
generate p p > z a [QCD]
add process p p > z a j [QCD]

And then I do showering with Pythia8 and decay z to b b~ using madspin. However, on launching the process, I meet with this error:

INFO: Using built-in libraries for PDFs
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_uux_za...
INFO: Compiling P0_ccx_za...
INFO: Compiling P0_ddx_za...
INFO: Compiling P0_ssx_za...
INFO: P0_ddx_za done.
INFO: Compiling P0_uxu_za...
INFO: P0_uux_za done.
INFO: Compiling P0_cxc_za...
INFO: P0_ccx_za done.
INFO: Compiling P0_dxd_za...
INFO: P0_ssx_za done.
INFO: Compiling P0_sxs_za...
INFO: P0_dxd_za done.
INFO: Compiling P0_gu_zau...
INFO: P0_cxc_za done.
INFO: Compiling P0_gc_zac...
INFO: P0_sxs_za done.
INFO: Compiling P0_gd_zad...
INFO: P0_uxu_za done.
INFO: Compiling P0_gs_zas...
STOP pz is not conserved (flag:CT692)
INFO: P0_gc_zac done.
INFO: Compiling P0_gux_zaux...
INFO: P0_gu_zau done.
INFO: Compiling P0_gcx_zacx...
INFO: P0_gd_zad done.
INFO: Compiling P0_gdx_zadx...
STOP pz is not conserved (flag:CT692)
INFO: P0_gcx_zacx done.
INFO: Compiling P0_gsx_zasx...
INFO: P0_gux_zaux done.
INFO: Compiling P0_ug_zau...
INFO: P0_gdx_zadx done.
INFO: Compiling P0_uux_zag...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
INFO: P0_gs_zas done.
INFO: Compiling P0_cg_zac...
INFO: P0_ug_zau done.
INFO: Compiling P0_ccx_zag...
INFO: P0_uux_zag done.
INFO: Compiling P0_dg_zad...
STOP pz is not conserved (flag:CT692)
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
INFO: P0_cg_zac done.
INFO: Compiling P0_ddx_zag...
INFO: P0_gsx_zasx done.
INFO: Compiling P0_sg_zas...
INFO: P0_ccx_zag done.
INFO: Compiling P0_ssx_zag...
STOP pz is not conserved (flag:CT692)
INFO: P0_dg_zad done.
INFO: Compiling P0_uxg_zaux...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
INFO: P0_sg_zas done.
INFO: Compiling P0_uxu_zag...
INFO: P0_ddx_zag done.
INFO: Compiling P0_cxg_zacx...
INFO: P0_ssx_zag done.
INFO: Compiling P0_cxc_zag...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
INFO: P0_uxg_zaux done.
INFO: Compiling P0_dxg_zadx...
STOP pz is not conserved (flag:CT692)
INFO: P0_cxg_zacx done.
INFO: Compiling P0_dxd_zag...
INFO: P0_uxu_zag done.
INFO: Compiling P0_sxg_zasx...
INFO: P0_cxc_zag done.
INFO: Compiling P0_sxs_zag...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_DENORMAL
STOP pz is not conserved (flag:CT692)
INFO: P0_dxg_zadx done.
STOP pz is not conserved (flag:CT692)
INFO: P0_sxg_zasx done.
INFO: P0_dxd_zag done.
INFO: P0_sxs_zag done.
INFO: Checking test output:
INFO: P0_uux_za
INFO: Result for test_ME:
INFO: Output of the failing test:
 Enter 0 to compute MC/MC(limit)
       1 to compute MC/ME(limit)
       2 to compute ME/ME(limit)
 Enter xi_i, y_ij to be used in coll/soft tests
  Enter -2 to generate them randomly
 Enter number of tests for soft and collinear limits
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
  ****************************************

       NNPDFDriver version 1.0.3
   Grid: NNPDF23nlo_as_0119_qed_mem0.grid
  ****************************************
 New value of alpha_s from PDF nn23nlo: 0.11899999999999999
 Give FKS configuration number ("0" loops over all)

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 1
 FKS partons are: i= 5 j= 1
 with PDGs: i= 21 j= 2

 Enter graph number (iconfig), '0' loops over all graphs
Using random seed offsets: 1 , 1 , 0
  with seed 32
 Ranmar initialization seeds 10730 9407
tau_min 1 1 : 0.91188E+02 -- 0.24119E+03
tau_min 2 1 : 0.91188E+02 -- 0.24119E+03
tau_min 3 1 : 0.91188E+02 -- 0.24119E+03
tau_min 4 1 : 0.91188E+02 -- 0.24119E+03

 Soft tests done for (Born) config 1
 Failures: 100
     Soft test 1 FAILED. Fraction of failures: 1.00

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 2
 FKS partons are: i= 5 j= 2
 with PDGs: i= 21 j= -2

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 99
     Soft test 2 FAILED. Fraction of failures: 0.99

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 3
 FKS partons are: i= 5 j= 1
 with PDGs: i= -2 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 4
     Soft test 3 PASSED. Fraction of failures: 0.04

 No collinear test for massive j_fks

 =================================================

 NEW FKS CONFIGURATION:
 FKS configuration number is 4
 FKS partons are: i= 5 j= 2
 with PDGs: i= 2 j= 21

 Enter graph number (iconfig), '0' loops over all graphs

 Soft tests done for (Born) config 1
 Failures: 1
     Soft test 4 PASSED. Fraction of failures: 0.01

 No collinear test for massive j_fks
Error detected in "launch auto "
write debug file /home/ad/Tish/my_computation/MG5_aMC_v2_6_7/pp2za_qcdnlo/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : Some tests failed, run cannot continue.
Please check that widths of final state particles (e.g. top) have been set to 0 in the param_card.dat.
quit
INFO:

The run-card is as follows:

#***********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat aMC@NLO *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a hash (#) are info or comments *
# *
# mind the format: value = variable ! comment *
# *
# Some of the values of variables can be list. These can either be *
# comma or space separated. *
# *
# To display additional parameter, you can use the command: *
# update to_full *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of LHE events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
  100000 = nevents ! Number of unweighted events requested
  -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
  -1 = nevt_job ! Max number of events per job in event generation.
                 ! (-1= no split).
#***********************************************************************
# Normalize the weights of LHE events such that they sum or average to *
# the total cross section *
#***********************************************************************
  average = event_norm ! valid settings: average, sum, bias
#***********************************************************************
# Number of points per itegration channel (ignored for aMC@NLO runs) *
#***********************************************************************
  0.01 = req_acc_fo ! Required accuracy (-1=ignored, and use the
                    ! number of points and iter. below)
# These numbers are ignored except if req_acc_FO is equal to -1
  5000 = npoints_fo_grid ! number of points to setup grids
  4 = niters_fo_grid ! number of iter. to setup grids
  10000 = npoints_fo ! number of points to compute Xsec
  6 = niters_fo ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
  0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
  1 = lpp1 ! beam 1 type (0 = no PDF)
  1 = lpp2 ! beam 2 type (0 = no PDF)
  6500.0 = ebeam1 ! beam 1 energy in GeV
  6500.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
  nn23nlo = pdlabel ! PDF set
  244600 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
              ! numbers for central PDF sets are allowed. Can be a list;
              ! PDF sets beyond the first are included via reweighting.
#***********************************************************************
# Include the NLO Monte Carlo subtr. terms for the following parton *
# shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
# WARNING: PYTHIA6PT works only for processes without FSR!!!! *
#***********************************************************************
  PYTHIA8 = parton_shower
  1.0 = shower_scale_factor ! multiply default shower starting
                                  ! scale by this factor
#***********************************************************************
# Renormalization and factorization scales *
# (Default functional form for the non-fixed scales is the sum of *
# the transverse masses divided by two of all final state particles *
# and partons. This can be changed in SubProcesses/set_scales.f or via *
# dynamical_scale_choice option) *
#***********************************************************************
  True = fixed_ren_scale ! if .true. use fixed ren scale
  True = fixed_fac_scale ! if .true. use fixed fac scale
  125.0 = mur_ref_fixed ! fixed ren reference scale
  125.0 = muf_ref_fixed ! fixed fact reference scale
  -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
           ! dynamical choices. Can be a list; scale choices beyond the
           ! first are included via reweighting
  1.0 = mur_over_ref ! ratio of current muR over reference muR
  1.0 = muf_over_ref ! ratio of current muF over reference muF
#***********************************************************************
# Reweight variables for scale dependence and PDF uncertainty *
#***********************************************************************
  1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
  1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
  True = reweight_scale ! Reweight to get scale variation using the
            ! rw_rscale and rw_fscale factors. Should be a list of
            ! booleans of equal length to dynamical_scale_choice to
            ! specify for which choice to include scale dependence.
  False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
            ! list booleans of equal length to lhaid to specify for
            ! which PDF set to include the uncertainties.
#***********************************************************************
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#***********************************************************************
  False = store_rwgt_info ! Store info for reweighting in LHE file
#***********************************************************************
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#***********************************************************************
  0 = ickkw
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
# written in the LHE event file *
#***********************************************************************
  15.0 = bwcutoff
#***********************************************************************
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
  0.7 = jetradius ! The radius parameter for the jet algorithm
  100.0 = ptj ! Min jet transverse momentum
  -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  0.0 = ptl ! Min lepton transverse momentum
  -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
  0.0 = drll ! Min distance between opposite sign lepton pairs
  0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
  0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
  30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
# all the other parameters are ignored. *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  150.0 = ptgmin ! Min photon transverse momentum
  -1.0 = etagamma ! Max photon abs(pseudo-rap)
  0.4 = r0gamma ! Radius of isolation code
  1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
  1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
  True = isoem ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# Cuts associated to MASSIVE particles identified by their PDG codes. *
# All cuts are applied to both particles and anti-particles, so use *
# POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
# {6:100, 25:200} for multiple particles *
#***********************************************************************
  {} = pt_min_pdg ! Min pT for a massive particle
  {} = pt_max_pdg ! Max pT for a massive particle
  {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
#***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
  0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#***********************************************************************

Could you tell what is going wrong here?

Regards,
Antara

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

Not sure to understand why you use such syntax.

> import model loop_qcd_qed_sm-full

1) Why do you use that model?
2) Why do you use the "-full" type of that model? This will likely set a lot of mass for light quark which is problematic for NLO computation

> define p = p b b~
> define j = j b b~

This should be automatic if you use a five flavor model. (which you do not).
If you want to do a five flavor computation you need to use a five flavor model
(while the "-full" is likely a three flavor model)

With the following syntax (4 flavor computation here)
> import model sm
> generate p p > z a [QCD]
> add process p p > z a j [QCD]
> output
> launch

I do not have any issue (at least so far)

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Antara Sarkar for more information if necessary.

To post a message you must log in.