MadGraph5_aMC@NLO

Command "launch auto " interrupted with error: Exception : Some of the run failed:

Asked by celik

Hello Dear Experts,

I am generating a process p p > z t t~ [QCD] however each time I am getting an error saying "launch auto interrepted with error.....". If i generate only 20 events, there is no issue popping up at all but when i run it with at least 100 events, then error message arises.

I have also no issue with generating LO processes. It happens only with zttbar NLO process and running MG on my linux computer. I asssume it has something to do with reweighting.
As I am newbie, I would appreciate any help. I am also pasting the content of the log file.

Thanks in advance,

Ali

Below is terminal output.
INFO: Events generated
reweight -from_cards
INFO: split the event file in bunch of 2500 events
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h35 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 1m 11s ]
INFO: Collect and combine the various output file.
Command "launch auto " interrupted with error:
Exception : Some of the run failed: Please read

The following lines are from log file:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.4 2018-11-09 *
#* *
#* 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 auto
Traceback (most recent call last):
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1501, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/amcatnlo_run_interface.py", line 1685, in do_launch
    self.exec_cmd('reweight -from_cards', postcmd=False)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/common_run_interface.py", line 2099, in do_reweight
    raise Exception, "Some of the run failed: Please read %s_%s_debug.log" % (f, self.run_tag)
Exception: Some of the run failed: Please read /home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/Events/run_03/events.lhe.gz_0.lhe_tag_1_debug.log
Value of current Options:
              text_editor : None
      notification_center : True
                    pjfry : None
       cluster_local_path : None
  default_unset_couplings : 99
       group_subprocesses : Auto
ignore_six_quark_processes : False
    loop_optimized_output : True
    cluster_status_update : (600, 30)
         fortran_compiler : None
               hepmc_path : None
                  collier : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
              auto_update : 7
             pythia8_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/pythia8
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : /home/ali/Desktop/MG5_aMC_v2_6_4/pythia-pgs
                  td_path : None
             delphes_path : /home/ali/Desktop/MG5_aMC_v2_6_4/Delphes
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/madanalysis5/madanalysis5
      exrootanalysis_path : /home/ali/Desktop/MG5_aMC_v2_6_4/ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : fastjet-config
                 run_mode : 2
              web_browser : /opt/google/chrome/google-chrome
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : None
         madanalysis_path : None
                   lhapdf : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lhapdf6/bin/lhapdf-config
             stdout_level : 20
                  nb_core : 4
            f2py_compiler : /usr/local/bin/f2py
                    ninja : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/MG5aMC_PY8_interface
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                 mg5_path : /home/ali/Desktop/MG5_aMC_v2_6_4
                  timeout : 900
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.4 2018-11-09 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_sm
generate p p > z t t~ [QCD]
output 31_ocak_zttbar_nlo
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR LOOP
###################################
Block loop
    1 9.118800e+01 # MU_R

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 80.419002 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.700000e+00 # ymb
    6 1.730000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 82 # gh
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
#***********************************************************************
# 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 *
#***********************************************************************
  200 = 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)
  7000.0 = ebeam1 ! beam 1 energy in GeV
  7000.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
  lhapdf = pdlabel ! PDF set
  267200 = 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!!!! *
#***********************************************************************
  HERWIG6 = 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 *
#***********************************************************************
  True = 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)
#***********************************************************************
# 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)
#***********************************************************************

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi Ali,
so, can you confirm that you get some printout on the screen saying that events have been generated and with a total cross section?
Is there any summary.txt file inside Events/run_XX ? can you paste the content?

INFO: Events generated
reweight -from_cards
INFO: split the event file in bunch of 2500 events
INFO: Idle: 1, Running: 0, Completed: 0 [ current time: 12h35 ]
INFO: Idle: 0, Running: 0, Completed: 1 [ 1m 11s ]
INFO: Collect and combine the various output file.
Command "launch auto " interrupted with error:
Exception : Some of the run failed: Please read
/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/Events/run_03/events.lhe.gz_0.lhe_tag_1_debug.log

What is the content of /home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/Events/run_03/events.lhe.gz_0.lhe_tag_1_debug.log ?

Thanks!

Marco

