MadGraph5_aMC@NLO

ValueError (driver.f)

Asked by Sihyun Jeon

Hi all,

Now I am having different kind of errors when using madspin.

Please help me once more.

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.6.0 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
Traceback (most recent call last):
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/amcatnlo_run_interface.py", line 1243, in do_launch
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HNtoDilepDijet_M-40/bin/internal/common_run_interface.py", line 3463, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1580, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1465, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/madgraph/interface/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/madgraph/various/misc.py", line 99, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/MadSpin/interface_madspin.py", line 582, in do_launch
    generate_all.run()
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/MadSpin/decay.py", line 2114, in run
    self.get_max_weight_from_event(decay_mapping)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/MadSpin/decay.py", line 3177, in get_max_weight_from_event
    weight = self.get_max_weight_from_fortran(decay['path'], event_map,numberps,self.options['BW_cut'])
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/MadSpin/decay.py", line 3314, in get_max_weight_from_fortran
    path, std_in)
  File "/Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/MadSpin/decay.py", line 3348, in loadfortran
    maxweight=float(external.stdout.readline())
ValueError: could not convert string to float: At line 1506 of file driver.f

Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
                  timeout : 60
       cluster_local_path : None
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HEPTools/lhapdf6/bin/lhapdf-config
             cluster_size : 100
           cluster_memory : None
                    pjfry : None
    cluster_status_update : (600, 30)
             cluster_time : None
            f2py_compiler : /Users/chasejeon/Library/Python/2.7/bin/f2py
                    ninja : /Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/HEPTools/lib
               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 : /Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : None
      exrootanalysis_path : None
         fortran_compiler : None
                  nb_core : 4
                  collier : None
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                    golem : None
             cpp_compiler : None
      notification_center : True
                 run_mode : 2
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 2017-08-16 *
#* *
#* 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 group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model /Users/chasejeon/Documents/Research/MG5_aMC_v2_6_0/models\
/SM_HeavyN_NLO
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~
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > n1 mu+ [QCD]
output HNtoDilepDijet_M-40 -nojpeg
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.733000e+02 # mt
      23 9.118760e+01 # mz
      25 1.257000e+02 # mh
      9900012 4.000000e+01 # mn1
      9900014 5.000000e+02 # mn2
      9900016 1.000000e+03 # mn3
      1 0.000000e+00 # d : 0.0
      2 0.000000e+00 # u : 0.0
      3 0.000000e+00 # s : 0.0
      4 0.000000e+00 # c : 0.0
      5 0.000000e+00 # b : 0.0
      11 0.000000e+00 # e- : 0.0
      12 0.000000e+00 # ve : 0.0
      13 0.000000e+00 # mu- : 0.0
      14 0.000000e+00 # vm : 0.0
      15 0.000000e+00 # ta- : 0.0
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 7.995123e+01 # w+ : cmath.sqrt(mz__exp__2/2. + cmath.sqrt(mz__exp__4/4. - (aew*cmath.pi*mz__exp__2)/(gf*sqrt__2)))
      9000002 9.118760e+01 # ghz : mz
      9000003 7.995123e+01 # ghwp : mw
      9000004 7.995123e+01 # ghwm : mw
###################################
## INFORMATION FOR NUMIXING
###################################
BLOCK NUMIXING #
      1 1.000000e-01 # ven1
      2 0.000000e+00 # ven2
      3 0.000000e+00 # ven3
      4 1.000000e-01 # vmun1
      5 0.000000e+00 # vmun2
      6 0.000000e+00 # vmun3
      7 0.000000e+00 # vtan1
      8 0.000000e+00 # vtan2
      9 0.000000e+00 # vtan3
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.279400e+02 # aewm1
      2 1.174560e-05 # gf
      3 1.184000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      6 1.733000e+02 # ymt
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 9900012 # n1
      1 0 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9900014
###################################
BLOCK QNUMBERS 9900014 # n2
      1 0 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
###################################
## INFORMATION FOR QNUMBERS 9900016
###################################
BLOCK QNUMBERS 9900016 # n3
      1 0 # 3 times electric charge
      2 2 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 0 # particle/antiparticle distinction (0=own anti)
