MadGraph5_aMC@NLO

P0_gg_httx/ajob1 1 all 0 0 launch ends with non zero status: 139. Stop all computation

Asked by Andrej Saibel

Dear all,

I'm sorry that I have to come back to https://answers.launchpad.net/mg5amcnlo/+question/679425

According to the recommendations in the question above, I'm moving my FO analysis to HwU from TopDrawer to have a more consistent set of uncertainties.

However, I've encountered two strange things:

1. After compiling the code, it crashes with the message in [1] during the setup of the grids. Unfortunately, the output of the log in [2] isn't very helpful for me.

2. I set explicitly mur_ref_fixed and muf_ref_fixed , but the block in the param_card with mu_r doesn't get changed (is still the one from the run before). Is this expected? You can find this information also in [2]

[1]

ype '0', 'auto', 'done' or just press enter when you are done.
INFO: will run in mode: NLO
set lhaid 303400
INFO: modify parameter lhaid of the run_card.dat to [303400]
set param_card mass 6 163
INFO: modify param_card information BLOCK mass with id (6,) set to 163.0
set param_card yukawa 6 163
INFO: modify param_card information BLOCK yukawa with id (6,) set to 163.0
set decay 6 0.0
INFO: modify param_card information BLOCK decay with id (6,) set to 0.0
set mur_ref_fixed 163
INFO: modify parameter mur_ref_fixed of the run_card.dat to 163.0
set muf_ref_fixed 163
INFO: modify parameter muf_ref_fixed of the run_card.dat to 163.0
set mur_over_ref 1
INFO: modify parameter mur_over_ref of the run_card.dat to 1.0
set muf_over_ref 1
INFO: modify parameter muf_over_ref of the run_card.dat to 1.0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. FO_analyse : FO_analyse_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, FO_analyse, enter path]
The answer to the previous question is not set in your input file
Use 0 value
INFO: Update the dependent parameter of the param_card.dat
INFO: Starting run
INFO: Compiling the code
INFO: Using built-in libraries for PDFs
INFO: Compiling source...
Gtk-Message: Failed to load module "canberra-gtk-module"
Gtk-Message: Failed to load module "canberra-gtk-module"
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_gg_httx...
INFO: Compiling P0_uux_httx...
INFO: Compiling P0_uxu_httx...
INFO: P0_uxu_httx done.
INFO: P0_uux_httx done.
INFO: P0_gg_httx done.
INFO: Checking test output:
INFO: P0_gg_httx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uux_httx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uxu_httx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: Starting run
INFO: Using 4 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
INFO: Setting up grids
/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ttH_Test_HwU/SubProcesses/P0_gg_httx/ajob1: line 34: 107244 Segmentation fault ../madevent_mintFO > log.txt < input_app.txt 2>&1
WARNING: program /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ttH_Test_HwU/SubProcesses/P0_gg_httx/ajob1 1 all 0 0 launch ends with non zero status: 139. Stop all computation
/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ttH_Test_HwU/SubProcesses/P0_uux_httx/ajob1: line 34: 107246 Segmentation fault ../madevent_mintFO > log.txt < input_app.txt 2>&1
/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ttH_Test_HwU/SubProcesses/P0_uxu_httx/ajob1: line 34: 107248 Terminated ../madevent_mintFO > log.txt < input_app.txt 2>&1
INFO: Idle: 0, Running: 0, Completed: 3 [ current time: 17h39 ]

[2]

launch auto
Traceback (most recent call last):
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/amcatnlo_run_interface.py", line 1679, in do_launch
    evt_file = self.run(mode, options)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/amcatnlo_run_interface.py", line 1843, in run
    self.run_all_jobs(jobs_to_run,integration_step)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/amcatnlo_run_interface.py", line 2116, in run_all_jobs
    self.wait_for_complete(run_type)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/interface/amcatnlo_run_interface.py", line 4657, in wait_for_complete
    self.cluster.wait(self.me_dir, update_status)
  File "/nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/madgraph/various/cluster.py", line 822, in wait
    raise Exception, self.fail_msg
