Warning: m_2^2 is negative in t-channel branching

Hi,

Did this happen at the MadSpin stage?
If yes, are you using the latest version of the code? (2.5.3)

If not, do you have a debug file created?

Cheers,

Olivier

> On 9 Mar 2017, at 16:38, Marc de Beurs <email address hidden> wrote:
>
> New question #545042 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/545042
>
> Dear experts,
>
> I am trying to run the following (single top) process:
>
> generate p p > p w+ b QED=3 QCD=0 EFT=0
>
> Where I use a model containing some EFT operators (turned off at this stage).
>
> Compilation and everything works fine, however during the generation step I receive the following errors:
>
> rm: cannot remove `ftn26': No such file or directory
> Command "generate_events run_01" interrupted with error:
> ValueError : could not convert string to float: Warning: m_2^2 is negative in t-channel branching 0.0000000000000000
>
> I have tried many things but I am not able to figure out what the problem is.
> At the process where I demand an intermediate top, the problem does not occur:
>
> generate p p > t > p w+ b QED=3 QCD=0 EFT=0
>
> My thought is that something in the parameter card might be set wrong (top width or something). Below you can find the full param_card that I am using.
>
> Thanks in advance!
> Kind regards, Marc
>
> ######################################################################
> ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
> ######################################################################
> ###################################
> ## INFORMATION FOR DIM6
> ###################################
> BLOCK DIM6 #
> 1 1.000000e+03 # lambda
> 3 1.000000e-10 # ctg
> 7 1.000000e-10 # cfq3
> 8 1.000000e-10 # ctw
> 12 1.725000e+02 # muprime
> ###################################
> ## INFORMATION FOR EFT
> ###################################
> BLOCK EFT #
> 1 -1.000000e-10 # c4f
> ###################################
> ## INFORMATION FOR LOOP
> ###################################
> BLOCK LOOP #
> 1 9.118800e+01 # mu_r
> ###################################
> ## INFORMATION FOR MASS
> ###################################
> BLOCK MASS #
> 6 1.725000e+02 # mt
> 15 1.777000e+00 # mta
> 23 9.118760e+01 # mz
> 25 1.250000e+02 # mh
> 34 1.000000e+08 # mwp
> 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
> 16 0.000000e+00 # vt : 0.0
> 21 0.000000e+00 # g : 0.0
> 22 0.000000e+00 # a : 0.0
> 24 7.982436e+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.982436e+01 # ghwp : mw
> 9000004 7.982436e+01 # ghwm : mw
> ###################################
> ## INFORMATION FOR SMINPUTS
> ###################################
> BLOCK SMINPUTS #
> 1 1.279000e+02 # aewm1
> 2 1.166370e-05 # gf
> 3 1.184000e-01 # as
> ###################################
> ## INFORMATION FOR YUKAWA
> ###################################
> BLOCK YUKAWA #
> 6 1.725000e+02 # ymt
> 15 1.777000e+00 # ymtau
> ###################################
> ## INFORMATION FOR DECAY
> ###################################
> BLOCK QNUMBERS 34 # wp+
> 1 3 # 3 times electric charge
> 2 3 # number of spin states (2s+1)
> 3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
> 4 1 # 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.485822e+00
> # BR NDA ID1 ID2 ...
> 1.000000e+00 2 24 5 # 1.48582192013
> #
> # 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.415902e+00
> # BR NDA ID1 ID2 ...
> 1.518058e-01 2 5 -5 # 0.366747999839
> 1.518058e-01 2 3 -3 # 0.366747999834
> 1.518058e-01 2 1 -1 # 0.366747999834
> 1.176194e-01 2 4 -4 # 0.284156909669
> 1.176194e-01 2 2 -2 # 0.284156909669
> 6.866355e-02 2 16 -16 # 0.1658844066
> 6.866355e-02 2 14 -14 # 0.1658844066
> 6.866355e-02 2 12 -12 # 0.1658844066
> 3.447711e-02 2 13 -13 # 0.0832933164354
> 3.447711e-02 2 11 -11 # 0.0832933164354
> 3.439891e-02 2 15 -15 # 0.0831044017136
> #
> # PDG Width
> DECAY 24 2.002840e+00
> # BR NDA ID1 ID2 ...
> 3.333609e-01 2 4 -3 # 0.667668389896
> 3.333609e-01 2 2 -1 # 0.667668389896
> 1.111203e-01 2 14 -13 # 0.222556129965
> 1.111203e-01 2 12 -11 # 0.222556129965
> 1.110377e-01 2 16 -15 # 0.222390692164
> #
> # PDG Width
> DECAY 25 1.166966e-03
> # BR NDA ID1 ID2 ...
> 2.217094e-01 2 15 -15 # 0.000258727280256
> 1.181869e-01 3 -2 1 24 # 0.000137919982104
> 1.181869e-01 3 -4 3 24 # 0.000137919982104
> 1.173428e-01 3 -24 -1 2 # 0.000136934997667
> 1.173428e-01 3 -24 -3 4 # 0.000136934997667
> 3.977795e-02 3 -24 -11 12 # 4.64194964053e-05
> 3.977795e-02 3 -24 -13 14 # 4.64194964053e-05
> 3.922138e-02 3 -16 15 24 # 4.57699970322e-05
> 3.917853e-02 3 -12 11 24 # 4.57199985137e-05
> 3.917853e-02 3 -14 13 24 # 4.57199985137e-05
> 3.892231e-02 3 -24 -15 16 # 4.5420997383e-05
> 1.079980e-02 3 -1 1 23 # 1.26029962047e-05
> 1.079980e-02 3 -3 3 23 # 1.26029962047e-05
> 1.077667e-02 3 -5 5 23 # 1.25760033624e-05
> 8.459116e-03 3 -2 2 23 # 9.87149676457e-06
> 8.459116e-03 3 -4 4 23 # 9.87149676457e-06
> 4.867896e-03 3 -12 12 23 # 5.68066636294e-06
> 4.867896e-03 3 -14 14 23 # 5.68066636294e-06
> 4.867896e-03 3 -16 16 23 # 5.68066636294e-06
> 2.458127e-03 3 -11 11 23 # 2.86854987486e-06
> 2.458127e-03 3 -13 13 23 # 2.86854987486e-06
> 2.360138e-03 3 -15 15 23 # 2.75419994817e-06
> #
> # PDG Width
> DECAY 34 0.000000e+00
> #
> # PDG Width
> DECAY 9000002 2.495200e+00
> #
> # PDG Width
> DECAY 9000003 2.085000e+00
> #
> # PDG Width
> DECAY 9000004 2.085000e+00
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi Olivier,

Thanks for your quick response!
The first warnings:

rm: cannot remove `ftn26': No such file or directory
rm: cannot remove `ftn26': No such file or directory
rm: cannot remove `ftn26': No such file or directory
rm: cannot remove `ftn26': No such file or directory

Happen during the madgraph even generation phase.

The final error:

Command "generate_events run_01" interrupted with error:
ValueError : could not convert string to float: Warning: m_2^2 is negative in t-channel branching 0.0000000000000000

Happens indeed during madspin.
I do not know if both problems are connected.

I am using the latest update of madgraph 2.5.3
Not sure how to check which version of madspin I am using, but I have installed through madgraph.

Here is my debug.log file:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.5.3 2017-03-09 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1414, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1369, in onecmd_orig
    return func(arg, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/madevent_interface.py", line 2119, in do_generate_events
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1441, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1369, in onecmd_orig
    return func(arg, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/common_run_interface.py", line 3388, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1556, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1441, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/interface/extended_cmd.py", line 1369, in onecmd_orig
    return func(arg, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadSpin/interface_madspin.py", line 574, in do_launch
    generate_all.run()
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadSpin/decay.py", line 2101, in run
    self.get_max_weight_from_event(decay_mapping)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadSpin/decay.py", line 3164, in get_max_weight_from_event
    weight = self.get_max_weight_from_fortran(decay['path'], event_map,numberps,self.options['BW_cut'])
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadSpin/decay.py", line 3303, in get_max_weight_from_fortran
    path, std_in)
  File "/data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadSpin/decay.py", line 3337, in loadfortran
    maxweight=float(external.stdout.readline())
ValueError: could not convert string to float: Warning: m_2^2 is negative in t-channel branching 0.0000000000000000

                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : None
             cluster_memory : None (user set)
               cluster_size : 100
              cluster_queue : None (user set)
                    nb_core : 32 (user set)
               cluster_time : 32 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/pythia-pgs (user set)
                    td_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/td (user set)
               delphes_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/Delphes (user set)
                thepeg_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/HEPTools/herwig++ (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/MadAnalysis (user set)
                     lhapdf : lhapdf-config
              f2py_compiler : None
                 hepmc_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/HEPTools/hepmc (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : /data/atlas/users/mdebeurs/hep/MG5_aMC_v2_4_3/ExRootAnalysis (user set)
                    timeout : 0 (user set)
               cpp_compiler : None
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.5.3 2017-03-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 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 TEFT_EW_bR_Wp-singletop
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
set complex_mass_scheme True
generate p p > p w+ b QED=3 QCD=0 EFT=0
output /data/atlas/users/mdebeurs/MG5/Test//halfchain-Test-SM-LO-100-T\
EFT_EW_bR_Wp-singletop/
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR DIM6
###################################
BLOCK DIM6 #
      1 1.000000e+03 # lambda
      3 1.000000e-10 # ctg
      7 1.000000e-10 # cfq3
      8 1.000000e-10 # ctw
      12 1.725000e+02 # muprime
###################################
## INFORMATION FOR EFT
###################################
BLOCK EFT #
      1 -1.000000e-10 # c4f
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.725000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118760e+01 # mz
      25 1.250000e+02 # mh
      34 1.000000e+08 # mwp
      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
      16 0.000000e+00 # vt : 0.0
      21 0.000000e+00 # g : 0.0
      22 0.000000e+00 # a : 0.0
      24 7.982436e+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.982436e+01 # ghwp : mw
      9000004 7.982436e+01 # ghwm : mw
###################################
## INFORMATION FOR SMINPUTS
###################################
BLOCK SMINPUTS #
      1 1.279000e+02 # aewm1
      2 1.166370e-05 # gf
      3 1.184000e-01 # as
###################################
## INFORMATION FOR YUKAWA
###################################
BLOCK YUKAWA #
      6 1.725000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
BLOCK QNUMBERS 34 # wp+
      1 3 # 3 times electric charge
      2 3 # number of spin states (2s+1)
      3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
      4 1 # 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.485822e+00
# BR NDA ID1 ID2 ...
   1.000000e+00 2 24 5 # 1.48582192013
#
# 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.415902e+00
# BR NDA ID1 ID2 ...
   1.518058e-01 2 5 -5 # 0.366747999839
   1.518058e-01 2 3 -3 # 0.366747999834
   1.518058e-01 2 1 -1 # 0.366747999834
   1.176194e-01 2 4 -4 # 0.284156909669
   1.176194e-01 2 2 -2 # 0.284156909669
   6.866355e-02 2 16 -16 # 0.1658844066
   6.866355e-02 2 14 -14 # 0.1658844066
   6.866355e-02 2 12 -12 # 0.1658844066
   3.447711e-02 2 13 -13 # 0.0832933164354
   3.447711e-02 2 11 -11 # 0.0832933164354
   3.439891e-02 2 15 -15 # 0.0831044017136
#
# PDG Width
DECAY 24 2.002840e+00
# BR NDA ID1 ID2 ...
   3.333609e-01 2 4 -3 # 0.667668389896
   3.333609e-01 2 2 -1 # 0.667668389896
   1.111203e-01 2 14 -13 # 0.222556129965
   1.111203e-01 2 12 -11 # 0.222556129965
   1.110377e-01 2 16 -15 # 0.222390692164
#
# PDG Width
DECAY 25 1.166966e-03
# BR NDA ID1 ID2 ...
   2.217094e-01 2 15 -15 # 0.000258727280256
   1.181869e-01 3 -2 1 24 # 0.000137919982104
   1.181869e-01 3 -4 3 24 # 0.000137919982104
   1.173428e-01 3 -24 -1 2 # 0.000136934997667
   1.173428e-01 3 -24 -3 4 # 0.000136934997667
   3.977795e-02 3 -24 -11 12 # 4.64194964053e-05
   3.977795e-02 3 -24 -13 14 # 4.64194964053e-05
   3.922138e-02 3 -16 15 24 # 4.57699970322e-05
   3.917853e-02 3 -12 11 24 # 4.57199985137e-05
   3.917853e-02 3 -14 13 24 # 4.57199985137e-05
   3.892231e-02 3 -24 -15 16 # 4.5420997383e-05
   1.079980e-02 3 -1 1 23 # 1.26029962047e-05
   1.079980e-02 3 -3 3 23 # 1.26029962047e-05
   1.077667e-02 3 -5 5 23 # 1.25760033624e-05
   8.459116e-03 3 -2 2 23 # 9.87149676457e-06
   8.459116e-03 3 -4 4 23 # 9.87149676457e-06
   4.867896e-03 3 -12 12 23 # 5.68066636294e-06
   4.867896e-03 3 -14 14 23 # 5.68066636294e-06
   4.867896e-03 3 -16 16 23 # 5.68066636294e-06
   2.458127e-03 3 -11 11 23 # 2.86854987486e-06
   2.458127e-03 3 -13 13 23 # 2.86854987486e-06
   2.360138e-03 3 -15 15 23 # 2.75419994817e-06
#
# PDG Width
DECAY 34 0.000000e+00
#
# PDG Width
DECAY 9000002 2.495200e+00
#
# PDG Width
DECAY 9000003 2.085000e+00
#
# PDG Width
DECAY 9000004 2.085000e+00
#*********************************************************************
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/conventions: *
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
#*********************************************************************
#
#*******************
# Running parameters
#*******************
#
#*********************************************************************
# Tag name for the run (one word) *
#*********************************************************************
  SM = run_tag ! name of the run
#*********************************************************************
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#*********************************************************************
  100 = nevents ! Number of unweighted events requested
  0 = iseed ! rnd seed (0=assigned automatically=default))
#*********************************************************************
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#*********************************************************************
     1 = lpp1 ! beam 1 type
     1 = lpp2 ! beam 2 type
     6500.0 = ebeam1 ! beam 1 total energy in GeV
     6500.0 = ebeam2 ! beam 2 total energy in GeV
#*********************************************************************
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#*********************************************************************
     0.0 = polbeam1 ! beam polarization for beam 1
     0.0 = polbeam2 ! beam polarization for beam 2
#*********************************************************************
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#*********************************************************************
     nn23lo1 = pdlabel ! PDF set
     230000 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#*********************************************************************
# Renormalization and factorization scales *
#*********************************************************************
 True = fixed_ren_scale ! if .true. use fixed ren scale
 True = fixed_fac_scale ! if .true. use fixed fac scale
 172.5 = scale ! fixed ren scale
 172.5 = dsqrt_q2fact1 ! fixed fact scale for pdf1
 172.5 = dsqrt_q2fact2 ! fixed fact scale for pdf2
 0 = dynamical_scale_choice ! Choose one of the preselected dynamical choices
 1.0 = scalefact ! scale factor for event-by-event scales
#*********************************************************************
# Type and output format
#*********************************************************************
  False = gridpack !True = setting up the grid pack
  -1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
  3.0 = lhe_version ! Change the way clustering information pass to shower.
  True = clusinfo ! include clustering tag in output
  average = event_norm ! average/sum. Normalization of the weight in the LHEF

#*********************************************************************
# Matching parameter (MLM only)
#*********************************************************************
 0 = ickkw ! 0 no matching, 1 MLM
 1.0 = alpsfact ! scale factor for QCD emission vx
 False = chcluster ! cluster only according to channel diag
 5 = asrwgtflavor ! highest quark flavor for a_s reweight
 False = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
                                   ! (turn off for VBF and single top processes)
 0.0 = xqcut ! minimum kt jet measure between partons
#*********************************************************************
#
#*********************************************************************
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#*********************************************************************
   0 = nhel ! using helicities importance sampling or not.
#*********************************************************************
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.ucl.ac.be/projects/madgraph/wiki/LOEventGenerationBias' *
#*********************************************************************
 None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
 {} = bias_parameters ! Specifies the parameters of the module.
#
#*******************************
# Parton level cuts definition *
#*******************************
#
#
#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  15.0 = bwcutoff ! (M+/-bwcutoff*Gamma)
#*********************************************************************
# Apply pt/E/eta/dr/mij/kt_durham cuts on decay products or not
# (note that etmiss/ptll/ptheavy/ht/sorted cuts always apply)
#*********************************************************************
   False = cut_decays ! Cut decay products
#*********************************************************************
# Standard Cuts *
#*********************************************************************
# Minimum and maximum pt's (for max, -1 means no cut) *
#*********************************************************************
 0.0 = ptj ! minimum pt for the jets
 0.0 = ptb ! minimum pt for the b
 0.0 = pta ! minimum pt for the photons
 0.0 = ptl ! minimum pt for the charged leptons
 0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
 0.0 = ptheavy ! minimum pt for one heavy final state
 -1.0 = ptjmax ! maximum pt for the jets
 -1.0 = ptbmax ! maximum pt for the b
 -1.0 = ptamax ! maximum pt for the photons
 -1.0 = ptlmax ! maximum pt for the charged leptons
 -1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#*********************************************************************
# Minimum and maximum E's (in the center of mass frame) *
#*********************************************************************
  0.0 = ej ! minimum E for the jets
  0.0 = eb ! minimum E for the b
  0.0 = ea ! minimum E for the photons
  0.0 = el ! minimum E for the charged leptons
  -1.0 = ejmax ! maximum E for the jets
 -1.0 = ebmax ! maximum E for the b
 -1.0 = eamax ! maximum E for the photons
 -1.0 = elmax ! maximum E for the charged leptons
#*********************************************************************
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#*********************************************************************
  -1.0 = etaj ! max rap for the jets
  -1.0 = etab ! max rap for the b
 -1.0 = etaa ! max rap for the photons
 -1.0 = etal ! max rap for the charged leptons
 0.0 = etajmin ! min rap for the jets
 0.0 = etabmin ! min rap for the b
 0.0 = etaamin ! min rap for the photons
 0.0 = etalmin ! main rap for the charged leptons
#*********************************************************************
# Minimum and maximum DeltaR distance *
#*********************************************************************
 0.0 = drjj ! min distance between jets
 0.0 = drbb ! min distance between b's
 0.0 = drll ! min distance between leptons
 0.0 = draa ! min distance between gammas
 0.0 = drbj ! min distance between b and jet
 0.0 = draj ! min distance between gamma and jet
 0.0 = drjl ! min distance between jet and lepton
 0.0 = drab ! min distance between gamma and b
 0.0 = drbl ! min distance between b and lepton
 0.0 = dral ! min distance between gamma and lepton
 -1.0 = drjjmax ! max distance between jets
 -1.0 = drbbmax ! max distance between b's
 -1.0 = drllmax ! max distance between leptons
 -1.0 = draamax ! max distance between gammas
 -1.0 = drbjmax ! max distance between b and jet
 -1.0 = drajmax ! max distance between gamma and jet
 -1.0 = drjlmax ! max distance between jet and lepton
 -1.0 = drabmax ! max distance between gamma and b
 -1.0 = drblmax ! max distance between b and lepton
 -1.0 = dralmax ! maxdistance between gamma and lepton
#*********************************************************************
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#*********************************************************************
 0.0 = mmjj ! min invariant mass of a jet pair
 0.0 = mmbb ! min invariant mass of a b pair
 0.0 = mmaa ! min invariant mass of gamma gamma pair
 0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 -1.0 = mmjjmax ! max invariant mass of a jet pair
 -1.0 = mmbbmax ! max invariant mass of a b pair
 -1.0 = mmaamax ! max invariant mass of gamma gamma pair
 -1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
 -1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#*********************************************************************
# Minimum and maximum pt for 4-momenta sum of leptons *
#*********************************************************************
 0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
 -1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#*********************************************************************
# Inclusive cuts *
#*********************************************************************
 0.0 = xptj ! minimum pt for at least one jet
 0.0 = xptb ! minimum pt for at least one b
 0.0 = xpta ! minimum pt for at least one photon
 0.0 = xptl ! minimum pt for at least one charged lepton
#*********************************************************************
# Control the pt's of the jets sorted by pt *
#*********************************************************************
 0.0 = ptj1min ! minimum pt for the leading jet in pt
 0.0 = ptj2min ! minimum pt for the second jet in pt
 0.0 = ptj3min ! minimum pt for the third jet in pt
 0.0 = ptj4min ! minimum pt for the fourth jet in pt
 -1.0 = ptj1max ! maximum pt for the leading jet in pt
 -1.0 = ptj2max ! maximum pt for the second jet in pt
 -1.0 = ptj3max ! maximum pt for the third jet in pt
 -1.0 = ptj4max ! maximum pt for the fourth jet in pt
 0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#*********************************************************************
# Control the pt's of leptons sorted by pt *
#*********************************************************************
 0.0 = ptl1min ! minimum pt for the leading lepton in pt
 0.0 = ptl2min ! minimum pt for the second lepton in pt
 0.0 = ptl3min ! minimum pt for the third lepton in pt
 0.0 = ptl4min ! minimum pt for the fourth lepton in pt
 -1.0 = ptl1max ! maximum pt for the leading lepton in pt
 -1.0 = ptl2max ! maximum pt for the second lepton in pt
 -1.0 = ptl3max ! maximum pt for the third lepton in pt
 -1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#*********************************************************************
# Control the Ht(k)=Sum of k leading jets *
#*********************************************************************
 0.0 = htjmin ! minimum jet HT=Sum(jet pt)
 -1.0 = htjmax ! maximum jet HT=Sum(jet pt)
 0.0 = ihtmin !inclusive Ht for all partons (including b)
 -1.0 = ihtmax !inclusive Ht for all partons (including b)
 0.0 = ht2min ! minimum Ht for the two leading jets
 0.0 = ht3min ! minimum Ht for the three leading jets
 0.0 = ht4min ! minimum Ht for the four leading jets
 -1.0 = ht2max ! maximum Ht for the two leading jets
 -1.0 = ht3max ! maximum Ht for the three leading jets
 -1.0 = ht4max ! maximum Ht for the four leading jets
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#***********************************************************************
 0.0 = ptgmin ! Min photon transverse momentum
 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)
#*********************************************************************
# WBF cuts *
#*********************************************************************
 0.0 = xetamin ! minimum rapidity for two jets in the WBF case
 0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#***********************************************************************
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#***********************************************************************
 -1.0 = ktdurham
 0.4 = dparameter
 -1.0 = ptlund
 1, 2, 3, 4, 5, 6, 21, 82 = pdgs_for_merging_cut ! PDGs for two cuts above
#*********************************************************************
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#*********************************************************************
 5 = maxjetflavor ! Maximum jet pdg code
#*********************************************************************
#
#*********************************************************************
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#*********************************************************************
   True = use_syst ! Enable systematics studies
#
#**************************************
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#**************************************
#
0.5 1 2 = sys_scalefact # factorization/renormalization scale factor
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_lo_as_0130_qed = sys_pdf # list of pdf sets
# MSTW2008nlo68cl.LHgrid 1 = sys_pdf
#

Hi,

The following fix should solve your problem.

Thanks,

Olivier

=== modified file 'MadSpin/src/driver.f'
--- MadSpin/src/driver.f 2016-11-04 15:29:19 +0000
+++ MadSpin/src/driver.f 2017-03-10 14:14:14 +0000
@@ -904,7 +904,7 @@
              else
                 ! might be negative because of numerical unstabilities
                 index_p2=itree(2,i-1)
- if (index_p2.gt.0) then
+ if (index_p2.ge.0) then
                    m2_tchan(i)=m(index_p2)
                 else
                    if (m(index_p2).gt.5d-2)then

 Hi Olivier,

Thanks for your response. I have changed the line what you suggested in the madspin driver.f.
But the same warning and error still occurs...

Do I maybe need to recompile (or make) madspin after this change?
Also I believe that this will not solve the warning that I get during the madgraph generation of events:

rm: cannot remove `ftn26': No such file or directory

Or doesn't this matter for the results?

Cheers Marc

On Friday, March 10, 2017 15:59 CET, Olivier Mattelaer <email address hidden> wrote:

> Your question #545042 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/545042
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> The following fix should solve your problem.
>
> Thanks,
>
>
> Olivier
>
>
> === modified file 'MadSpin/src/driver.f'
> --- MadSpin/src/driver.f 2016-11-04 15:29:19 +0000
> +++ MadSpin/src/driver.f 2017-03-10 14:14:14 +0000
> @@ -904,7 +904,7 @@
> else
> ! might be negative because of numerical unstabilities
> index_p2=itree(2,i-1)
> - if (index_p2.gt.0) then
> + if (index_p2.ge.0) then
> m2_tchan(i)=m(index_p2)
> else
> if (m(index_p2).gt.5d-2)then
>
> -- If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/mg5amcnlo/+question/545042/+confirm?answer_id=2
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/545042
>
> You received this question notification because you asked the question.

Dear Olivier,

Yes indeed, I had turned all the cuts off. Setting a cut on the jet pt solved my problem.
Thanks so much for your help!

Cheers Marc