Revision history for this message
celik (watzinki) said : 		#2

Hello Marco,
Thanks for your prompt reply indeed.
Events are generated with cross section. Dumping the content of the summary txt file and .log file respectively:

 --------------------------------------------------------------
      Summary:
      Process p p > z t t~ [QCD]
      Run at p-p collider (7000.0 + 7000.0 GeV)
      Number of events generated: 200
      Total cross section: 8.736e-01 +- 3.7e-03 pb
   --------------------------------------------------------------
Content of .log file:

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.6.4 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 *
#* *
#************************************************************
reweight /home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/Events/run_03/events.lhe.gz_0.lhe -from_cards --multicore=create
Traceback (most recent call last):
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/extended_cmd.py", line 1501, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/common_run_interface.py", line 2140, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1643, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 546, in do_launch
    weight = self.calculate_weight(event)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 954, in calculate_weight
    return self.calculate_nlo_weight(event)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 1028, in calculate_nlo_weight
    w_origV = self.calculate_matrix_element(cevent, 'V0', scale2=scale2)
  File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 1155, in calculate_matrix_element
    orig_order, Pdir, hel_dict = self.id_to_path[tag]
KeyError: (((-3, 3), (-6, 6, 23)), 'V')
Value of current Options:
              text_editor : None
      notification_center : True
                    pjfry : None
       cluster_local_path : None
    cluster_status_update : (600, 30)
               hepmc_path : None
          pythia-pgs_path : /home/ali/Desktop/MG5_aMC_v2_6_4/pythia-pgs
              thepeg_path : None
        madanalysis5_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/madanalysis5/madanalysis5
                      OLP : MadLoop
                 applgrid : applgrid-config
                 run_mode : 2
        cluster_temp_path : None
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lhapdf6/bin/lhapdf-config
      output_dependencies : external
            f2py_compiler : /usr/local/bin/f2py
                    ninja : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
   automatic_html_opening : True
       cluster_retry_wait : 300
      exrootanalysis_path : /home/ali/Desktop/MG5_aMC_v2_6_4/ExRootAnalysis
                  timeout : 900
                  nb_core : 4
         fortran_compiler : None
                  collier : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
             pythia8_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/pythia8
                hwpp_path : None
                    golem : None
                  td_path : None
             delphes_path : /home/ali/Desktop/MG5_aMC_v2_6_4/Delphes
              auto_update : 7
             cluster_type : condor
               eps_viewer : None
                  fastjet : fastjet-config
              web_browser : /opt/google/chrome/google-chrome
             cluster_size : 100
           cluster_memory : None
                  amcfast : amcfast-config
             stdout_level : None
             cluster_time : None
mg5amc_py8_interface_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/MG5aMC_PY8_interface
         cluster_nb_retry : 1
                 mg5_path : /home/ali/Desktop/MG5_aMC_v2_6_4
             syscalc_path : None
             cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.4 2018-11-09 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#************************************************************
set default_unset_couplings 99
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model loop_sm
generate p p > z t t~ [QCD]
output 31_ocak_zttbar_nlo
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

###################################
## INFORMATION FOR LOOP
###################################
Block loop
    1 9.118800e+01 # MU_R

###################################
## INFORMATION FOR MASS
###################################
Block mass
    5 4.700000e+00 # MB
    6 1.730000e+02 # MT
   15 1.777000e+00 # MTA
   23 9.118800e+01 # MZ
   25 1.250000e+02 # MH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
  1 0.000000 # d : 0.0
  2 0.000000 # u : 0.0
  3 0.000000 # s : 0.0
  4 0.000000 # c : 0.0
  11 0.000000 # e- : 0.0
  12 0.000000 # ve : 0.0
  13 0.000000 # mu- : 0.0
  14 0.000000 # vm : 0.0
  16 0.000000 # vt : 0.0
  21 0.000000 # g : 0.0
  22 0.000000 # a : 0.0
  24 80.419002 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))

###################################
## INFORMATION FOR SMINPUTS
###################################
Block sminputs
    1 1.325070e+02 # aEWM1
    2 1.166390e-05 # Gf
    3 1.180000e-01 # aS

