MadGraph5_aMC@NLO

Single top t channel with aMC@NLO

Asked by Benedikt Maier

Dear experts,

sorry for opening another thread, but I'm not sure if it is related to the one I posted recently (s channel single top, which I'm still working on).

I just checked out the latest official

bzr branch lp:madgraph5

and am generating

generate p p > t~ b j $$ w+ w- [QCD]

I am modifying the setscales.f in the sense that I use a certain dynamical scale, a la:

tmp = 4.7**2 + pt(pp(0,nexternal-1))**2
        tmp=4*sqrt(tmp)
        temp_scale_id='4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mur2'

With the MG5 version 2_0_0 this worked and I generated the process successfully.

However, now the calculation crashes with the following error (I put the logfile at the end)

##########################################################

 INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Generating events without running the shower.
INFO: Setting up grid
INFO: Idle: 0, Running: 8, Completed: 0 [ current time: 14h05 ]
INFO: Idle: 0, Running: 7, Completed: 1 [ 46.5s ]
INFO: Idle: 0, Running: 6, Completed: 2 [ 47.6s ]
INFO: Idle: 0, Running: 5, Completed: 3 [ 53.4s ]
INFO: Idle: 0, Running: 4, Completed: 4 [ 1m 32s ]
INFO: Idle: 0, Running: 3, Completed: 5 [ 1m 33s ]
INFO: Idle: 0, Running: 2, Completed: 6 [ 1m 52s ]
INFO: Idle: 0, Running: 1, Completed: 7 [ 4m 59s ]
INFO: Idle: 0, Running: 0, Completed: 8 [ 4m 59s ]
INFO: Determining the number of unweighted events per channel

      Intermediate results:
      Random seed: 15
      Total cross-section: 6.168e+01 +- 1.3e+00 pb
      Total abs(cross-section): 1.899e+02 +- 1.8e+00 pb

