Trying untested ttH+2Jets process in aMC@NLO

Asked by Marco A. Harrendorf on 2014-05-19

Hello Generator-Experts,

at the moment we are working on FxFx-Merging for the ttH process.
For that reason I was also trying to generate the untested p p > t t~ h j j [QCD] process.

After I increased the RAM up to 32 Gb, I can create the working folder by using the command "output".
Just for your information, if you are interested:
During the "output" command the RAM demand goes up to 22 Gb. Finally when aMC@NLO "Generate Webpages" the RAM demand goes up to 30 Gb. I already switched off the JPEG generation.
So maybe in a future release it would be also nice, if one can switch off the webpage generation also.

When I try to launch the generated ttH+2Jets process I get the following output, before aMC@NLO crashes:
--------------------------------------------------------------------------------------------------------------------------------------

INFO: P0_uxu_httxccx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_gu_httxug
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_uux_httxgg
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_gux_httxuxg
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_gg_httxuux
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
INFO: P0_gg_httxgg
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
Command "launch auto " interrupted with error:
IOError : [Errno 2] No such file or directory: '/local/scratch/harrendorf/tthjj_test2/SubProcesses/P0_gg_httxgg/check_poles.log'
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/local/scratch/harrendorf/tthjj_test2/run_01_tag_1_debug.log'.
Please attach this file to your report.

-------------------------------------------------------------------------------------------------------------------------------

The mentioned "run_01_tag_1_debug.log" log file looks like this:
Comment: We want to use a new FxFx-Merging method with Pythia8 in the CMSSW Framework. Because of this 3 = ickkw is set in the run card.

-------------------------------------------------------------------------------------------------------------------------------

#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.1.1 2014-03-31 *
#* *
#* 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 "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/extended_cmd.py", line 872, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/extended_cmd.py", line 865, in onecmd_orig
    return func(arg, **opt)
  File "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 1184, in do_launch
    self.compile(mode, options)
  File "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 3283, in compile
    self.check_tests(test, this_dir)
  File "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 3297, in check_tests
    return self.parse_check_poles_log(pjoin(dir, '%s.log' % test))
  File "/usr/users/harrendorf/mg5/MG5_aMC_v2_1_1/madgraph/interface/amcatnlo_run_interface.py", line 3316, in parse_check_poles_log
    content = open(log).read()
IOError: [Errno 2] No such file or directory: '/local/scratch/harrendorf/tthjj_test2/SubProcesses/P0_gg_httxgg/check_poles.log'
Value of current Options:
              text_editor : emacs
              web_browser : None
        cluster_temp_path : None
       group_subprocesses : False
                  timeout : 0
            cluster_queue : None
         madanalysis_path : None
                   lhapdf : lhapdf-config
ignore_six_quark_processes : False
             stdout_level : 20
    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
      output_dependencies : external
          pythia-pgs_path : None
                 mg5_path : /usr/users/harrendorf/mg5/MG5_aMC_v2_1_1
                  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 : None
                    gauge : unitary
      complex_mass_scheme : False
                  nb_core : 7
                 run_mode : 2
#************************************************************
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.1.1 2014-03-31 *
#* *
#* 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 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~
import model loop_sm-no_b_mass
define p = p b b~
define j = j b b~
generate p p > t t~ h j j [QCD]
output tthjj_test2 -nojpeg
######################################################################
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
######################################################################
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
######################################################################

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

###################################
## INFORMATION FOR MASS
###################################
Block mass
    6 1.725000e+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
  5 0.000000 # b : 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
    6 1.725000e+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 *
#***********************************************************************
 100000 = 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
 5000 = npoints_FO_grid ! number of points to setup grids
 4 = niters_FO_grid ! number of iter. to setup grids
 10000 = npoints_FO ! number of points to compute Xsec
 6 = niters_FO ! number of iter. to compute Xsec
#***********************************************************************
# Random number seed *
#***********************************************************************
     0 = iseed ! rnd seed (0=assigned automatically=default))
#***********************************************************************
# Collider type and energy *
#***********************************************************************
    1 = lpp1 ! beam 1 type (0 = no PDF)
    1 = lpp2 ! beam 2 type (0 = no PDF)
 4000 = ebeam1 ! beam 1 energy in GeV
 4000 = ebeam2 ! beam 2 energy in GeV
#***********************************************************************
# PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
#***********************************************************************
 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
 91.188 = muR_ref_fixed ! fixed ren reference scale
 91.188 = muF1_ref_fixed ! fixed fact reference scale for pdf1
 91.188 = 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
 91.188 = 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 *
#***********************************************************************
 .true. = reweight_scale ! reweight to get scale dependence
  0.5 = rw_Rscale_down ! lower bound for ren scale variations
  2.0 = rw_Rscale_up ! upper bound for ren scale variations
  0.5 = rw_Fscale_down ! lower bound for fact scale variations
  2.0 = rw_Fscale_up ! upper bound for fact scale variations
 .false. = reweight_PDF ! reweight to get PDF uncertainty
  21101 = PDF_set_min ! First of the error PDF sets
  21140 = 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. *
#***********************************************************************
 3 = ickkw ! 0 no merging, 3 FxFx merging
#***********************************************************************
#
#***********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) *
#***********************************************************************
 15 = 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 = 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 a jet when applying cuts.*
# At least all massless quarks of the model should be included here. *
#***********************************************************************
 5 = maxjetflavor
#***********************************************************************

---------------------------------------------------------------------------------------------------------------------------------

I already checked, that the "SubProcesses/P0_gg_httxgg/check_poles.log" log file is really missing.