###################################
## INFORMATION FOR YUKAWA
###################################
Block yukawa
    5 4.700000e+00 # ymb
    6 1.730000e+02 # ymt
   15 1.777000e+00 # ymtau

###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 1.491500e+00 # WT
DECAY 23 2.441404e+00 # WZ
DECAY 24 2.047600e+00 # WW
DECAY 25 6.382339e-03 # WH
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # ta- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
#===========================================================
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#===========================================================

Block QNUMBERS 82 # gh
        1 0 # 3 times electric charge
        2 1 # number of spin states (2S+1)
        3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
        4 1 # Particle/Antiparticle distinction (0=own anti)
#***********************************************************************
# 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 *
#***********************************************************************
  200 = 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)
  7000.0 = ebeam1 ! beam 1 energy in GeV
  7000.0 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
  lhapdf = pdlabel ! PDF set
  267200 = 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!!!! *
#***********************************************************************
  HERWIG6 = 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 *
#***********************************************************************
  True = 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)
#***********************************************************************
# 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)
#***********************************************************************

Revision history for this message
marco zaro (marco-zaro) said : 		#3

Hi Ali,
so, the error occurs in the reweight stage. I do not know if this is really needed, as one does already the scale/PDF reweighting before the printout of the summary.
Olivier, can you comment on this?
cheers,

Marco