Exception: program /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ttH_Test_HwU/SubProcesses/P0_gg_httx/ajob1 1 all 0 0 launch ends with non zero status: 139. Stop all computation
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 : /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/HEPTools/lib
              auto_update : 0
             pythia8_path : None
                hwpp_path : None
low_mem_multicore_nlo_generation : False
                    golem : None
          pythia-pgs_path : None
                  td_path : None
             delphes_path : None
              thepeg_path : None
             cluster_type : condor
        madanalysis5_path : /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/HEPTools/madanalysis5/madanalysis5
      exrootanalysis_path : /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/ExRootAnalysis
                      OLP : MadLoop
                 applgrid : applgrid-config
               eps_viewer : None
                  fastjet : None
                 run_mode : 2
              web_browser : None
   automatic_html_opening : False
        cluster_temp_path : None
             cluster_size : 100
            cluster_queue : None
             syscalc_path : None
         madanalysis_path : None
                   lhapdf : /nfs/dust/cms/user/asaibel/MCGenerators/LHAPDF-6.2.1_Installation/bin/lhapdf-config
             stdout_level : 20
                  nb_core : 4
            f2py_compiler : None
                    ninja : /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5/HEPTools/lib
                  amcfast : amcfast-config
       cluster_retry_wait : 300
      output_dependencies : external
           crash_on_error : False
mg5amc_py8_interface_path : None
         loop_color_flows : False
                  samurai : None
         cluster_nb_retry : 1
                 mg5_path : /nfs/dust/cms/user/asaibel/MCGenerators/MG5aMCatNLOForttH/MG5_aMC_v2_6_5
                  timeout : 60
                    gauge : unitary
      complex_mass_scheme : False
             cpp_compiler : None
   max_npoint_for_channel : 0
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.5 2018-02-03 *
#* *
#* 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 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~
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 sm
define p = g u c d s u~ c~ d~ s~ b b~
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
define j = p
import model loop_sm-no_b_mass
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5
define j = p
define p = p b b~
define j = j b b~
generate p p > t t~ h [QCD]
output ttH_QCD_basic
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
######################################################################
###################################
## INFORMATION FOR LOOP
###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################
## INFORMATION FOR MASS
###################################
BLOCK MASS #
      6 1.630000e+02 # mt
      15 1.777000e+00 # mta
      23 9.118800e+01 # mz
      25 1.250000e+02 # mh
      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 8.041900e+01 # 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 #
      6 1.630000e+02 # ymt
      15 1.777000e+00 # ymtau
###################################
## INFORMATION FOR DECAY
###################################
DECAY 6 0.000000e+00 # wt
DECAY 23 2.441404e+00 # wz
DECAY 24 2.047600e+00 # ww
DECAY 25 6.382339e-03 # wh
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
###################################
## INFORMATION FOR QNUMBERS 82
###################################
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 *
#***********************************************************************
  20000000 = 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.0001 = 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
  1000000 = npoints_fo_grid ! number of points to setup grids
  9 = niters_fo_grid ! number of iter. to setup grids
  2000000 = npoints_fo ! number of points to compute Xsec
  9 = 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. *
#***********************************************************************
  nn23nlo = pdlabel ! PDF set
  303400 = 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) *
#***********************************************************************
  True = fixed_ren_scale ! if .true. use fixed ren scale
  True = fixed_fac_scale ! if .true. use fixed fac scale
  163.0 = mur_ref_fixed ! fixed ren reference scale
  163.0 = muf_ref_fixed ! fixed fact reference scale
  0 = 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
  False = 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.
  False = reweight_pdf ! Reweight to get PDF uncertainty. Should be a
            ! list booleans of equal length to lhaid to specify for
            ! which PDF set to include the uncertainties.