#
#*************************
# Decay widths *
#*************************
#
# PDG Width
DECAY 1 0.000000e+00
#
# PDG Width
DECAY 2 0.000000e+00
#
# PDG Width
DECAY 3 0.000000e+00
#
# PDG Width
DECAY 4 0.000000e+00
#
# PDG Width
DECAY 5 0.000000e+00
#
# PDG Width
DECAY 6 1.350000e+00
#
# PDG Width
DECAY 11 0.000000e+00
#
# PDG Width
DECAY 12 0.000000e+00
#
# PDG Width
DECAY 13 0.000000e+00
#
# PDG Width
DECAY 14 0.000000e+00
#
# PDG Width
DECAY 15 0.000000e+00
#
# PDG Width
DECAY 16 0.000000e+00
#
# PDG Width
DECAY 21 0.000000e+00
#
# PDG Width
DECAY 22 0.000000e+00
#
# PDG Width
DECAY 23 2.495200e+00
#
# PDG Width
DECAY 24 2.085000e+00
#
# PDG Width
DECAY 25 4.170000e-03
#
# PDG Width
DECAY 9000002 2.495200e+00
#
# PDG Width
DECAY 9000003 2.085000e+00
#
# PDG Width
DECAY 9000004 2.085000e+00
#
# PDG Width
DECAY 9900012 1.325116e-06
# BR NDA ID1 ID2 ...
   6.393778e-02 3 -13 -2 1 # 8.47249752825e-08
   6.393778e-02 3 -13 -4 3 # 8.47249752825e-08
   6.342462e-02 3 -11 -2 1 # 8.40449787559e-08
   6.342462e-02 3 -11 -4 3 # 8.40449787559e-08
   6.336425e-02 3 -1 2 13 # 8.3964981503e-08
   6.336425e-02 3 -3 4 13 # 8.3964981503e-08
   6.255300e-02 3 -1 2 11 # 8.2889981148e-08
   6.255300e-02 3 -3 4 11 # 8.2889981148e-08
   2.110456e-02 3 -13 11 14 # 2.7965990129e-08
   2.110456e-02 3 -15 11 16 # 2.7965990129e-08
   2.102080e-02 3 -13 -12 11 # 2.78549984128e-08
   2.102080e-02 3 -16 -13 15 # 2.78549984128e-08
   2.095099e-02 3 -11 12 13 # 2.77624920648e-08
   2.095099e-02 3 -15 13 16 # 2.77624920648e-08
   2.088534e-02 3 -14 -11 13 # 2.76754981994e-08
   2.088534e-02 3 -16 -11 15 # 2.76754981994e-08
   1.256267e-02 3 -11 11 12 # 1.66469950197e-08
   1.252645e-02 3 -12 -11 11 # 1.65989993182e-08
   1.251060e-02 3 -14 -13 13 # 1.65779962296e-08
   1.244193e-02 3 -13 13 14 # 1.64870005139e-08
   1.133410e-02 3 -1 1 14 # 1.50189972556e-08
   1.133410e-02 3 -3 3 14 # 1.50189972556e-08
   1.133410e-02 3 -5 5 14 # 1.50189972556e-08
   1.120682e-02 3 -12 -1 1 # 1.48503364911e-08
   1.120682e-02 3 -12 -3 3 # 1.48503364911e-08
   1.120682e-02 3 -12 -5 5 # 1.48503364911e-08
   1.114971e-02 3 -14 -1 1 # 1.47746591164e-08
   1.114971e-02 3 -14 -3 3 # 1.47746591164e-08
   1.114971e-02 3 -14 -5 5 # 1.47746591164e-08
   1.112053e-02 3 -1 1 12 # 1.47359922315e-08
   1.112053e-02 3 -3 3 12 # 1.47359922315e-08
   1.112053e-02 3 -5 5 12 # 1.47359922315e-08
   8.759986e-03 3 -14 -2 2 # 1.16079976084e-08
   8.759986e-03 3 -14 -4 4 # 1.16079976084e-08
   8.729045e-03 3 -12 -2 2 # 1.15669971942e-08
   8.729045e-03 3 -12 -4 4 # 1.15669971942e-08
   8.719989e-03 3 -2 2 14 # 1.15549969437e-08
   8.719989e-03 3 -4 4 14 # 1.15549969437e-08
   8.668296e-03 3 -2 2 12 # 1.14864977223e-08
   8.668296e-03 3 -4 4 12 # 1.14864977223e-08
   5.150491e-03 3 -14 -12 12 # 6.82499803196e-09
   5.150491e-03 3 -16 -14 16 # 6.82499803196e-09
   5.103326e-03 3 -12 12 14 # 6.76249893582e-09
   5.103326e-03 3 -16 14 16 # 6.76249893582e-09
   5.093893e-03 3 -14 12 14 # 6.74999911659e-09
   5.093893e-03 3 -16 12 16 # 6.74999911659e-09
   5.065593e-03 3 -14 -12 14 # 6.71249833379e-09
   5.065593e-03 3 -16 -12 16 # 6.71249833379e-09
   2.577208e-03 3 -14 -14 14 # 3.41509955613e-09
   2.551851e-03 3 -13 12 13 # 3.38149858972e-09
   2.551851e-03 3 -15 12 15 # 3.38149858972e-09
   2.546456e-03 3 -11 11 14 # 3.3743495889e-09
   2.546456e-03 3 -15 14 15 # 3.3743495889e-09
   2.538645e-03 3 -14 -11 11 # 3.36399910782e-09
   2.538645e-03 3 -15 -14 15 # 3.36399910782e-09
   2.530495e-03 3 -14 14 14 # 3.35319941242e-09
   2.524835e-03 3 -13 -12 13 # 3.34569925586e-09
   2.524835e-03 3 -15 -12 15 # 3.34569925586e-09
   2.518043e-03 3 -12 12 12 # 3.33669906799e-09
   2.511629e-03 3 -12 -12 12 # 3.32819977396e-09