> On 31 Jan 2019, at 11:12, ali <email address hidden> wrote:
>
> Question #678222 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/678222
>
> Status: Needs information => Open
>
> ali gave more information on the question:
> Hello Marco,
> Thanks for your prompt reply indeed.
> Events are generated with cross section. Dumping the content of the summary txt file and .log file respectively:
>
>
> --------------------------------------------------------------
> Summary:
> Process p p > z t t~ [QCD]
> Run at p-p collider (7000.0 + 7000.0 GeV)
> Number of events generated: 200
> Total cross section: 8.736e-01 +- 3.7e-03 pb
> --------------------------------------------------------------
> Content of .log file:
>
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 5.2.6.4 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 *
> #* *
> #************************************************************
> reweight /home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/Events/run_03/events.lhe.gz_0.lhe -from_cards --multicore=create
> Traceback (most recent call last):
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/extended_cmd.py", line 1501, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/bin/31_ocak_zttbar_nlo/bin/internal/common_run_interface.py", line 2140, in do_reweight
> reweight_cmd.import_command_file(path)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1643, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/various/misc.py", line 100, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 546, in do_launch
> weight = self.calculate_weight(event)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 954, in calculate_weight
> return self.calculate_nlo_weight(event)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 1028, in calculate_nlo_weight
> w_origV = self.calculate_matrix_element(cevent, 'V0', scale2=scale2)
> File "/home/ali/Desktop/MG5_aMC_v2_6_4/madgraph/interface/reweight_interface.py", line 1155, in calculate_matrix_element
> orig_order, Pdir, hel_dict = self.id_to_path[tag]
> KeyError: (((-3, 3), (-6, 6, 23)), 'V')
> Value of current Options:
> text_editor : None
> notification_center : True
> pjfry : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> hepmc_path : None
> pythia-pgs_path : /home/ali/Desktop/MG5_aMC_v2_6_4/pythia-pgs
> thepeg_path : None
> madanalysis5_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/madanalysis5/madanalysis5
> OLP : MadLoop
> applgrid : applgrid-config
> run_mode : 2
> cluster_temp_path : None
> cluster_queue : None
> madanalysis_path : None
> lhapdf : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lhapdf6/bin/lhapdf-config
> output_dependencies : external
> f2py_compiler : /usr/local/bin/f2py
> ninja : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
> automatic_html_opening : True
> cluster_retry_wait : 300
> exrootanalysis_path : /home/ali/Desktop/MG5_aMC_v2_6_4/ExRootAnalysis
> timeout : 900
> nb_core : 4
> fortran_compiler : None
> collier : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/lib
> pythia8_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/pythia8
> hwpp_path : None
> golem : None
> td_path : None
> delphes_path : /home/ali/Desktop/MG5_aMC_v2_6_4/Delphes
> auto_update : 7
> cluster_type : condor
> eps_viewer : None
> fastjet : fastjet-config
> web_browser : /opt/google/chrome/google-chrome
> cluster_size : 100
> cluster_memory : None
> amcfast : amcfast-config
> stdout_level : None
> cluster_time : None
> mg5amc_py8_interface_path : /home/ali/Desktop/MG5_aMC_v2_6_4/HEPTools/MG5aMC_PY8_interface
> cluster_nb_retry : 1
> mg5_path : /home/ali/Desktop/MG5_aMC_v2_6_4
> syscalc_path : None
> cpp_compiler : None
> #************************************************************
> #* MadGraph5_aMC@NLO *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.4 2018-11-09 *
> #* *
> #* The MadGraph5_aMC@NLO Development Team - Find us at *
> #* https://server06.fynu.ucl.ac.be/projects/madgraph *
> #* *
> #************************************************************
> #* *
> #* Command File for MadGraph5_aMC@NLO *
> #* *
> #* run as ./bin/mg5_aMC filename *
> #* *
> #************************************************************
> set default_unset_couplings 99
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set low_mem_multicore_nlo_generation False
> set loop_color_flows False
> set gauge unitary
> set complex_mass_scheme False
> set max_npoint_for_channel 0
> import model sm
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model loop_sm
> generate p p > z t t~ [QCD]
> output 31_ocak_zttbar_nlo
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> ######################################################################
> ## ##
> ## Width set on Auto will be computed following the information ##
> ## present in the decay.py files of the model. ##
> ## See arXiv:1402.1178 for more details. ##
> ## ##
> ######################################################################
>
> ###################################
> ## INFORMATION FOR LOOP
> ###################################
> Block loop
> 1 9.118800e+01 # MU_R
>
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> Block mass
> 5 4.700000e+00 # MB
> 6 1.730000e+02 # MT
> 15 1.777000e+00 # MTA
> 23 9.118800e+01 # MZ
> 25 1.250000e+02 # MH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> 1 0.000000 # d : 0.0
> 2 0.000000 # u : 0.0
> 3 0.000000 # s : 0.0
> 4 0.000000 # c : 0.0
> 11 0.000000 # e- : 0.0
> 12 0.000000 # ve : 0.0
> 13 0.000000 # mu- : 0.0
> 14 0.000000 # vm : 0.0
> 16 0.000000 # vt : 0.0
> 21 0.000000 # g : 0.0
> 22 0.000000 # a : 0.0
> 24 80.419002 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
>
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> Block sminputs
> 1 1.325070e+02 # aEWM1
> 2 1.166390e-05 # Gf
> 3 1.180000e-01 # aS
>
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> Block yukawa
> 5 4.700000e+00 # ymb
> 6 1.730000e+02 # ymt
> 15 1.777000e+00 # ymtau
>
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> DECAY 6 1.491500e+00 # WT
> DECAY 23 2.441404e+00 # WZ
> DECAY 24 2.047600e+00 # WW
> DECAY 25 6.382339e-03 # WH
> ## Dependent parameters, given by model restrictions.
> ## Those values should be edited following the
> ## analytical expression. MG5 ignores those values
> ## but they are important for interfacing the output of MG5
> ## to external program such as Pythia.
> DECAY 1 0.000000 # d : 0.0
> DECAY 2 0.000000 # u : 0.0
> DECAY 3 0.000000 # s : 0.0
> DECAY 4 0.000000 # c : 0.0
> DECAY 5 0.000000 # b : 0.0
> DECAY 11 0.000000 # e- : 0.0
> DECAY 12 0.000000 # ve : 0.0
> DECAY 13 0.000000 # mu- : 0.0
> DECAY 14 0.000000 # vm : 0.0
> DECAY 15 0.000000 # ta- : 0.0
> DECAY 16 0.000000 # vt : 0.0
> DECAY 21 0.000000 # g : 0.0
> DECAY 22 0.000000 # a : 0.0
> #===========================================================
> # QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
> #===========================================================
>
> Block QNUMBERS 82 # gh
> 1 0 # 3 times electric charge
> 2 1 # number of spin states (2S+1)
> 3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # Particle/Antiparticle distinction (0=own anti)
> #***********************************************************************
> # 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 *
> #***********************************************************************
> 200 = 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)
> 7000.0 = ebeam1 ! beam 1 energy in GeV
> 7000.0 = ebeam2 ! beam 2 energy in GeV
> #***********************************************************************
> # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> #***********************************************************************
> lhapdf = pdlabel ! PDF set
> 267200 = 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!!!! *
> #***********************************************************************
> HERWIG6 = 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 *
> #***********************************************************************
> True = 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)
> #***********************************************************************
> # 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)
> #***********************************************************************
>
> --
> You received this question notification because you are assigned to this
> question.

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