#***********************************************************************
# Store reweight information in the LHE file for off-line model- *
# parameter reweighting at NLO+PS accuracy *
#***********************************************************************
  False = store_rwgt_info ! Store info for reweighting in LHE file
#***********************************************************************
# ickkw parameter: *
# 0: No merging *
# 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
# level. After showering an MLM-type merging should be applied as *
# well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
# 4: UNLOPS merging (with pythia8 only). No interface from within *
# MG5_aMC available, but available in Pythia8. *
# -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
#***********************************************************************
  0 = ickkw
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
# written in the LHE event file *
#***********************************************************************
  15.0 = bwcutoff
#***********************************************************************
# Cuts on the jets. Jet clustering is performed by FastJet. *
# - When matching to a parton shower, these generation cuts should be *
# considerably softer than the analysis cuts. *
# - More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  1.0 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
  0.7 = jetradius ! The radius parameter for the jet algorithm
  10.0 = ptj ! Min jet transverse momentum
  -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
#***********************************************************************
# Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  0.0 = ptl ! Min lepton transverse momentum
  -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
  0.0 = drll ! Min distance between opposite sign lepton pairs
  0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
  0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
  30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
#***********************************************************************
# Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
# all the other parameters are ignored. *
# More specific cuts can be specified in SubProcesses/cuts.f *
#***********************************************************************
  20.0 = ptgmin ! Min photon transverse momentum
  -1.0 = etagamma ! Max photon abs(pseudo-rap)
  0.4 = r0gamma ! Radius of isolation code
  1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
  1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
  True = isoem ! isolate photons from EM energy (photons and leptons)
#***********************************************************************
# 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:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Valentin Hirschi Edit question
Last query:
Last reply:
Revision history for this message
Andrej Saibel (asaibel) said : 		#1

Adding some more info:

The following is my HwU card:

c
c This file contains the default histograms for fixed order runs: it
c only plots the total rate as an example. It can be used as a template
c to make distributions for other observables.
c
c This uses the HwU package and generates histograms in the HwU/GnuPlot
c format. This format is human-readable. After running, the histograms
c can be found in the Events/run_XX/ directory.
c
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      subroutine analysis_begin(nwgt,weights_info)
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c This subroutine is called once at the start of each run. Here the
c histograms should be declared.
c
c Declare the histograms using 'HwU_book'.
c o) The first argument is an integer that labels the histogram. In
c the analysis_end and analysis_fill subroutines this label is used
c to keep track of the histogram. The label should be a number
c starting at 1 and be increased for each plot.
c o) The second argument is a string that will apear above the
c histogram. Do not use brackets "(" or ")" inside this string.
c o) The third, forth and fifth arguments are the number of bis, the
c lower edge of the first bin and the upper edge of the last
c bin.
c o) When including scale and/or PDF uncertainties, declare a
c histogram for each weight, and compute the uncertainties from the
c final set of histograms
c
      implicit none
c When including scale and/or PDF uncertainties the total number of
c weights considered is nwgt
      integer nwgt
c In the weights_info, there is an text string that explains what each
c weight will mean. The size of this array of strings is equal to nwgt.
      character*(*) weights_info(*)
c Initialize the histogramming package (HwU). Pass the number of
c weights and the information on the weights:
      call HwU_inithist(nwgt,weights_info)
c declare (i.e. book) the histograms
      call HwU_book(1,'total rate ', 5,0.5d0,5.5d0)
      call HwU_book(2,'total rate Born ', 5,0.5d0,5.5d0)
      call HwU_book(1001,'pt H',10.0d0,0.0d0,750d0)
      call HwU_book(1002,'pt anti',10.0d0,0d0, 750d0)
      call HwU_book(1003,'inv mass',20d0,440d0,2000d0)
      call HwU_book(1004,'y anti',0.20d0,-3.0d0,3.0d0)
      call HwU_book(1005,'y H',0.20d0,-3.0d0,3.0d0)
      call HwU_book(1006,'tt inv',20.0d0,340d0,1200d0)
      call HwU_book(1007,'zdist',0.05d0,0.0d0,1.1d0)
      call HwU_book(1008,'mh inv',10.0d0,120.0d0,180.0d0)
      call HwU_book(1009,'mt inv',2.0d0,120.0d0,180.0d0)
      call HwU_book(1010,'pt top',10.0d0,0d0,750d0)
      call HwU_book(1011,'y top',0.20d0,-3.0d0,3.0d0)