#
# PDG Width
DECAY 9900014 1.500000e+00
#
# PDG Width
DECAY 9900016 1.230000e+01
#***********************************************************************
# 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 *
#***********************************************************************
#
#*******************
# 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 *
#***********************************************************************
  100 = nevents ! Number of unweighted events requested
  0.001 = 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 ! 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 *
#***********************************************************************
  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
  25000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#***********************************************************************
# 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 of all final state particles and partons. This *
# can be changed in SubProcesses/set_scales.f) *
#***********************************************************************
  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 = muF1_ref_fixed ! fixed fact reference scale for pdf1
  91.118 = muF2_ref_fixed ! fixed fact reference scale for pdf2
  -1 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
#***********************************************************************
# Renormalization and factorization scales (advanced and NLO options) *
#***********************************************************************
  False = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
  91.118 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
  1.0 = muR_over_ref ! ratio of current muR over reference muR
  1.0 = muF1_over_ref ! ratio of current muF1 over reference muF1
  1.0 = muF2_over_ref ! ratio of current muF2 over reference muF2
  1.0 = QES_over_ref ! ratio of current QES over reference QES
#***********************************************************************
# Reweight flags to get scale dependence and PDF uncertainty *
# For scale dependence: factor rw_scale_up/down around central scale *
# For PDF uncertainty: use LHAPDF with supported set *
#***********************************************************************
  True = reweight_scale ! reweight to get scale dependence
  0.5 = rw_Rscale_down ! lower bound for ren scale variations
  2.0 = rw_Rscale_up ! upper bound for ren scale variations
  0.5 = rw_Fscale_down ! lower bound for fact scale variations
  2.0 = rw_Fscale_up ! upper bound for fact scale variations
  True = reweight_PDF ! reweight to get PDF uncertainty
  292201 = PDF_set_min ! First of the error PDF sets
  292302 = PDF_set_max ! Last of the error PDF sets
#***********************************************************************
# 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) *
#***********************************************************************
  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
  15.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 gen cuts can be specified in SubProcesses/cuts.f) *
#***********************************************************************
  4.0 = ptl ! Min lepton transverse momentum
  3.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
  4.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 *
#***********************************************************************
  0.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)
#***********************************************************************
# Maximal PDG code for quark to be considered a jet when applying cuts.*
# At least all massless quarks of the model should be included here. *
#***********************************************************************
  5 = maxjetflavor
 #***********************************************************************
# For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
#***********************************************************************
  0 = iappl ! aMCfast switch (0=OFF, 1=prepare APPLgrids, 2=fill grids)
#***********************************************************************

Question information

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

Hi,

Looks like that you are not running in "onshell" mode this time. Is this on purpose?
if it is, could you copy/paste your madspin_card? I will try to reproduce the problem.

Cheers,

Olivier

Revision history for this message
Sihyun Jeon (shjeon) said : 		#2

Ah yes sorry..
I made a mistake and didn't have a problem then.

It now works great!
Thank you for all the helps

Revision history for this message
ChristopherClarkMucend (christopherclarkmucend) said : 		#3

Cool topic on programming. As for me, such resources are very helpful in doing homework. Everyone knows that there are many resources to help with homework in humanitarian subjects, such as https://cheetahpapers.com/research-paper/pay/ .But few resources that help to do homework in programming. This resource is great for helping you not only do your programming homework and generally develop your knowledge.