Problem generating e+e->jj @NLO

Asked by Javier Fernandez

Hi,

I'm trying to generate process e+e- > j j @ NLO, the program crashes when doing event reweight, it looks like if the events were not generated[1]

I attach datacards[2], and debug.log [3].

Any help is very apreciated.
Thank you very much in advance,
Javi

[1]
      Intermediate results:
      Random seed: 33
      Total cross section: 3.226e+04 +- 3.5e+01 pb
      Total abs(cross section): 3.811e+04 +- 3.3e+01 pb

INFO: Generating events
INFO: Idle: 0, Running: 1, Completed: 15 [ current time: 12h07 ]
INFO: Idle: 0, Running: 0, Completed: 16 [ 1s ]
INFO: Doing reweight
INFO: Idle: 0, Running: 0, Completed: 16 [ current time: 12h07 ]
tail: cannot open `/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aMCatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/SubProcesses/P0_epem_uux/GF2_1/events.lhe.rwgt' for reading: No such file or directory
Error detected in "launch -n pilotrun"
write debug file /afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aMCatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/pilotrun_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
aMCatNLOError : An error occurred during reweight. Check the'reweight_xsec_events.output' files inside the 'SubProcesses/P*/G*/ directories for details

[2] proc.card
import model loop_sm-ckm

generate e+ e- > j j [QCD] @0

output eejjNLObmass -nojpeg

run.card
#***********************************************************************
# 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. *
#***********************************************************************
#
#*******************
# 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 *
#***********************************************************************
 300 = nevents ! Number of unweighted events requested
 0.001 = req_acc ! Required accuracy (-1=auto determined from nevents)
 20 = 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 ! average or sum
#***********************************************************************
# 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 = lpp1 ! beam 1 type (0 = no PDF)
 0 = lpp2 ! beam 2 type (0 = no PDF)
 45.6 = ebeam1 ! beam 1 energy in GeV
 45.6 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 lhapdf = pdlabel ! PDF set
 292200, 90000, 90200, 90400, 13100, 13164, 13166, 25200, 25260, 11000, 11067, 1
1069 = 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) *
#***********************************************************************
 False = fixed_ren_scale ! if .true. use fixed ren scale
 False = fixed_fac_scale ! if .true. use fixed fac scale
 91.118 = muR_ref_fixed ! fixed ren reference scale
 91.118 = 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.
 True = 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
 10.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 *
#***********************************************************************
 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)
#***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
#***********************************************************************

[3] #************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.4.2 20xx-xx-xx *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* and *
#* http://amcatnlo.cern.ch *
#* *
#************************************************************
#* *
#* Command File for aMCatNLO *
#* *
#* run as ./bin/aMCatNLO.py filename *
#* *
#************************************************************
launch -n pilotrun
Traceback (most recent call last):
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/extended_
cmd.py", line 1009, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/extended_
cmd.py", line 964, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/amcatnlo_
run_interface.py", line 1205, in do_launch
    evt_file = self.run(mode, options)
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/amcatnlo_
run_interface.py", line 1445, in run
    return self.reweight_and_collect_events(options, mode, nevents, event_norm)
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/amcatnlo_
run_interface.py", line 2781, in reweight_and_collect_events
    scale_pdf_info = self.run_reweight(options['reweightonly'])
  File "/afs/cern.ch/work/j/jfernan/public/ttGen/genproductions/bin/MadGraph5_aM
CatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/processtmp/bin/internal/amcatnlo_
run_interface.py", line 3598, in run_reweight
    '\'SubProcesses/P*/G*/ directories for details')
aMCatNLOError: An error occurred during reweight. Check the'reweight_xsec_events
.output' files inside the 'SubProcesses/P*/G*/ directories for details
Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
                  timeout : 60
       cluster_local_path : /cvmfs/cp3.uclouvain.be/madgraph/
            cluster_queue : local
         madanalysis_path : None
                   lhapdf : /cvmfs/cms.cern.ch/slc6_amd64_gcc481/ext
ernal/lhapdf/6.1.6-ddibom/share/LHAPDF//../../bin/lhapdf-config
             cluster_size : 100
           cluster_memory : None
                    pjfry : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : None
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : None
          pythia-pgs_path : None
                 mg5_path : /afs/cern.ch/work/j/jfernan/public/ttGen
/genproductions/bin/MadGraph5_aMCatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/M
G5_aMC_v2_4_2
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
         cluster_nb_retry : 1
         fortran_compiler : None
                  nb_core : 8
              auto_update : 0
      exrootanalysis_path : None
               eps_viewer : None
             syscalc_path : /afs/cern.ch/work/j/jfernan/public/ttGen
/genproductions/bin/MadGraph5_aMCatNLO/eejjNLObmass/eejjNLObmass_gridpack/work/M
G5_aMC_v2_4_2/SysCalc
                  fastjet : None
                    golem : None
             cpp_compiler : None
      notification_center : True
                 run_mode : 2

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Rikkert Frederix Edit question
Solved by:
Javier Fernandez
Solved:
Last query:
Last reply:
Revision history for this message
Rikkert Frederix (frederix) said :
#1

Dear Javi,

It seems like you are trying to reweight with different PDF sets. From you run_card:

 lhapdf = pdlabel ! PDF set
 292200, 90000, 90200, 90400, 13100, 13164, 13166, 25200, 25260, 11000, 11067, 11069 = 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.

Of course, this doesn't really make sense for lepton-lepton collisions. I would have thought that the code recognised this behaviour and automatically skip it, but this obviously is not happening. If you replace it by

 lhapdf = pdlabel
 292200 = lhaid

the reweighting code should run smoothly. I'll put a fix in the development version of the MG5_aMC code to automatically detect such input and make sure that the code runs even with the long list of PDF sets.

best,
Rikkert

Revision history for this message
Rikkert Frederix (frederix) said :
#2

Addendum:

similarly, you should change

True = reweight_PDF

to

False = reweight_PDF

Best regards,
Rikkert

Revision history for this message
Javier Fernandez (javier-fernandez-h) said :
#3

Thanks Rikkert, this solved the issue. Naively I had thought that with e+e- input the PDF stuff was ignored since the following message is thrown:

INFO: Lepton-Lepton collision: Ignoring 'pdlabel' and 'lhaid' in the run_card.

Thank you
Javi