c call HwU_book(1012,'pt ttbarH',10d0,0.0d0,750.0d0)
      return
      end

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      subroutine analysis_end(dummy)
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c This subroutine is called once at the end of the run. Here the
c histograms are written to disk. Note that this is done for each
c integration channel separately. There is an external script that will
c read the HwU data files in each of the integration channels and
c combines them by summing all the bins in a final single HwU data file
c to be put in the Events/run_XX directory, together with a gnuplot
c file to convert them to a postscript histogram file.
      implicit none
      double precision dummy
      call HwU_write_file
      return
      end

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
      subroutine analysis_fill(p,istatus,ipdg,wgts,ibody)
cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
c This subroutine is called for each n-body and (n+1)-body configuration
c that passes the generation cuts. Here the histrograms are filled.
      implicit none
c This includes the 'nexternal' parameter that labels the number of
c particles in the (n+1)-body process
      include 'nexternal.inc'
c This is an array which is '-1' for initial state and '1' for final
c state particles
      integer istatus(nexternal)
c This is an array with (simplified) PDG codes for the particles. Note
c that channels that are combined (i.e. they have the same matrix
c elements) are given only 1 set of PDG codes. This means, e.g., that
c when using a 5-flavour scheme calculation (massless b quark), no
c b-tagging can be applied.
      integer iPDG(nexternal)
c The array of the momenta and masses of the initial and final state
c particles in the lab frame. The format is "E, px, py, pz, mass", while
c the second dimension loops over the particles in the process. Note
c that these are the (n+1)-body particles; for the n-body there is one
c momenta equal to all zero's (this is not necessarily the last particle
c in the list). If one uses IR-safe obserables only, there should be no
c difficulty in using this.
      double precision p(0:4,nexternal)
c The weight of the current phase-space point is wgts(1). If scale
c and/or PDF uncertainties are included through reweighting, the rest of
c the array contains the list of weights in the same order as described
c by the weigths_info strings in analysis_begin
      double precision wgts(*)
c The ibody variable is:
c ibody=1 : (n+1)-body contribution
c ibody=2 : n-body contribution (excluding the Born)
c ibody=3 : Born contribution
c The histograms need to be filled for all these contribution to get the
c physics NLO results. (Note that the adaptive phase-space integration
c is optimized using the sum of the contributions, therefore plotting
c them separately might lead to larger than expected statistical
c fluctuations).
      integer ibody
c local variable
      double precision var,pt_top,pt_5,tt_inv,inv_mass,pt_H,y_H,pt_j,y_5,y_4,pt_4,mh_inv,mt_inv,zdist,zden
      double precision pi
      PARAMETER (PI=3.14159265358979312D0)
c
c Fill the histograms here using a call to the HwU_fill()
c subroutine. The first argument is the histogram label, the second is
c the numerical value of the variable to plot for the current
c phase-space point and the final argument is the weight of the current
c phase-space point.
      var=1d0