Do you have any idea, what I could try to get the ttH+2Jet process running or is it hopeless?

Thank you in advance.

Kind regards,

Marco

Question information

Language:
English Edit question
Status:
Open
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Rikkert Frederix Edit question
Last query:
2014-06-16
Last reply:
2014-06-05
Rikkert Frederix (frederix) said : #1

Dear Marco,

Is the SubProcesses/P0_gg_httxgg/check_poles executable there?

If so, can you run ./check_poles and follow the instructions written on the screen. Give it 20 PS points to test with the default precision.

If not, can you try to compile it? Run
export madloop=true
make check_poles
(If it complains about fastjet not being there, you have to set --in the beginning of the makefile:
fastjet_config=<PATH_TO_FASTJET>/bin/fastjet-config
and then run it?

Let me know if this gives any useful information.

Best regards,
Rikkert

Hello Rikkert,

sorry for the long delay.

The executable check_poles was not in the mentioned directory.

When I try to compile check_poles I get the following output:

-----------------------------------------------
f77 check_poles.f -o check_poles
   MAIN driver:
Error on line 13 of genps.inc: unbalanced parentheses, statement skipped
Error processing common blocks before line 61: Declaration error for fixed_renfixed_fac_scale: attempt to use undefined variable
Error processing common blocks before line 61: Declaration error for hmult: attempt to use undefined variable
Error processing common blocks before line 61: Declaration error for rw_fscale_down: common alignment
Error processing common blocks before line 61: Declaration error for rw_fscale_up: common alignment
Error processing common blocks before line 61: Declaration error for rw_rsdo_rwgt_pdf: attempt to use undefined variable
Error on line 89: syntax error
Error on line 90: misplaced ENDDO
Warning on line 237: RETURN statement in main or block data
Error on line 238: DO loop or BLOCK IF not closed
Warning on line 238: local variable ng never used
Warning on line 238: local variable hm never used
Warning on line 238: local variable rw_ never used
Warning on line 238: local variable do_rwgt_pdf never used
Warning on line 238: local variable pswgt never used
Warning on line 238: local variable maxdim never used
Warning on line 238: local variable rw_rscale_down never used
Warning on line 238: local variable fixed_fac_scale never used
Warning on line 238: local variable fixed_ren_scale never used
   rambo:
   rans:
   outfun:
/usr/bin/f77: aborting compilation
make: *** [check_poles] Fehler 25

---------------------------------------------------

The genps.inc file looks like:
----------------------------------------------------
*************************************************************************
c Parameters used by genps
c*************************************************************************
c*************************************************************************
c Parameters used by genps and dsample, you must recompile
c dsample if you change anything below
c*************************************************************************
      include 'maxparticles.inc'
      integer max_tree , max_dim
      parameter (max_tree=200, max_dim=max_tree*(max_branch+1))
      integer ng , maxdim , maxinvar , maxconfigs
c parameter (ng = 96, maxdim = 25, maxinvar= 50 , maxconfigs=10)
      parameter (ng = 96, maxdim = 3*(max_particles-2)-1, maxinvar= 4*max_particles, maxconfigs=10)

      include 'maxconfigs.inc'

      double precision xgmin, xgmax
      parameter (xgmin=-1d0, xgmax=1d0)

      integer maxproc, maxamps
      parameter (maxproc=100, maxamps=7500)

c integer maxproc
c parameter (maxproc=100)

      integer maxevents !Requires about 1K/event
      parameter (maxevents=100000) !Maximum # events to write to disk

c*************************************************************************
c Parameters used for parrallel running
c*************************************************************************
      integer max_host ,maxplace ,maxpoints ,maxans
      parameter (max_host=9,maxplace=9,maxpoints=10,maxans=5)
      integer maxprb
      parameter (maxprb = maxconfigs*maxplace*maxpoints)
      integer maxfprb
      parameter (maxfprb = maxinvar*maxplace*maxpoints)

----------------------------------------------------

Best regards,

Marco

Rikkert Frederix (frederix) said : #3

It looks like you are not compiling with the same compiler as you have used for the rest of the code. Please use the same.

Best,
Rikkert

Hello Rikkert,

I made some checks and tests:

* I used the same compiler as before and my settings work for more basic
processes like ttH.
* The check_poles executable is build for none of the subprocesses.
* There is also no makefile in the subprocesses folders.
* Creating the process new from scratch does not solve my problems.

* I copied the makefile from the basic ttH process to the ttH+2Jets
subprocess folders and now I can build the check_poles executable.

Is there any reason, why the makefiles are not written in the subprocess
folders?
Can this behaviour caused by switching off the JPEG generation?

Best regards,

Marco

--
Marco A. Harrendorf <email address hidden>
KIT - IEKP

Am Donnerstag, den 05.06.2014, 13:56 +0000 schrieb Rikkert Frederix:
> Your question #248891 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/248891
>
> Status: Open => Needs information
>
> Rikkert Frederix requested more information:
> It looks like you are not compiling with the same compiler as you have
> used for the rest of the code. Please use the same.
>
> Best,
> Rikkert
>

Rikkert Frederix (frederix) said : #5

Hi Marco,

I have no idea what's going on. I cannot reproduce the non-creation of the makefiles. It does not depend on the JPEG generation.

However, given that the process is not created correctly (missing makefiles), I cannot guarantee that the rest of the process is correct.

Best,
Rikkert

Can you help with this problem?

Provide an answer of your own, or ask Marco A. Harrendorf for more information if necessary.

To post a message you must log in.