fortran runtime erro: Bad integer for item 15 in list input

Hi,

Thanks for reporting this, will look at this, could you indicate what is your fortran compiler version?
It is likely compiler specific and this will help me to figure out what the issue is.

Cheers,

Olivier

GNU Fortran (GCC) 4.8.5.

On Thu, Jul 13, 2023 at 5:20 PM Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #707272 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/707272
>
> Status: Open => Needs information
>
> Olivier Mattelaer requested more information:
> Hi,
>
> Thanks for reporting this, will look at this, could you indicate what is
> your fortran compiler version?
> It is likely compiler specific and this will help me to figure out what
> the issue is.
>
> Cheers,
>
> Olivier
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https://answers.launchpad.net/mg5amcnlo/+question/707272
>
> You received this question notification because you asked the question.
>

Hi.
When I make the reweight_scale False, then the process run and doesn't give any Fortran runtime error. But I am not prety much sure about this changes, either it cause some problems in the lhe file or not.
Here is the run_card.dat with the above given changes:
#***********************************************************************
# 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 *
#***********************************************************************
  10000 = 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).
#***********************************************************************
# Output format
#***********************************************************************
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  average = event_norm ! average/sum/bias. Normalization of the weight in the LHEF
#***********************************************************************
# 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 *
# 0 = no PDF *
# 1/-1 = proton/antiproton *
# 3/-3 = electron/positron with ISR/Beamstrahlung; *
# 4/-4 = muon/antimuon with ISR/Beamstrahlung; *
#***********************************************************************
  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. *
#***********************************************************************
  lhapdf = 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.
  0 = pdfscheme ! the scheme of the input PDFs. 0->MSbar; 1->DIS
             ! 2->eta (leptonic); 3->beta (leptonic)
      ! 4->mixed (leptonic); 6->delta (leptonic)
                    ! if pdlabel==emela, this is set automatically
#***********************************************************************
# The following block is specific to lepton collisions (lpp=+-3) *
#***********************************************************************
  True = photons_from_lepton ! whether to include or not photons from
                              ! lepton ISR
#***********************************************************************
# 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
  False = mcatnlo_delta ! use MC@NLO-Delta matching, arXiv:2002.12716
                        ! (only with Pythia8309 or later)
#***********************************************************************
# 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
  344.0 = mur_ref_fixed ! fixed ren reference scale
  344.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
  False = 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
#***********************************************************************
# Customization of the code. List of files containing user hook function
#***********************************************************************
   = custom_fcts ! List of files containing user hook function
#***********************************************************************
# 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. *
# - If gamma_is_j, photons are also clustered with jets. *
# Otherwise, they will be treated as tagged particles and photon *
# isolation will be applied. Note that photons in the real emission *
# will always be clustered with QCD partons. *
# - 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
  10.0 = ptj ! Min jet transverse momentum
  -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
  False = gamma_is_j ! Wether to cluster photons as jets or not
#***********************************************************************
# 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
#***********************************************************************
# Fermion-photon recombination parameters *
# If Rphreco=0, no recombination is performed *
#***********************************************************************
  0.1 = rphreco ! Minimum fermion-photon distance for recombination
  -1.0 = etaphreco ! Maximum abs(pseudo-rap) for photons to be recombined (a value .lt.0 means no cut)
  False = lepphreco ! Recombine photons and leptons together
  False = quarkphreco ! Recombine photons and quarks together
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# Not applied if gamma_is_j *
# When ptgmin=0, all the other parameters are ignored *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  20.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
#***********************************************************************
# Use PineAPPL to generate PDF-independent fast-interpolation grid *
# (https://zenodo.org/record/3992765#.X2EWy5MzbVo) *
#***********************************************************************
  False = pineappl ! PineAPPL switch
#***********************************************************************
# Folding parameters for S-events to reduce the number of negatively *
# weighted events. Allowed values are 1, 2, 4 or 8 for each of the *
# three variables. Typically, folding in xi_i or y_ij results in the *
# largest reduction of negatively weighted events. (arXiv:2002.12716) *
#***********************************************************************
  1, 1, 1 = folding ! correspond to folding in xi_i, y_ij, and phi_i
#***********************************************************************

So I do reproduce the error in version 3.5.1 but not in our Long Term Stable version (I'm using the version 2.9.16 which is not yet released, but 2.9.15 should be the same).
So another option is to move back to our Long Term Stable release.

@Rik: Do you think that this might be related to MCatNLO-Delta?

Cheers,

Olivier

Thanks Olivier Mattelaer, that solved my question.

Thank you very much.
I installed the LTS version(2.9.15) which works fine, and generates the NLO
processes without any Fortran runtime error.

On Wed, Jul 19, 2023 at 10:20 PM Muhammad sahil <
<email address hidden>> wrote:

> Your question #707272 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/707272
>
> Status: Answered => Solved
>
> You confirmed that the question is solved:
> Thanks Olivier Mattelaer, that solved my question.
>
> --
> You received this question notification because you asked the question.
>

Hi both,

The bug in 3.5.0 with the reweighting has nothing to do with the Delta implementation. Rather, it has been there from the very first v3 of the code: there is a space missing between two integers, and, therefore, if the latter of the two becomes to large, it merges with the first one in the line. Hence, the fix is a rather easy one. At line 3793 of SubProocesses/fks_singular.f replace

30 format(i15,i2,6(1x,d14.8),6(1x,i2),1x,i8,1x,d18.12,1x,d18.12)

by

30 format(i15,1x,i2,6(1x,d14.8),6(1x,i2),1x,i8,1x,d18.12,1x,d18.12)

Note the extra 1x in the 2nd position.

Olivier, which git branch should this be committed to?

best,
Rikkert

Thanks Rik,

For the branch it should be in the 3.5.2 branch.

Thanks,

Olivier

> On 20 Jul 2023, at 20:45, Rikkert Frederix <email address hidden> wrote:
>
> Question #707272 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/707272
>
> Rikkert Frederix posted a new comment:
> Hi both,
>
> The bug in 3.5.0 with the reweighting has nothing to do with the Delta
> implementation. Rather, it has been there from the very first v3 of the
> code: there is a space missing between two integers, and, therefore, if
> the latter of the two becomes to large, it merges with the first one in
> the line. Hence, the fix is a rather easy one. At line 3793 of
> SubProocesses/fks_singular.f replace
>
> 30 format(i15,i2,6(1x,d14.8),6(1x,i2),1x,i8,1x,d18.12,1x,d18.12)
>
> by
>
> 30 format(i15,1x,i2,6(1x,d14.8),6(1x,i2),1x,i8,1x,d18.12,1x,d18.12)
>
> Note the extra 1x in the 2nd position.
>
> Olivier, which git branch should this be committed to?
>
> best,
> Rikkert
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.