c always fill the total rate
      call HwU_fill(1,var,wgts)

      pt_5 = dsqrt(p(1,5)**2+p(2,5)**2)
      inv_mass = dsqrt((p(0,3)+p(0,4)+p(0,5))**2-(p(1,3)+p(1,4)+p(1,5))**2-(p(2,3)+p(2,4)+p(2,5))**2-(p(3,3)+p(3,4)+p(3,5))**2)
      tt_inv = dsqrt((p(0,4)+p(0,5))**2-(p(1,4)+p(1,5))**2-(p(2,4)+p(2,5))**2-(p(3,4)+p(3,5))**2)
      pt_H = dsqrt(p(1,3)**2+p(2,3)**2)
      y_H = 1.0d0/2.0d0*Log((p(0,3)+p(3,3))/(p(0,3)-p(3,3)))
      pt_j = dsqrt(p(1,6)**2+p(2,6)**2)
      y_5 = 1.0d0/2.0d0*Log((p(0,5)+p(3,5))/(p(0,5)-p(3,5)))
      mh_inv = dsqrt(p(0,3)**2-p(1,3)**2-p(2,3)**2-p(3,3)**2)
      mt_inv = dsqrt(p(0,5)**2-p(1,5)**2-p(2,5)**2-p(3,5)**2)
      pt_4 = dsqrt(p(1,4)**2+p(2,4)**2)
      y_4 = 1.0d0/2.0d0*Log((p(0,4)+p(3,4))/(p(0,4)-p(3,4)))
      zdist = (p(0,3)+p(0,4)+p(0,5))**2-(p(1,3)+p(1,4)+p(1,5))**2-(p(2,3)+p(2,4)+p(2,5))**2-(p(3,3)+p(3,4)+p(3,5))**2
      zden = (p(0,3)+p(0,4)+p(0,5)+p(0,6))**2-(p(1,3)+p(1,4)+p(1,5)+p(1,6))**2
      zden= zden-(p(2,3)+p(2,4)+p(2,5)+p(2,6))**2-(p(3,3)+p(3,4)+p(3,5)+p(3,6))**2
      zdist = zdist/zden

      call HwU_fill(1001,pt_H,wgts)
      call HwU_fill(1002,pt_5,wgts)
      call HwU_fill(1003,inv_mass,wgts)
      call HwU_fill(1004,y_5,wgts)
c call HwU_fill(1005,y_H,wgts/0.25d0)
      call HwU_fill(1005,y_H,wgts)
      call HwU_fill(1006,tt_inv,wgts)
      call HwU_fill(1007,zdist,wgts)
      call HwU_fill(1008,mh_inv,wgts)
      call HwU_fill(1009,mt_inv,wgts)
      call HwU_fill(1010,pt_4,wgts)
      call HwU_fill(1011,y_4,wgts)
c call HwU_fill(1012,pt_ttbarH,wgts/10.0d0)

c only fill the total rate for the Born when ibody=3
      if (ibody.eq.3) call HwU_fill(2,var,wgts)
      return
      end

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

error 139 is a segmentation fault.
One advise is to modify Source/make_opts to increase the compiler check flag, this should help you to find at which line of your HWU code, you made a mistake.

you can for example uncomment the following line of that file.
  #FFLAGS+= -g -fbounds-check -ffpe-trap=invalid,zero,overflow,underflow,denormal -Wall

The critical flags being -fbounds-check

Cheers,

Olivier

Revision history for this message
Andrej Saibel (asaibel) said : 		#3

Dear Olivier,

thanks a lot for your answer.
What about the second issue? is this setting irrelevant?

2. I set explicitly mur_ref_fixed and muf_ref_fixed , but the block in the param_card with mu_r doesn't get changed (is still the one from the run before). Is this expected? You can find this information also in [2]

my input:

order=LO
fixed_order=ON
set lhaid 303400
set param_card mass 6 163
set param_card yukawa 6 163
set decay 6 0.0
set mur_ref_fixed 163
set muf_ref_fixed 163
set mur_over_ref 1
set muf_over_ref 1

I still get the following in the param card:

###################################
BLOCK LOOP #
      1 9.118800e+01 # mu_r
###################################

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

HI,

I think that in such cases, the value inside the param_card is irrelevant (like the value of alpha_s).
But Valentin should know more about that and will comment if needed.

Cheers,

Olivier

Can you help with this problem?

Provide an answer of your own, or ask Andrej Saibel for more information if necessary.

To post a message you must log in.