Hi,

What kind of reweighting are you doing?
Can you pass your reweight_card.dat?

Note that such function is NOT for scale/PDF but for generating additional weights for a BSM model.

Cheers,

Olivier

Revision history for this message
celik (watzinki) said : 		#5

Hi Olivier,
I use the default reweight card without any modifications in it. I just turned on reweight="NLO " flag along with Herwig, Madanalysis etc on the prompt screen right before the launch command.

#*************************************************************************
# Reweight Module *
# Matrix-Element reweighting at LO/NLO *
# Mattelaer Olivier arxiv:1607.00763 *
#*************************************************************************
#
# Note:
# 1) the value of alpha_s will be used from the event so the value in
# the param_card is not taken into account.
# 2) It is (in general) dangerous/wrong to change parameters by a large
# amount, if this changes the shape of the matrix elements a lot.
# (For example, changing a particle's mass by much more than its
# width leads to very inaccurate result). In such a case, separate
# event generation runs are needed.
#
#************************************************************************
# ENTER YOUR COMMANDS BELOW.
#************************************************************************

 change mode NLO
                   # has no effect since only "LO" mode is allowed.

launch
# SPECIFY A PATH OR USE THE SET COMMAND LIKE THIS:
# set sminputs 1 130 # modify 1/alpha_EW

#************************************************************************
# Manual: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/Reweight
#************************************************************************
#
# Example of (standard) code for the computation of two weights:
#
# launch ! tag to start the computation of the first weight
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# set BLOCKNAME ID VALUE ! rule to modify the current param_card
# launch ! start to compute a second weight
# /home/Cards/param_card_2.dat ! you can also enter a path to a valid card
#
# Note: The command to specify the parameter are AFTER the associated "launch"
#
# Possible options:
# You can enter one of the following lines to customize the reweighting
# procedure. These need to be given before the 'launch' command.
#
# change model NAME : use another model for the matrix-elements to reweight
# with. In this case you need to provide the path to a correct
# param_card for the new model; you cannot modify the original one
# with the 'set' command.
# change process DEF [--add]: change the process by which you reweight.
# The initial and final state particles of the new process should
# be exactly identical to the ones in the original process.
# change helicity False: perform the reweighting by helicity summed
# matrix-elements even if the events have been written with a
# single helicity state.
# change mode XXX: change the type of reweighting performed.
# allowed values: LO, NLO, LO+NLO
# - This command has no effect for reweighting an .lhe event file with LO accuracy.
# In that case LO mode is always used (whatever entry is set).
# - When the .lhe file reweighted is at NLO accuracy, then all modes are allowed.
# * "LO" is an approximate leading order method
# * "NLO" is the NLO accurate method
# * "LO+NLO" runs both
# - "NLO" and "LO+NLO" modes requires 'store_rwgt_info' equals True (run_card.dat)
# If the reweighting is done at generation level this parameter will
# automatically be set on True.
#************************************************************************

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

Hi,

Then this is the normal behavior. They are not good default for the re-weighitng module.
The only good default is to have reweight=OFF

Setting it to NLO but not specify any model/new benchmark point is just a waste of cpu (since the weigths will always be one)
so I make the code to crash in this case...

Cheers,

Olivier

Revision history for this message
celik (watzinki) said : 		#7

Thanks Olivier!
One more question. Does setting reweight=off omit pdf reweighting? I mean, setting it "OFF" would also effect/ turn off applying pdf reweighting? I am using lhapdf and made required modifications in the run_card.dat .

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

systematics computation is something completly different. So even with reweight=OFF, you will have it.

Cheers,

Olivier

Revision history for this message
celik (watzinki) said : 		#9

Thanks Olivier Mattelaer, that solved my question.