INFO: Computing upper envelope
INFO: Idle: 0, Running: 8, Completed: 0 [ current time: 14h10 ]
INFO: Idle: 0, Running: 7, Completed: 1 [ 1m 36s ]
INFO: Idle: 0, Running: 6, Completed: 2 [ 1m 45s ]
INFO: Idle: 0, Running: 5, Completed: 3 [ 1m 47s ]
INFO: Idle: 0, Running: 4, Completed: 4 [ 4m 39s ]
/afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/ajob1: line 13: 7225 Floating point exception../madevent_mintMC > log.txt < input_app.txt 2>&1
WARNING: program /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/ajob1 2 F 1 launch ends with non zero status: 136. Stop all computation
WARNING: Last 15 lines of logfile /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/*/log.txt:
 ERROR in OneLOop dilog2_r: r1,r2 = .7036300903284785162E-16, .7964685568044487834E-15, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .5231995004828015039E-15, .1353689966652817739E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .4484567146995441744E-15, .1160305685702415262E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .4484567146995441744E-15, .1160305685702415262E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .5979422862660588335E-15, .1547074247603220217E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .7512099217322685834E-15, .1434639658806850089E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .6009679373858149062E-15, .1147711727045479992E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .6009679373858149062E-15, .1147711727045479992E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .4549226662644539214E-15, .9631071934392942122E-15, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .6065635550192718623E-15, .1284142924585725222E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .6065635550192718623E-15, .1284142924585725222E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .1454738538845453516E-16, .1351023166478349429E-14, returning 0
 ERROR in OneLOop dilog2_r: r1,r2 = .7472302259718281443E-16, .1435535743846529388E-14, returning 0
Time in seconds: 1672

INFO: remove job currently running
INFO: remove job currently running
INFO: remove job currently running
Command "launch auto " interrupted with error:
Exception : program /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/ajob1 2 F 1 launch ends with non zero status: 136. Stop all computation
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/run_01_tag_1_debug.log'.
Please attach this file to your report.
INFO:

quit
INFO:

MG5_aMC>exit

[1]+ Done emacs /afs/cern.ch/work/b/bmaier/public/test/SubProcesses/setscales.f
[lxplus418] /afs/cern.ch/user/b/bmaier/private/madgraph5 > cat /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/run_01_tag_1_debug.log
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.0.1 2014-01-20 *
#* *
#* 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 "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/extended_cmd.py", line 841, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/extended_cmd.py", line 834, in onecmd_orig
    return func(arg, **opt)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/amcatnlo_run_interface.py", line 1031, in do_launch
    evt_file = self.run(mode, options)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/amcatnlo_run_interface.py", line 1344, in run
    self.run_all(job_dict, [['2', 'F', '%d' % i]], status, split_jobs = split)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/amcatnlo_run_interface.py", line 2624, in run_all
    self.wait_for_complete(run_type)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/interface/amcatnlo_run_interface.py", line 2592, in wait_for_complete
    self.cluster.wait(self.me_dir, update_status)
  File "/afs/cern.ch/user/b/bmaier/private/madgraph5/madgraph/various/cluster.py", line 544, in wait
    raise Exception, msg
Exception: program /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/ajob1 2 F 1 launch ends with non zero status: 136. Stop all computation
Value of current Options:
              text_editor : None
              web_browser : None
        cluster_temp_path : None
       group_subprocesses : False
                  timeout : 60
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : /afs/cern.ch/user/b/bmaier/private/LHAPDF_589/bin/lhapdf-config
ignore_six_quark_processes : False
                 mg5_path : /afs/cern.ch/user/b/bmaier/private/madgraph5
    loop_optimized_output : True
    cluster_status_update : (600, 30)
               hepmc_path : None
             pythia8_path : None
                hwpp_path : None
   automatic_html_opening : False
       cluster_retry_wait : 300
             stdout_level : 10
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
      exrootanalysis_path : None
         fortran_compiler : None
                      OLP : MadLoop
              auto_update : 7
         cluster_nb_retry : 1
               eps_viewer : None
             syscalc_path : None
                  fastjet : /afs/cern.ch/user/b/bmaier/private/fastjet-3.0.1/fastjet-config
                    gauge : unitary
      complex_mass_scheme : False
                  nb_core : 16
                 run_mode : 2
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.0.1 2014-01-20 *
#* *
#* 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 gauge unitary
set complex_mass_scheme False
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~
set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set gauge unitary
set complex_mass_scheme False
import model loop_sm
define j = g u c d s u~ c~ d~ s~
define p = j
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > t~ b j $$ w+ w- [QCD]
output /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS
######################################################################
## 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. By default, ##
## this is only 1->2 decay modes. ##
## ##
######################################################################

###################################
## 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)))
  82 0.000000 # gh : 0.0

###################################
## 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
DECAY 82 0.000000 # gh : 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 *
#***********************************************************************
#
#*******************
# Running parameters
#*******************
#
#***********************************************************************
# Tag name for the run (one word) *
#***********************************************************************
  tag_1 = run_tag ! name of the run
#***********************************************************************
# Number of events (and their normalization) and the required *
# (relative) accuracy on the Xsec. *
# These values are ignored for fixed order runs *
#***********************************************************************
500000 = nevents ! Number of unweighted events requested
    -1 = req_acc ! Required accuracy (-1=auto determined from nevents)
    -1 = nevt_job! Max number of events per job in event generation.
                 ! (-1= no split).
average = event_norm ! Normalize events to sum or average to the X sect.
#***********************************************************************
# 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
 500 = npoints_FO_grid ! number of points to setup grids
 4 = niters_FO_grid ! number of iter. to setup grids
 20 = npoints_FO ! number of points to compute Xsec
 6 = niters_FO ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
     15 = 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 = ebeam1 ! beam 1 energy in GeV
 6500 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 lhapdf = pdlabel ! PDF set
 10981 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
# cteq6_m = pdlabel ! PDF set
# 21100 = 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
#***********************************************************************
# 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) *
#***********************************************************************
 F = fixed_ren_scale ! if .true. use fixed ren scale
 F = fixed_fac_scale ! if .true. use fixed fac scale
 172.5 = muR_ref_fixed ! fixed ren reference scale
 172.5 = muF1_ref_fixed ! fixed fact reference scale for pdf1
 172.5 = muF2_ref_fixed ! fixed fact reference scale for pdf2
#***********************************************************************
# Renormalization and factorization scales (advanced and NLO options) *
#***********************************************************************
 F = fixed_QES_scale ! if .true. use fixed Ellis-Sexton scale
 172.5 = QES_ref_fixed ! fixed Ellis-Sexton reference scale
 1 = muR_over_ref ! ratio of current muR over reference muR
 1 = muF1_over_ref ! ratio of current muF1 over reference muF1
 1 = muF2_over_ref ! ratio of current muF2 over reference muF2
 1 = 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 *
#***********************************************************************
 .false. = 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
 .false. = reweight_PDF ! reweight to get PDF uncertainty
  10901 = PDF_set_min ! First of the error PDF sets
  10952 = PDF_set_max ! Last of the error PDF sets
#***********************************************************************
# 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 more details. *
#***********************************************************************
 0 = ickkw ! 0 no merging, 3 FxFx merging
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) *
#***********************************************************************
 15 = bwcutoff
#***********************************************************************
# Cuts on the jets *
# 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 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
 0.7 = jetradius ! The radius parameter for the jet algorithm
  10 = ptj ! Min jet transverse momentum
  -1 = 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) *
#***********************************************************************
   0 = ptl ! Min lepton transverse momentum
  -1 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
   0 = drll ! Min distance between opposite sign lepton pairs
  30 = mll ! Min inv. mass of all oppositely charged lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
#***********************************************************************
  20 = ptgmin ! Min photon transverse momentum
  -1 = 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 as a jet *
#***********************************************************************
 5 = maxjetflavor
#***********************************************************************

###############################################
Do you have any idea what is the reason for that?

Thank you,

Benedikt

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Last query:
Last reply:
Revision history for this message
Benedikt Maier (bmaier) said : 		#1

Dear all,

I hope this thread becomes not too confusing, but maybe this is of any importance:

As mentioned in the post above, I receive the above error when I use a modified setscales.f (although this worked with MG5_v2_0_0).

I now tried to generate it with the default setscales.f, otherwise same settings, and this generation run worked! However, the debug information gives me warnings about O(300 000) "errors":

INFO: Generating events
INFO: Idle: 0, Running: 8, Completed: 0 [ current time: 19h36 ]
INFO: Idle: 0, Running: 7, Completed: 1 [ 1m 48s ]
INFO: Idle: 0, Running: 6, Completed: 2 [ 2m 36s ]
INFO: Idle: 0, Running: 5, Completed: 3 [ 32m 58s ]
INFO: Idle: 0, Running: 4, Completed: 4 [ 55m 34s ]
INFO: Idle: 0, Running: 3, Completed: 5 [ 4h 33m ]
INFO: Idle: 0, Running: 2, Completed: 6 [ 4h 35m ]
INFO: Idle: 0, Running: 1, Completed: 7 [ 8h 0m ]
INFO: Idle: 0, Running: 0, Completed: 8 [ 8h 5m ]
INFO: Collecting events
DEBUG:
  Number of loop ME evaluations (by MadLoop): 2213074
    Stability unknown: 0
    Stable PS point: 2213068
    Unstable PS point (and rescued): 6
    Unstable PS point (and not rescued): 0
    Only double precision used: 2213068
    Quadruple precision used: 6
    Initialization phase-space points: 24

  Statistics on virtual integration optimization :
    Maximum virt fraction computed 0.045 (P0_uxg_btxdx/GF1)
    Minimum virt fraction computed 0.006 (P0_uxg_btxdx/GF2)
    Average virt fraction computed 0.018
  Stats below exclude negligible channels (2 excluded out of 8)
    Maximum virt ratio used 20.25 (P0_gd_btxu/GF2)
    Maximum virt ratio found from grids 27.20 (P0_gd_btxu/GF2)
    Max. MC err. on virt ratio from grids 2.1 % (P0_uxg_btxdx/GF1)
    Maximum MC error on abs virt 13.3 % (P0_dg_btxu/GF1)

  Inclusive timing profile :
    Overall slowest channel 5:18:28 (P0_dg_btxu/GF1 [step 2])
    Average channel running time 1:01:15
    Aggregated total running time 1 day, 8:40:14
  Timing profile for <PDF_Engine> :
    Overall fraction of time 25.707 %
    Largest fraction of time 30.679 % (P0_gux_btxdx/GF1 [step 1])
    Smallest fraction of time 19.789 % (P0_gd_btxu/GF2 [step 0])
  Timing profile for <clustering> :
    Overall fraction of time 18.061 %
    Largest fraction of time 22.184 % (P0_gux_btxdx/GF2 [step 1])
    Smallest fraction of time 14.061 % (P0_gd_btxu/GF2 [step 0])
  Timing profile for <other_tasks> :
    Overall fraction of time 17.972 %
    Largest fraction of time 21.404 % (P0_dg_btxu/GF1 [step 1])
    Smallest fraction of time 13.747 % (P0_gd_btxu/GF2 [step 0])
  Timing profile for <IS_evaluation> :
    Overall fraction of time 15.374 %
    Largest fraction of time 17.525 % (P0_uxg_btxdx/GF2 [step 1])
    Smallest fraction of time 11.857 % (P0_gd_btxu/GF2 [step 0])
  Timing profile for <OneLoop_Engine> :
    Overall fraction of time 13.880 %
    Largest fraction of time 32.955 % (P0_gd_btxu/GF2 [step 0])
    Smallest fraction of time 1.410 % (P0_uxg_btxdx/GF2 [step 1])
  Timing profile for <Reals_evaluation> :
    Overall fraction of time 6.896 %
    Largest fraction of time 8.607 % (P0_gd_btxu/GF1 [step 1])
    Smallest fraction of time 5.800 % (P0_uxg_btxdx/GF1)
  Timing profile for <PS_Generation> :
    Overall fraction of time 2.108 %
    Largest fraction of time 2.471 % (P0_uxg_btxdx/GF2 [step 1])
    Smallest fraction of time 1.574 % (P0_gd_btxu/GF2 [step 0])
      WARNING:: A total of 335156 errors have been found in the following log files:
       > 82136 errors in tchan_antitop_2_NLO_4FS/SubProcesses/P0_gd_btxu/GF1/log.txt
       > 278 errors in tchan_antitop_2_NLO_4FS/SubProcesses/P0_gd_btxu/GF1/log_MINT0.txt
       > 1094 errors in tchan_antitop_2_NLO_4FS/SubProcesses/P0_gd_btxu/GF1/log_MINT1.txt
      And another 251648 errors in 29 other log files

INFO:
      Summary:
      Process p p > t~ b j $$ w+ w- [QCD]
      Run at p-p collider (6500 + 6500 GeV)
      Total cross-section: 6.956e+01 +- 2.8e-01 pb
      Number of events generated: 500000
      Parton shower to be used: PYTHIA8
      Fraction of negative weights: 0.37
      Total running time : 11h 43m

So I think it is basically about the question: What has changed from v2_0_0 to v2_0_1 that my modified setscales.f does not work anymore? Or what can I do in order to make it work again (we want to use the dynamical scale of the b quark, as proposed in one of R. Frederix' paper)?
And do I need to be concerned about 335156 errors :)?

Thank you very much!
Benedikt Maier

Revision history for this message
Valentin Hirschi (valentin-hirschi) said : 		#2

Hi Benedikt,

Firstly, in you scale setting, in the expression

tmp = 4.7**2 + pt(pp(0,nexternal-1))**2
        tmp=4*sqrt(tmp)
        temp_scale_id='4*Sqrt[m(b)**2 + pT(b)**2], b=b quark mur2'

the momentum pp(0,nexternal-1) seems to pick the top quark and not the bottom one, according to the order set by aMC@NLO and that you can see from the P* folder name:

/afs/[...]/SubProcesses/P0_gux_btxdx/ajob1

Secondly, I see that in the run with your custom scale, the job stopped because of a floating point expression.
Could you make sure that the file

[...]/SubProcesses/trapfpe.c

is empty. In the development version this file contains special code to trap floating point exception which can also be triggered on harmless ones. If it is not, please empty it, recompile an try again.

It might also be instructive to run the 'job' which crashed with a floating point exception independently and locally on the front end cluster. This gives you the luxury of running it with gdb so as to get pointed to the origin of the FPE. I.e.:

cd /afs/cern.ch/work/b/bmaier/public/lhe_files_CSA14/tchan_antitop_NLO_4FS/SubProcesses/P0_gux_btxdx/

then

gdb ajob1 2 F 1

followed by

run

This might help to understand what went wrong.

Finally, concerning the many Errors, this is fixed in the current development version which should be released soon. They are harmless and you can disregard them for now.

Let me know if any of the above helps, and if not I'll try to reproduce your issue to investigate it further.

Revision history for this message
Valentin Hirschi (valentin-hirschi) said : 		#3

Closing the ticket.

Can you help with this problem?

Provide an answer of your own, or ask Benedikt Maier for more information if necessary.

To post a message you must log in.