MadSpin message: CRITICAL: over_weight

Asked by Marc de Beurs on 2018-08-22

Dear Madspin experts,

When I am running a single top sample, (t-channel or s-channel), and I want the tops to decay lepontically (t > w+ b, w+ > l+ vl) in madspin, I sometimes run into the critical message of:

CRITICAL: over_weight: 1.24394983756 ('tx_bxwm_wm_tamvtx',), occurence: 0.0492610837438%, occurence_channel: 0.392156862745% production_tag:P4_us_txbuc [((1, 4), (-6, 2, 4, 5))], decay:P4_us_txbuc_tx_bxwm_wm_tamvtx [('tx_bxwm_wm_tamvtx',)], BW_cut: 13.097 [decay.py at line 2297]

Is this something to worry about?

I am running MG5 version 2.6.2 and below you can find the std output which also contains the cards (I am running it inside the atlas framework, so you can ignore all the messages that are unknown to you):

Cheers Marc

11:34:57 Tue Aug 21 11:34:57 CEST 2018
11:34:57 Preloading tcmalloc_minimal.so
11:34:57 Preloading /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/IntelSoftware/linux/x86_64/xe2013/composer_xe_2013.3.163/compiler/lib/intel64/libintlc.so.5:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/IntelSoftware/linux/x86_64/xe2013/composer_xe_2013.3.163/compiler/lib/intel64/libimf.so
11:34:57 Py:Athena INFO including file "AthenaCommon/Preparation.py"
11:34:57 Py:Athena INFO using release [MCProd-19.2.5.33.4] [x86_64-slc6-gcc47-opt] [19.2.X.Y.Z-VAL-MCProd/rel_0] -- built on [2018-05-23 11:02]
11:34:57 Py:Athena INFO including file "AthenaCommon/Bootstrap.py"
11:34:57 Py:Athena INFO including file "AthenaCommon/Atlas.UnixStandardJob.py"
11:34:59 Py:Athena INFO including file "AthenaCommon/Execution.py"
11:34:59 Py:Athena INFO including file "runargs.generate.py"
11:34:59 Py:Athena INFO including file "EvgenJobTransforms/skeleton.GENtoEVGEN.py"
11:34:59 Py:Athena INFO including file "AthenaCommon/Atlas_Gen.UnixStandardJob.py"
11:34:59 Py:Athena INFO including file "PartPropSvc/PartPropSvc.py"
11:35:00 Py:Athena INFO including file "MC15.412049.aMcAtNloHerwig7EG_H7UE_singletop_schan_lept.py"
11:35:00 Py:Athena INFO including file "./MadGraphControl_schan_NLO_Herwig7_A14.py"
11:35:00 Py:MadGraphUtils INFO Assuming that proc_card_mg5.dat is a well-formatted process card
11:35:00 Py:MadGraphUtils INFO Started process generation at Tue Aug 21 11:35:00 2018
11:35:04 [?1034h************************************************************
11:35:04 * *
11:35:04 * W E L C O M E to *
11:35:04 * M A D G R A P H 5 _ a M C @ N L O *
11:35:04 * *
11:35:04 * *
11:35:04 * * * *
11:35:04 * * * * * *
11:35:04 * * * * * 5 * * * * *
11:35:04 * * * * * *
11:35:04 * * * *
11:35:04 * *
11:35:04 * VERSION 2.6.2 2018-04-29 *
11:35:04 * *
11:35:04 * The MadGraph5_aMC@NLO Development Team - Find us at *
11:35:04 * https://server06.fynu.ucl.ac.be/projects/madgraph *
11:35:04 * and *
11:35:04 * http://amcatnlo.web.cern.ch/amcatnlo/ *
11:35:04 * *
11:35:04 * Type 'help' for in-line help. *
11:35:04 * Type 'tutorial' to learn how MG5 works *
11:35:04 * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
11:35:04 * Type 'tutorial MadLoop' to learn how MadLoop works *
11:35:04 * *
11:35:04 ************************************************************
11:35:04 load MG5 configuration from ../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
11:35:04 set collier to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/collier/1.1/x86_64-slc6-gcc47-opt
11:35:04 fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
11:35:04 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
11:35:04 MG5_aMC> set fastjet /PATH/TO/fastjet-config
11:35:04
11:35:04 set lhapdf to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.4/x86_64-slc6-gcc47-opt/bin/lhapdf-config
11:35:04 set ninja to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib
11:35:04 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
11:35:04 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
11:35:04 import /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/proc_card_mg5.dat
11:35:04 The import format was not given, so we guess it as command
11:35:04 import model loop_sm-no_b_mass
11:35:05 INFO: reload from .py file
11:35:05 INFO: load particles
11:35:05 INFO: load vertices
11:35:05 INFO: Restrict model loop_sm with file ../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
11:35:05 INFO: Run "set stdout_level DEBUG" before import for more information.
11:35:05 INFO: Change particles name to pass to MG5 convention
11:35:05 Defined multiparticle p = g u c d s u~ c~ d~ s~
11:35:05 Defined multiparticle j = g u c d s u~ c~ d~ s~
11:35:05 Defined multiparticle l+ = e+ mu+
11:35:05 Defined multiparticle l- = e- mu-
11:35:05 Defined multiparticle vl = ve vm vt
11:35:05 Defined multiparticle vl~ = ve~ vm~ vt~
11:35:05 Pass the definition of 'j' and 'p' to 5 flavour scheme.
11:35:05 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
11:35:05 define l+ = l+ ta+
11:35:05 Defined multiparticle l+ = e+ mu+ ta+
11:35:05 define l- = l- ta-
11:35:05 Defined multiparticle l- = e- mu- ta-
11:35:05 generate p p > w+ > t b~ [QCD]
11:35:05 WARNING: Process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > w+> t b~ WEIGHTED=4 [ all = QCD ] can have real emission processes which are not finite.
11:35:05 To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
11:35:05 g, b, c, s, u, d, d~, u~, s~, c~, b~ 
11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: u d~ > w+ > t b~ [ all = QCD ] (1 / 4)
11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: c s~ > w+ > t b~ [ all = QCD ] (2 / 4)
11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: d~ u > w+ > t b~ [ all = QCD ] (3 / 4)
11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: s~ c > w+ > t b~ [ all = QCD ] (4 / 4)
11:35:05 INFO: Generating virtual matrix elements using MadLoop:
11:35:05 INFO: Generating virtual matrix element with MadLoop for process: u d~ > w+ > t b~ QED<=2 [ all = QCD ] (1 / 4)
11:35:05 INFO: Generating virtual matrix elements using MadLoop:
11:35:05 INFO: Generating virtual matrix element with MadLoop for process: c s~ > w+ > t b~ QED<=2 [ all = QCD ] (2 / 4)
11:35:06 INFO: Generating virtual matrix elements using MadLoop:
11:35:06 INFO: Generating virtual matrix element with MadLoop for process: d~ u > w+ > t b~ QED<=2 [ all = QCD ] (3 / 4)
11:35:06 INFO: Generating virtual matrix elements using MadLoop:
11:35:06 INFO: Generating virtual matrix element with MadLoop for process: s~ c > w+ > t b~ QED<=2 [ all = QCD ] (4 / 4)
11:35:06 INFO: Generated 4 subprocesses with 32 real emission diagrams, 4 born diagrams and 8 virtual diagrams
11:35:06 add process p p > w- > t~ b [QCD]
11:35:06 WARNING: Process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > w-> t~ b WEIGHTED=4 [ all = QCD ] can have real emission processes which are not finite.
11:35:06 To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
11:35:06 g, b, c, s, u, d, d~, u~, s~, c~, b~ 
11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: d u~ > w- > t~ b [ all = QCD ] (1 / 4)
11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: s c~ > w- > t~ b [ all = QCD ] (2 / 4)
11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: u~ d > w- > t~ b [ all = QCD ] (3 / 4)
11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: c~ s > w- > t~ b [ all = QCD ] (4 / 4)
11:35:06 INFO: Generating virtual matrix elements using MadLoop:
11:35:06 INFO: Generating virtual matrix element with MadLoop for process: d u~ > w- > t~ b QED<=2 [ all = QCD ] (1 / 4)
11:35:06 INFO: Generating virtual matrix elements using MadLoop:
11:35:06 INFO: Generating virtual matrix element with MadLoop for process: s c~ > w- > t~ b QED<=2 [ all = QCD ] (2 / 4)
11:35:06 INFO: Generating virtual matrix elements using MadLoop:
11:35:06 INFO: Generating virtual matrix element with MadLoop for process: u~ d > w- > t~ b QED<=2 [ all = QCD ] (3 / 4)
11:35:07 INFO: Generating virtual matrix elements using MadLoop:
11:35:07 INFO: Generating virtual matrix element with MadLoop for process: c~ s > w- > t~ b QED<=2 [ all = QCD ] (4 / 4)
11:35:07 INFO: Generated 4 subprocesses with 32 real emission diagrams, 4 born diagrams and 8 virtual diagrams
11:35:07 output -f
11:35:07 INFO: Writing out the aMC@NLO code, using optimized Loops
11:35:07 INFO: initialize a new directory: PROCNLO_loop_sm-no_b_mass_0
11:35:11 INFO: remove old information in PROCNLO_loop_sm-no_b_mass_0
11:35:18 INFO: When using the 'output_dependencies=internal' MG5_aMC option, the (optional) reduction library ninja cannot be employed because it is not distributed with the MG5_aMC code so that it cannot be copied locally.
11:35:18 INFO: When using the 'output_dependencies=internal' MG5_aMC option, the (optional) reduction library collier cannot be employed because it is not distributed with the MG5_aMC code so that it cannot be copied locally.
11:35:18 INFO: Generating real emission matrix-elements...
11:35:19 INFO: Generating Helas calls for FKS process: u d~ > w+ > t b~ [ all = QCD ] (1 / 8)
11:35:19 INFO: Processing color information for process: u d~ > w+ > t b~ [ all = QCD ]
11:35:19 INFO: Generating Helas calls for FKS process: c s~ > w+ > t b~ [ all = QCD ] (2 / 8)
11:35:19 INFO: Generating Helas calls for FKS process: d~ u > w+ > t b~ [ all = QCD ] (3 / 8)
11:35:20 INFO: Processing color information for process: d~ u > w+ > t b~ [ all = QCD ]
11:35:20 INFO: Generating Helas calls for FKS process: s~ c > w+ > t b~ [ all = QCD ] (4 / 8)
11:35:20 INFO: Generating Helas calls for FKS process: d u~ > w- > t~ b [ all = QCD ] (5 / 8)
11:35:20 INFO: Processing color information for process: d u~ > w- > t~ b [ all = QCD ]
11:35:20 INFO: Generating Helas calls for FKS process: s c~ > w- > t~ b [ all = QCD ] (6 / 8)
11:35:21 INFO: Generating Helas calls for FKS process: u~ d > w- > t~ b [ all = QCD ] (7 / 8)
11:35:21 INFO: Processing color information for process: u~ d > w- > t~ b [ all = QCD ]
11:35:21 INFO: Generating Helas calls for FKS process: c~ s > w- > t~ b [ all = QCD ] (8 / 8)
11:35:22 INFO: ... Done
11:35:22 INFO: Writing files in P0_udx_wp_tbx (1 / 4)
11:35:22 INFO: Creating files in directory V0_udx_wp_tbx
11:35:22 ALOHA: aloha creates FFV1 set of routines with options: L2,P0
11:35:22 INFO: Computing diagram color coefficients
11:35:22 INFO: Drawing loop Feynman diagrams for Process: u d~ > w+ > t b~ QED<=2 [ all = QCD ]
11:35:22 INFO: Generating born Feynman diagrams for Process: u d~ > w+ > t b~ QED<=2 [ all = QCD ]
11:35:23 INFO: Writing files in P0_dxu_wp_tbx (2 / 4)
11:35:23 INFO: Creating files in directory V0_dxu_wp_tbx
11:35:23 INFO: Computing diagram color coefficients
11:35:23 INFO: Drawing loop Feynman diagrams for Process: d~ u > w+ > t b~ QED<=2 [ all = QCD ]
11:35:23 INFO: Generating born Feynman diagrams for Process: d~ u > w+ > t b~ QED<=2 [ all = QCD ]
11:35:24 INFO: Writing files in P0_dux_wm_txb (3 / 4)
11:35:24 INFO: Creating files in directory V0_dux_wm_txb
11:35:24 INFO: Computing diagram color coefficients
11:35:24 INFO: Drawing loop Feynman diagrams for Process: d u~ > w- > t~ b QED<=2 [ all = QCD ]
11:35:24 INFO: Generating born Feynman diagrams for Process: d u~ > w- > t~ b QED<=2 [ all = QCD ]
11:35:25 INFO: Writing files in P0_uxd_wm_txb (4 / 4)
11:35:25 INFO: Creating files in directory V0_uxd_wm_txb
11:35:25 INFO: Computing diagram color coefficients
11:35:25 INFO: Drawing loop Feynman diagrams for Process: u~ d > w- > t~ b QED<=2 [ all = QCD ]
11:35:25 INFO: Generating born Feynman diagrams for Process: u~ d > w- > t~ b QED<=2 [ all = QCD ]
11:35:26 History written to /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/proc_card_mg5.dat
11:35:26 ALOHA: aloha creates FFV2 routines
11:35:26 ALOHA: aloha creates FFV1 routines
11:35:27 save configuration file to /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:35:27 INFO: Use Fortran compiler gfortran
11:35:27 INFO: Use c++ compiler g++
11:35:27 INFO: Generate jpeg diagrams
11:35:35 INFO: Generate web pages
11:35:44 Type "launch" to generate events from this process, or see
11:35:44 /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/README
11:35:44 Run "open index.html" to see more information about this process.
11:35:44 quit
11:35:44 Py:MadGraphUtils INFO Finished process generation at Tue Aug 21 11:35:44 2018
11:35:44 Py:MadGraphUtils WARNING Path for option fastjet not found in original config
11:35:44 Py:MadGraphUtils INFO Modifying config paths to avoid use of afs:
11:35:44 Py:MadGraphUtils INFO {'ninja': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib', 'lhapdf': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//lhapdf/6.1.5/x86_64-slc6-gcc47-opt/bin/lhapdf-config', 'syscalc_path': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//syscalc/1.1.4/x86_64-slc6-gcc47-opt', 'collier': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//collier/1.1/x86_64-slc6-gcc47-opt'}
11:35:44 Py:MadGraphUtils INFO Copying default run_card.dat from PROCNLO_loop_sm-no_b_mass_0/Cards/run_card.dat
11:35:44 Py:MadGraphUtils INFO Found old card at run_card.tmp.dat
11:35:44 Py:MadGraphUtils INFO Found old param card at PROCNLO_loop_sm-no_b_mass_0/Cards/param_card.dat
11:35:44 Py:MadGraphUtils INFO proc_card:
11:35:44
11:35:44 import model loop_sm-no_b_mass
11:35:44 define l+ = l+ ta+
11:35:44 define l- = l- ta-
11:35:44 generate p p > w+ > t b~ [QCD]
11:35:44 add process p p > w- > t~ b [QCD]
11:35:44 output -f
11:35:44 Py:MadGraphUtils INFO run_card:
11:35:44 #***********************************************************************
11:35:44 # MadGraph5_aMC@NLO *
11:35:44 # *
11:35:44 # run_card.dat aMC@NLO *
11:35:44 # *
11:35:44 # This file is used to set the parameters of the run. *
11:35:44 # *
11:35:44 # Some notation/conventions: *
11:35:44 # *
11:35:44 # Lines starting with a hash (#) are info or comments *
11:35:44 # *
11:35:44 # mind the format: value = variable ! comment *
11:35:44 # *
11:35:44 # Some of the values of variables can be list. These can either be *
11:35:44 # comma or space separated. *
11:35:44 # *
11:35:44 # To display additional parameter, you can use the command: *
11:35:44 # update to_full *
11:35:44 #***********************************************************************
11:35:44 #
11:35:44 #*******************
11:35:44 # Running parameters
11:35:44 #*******************
11:35:44 #
11:35:44 #***********************************************************************
11:35:44 # Tag name for the run (one word) *
11:35:44 #***********************************************************************
11:35:44 tag_1 = run_tag ! name of the run
11:35:44 #***********************************************************************
11:35:44 # Number of LHE events (and their normalization) and the required *
11:35:44 # (relative) accuracy on the Xsec. *
11:35:44 # These values are ignored for fixed order runs *
11:35:44 #***********************************************************************
11:35:44 5500 = nevents ! Number of unweighted events requested
11:35:44 -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
11:35:44 -1 = nevt_job! Max number of events per job in event generation.
11:35:44 ! (-1= no split).
11:35:44 #***********************************************************************
11:35:44 # Normalize the weights of LHE events such that they sum or average to *
11:35:44 # the total cross section *
11:35:44 #***********************************************************************
11:35:44 average = event_norm ! valid settings: average, sum, bias
11:35:44 #***********************************************************************
11:35:44 # Number of points per itegration channel (ignored for aMC@NLO runs) *
11:35:44 #***********************************************************************
11:35:44 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
11:35:44 ! number of points and iter. below)
11:35:44 # These numbers are ignored except if req_acc_FO is equal to -1
11:35:44 5000 = npoints_FO_grid ! number of points to setup grids
11:35:44 4 = niters_FO_grid ! number of iter. to setup grids
11:35:44 10000 = npoints_FO ! number of points to compute Xsec
11:35:44 6 = niters_FO ! number of iter. to compute Xsec
11:35:44 #***********************************************************************
11:35:44 # Random number seed *
11:35:44 #***********************************************************************
11:35:44 1 = iseed ! rnd seed (0=assigned automatically=default))
11:35:44 #***********************************************************************
11:35:44 # Collider type and energy *
11:35:44 #***********************************************************************
11:35:44 1 = lpp1 ! beam 1 type (0 = no PDF)
11:35:44 1 = lpp2 ! beam 2 type (0 = no PDF)
11:35:44 6500 = ebeam1 ! beam 1 energy in GeV
11:35:44 6500 = ebeam2 ! beam 2 energy in GeV
11:35:44 #***********************************************************************
11:35:44 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
11:35:44 #***********************************************************************
11:35:44 lhapdf = pdlabel ! PDF set
11:35:44 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
11:35:44 ! numbers for central PDF sets are allowed. Can be a list;
11:35:44 ! PDF sets beyond the first are included via reweighting.
11:35:44 #***********************************************************************
11:35:44 # Include the NLO Monte Carlo subtr. terms for the following parton *
11:35:44 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
11:35:44 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
11:35:44 #***********************************************************************
11:35:44 HERWIGPP = parton_shower
11:35:44 1.0 = shower_scale_factor ! multiply default shower starting
11:35:44 ! scale by this factor
11:35:44 #***********************************************************************
11:35:44 # Renormalization and factorization scales *
11:35:44 # (Default functional form for the non-fixed scales is the sum of *
11:35:44 # the transverse masses divided by two of all final state particles *
11:35:44 # and partons. This can be changed in SubProcesses/set_scales.f or via *
11:35:44 # dynamical_scale_choice option) *
11:35:44 #***********************************************************************
11:35:44 True = fixed_ren_scale ! if .true. use fixed ren scale
11:35:44 True = fixed_fac_scale ! if .true. use fixed fac scale
11:35:44 172.5 = muR_ref_fixed ! fixed ren reference scale
11:35:44 172.5 = muF_ref_fixed ! fixed fact reference scale
11:35:44 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
11:35:44 ! dynamical choices. Can be a list; scale choices beyond the
11:35:44 ! first are included via reweighting
11:35:44 1.0 = muR_over_ref ! ratio of current muR over reference muR
11:35:44 1.0 = muF_over_ref ! ratio of current muF over reference muF
11:35:44 #***********************************************************************
11:35:44 # Reweight variables for scale dependence and PDF uncertainty *
11:35:44 #***********************************************************************
11:35:44 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
11:35:44 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
11:35:44 False = reweight_scale ! Reweight to get scale variation using the
11:35:44 ! rw_rscale and rw_fscale factors. Should be a list of
11:35:44 ! booleans of equal length to dynamical_scale_choice to
11:35:44 ! specify for which choice to include scale dependence.
11:35:44 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
11:35:44 ! list booleans of equal length to lhaid to specify for
11:35:44 ! which PDF set to include the uncertainties.
11:35:44 #***********************************************************************
11:35:44 # Store reweight information in the LHE file for off-line model- *
11:35:44 # parameter reweighting at NLO+PS accuracy *
11:35:44 #***********************************************************************
11:35:44 False = store_rwgt_info ! Store info for reweighting in LHE file
11:35:44 #***********************************************************************
11:35:44 # ickkw parameter: *
11:35:44 # 0: No merging *
11:35:44 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
11:35:44 # level. After showering an MLM-type merging should be applied as *
11:35:44 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
11:35:44 # 4: UNLOPS merging (with pythia8 only). No interface from within *
11:35:44 # MG5_aMC available, but available in Pythia8. *
11:35:44 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
11:35:44 #***********************************************************************
11:35:44 0 = ickkw
11:35:44 #***********************************************************************
11:35:44 #
11:35:44 #***********************************************************************
11:35:44 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
11:35:44 # written in the LHE event file *
11:35:44 #***********************************************************************
11:35:44 50 = bwcutoff
11:35:44 #***********************************************************************
11:35:44 # Cuts on the jets. Jet clustering is performed by FastJet. *
11:35:44 # - When matching to a parton shower, these generation cuts should be *
11:35:44 # considerably softer than the analysis cuts. *
11:35:44 # - More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
11:35:44 0.4 = jetradius ! The radius parameter for the jet algorithm
11:35:44 0.1 = ptj ! Min jet transverse momentum
11:35:44 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:44 #***********************************************************************
11:35:44 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 0.0 = ptl ! Min lepton transverse momentum
11:35:44 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:44 0.0 = drll ! Min distance between opposite sign lepton pairs
11:35:44 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
11:35:44 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
11:35:44 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
11:35:44 #***********************************************************************
11:35:44 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
11:35:44 # all the other parameters are ignored. *
11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 20.0 = ptgmin ! Min photon transverse momentum
11:35:44 -1.0 = etagamma ! Max photon abs(pseudo-rap)
11:35:44 0.4 = R0gamma ! Radius of isolation code
11:35:44 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
11:35:44 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
11:35:44 True = isoEM ! isolate photons from EM energy (photons and leptons)
11:35:44 #***********************************************************************
11:35:44 # Cuts associated to MASSIVE particles identified by their PDG codes. *
11:35:44 # All cuts are applied to both particles and anti-particles, so use *
11:35:44 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
11:35:44 # {6:100, 25:200} for multiple particles *
11:35:44 #***********************************************************************
11:35:44 {} = pt_min_pdg ! Min pT for a massive particle
11:35:44 {} = pt_max_pdg ! Max pT for a massive particle
11:35:44 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
11:35:44 #***********************************************************************
11:35:44 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
11:35:44 #***********************************************************************
11:35:44 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
11:35:44 #***********************************************************************
11:35:44 Py:MadGraphUtils INFO param_card:
11:35:44 ######################################################################
11:35:44 ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
11:35:44 ######################################################################
11:35:44 ## ##
11:35:44 ## Width set on Auto will be computed following the information ##
11:35:44 ## present in the decay.py files of the model. ##
11:35:44 ## See arXiv:1402.1178 for more details. ##
11:35:44 ## ##
11:35:44 ######################################################################
11:35:44
11:35:44 ###################################
11:35:44 ## INFORMATION FOR LOOP
11:35:44 ###################################
11:35:44 Block loop
11:35:44 1 9.118800e+01 # MU_R
11:35:44
11:35:44 ###################################
11:35:44 ## INFORMATION FOR MASS
11:35:44 ###################################
11:35:44 Block mass
11:35:44 6 1.725000e+02 # MT
11:35:44 15 1.776820e+00 # MTA
11:35:44 23 9.118760e+01 # MZ
11:35:44 25 1.250000e+02 # MH
11:35:44 ## Dependent parameters, given by model restrictions.
11:35:44 ## Those values should be edited following the
11:35:44 ## analytical expression. MG5 ignores those values
11:35:44 ## but they are important for interfacing the output of MG5
11:35:44 ## to external program such as Pythia.
11:35:44 1 0.000000 # d : 0.0
11:35:44 2 0.000000 # u : 0.0
11:35:44 3 0.000000 # s : 0.0
11:35:44 4 0.000000 # c : 0.0
11:35:44 5 0.000000 # b : 0.0
11:35:44 11 0.000000 # e- : 0.0
11:35:44 12 0.000000 # ve : 0.0
11:35:44 13 0.000000 # mu- : 0.0
11:35:44 14 0.000000 # vm : 0.0
11:35:44 16 0.000000 # vt : 0.0
11:35:44 21 0.000000 # g : 0.0
11:35:44 22 0.000000 # a : 0.0
11:35:44 24 80.399000 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
11:35:44
11:35:44 ###################################
11:35:44 ## INFORMATION FOR SMINPUTS
11:35:44 ###################################
11:35:44 Block sminputs
11:35:44 1 1.323489e+02 # aEWM1
11:35:44 2 1.166370e-05 # Gf
11:35:44 3 1.184000e-01 # aS
11:35:44
11:35:44 ###################################
11:35:44 ## INFORMATION FOR YUKAWA
11:35:44 ###################################
11:35:44 Block yukawa
11:35:44 6 1.725000e+02 # ymt
11:35:44 15 1.777000e+00 # ymtau
11:35:44
11:35:44 ###################################
11:35:44
11:35:44 #*************************
11:35:44 # Decay widths *
11:35:44 #*************************
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 1 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 2 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 3 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 4 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 5 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 6 1.320000e+00
11:35:44 # BR NDA ID1 ID2 ...
11:35:44 1.000000e+00 2 5 24 # 1.32
11:35:44 DECAY -6 1.320000e+00
11:35:44 # BR NDA ID1 ID2 ...
11:35:44 1.000000e+00 2 -5 -24 # 1.32
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 11 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 12 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 13 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 14 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 15 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 16 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 21 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 22 0.000000e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 23 2.495200e+00
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY -24 2.085000e+00
11:35:44 # BR NDA ID1 ID2 ...
11:35:44 3.3370000e-01 2 1 -2
11:35:44 3.337000e-01 2 3 -4
11:35:44 1.082000e-01 2 11 -12
11:35:44 1.082000e-01 2 13 -14
11:35:44 1.082000e-01 2 15 -16
11:35:44 DECAY 24 2.085000e+00
11:35:44 # BR NDA ID1 ID2 ...
11:35:44 3.337000e-01 2 -1 2
11:35:44 3.337000e-01 2 -3 4
11:35:44 1.082000e-01 2 -11 12
11:35:44 1.082000e-01 2 -13 14
11:35:44 1.082000e-01 2 -15 16
11:35:44 #
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 25 6.382339e-03
11:35:44 #
11:35:44 # PDG Width
11:35:44 DECAY 82 0.000000e+00
11:35:44 Py:MadGraphUtils INFO Found MadGraph/MadGraph5_aMC@NLO version 2.6.2
11:35:44 Py:MadGraphUtils INFO Started generating gridpack at Tue Aug 21 11:35:44 2018
11:35:44 Py:MadGraphUtils WARNING >>>>>> THIS KIND OF JOB SHOULD ONLY BE RUN LOCALLY - NOT IN GRID JOBS <<<<<<
11:35:44 Py:MadGraphUtils INFO Setting nevents = 1.
11:35:44 Py:MadGraphUtils INFO Setting req_acc = 0.001000
11:35:44 Py:MadGraphUtils INFO run_card.dat: run_card.dat
11:35:44 #***********************************************************************
11:35:44 # MadGraph5_aMC@NLO *
11:35:44 # *
11:35:44 # run_card.dat aMC@NLO *
11:35:44 # *
11:35:44 # This file is used to set the parameters of the run. *
11:35:44 # *
11:35:44 # Some notation/conventions: *
11:35:44 # *
11:35:44 # Lines starting with a hash (#) are info or comments *
11:35:44 # *
11:35:44 # mind the format: value = variable ! comment *
11:35:44 # *
11:35:44 # Some of the values of variables can be list. These can either be *
11:35:44 # comma or space separated. *
11:35:44 # *
11:35:44 # To display additional parameter, you can use the command: *
11:35:44 # update to_full *
11:35:44 #***********************************************************************
11:35:44 #
11:35:44 #*******************
11:35:44 # Running parameters
11:35:44 #*******************
11:35:44 #
11:35:44 #***********************************************************************
11:35:44 # Tag name for the run (one word) *
11:35:44 #***********************************************************************
11:35:44 tag_1 = run_tag ! name of the run
11:35:44 #***********************************************************************
11:35:44 # Number of LHE events (and their normalization) and the required *
11:35:44 # (relative) accuracy on the Xsec. *
11:35:44 # These values are ignored for fixed order runs *
11:35:44 #***********************************************************************
11:35:44 1 = nevents ! Number of unweighted events requested
11:35:44 0.001000 = req_acc ! Required accuracy (-1=auto determined from nevents)
11:35:44 -1 = nevt_job! Max number of events per job in event generation.
11:35:44 ! (-1= no split).
11:35:44 #***********************************************************************
11:35:44 # Normalize the weights of LHE events such that they sum or average to *
11:35:44 # the total cross section *
11:35:44 #***********************************************************************
11:35:44 average = event_norm ! valid settings: average, sum, bias
11:35:44 #***********************************************************************
11:35:44 # Number of points per itegration channel (ignored for aMC@NLO runs) *
11:35:44 #***********************************************************************
11:35:44 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
11:35:44 ! number of points and iter. below)
11:35:44 # These numbers are ignored except if req_acc_FO is equal to -1
11:35:44 5000 = npoints_FO_grid ! number of points to setup grids
11:35:44 4 = niters_FO_grid ! number of iter. to setup grids
11:35:44 10000 = npoints_FO ! number of points to compute Xsec
11:35:44 6 = niters_FO ! number of iter. to compute Xsec
11:35:44 #***********************************************************************
11:35:44 # Random number seed *
11:35:44 #***********************************************************************
11:35:44 1 = iseed ! rnd seed (0=assigned automatically=default))
11:35:44 #***********************************************************************
11:35:44 # Collider type and energy *
11:35:44 #***********************************************************************
11:35:44 1 = lpp1 ! beam 1 type (0 = no PDF)
11:35:44 1 = lpp2 ! beam 2 type (0 = no PDF)
11:35:44 6500 = ebeam1 ! beam 1 energy in GeV
11:35:44 6500 = ebeam2 ! beam 2 energy in GeV
11:35:44 #***********************************************************************
11:35:44 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
11:35:44 #***********************************************************************
11:35:44 lhapdf = pdlabel ! PDF set
11:35:44 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
11:35:44 ! numbers for central PDF sets are allowed. Can be a list;
11:35:44 ! PDF sets beyond the first are included via reweighting.
11:35:44 #***********************************************************************
11:35:44 # Include the NLO Monte Carlo subtr. terms for the following parton *
11:35:44 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
11:35:44 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
11:35:44 #***********************************************************************
11:35:44 HERWIGPP = parton_shower
11:35:44 1.0 = shower_scale_factor ! multiply default shower starting
11:35:44 ! scale by this factor
11:35:44 #***********************************************************************
11:35:44 # Renormalization and factorization scales *
11:35:44 # (Default functional form for the non-fixed scales is the sum of *
11:35:44 # the transverse masses divided by two of all final state particles *
11:35:44 # and partons. This can be changed in SubProcesses/set_scales.f or via *
11:35:44 # dynamical_scale_choice option) *
11:35:44 #***********************************************************************
11:35:44 True = fixed_ren_scale ! if .true. use fixed ren scale
11:35:44 True = fixed_fac_scale ! if .true. use fixed fac scale
11:35:44 172.5 = muR_ref_fixed ! fixed ren reference scale
11:35:44 172.5 = muF_ref_fixed ! fixed fact reference scale
11:35:44 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
11:35:44 ! dynamical choices. Can be a list; scale choices beyond the
11:35:44 ! first are included via reweighting
11:35:44 1.0 = muR_over_ref ! ratio of current muR over reference muR
11:35:44 1.0 = muF_over_ref ! ratio of current muF over reference muF
11:35:44 #***********************************************************************
11:35:44 # Reweight variables for scale dependence and PDF uncertainty *
11:35:44 #***********************************************************************
11:35:44 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
11:35:44 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
11:35:44 False = reweight_scale ! Reweight to get scale variation using the
11:35:44 ! rw_rscale and rw_fscale factors. Should be a list of
11:35:44 ! booleans of equal length to dynamical_scale_choice to
11:35:44 ! specify for which choice to include scale dependence.
11:35:44 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
11:35:44 ! list booleans of equal length to lhaid to specify for
11:35:44 ! which PDF set to include the uncertainties.
11:35:44 #***********************************************************************
11:35:44 # Store reweight information in the LHE file for off-line model- *
11:35:44 # parameter reweighting at NLO+PS accuracy *
11:35:44 #***********************************************************************
11:35:44 False = store_rwgt_info ! Store info for reweighting in LHE file
11:35:44 #***********************************************************************
11:35:44 # ickkw parameter: *
11:35:44 # 0: No merging *
11:35:44 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
11:35:44 # level. After showering an MLM-type merging should be applied as *
11:35:44 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
11:35:44 # 4: UNLOPS merging (with pythia8 only). No interface from within *
11:35:44 # MG5_aMC available, but available in Pythia8. *
11:35:44 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
11:35:44 #***********************************************************************
11:35:44 0 = ickkw
11:35:44 #***********************************************************************
11:35:44 #
11:35:44 #***********************************************************************
11:35:44 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
11:35:44 # written in the LHE event file *
11:35:44 #***********************************************************************
11:35:44 50 = bwcutoff
11:35:44 #***********************************************************************
11:35:44 # Cuts on the jets. Jet clustering is performed by FastJet. *
11:35:44 # - When matching to a parton shower, these generation cuts should be *
11:35:44 # considerably softer than the analysis cuts. *
11:35:44 # - More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
11:35:44 0.4 = jetradius ! The radius parameter for the jet algorithm
11:35:44 0.1 = ptj ! Min jet transverse momentum
11:35:44 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:44 #***********************************************************************
11:35:44 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 0.0 = ptl ! Min lepton transverse momentum
11:35:44 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:44 0.0 = drll ! Min distance between opposite sign lepton pairs
11:35:44 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
11:35:44 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
11:35:44 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
11:35:44 #***********************************************************************
11:35:44 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
11:35:44 # all the other parameters are ignored. *
11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:44 #***********************************************************************
11:35:44 20.0 = ptgmin ! Min photon transverse momentum
11:35:44 -1.0 = etagamma ! Max photon abs(pseudo-rap)
11:35:44 0.4 = R0gamma ! Radius of isolation code
11:35:44 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
11:35:44 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
11:35:44 True = isoEM ! isolate photons from EM energy (photons and leptons)
11:35:44 #***********************************************************************
11:35:44 # Cuts associated to MASSIVE particles identified by their PDG codes. *
11:35:44 # All cuts are applied to both particles and anti-particles, so use *
11:35:44 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
11:35:44 # {6:100, 25:200} for multiple particles *
11:35:44 #***********************************************************************
11:35:44 {} = pt_min_pdg ! Min pT for a massive particle
11:35:44 {} = pt_max_pdg ! Max pT for a massive particle
11:35:44 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
11:35:44 #***********************************************************************
11:35:44 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
11:35:44 #***********************************************************************
11:35:44 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
11:35:44 #***********************************************************************
11:35:44 Py:MadGraphUtils INFO Run schan will be performed in mode 2 with 8 jobs in parallel.
11:35:44 Py:MadGraphUtils INFO Decaying with MadSpin. Moving card (madspin_card.dat) into place.
11:35:44 Py:MadGraphUtils INFO creating local LHAPDF dir: MGC_LHAPDF/
11:35:44 Py:MadGraphUtils INFO Found LHAPDF ID=260000, name=NNPDF30_nlo_as_0118!
11:35:44 Py:MadGraphUtils INFO linking /cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current//NNPDF30_nlo_as_0118 --> MGC_LHAPDF/NNPDF30_nlo_as_0118
11:35:44 Py:MadGraphUtils INFO linking /cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current//pdfsets.index --> MGC_LHAPDF/pdfsets.index
11:35:44 Py:MadGraphUtils INFO linking /cvmfs/atlas.cern.ch/repo/sw/Generators/lhapdfsets/current//lhapdf.conf --> MGC_LHAPDF/lhapdf.conf
11:35:44 Py:MadGraphUtils INFO Path to LHAPDF install dir:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt/share/LHAPDF
11:35:44 Py:MadGraphUtils INFO Path to LHAPDF data dir: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/MGC_LHAPDF
11:35:46 Willing to acquire file lhapdf-config from
11:35:46 /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/share/lhapdf-config
11:35:46 copy lhapdf-config
11:35:48 Py:MadGraphUtils INFO lhapdf-config --version: 6.1.5
11:35:49 Py:MadGraphUtils INFO lhapdf-config --prefix: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt
11:35:49 Py:MadGraphUtils INFO lhapdf-config --libdir: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt/lib
11:35:50 Py:MadGraphUtils INFO lhapdf-config --datadir: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/MGC_LHAPDF
11:35:50 Py:MadGraphUtils INFO lhapdf-config --pdfsets-path:
11:35:50 Py:MadGraphUtils INFO Setting lhapdf = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/lhapdf-config in PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:35:50 Py:MadGraphUtils INFO Creating links for LHAPDF
11:35:50 Py:MadGraphUtils INFO Available PDFs are:
11:35:50 Py:MadGraphUtils INFO ['NNPDF30_nlo_as_0118', 'lhapdf.conf', 'pdfsets.index']
11:35:50 Py:MadGraphUtils INFO For your information, the libraries available are (should include LHAPDF):
11:35:50 Py:MadGraphUtils INFO ['PDFsets', 'Pdfdata', 'libcts.a', 'libiregi.a', 'mpmodule.mod']
11:35:50 Py:MadGraphUtils INFO Path to fastjet install dir:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc47-opt
11:35:52 Willing to acquire file fastjet-config from
11:35:52 /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/share/fastjet-config
11:35:52 copy fastjet-config
11:35:52 Py:MadGraphUtils INFO fastjet-config --version: 3.0.3
11:35:52 Py:MadGraphUtils INFO fastjet-config --prefix: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc47-opt
11:35:52 Py:MadGraphUtils INFO Setting fastjet = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/fastjet-config in PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:35:52 Py:MadGraphUtils INFO Now I will hack the make files a bit. Apologies, but there seems to be no good way around this.
11:35:52 Py:MadGraphUtils INFO Configuring the fancy gfortran compiler instead of g77 / f77
11:35:52 Py:MadGraphUtils INFO Make file hacking complete.
11:35:52 Py:MadGraphUtils INFO proc_card:
11:35:52
11:35:52 import model loop_sm-no_b_mass
11:35:52 define l+ = l+ ta+
11:35:52 define l- = l- ta-
11:35:52 generate p p > w+ > t b~ [QCD]
11:35:52 add process p p > w- > t~ b [QCD]
11:35:52 output -f
11:35:52 Py:MadGraphUtils INFO run_card:
11:35:52 #***********************************************************************
11:35:52 # MadGraph5_aMC@NLO *
11:35:52 # *
11:35:52 # run_card.dat aMC@NLO *
11:35:52 # *
11:35:52 # This file is used to set the parameters of the run. *
11:35:52 # *
11:35:52 # Some notation/conventions: *
11:35:52 # *
11:35:52 # Lines starting with a hash (#) are info or comments *
11:35:52 # *
11:35:52 # mind the format: value = variable ! comment *
11:35:52 # *
11:35:52 # Some of the values of variables can be list. These can either be *
11:35:52 # comma or space separated. *
11:35:52 # *
11:35:52 # To display additional parameter, you can use the command: *
11:35:52 # update to_full *
11:35:52 #***********************************************************************
11:35:52 #
11:35:52 #*******************
11:35:52 # Running parameters
11:35:52 #*******************
11:35:52 #
11:35:52 #***********************************************************************
11:35:52 # Tag name for the run (one word) *
11:35:52 #***********************************************************************
11:35:52 tag_1 = run_tag ! name of the run
11:35:52 #***********************************************************************
11:35:52 # Number of LHE events (and their normalization) and the required *
11:35:52 # (relative) accuracy on the Xsec. *
11:35:52 # These values are ignored for fixed order runs *
11:35:52 #***********************************************************************
11:35:52 1 = nevents ! Number of unweighted events requested
11:35:52 0.001000 = req_acc ! Required accuracy (-1=auto determined from nevents)
11:35:52 -1 = nevt_job! Max number of events per job in event generation.
11:35:52 ! (-1= no split).
11:35:52 #***********************************************************************
11:35:52 # Normalize the weights of LHE events such that they sum or average to *
11:35:52 # the total cross section *
11:35:52 #***********************************************************************
11:35:52 average = event_norm ! valid settings: average, sum, bias
11:35:52 #***********************************************************************
11:35:52 # Number of points per itegration channel (ignored for aMC@NLO runs) *
11:35:52 #***********************************************************************
11:35:52 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
11:35:52 ! number of points and iter. below)
11:35:52 # These numbers are ignored except if req_acc_FO is equal to -1
11:35:52 5000 = npoints_FO_grid ! number of points to setup grids
11:35:52 4 = niters_FO_grid ! number of iter. to setup grids
11:35:52 10000 = npoints_FO ! number of points to compute Xsec
11:35:52 6 = niters_FO ! number of iter. to compute Xsec
11:35:52 #***********************************************************************
11:35:52 # Random number seed *
11:35:52 #***********************************************************************
11:35:52 1 = iseed ! rnd seed (0=assigned automatically=default))
11:35:52 #***********************************************************************
11:35:52 # Collider type and energy *
11:35:52 #***********************************************************************
11:35:52 1 = lpp1 ! beam 1 type (0 = no PDF)
11:35:52 1 = lpp2 ! beam 2 type (0 = no PDF)
11:35:52 6500 = ebeam1 ! beam 1 energy in GeV
11:35:52 6500 = ebeam2 ! beam 2 energy in GeV
11:35:52 #***********************************************************************
11:35:52 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
11:35:52 #***********************************************************************
11:35:52 lhapdf = pdlabel ! PDF set
11:35:52 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
11:35:52 ! numbers for central PDF sets are allowed. Can be a list;
11:35:52 ! PDF sets beyond the first are included via reweighting.
11:35:52 #***********************************************************************
11:35:52 # Include the NLO Monte Carlo subtr. terms for the following parton *
11:35:52 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
11:35:52 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
11:35:52 #***********************************************************************
11:35:52 HERWIGPP = parton_shower
11:35:52 1.0 = shower_scale_factor ! multiply default shower starting
11:35:52 ! scale by this factor
11:35:52 #***********************************************************************
11:35:52 # Renormalization and factorization scales *
11:35:52 # (Default functional form for the non-fixed scales is the sum of *
11:35:52 # the transverse masses divided by two of all final state particles *
11:35:52 # and partons. This can be changed in SubProcesses/set_scales.f or via *
11:35:52 # dynamical_scale_choice option) *
11:35:52 #***********************************************************************
11:35:52 True = fixed_ren_scale ! if .true. use fixed ren scale
11:35:52 True = fixed_fac_scale ! if .true. use fixed fac scale
11:35:52 172.5 = muR_ref_fixed ! fixed ren reference scale
11:35:52 172.5 = muF_ref_fixed ! fixed fact reference scale
11:35:52 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
11:35:52 ! dynamical choices. Can be a list; scale choices beyond the
11:35:52 ! first are included via reweighting
11:35:52 1.0 = muR_over_ref ! ratio of current muR over reference muR
11:35:52 1.0 = muF_over_ref ! ratio of current muF over reference muF
11:35:52 #***********************************************************************
11:35:52 # Reweight variables for scale dependence and PDF uncertainty *
11:35:52 #***********************************************************************
11:35:52 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
11:35:52 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
11:35:52 False = reweight_scale ! Reweight to get scale variation using the
11:35:52 ! rw_rscale and rw_fscale factors. Should be a list of
11:35:52 ! booleans of equal length to dynamical_scale_choice to
11:35:52 ! specify for which choice to include scale dependence.
11:35:52 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
11:35:52 ! list booleans of equal length to lhaid to specify for
11:35:52 ! which PDF set to include the uncertainties.
11:35:52 #***********************************************************************
11:35:52 # Store reweight information in the LHE file for off-line model- *
11:35:52 # parameter reweighting at NLO+PS accuracy *
11:35:52 #***********************************************************************
11:35:52 False = store_rwgt_info ! Store info for reweighting in LHE file
11:35:52 #***********************************************************************
11:35:52 # ickkw parameter: *
11:35:52 # 0: No merging *
11:35:52 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
11:35:52 # level. After showering an MLM-type merging should be applied as *
11:35:52 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
11:35:52 # 4: UNLOPS merging (with pythia8 only). No interface from within *
11:35:52 # MG5_aMC available, but available in Pythia8. *
11:35:52 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
11:35:52 #***********************************************************************
11:35:52 0 = ickkw
11:35:52 #***********************************************************************
11:35:52 #
11:35:52 #***********************************************************************
11:35:52 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
11:35:52 # written in the LHE event file *
11:35:52 #***********************************************************************
11:35:52 50 = bwcutoff
11:35:52 #***********************************************************************
11:35:52 # Cuts on the jets. Jet clustering is performed by FastJet. *
11:35:52 # - When matching to a parton shower, these generation cuts should be *
11:35:52 # considerably softer than the analysis cuts. *
11:35:52 # - More specific cuts can be specified in SubProcesses/cuts.f *
11:35:52 #***********************************************************************
11:35:52 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
11:35:52 0.4 = jetradius ! The radius parameter for the jet algorithm
11:35:52 0.1 = ptj ! Min jet transverse momentum
11:35:52 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:52 #***********************************************************************
11:35:52 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
11:35:52 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:52 #***********************************************************************
11:35:52 0.0 = ptl ! Min lepton transverse momentum
11:35:52 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
11:35:52 0.0 = drll ! Min distance between opposite sign lepton pairs
11:35:52 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
11:35:52 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
11:35:52 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
11:35:52 #***********************************************************************
11:35:52 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
11:35:52 # all the other parameters are ignored. *
11:35:52 # More specific cuts can be specified in SubProcesses/cuts.f *
11:35:52 #***********************************************************************
11:35:52 20.0 = ptgmin ! Min photon transverse momentum
11:35:52 -1.0 = etagamma ! Max photon abs(pseudo-rap)
11:35:52 0.4 = R0gamma ! Radius of isolation code
11:35:52 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
11:35:52 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
11:35:52 True = isoEM ! isolate photons from EM energy (photons and leptons)
11:35:52 #***********************************************************************
11:35:52 # Cuts associated to MASSIVE particles identified by their PDG codes. *
11:35:52 # All cuts are applied to both particles and anti-particles, so use *
11:35:52 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
11:35:52 # {6:100, 25:200} for multiple particles *
11:35:52 #***********************************************************************
11:35:52 {} = pt_min_pdg ! Min pT for a massive particle
11:35:52 {} = pt_max_pdg ! Max pT for a massive particle
11:35:52 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
11:35:52 #***********************************************************************
11:35:52 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
11:35:52 #***********************************************************************
11:35:52 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
11:35:52 #***********************************************************************
11:35:52 Py:MadGraphUtils INFO param_card:
11:35:52 ######################################################################
11:35:52 ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
11:35:52 ######################################################################
11:35:52 ## ##
11:35:52 ## Width set on Auto will be computed following the information ##
11:35:52 ## present in the decay.py files of the model. ##
11:35:52 ## See arXiv:1402.1178 for more details. ##
11:35:52 ## ##
11:35:52 ######################################################################
11:35:52
11:35:52 ###################################
11:35:52 ## INFORMATION FOR LOOP
11:35:52 ###################################
11:35:52 Block loop
11:35:52 1 9.118800e+01 # MU_R
11:35:52
11:35:52 ###################################
11:35:52 ## INFORMATION FOR MASS
11:35:52 ###################################
11:35:52 Block mass
11:35:52 6 1.725000e+02 # MT
11:35:52 15 1.776820e+00 # MTA
11:35:52 23 9.118760e+01 # MZ
11:35:52 25 1.250000e+02 # MH
11:35:52 ## Dependent parameters, given by model restrictions.
11:35:52 ## Those values should be edited following the
11:35:52 ## analytical expression. MG5 ignores those values
11:35:52 ## but they are important for interfacing the output of MG5
11:35:52 ## to external program such as Pythia.
11:35:52 1 0.000000 # d : 0.0
11:35:52 2 0.000000 # u : 0.0
11:35:52 3 0.000000 # s : 0.0
11:35:52 4 0.000000 # c : 0.0
11:35:52 5 0.000000 # b : 0.0
11:35:52 11 0.000000 # e- : 0.0
11:35:52 12 0.000000 # ve : 0.0
11:35:52 13 0.000000 # mu- : 0.0
11:35:52 14 0.000000 # vm : 0.0
11:35:52 16 0.000000 # vt : 0.0
11:35:52 21 0.000000 # g : 0.0
11:35:52 22 0.000000 # a : 0.0
11:35:52 24 80.399000 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
11:35:52
11:35:52 ###################################
11:35:52 ## INFORMATION FOR SMINPUTS
11:35:52 ###################################
11:35:52 Block sminputs
11:35:52 1 1.323489e+02 # aEWM1
11:35:52 2 1.166370e-05 # Gf
11:35:52 3 1.184000e-01 # aS
11:35:52
11:35:52 ###################################
11:35:52 ## INFORMATION FOR YUKAWA
11:35:52 ###################################
11:35:52 Block yukawa
11:35:52 6 1.725000e+02 # ymt
11:35:52 15 1.777000e+00 # ymtau
11:35:52
11:35:52 ###################################
11:35:52
11:35:52 #*************************
11:35:52 # Decay widths *
11:35:52 #*************************
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 1 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 2 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 3 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 4 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 5 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 6 1.320000e+00
11:35:52 # BR NDA ID1 ID2 ...
11:35:52 1.000000e+00 2 5 24 # 1.32
11:35:52 DECAY -6 1.320000e+00
11:35:52 # BR NDA ID1 ID2 ...
11:35:52 1.000000e+00 2 -5 -24 # 1.32
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 11 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 12 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 13 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 14 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 15 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 16 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 21 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 22 0.000000e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 23 2.495200e+00
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY -24 2.085000e+00
11:35:52 # BR NDA ID1 ID2 ...
11:35:52 3.3370000e-01 2 1 -2
11:35:52 3.337000e-01 2 3 -4
11:35:52 1.082000e-01 2 11 -12
11:35:52 1.082000e-01 2 13 -14
11:35:52 1.082000e-01 2 15 -16
11:35:52 DECAY 24 2.085000e+00
11:35:52 # BR NDA ID1 ID2 ...
11:35:52 3.337000e-01 2 -1 2
11:35:52 3.337000e-01 2 -3 4
11:35:52 1.082000e-01 2 -11 12
11:35:52 1.082000e-01 2 -13 14
11:35:52 1.082000e-01 2 -15 16
11:35:52 #
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 25 6.382339e-03
11:35:52 #
11:35:52 # PDG Width
11:35:52 DECAY 82 0.000000e+00
11:35:52 Py:MadGraphUtils INFO Setting automatic_html_opening = False
11:35:52 Py:MadGraphUtils INFO Setting run_mode = 2
11:35:52 Py:MadGraphUtils INFO Setting number of cores = 8 in Cards/amcatnlo_configuration.txt
11:35:52 Py:MadGraphUtils INFO amcatnlo_configuration.txt: Cards/amcatnlo_configuration.txt
11:35:52 ################################################################################
11:35:52 #
11:35:52 # Copyright (c) 2009 The MadGraph5_aMC@NLO Development team and Contributors
11:35:52 #
11:35:52 # This file is a part of the MadGraph5_aMC@NLO project, an application which
11:35:52 # automatically generates Feynman diagrams and matrix elements for arbitrary
11:35:52 # high-energy processes in the Standard Model and beyond.
11:35:52 #
11:35:52 # It is subject to the MadGraph5_aMC@NLO license which should accompany this
11:35:52 # distribution.
11:35:52 #
11:35:52 # For more information, visit madgraph.phys.ucl.ac.be and amcatnlo.web.cern.ch
11:35:52 #
11:35:52 ################################################################################
11:35:52 #
11:35:52 # This File contains some configuration variable for MadGraph/MadEvent
11:35:52 #
11:35:52 # Line starting by #! are comment and should remain commented
11:35:52 # Line starting with # should be uncommented if you want to modify the default
11:35:52 # value.
11:35:52 # Current value for all options can seen by typing "display options"
11:35:52 # after either ./bin/mg5_aMC or ./bin/madevent
11:35:52 #
11:35:52 # You can place this files in ~/.mg5/mg5_configuration.txt if you have more than
11:35:52 # one version of MG5.
11:35:52 #
11:35:52 ################################################################################
11:35:52 #! Prefered Fortran Compiler
11:35:52 #! If None: try to find g77 or gfortran on the system
11:35:52 #!
11:35:52 # fortran_compiler = None
11:35:52 # f2py_compiler = None
11:35:52
11:35:52
11:35:52 #! Prefered C++ Compiler
11:35:52 #! If None: try to find g++ or clang on the system
11:35:52 #!
11:35:52 # cpp_compiler = None
11:35:52
11:35:52 #! Prefered Text Editor
11:35:52 #! Default: use the shell default Editor
11:35:52 #! or try to find one available on the system
11:35:52 #! Be careful: Only shell based editor are allowed
11:35:52 # text_editor = None
11:35:52
11:35:52 #! Prefered WebBrower
11:35:52 #! If None: try to find one available on the system
11:35:52 # web_browser = None
11:35:52
11:35:52 #! Prefered PS viewer
11:35:52 #! If None: try to find one available on the system
11:35:52 # eps_viewer = None
11:35:52
11:35:52 #! Time allowed to answer question (if no answer takes default value)
11:35:52 #! 0: No time limit
11:35:52 # timeout = 60
11:35:52
11:35:52 #! Pythia8 path.
11:35:52 #! Defines the path to the pythia8 installation directory (i.e. the
11:35:52 #! on containing the lib, bin and include directories) .
11:35:52 #! If using a relative path, that starts from the mg5 directory
11:35:52 # pythia8_path = ./HEPTools/pythia8
11:35:52
11:35:52 #! MG5aMC_PY8_interface path
11:35:52 #! Defines the path of the C++ driver file that is used by MG5_aMC to
11:35:52 #! steer the Pythia8 shower.
11:35:52 #! Can be installed directly from within MG5_aMC with the following command:
11:35:52 #! MG5_aMC> install mg5amc_py8_interface
11:35:52 # mg5amc_py8_interface_path = ./HEPTools/MG5aMC_PY8_interface
11:35:52
11:35:52 #! Herwig++/Herwig7 paths
11:35:52 #! specify here the paths also to HepMC ant ThePEG
11:35:52 #! define the path to the herwig++, thepeg and hepmc directories.
11:35:52 #! paths can be absolute or relative from mg5 directory
11:35:52 #! WARNING: if Herwig7 has been installed with the bootstrap script,
11:35:52 #! then please set thepeg_path and hepmc_path to the same value as
11:35:52 #! hwpp_path
11:35:52 # hwpp_path =
11:35:52 # thepeg_path =
11:35:52 # hepmc_path =
11:35:52
11:35:52 #! Control when MG5 checks if he is up-to-date.
11:35:52 #! Enter the number of day between two check (0 means never)
11:35:52 #! A question is always asked before any update
11:35:52 # auto_update = 7
11:35:52
11:35:52 ################################################################################
11:35:52 # INFO FOR MADEVENT / aMC@NLO
11:35:52 ################################################################################
11:35:52 # If this file is in a MADEVENT Template. 'main directory' is the directory
11:35:52 # containing the SubProcesses directory. Otherwise this is the MadGraph5_aMC@NLO main
11:35:52 # directory (containing the directories madgraph and Template)
11:35:52
11:35:52 #! Allow/Forbid the automatic opening of the web browser (on the status page)
11:35:52 #! when launching MadEvent [True/False]
11:35:52 automatic_html_opening = False
11:35:52 #! allow notification of finished job in the notification center (Mac Only)
11:35:52 # notification_center = True
11:35:52
11:35:52
11:35:52 #! Default Running mode
11:35:52 #! 0: single machine/ 1: cluster / 2: multicore
11:35:52 run_mode = 2
11:35:52
11:35:52 #! Cluster Type [pbs|sge|condor|lsf|ge|slurm|htcaas|htcaas2] Use for cluster run only
11:35:52 #! And cluster queue (or partition for slurm)
11:35:52 #! And size of the cluster (some part of the code can adapt splitting accordingly)
11:35:52 # cluster_type = condor
11:35:52 # cluster_queue = madgraph
11:35:52 # cluster_size = 150
11:35:52
11:35:52 #! Path to a node directory to avoid direct writing on the central disk
11:35:52 #! Note that condor clusters avoid direct writing by default (therefore this
11:35:52 #! options does not affect condor clusters)
11:35:52 # cluster_temp_path = None
11:35:52
11:35:52 #! path to a node directory where local file can be found (typically pdf)
11:35:52 #! to avoid to send them to the node (if cluster_temp_path is on True or condor)
11:35:52 # cluster_local_path = None # example: /cvmfs/cp3.uclouvain.be/madgraph/
11:35:52
11:35:52 #! Cluster waiting time for status update
11:35:52 #! First number is when the number of waiting job is higher than the number
11:35:52 #! of running one (time in second). The second number is in the second case.
11:35:52 # cluster_status_update = 600 30
11:35:52
11:35:52 #! How to deal with failed submission (can occurs on cluster mode)
11:35:52 #! 0: crash, -1: print error, hangs the program up to manual instructions, N(>0) retry up to N times.
11:35:52 # cluster_nb_retry = 1
11:35:52
11:35:52 #! How much time to wait for the output file before resubmission/crash (filesystem can be very slow)
11:35:52 # cluster_retry_wait = 300
11:35:52
11:35:52 #! Nb_core to use (None = all) This is use only for multicore run
11:35:52 #! This correspond also to the number core used for code compilation for cluster mode
11:35:52 nb_core = 8
11:35:52
11:35:52 #! Pythia-PGS Package
11:35:52 #! relative path start from main directory
11:35:52 # pythia-pgs_path = ./pythia-pgs
11:35:52
11:35:52 #! Delphes Package
11:35:52 #! relative path start from main directory
11:35:52 # delphes_path = ./Delphes
11:35:52
11:35:52 #! MadAnalysis4 fortran-based package [for basic analysis]
11:35:52 #! relative path start from main directory
11:35:52 # madanalysis_path = ./MadAnalysis
11:35:52
11:35:52 #! MadAnalysis5 python-based Package [For advanced analysis]
11:35:52 #! relative path start from main directory
11:35:52 # madanalysis5_path = ./HEPTools/madanalysis5/madanalysis5
11:35:52
11:35:52 #! ExRootAnalysis Package
11:35:52 #! relative path start from main directory
11:35:52 # exrootanalysis_path = ./ExRootAnalysis
11:35:52
11:35:52 #! TOPDRAWER PATH
11:35:52 #! Path to the directory containing td executables
11:35:52 #! relative path start from main directory
11:35:52 # td_path = ./td
11:35:52
11:35:52 #! lhapdf-config
11:35:52 #! If None: try to find one available on the system
11:35:52 lhapdf = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/lhapdf-config
11:35:52
11:35:52 #! fastjet-config
11:35:52 #! If None: try to find one available on the system
11:35:52 fastjet = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/fastjet-config
11:35:52
11:35:52 #! MCatNLO-utilities
11:35:52 #! relative path starting from main directory
11:35:52 # MCatNLO-utilities_path = ./MCatNLO-utilities
11:35:52
11:35:52 #! Set what OLP to use for the loop ME generation
11:35:52 # OLP = MadLoop
11:35:52
11:35:52 #! Set the PJFRy++ directory containing pjfry's library
11:35:52 #! if auto: try to find it automatically on the system (default)
11:35:52 #! if '' or None: disabling pjfry
11:35:52 #! if pjfry=/PATH/TO/pjfry/lib: use that specific installation path for PJFry++
11:35:52 # pjfry = auto
11:35:52
11:35:52 #! Set the Golem95 directory containing golem's library
11:35:52 #! It only supports version higher than 1.3.0
11:35:52 #! if auto: try to find it automatically on the system (default)
11:35:52 #! if '' or None: disabling Golem95
11:35:52 #! if golem=/PATH/TO/golem/lib: use that speficif installation path for Golem95
11:35:52 # golem = auto
11:35:52
11:35:52 #! Set the samurai directory containing samurai's library
11:35:52 #! It only supports version higher than 2.0.0
11:35:52 #! if auto: try to find it automatically on the system (default)
11:35:52 #! if '' or None: disabling samurai
11:35:52 #! if samurai=/PATH/TO/samurai/lib: use that specific installation path for samurai
11:35:52 # samurai = None
11:35:52
11:35:52 #! Set the Ninja directory containing ninja's library
11:35:52 #! if '' or None: disabling ninja
11:35:52 #! if ninja=/PATH/TO/ninja/lib: use that specific installation path for ninja
11:35:52 ninja = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib
11:35:52
11:35:52 #! Set the COLLIER directory containing COLLIER's library
11:35:52 #! if '' or None: disabling COLLIER
11:35:52 #! if ninja=/PATH/TO/ninja/lib: use that specific installation path for COLLIER
11:35:52 # Note that it is necessary that you have generated a static library for COLLIER
11:35:52 collier = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//collier/1.1/x86_64-slc6-gcc47-opt
11:35:52
11:35:52 #! Set how MadLoop dependencies (such as CutTools) should be handled
11:35:52 #! > external : ML5 places a link to the MG5_aMC-wide libraries
11:35:52 #! > internal : ML5 copies all dependencies in the output so that it is independent
11:35:52 #! > environment_paths : ML5 searches for the dependencies in your environment path
11:35:52 output_dependencies = internal
11:35:52
11:35:52
11:35:52 #! SysCalc PATH
11:35:52 #! Path to the directory containing syscalc executables
11:35:52 #! relative path start from main directory
11:35:52 syscalc_path = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//syscalc/1.1.4/x86_64-slc6-gcc47-opt
11:35:52
11:35:52 #! Absolute paths to config scripts in the bin directories for APPLgrid and aMCFast.
11:35:52 # applgrid = applgrid-config
11:35:52 # amcfast = amcfast-config
11:35:52
11:35:52
11:35:52
11:35:52 # MG5 MAIN DIRECTORY
11:35:52 mg5_path = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt
11:35:52 Py:MadGraphUtils INFO Removing Cards/shower_card.dat to ensure we get parton level events only
11:35:52 Py:MadGraphUtils INFO "rw_rscale" = 1.0, 2.0, 0.5
11:35:52 Py:MadGraphUtils INFO "ptl" = 0.0
11:35:52 Py:MadGraphUtils INFO "xn" = 1.0
11:35:52 Py:MadGraphUtils INFO "ptj" = 0.1
11:35:52 Py:MadGraphUtils INFO "niters_FO" = 6
11:35:52 Py:MadGraphUtils INFO "nevt_job" = -1
11:35:52 Py:MadGraphUtils INFO "drll" = 0.0
11:35:52 Py:MadGraphUtils INFO "fixed_ren_scale" = True
11:35:52 Py:MadGraphUtils INFO "npoints_FO_grid" = 5000
11:35:52 Py:MadGraphUtils INFO "iseed" = 1
11:35:52 Py:MadGraphUtils INFO "muR_ref_fixed" = 172.5
11:35:52 Py:MadGraphUtils INFO "dynamical_scale_choice" = -1
11:35:52 Py:MadGraphUtils INFO "R0gamma" = 0.4
11:35:52 Py:MadGraphUtils INFO "reweight_scale" = False
11:35:52 Py:MadGraphUtils INFO "store_rwgt_info" = False
11:35:52 Py:MadGraphUtils INFO "lpp1" = 1
11:35:52 Py:MadGraphUtils INFO "nevents" = 1
11:35:52 Py:MadGraphUtils INFO "pt_max_pdg" = {}
11:35:52 Py:MadGraphUtils INFO "ickkw" = 0
11:35:52 Py:MadGraphUtils INFO "iappl" = 0
11:35:52 Py:MadGraphUtils INFO "etal" = -1.0
11:35:52 Py:MadGraphUtils INFO "epsgamma" = 1.0
11:35:52 Py:MadGraphUtils INFO "req_acc_FO" = 0.01
11:35:52 Py:MadGraphUtils INFO "pt_min_pdg" = {}
11:35:52 Py:MadGraphUtils INFO "etagamma" = -1.0
11:35:52 Py:MadGraphUtils INFO "drll_sf" = 0.0
11:35:52 Py:MadGraphUtils INFO "muR_over_ref" = 1.0
11:35:52 Py:MadGraphUtils INFO "pdlabel" = lhapdf
11:35:52 Py:MadGraphUtils INFO "parton_shower" = HERWIGPP
11:35:52 Py:MadGraphUtils INFO "event_norm" = average
11:35:52 Py:MadGraphUtils INFO "mll_sf" = 30.0
11:35:52 Py:MadGraphUtils INFO "muF_ref_fixed" = 172.5
11:35:52 Py:MadGraphUtils INFO "bwcutoff" = 50
11:35:52 Py:MadGraphUtils INFO "ptgmin" = 20.0
11:35:52 Py:MadGraphUtils INFO "rw_fscale" = 1.0, 2.0, 0.5
11:35:52 Py:MadGraphUtils INFO "mxx_min_pdg" = {}
11:35:52 Py:MadGraphUtils INFO "npoints_FO" = 10000
11:35:52 Py:MadGraphUtils INFO "reweight_PDF" = False
11:35:52 Py:MadGraphUtils INFO "muF_over_ref" = 1.0
11:35:52 Py:MadGraphUtils INFO "isoEM" = True
11:35:52 Py:MadGraphUtils INFO "mll" = 0.0
11:35:52 Py:MadGraphUtils INFO "fixed_fac_scale" = True
11:35:52 Py:MadGraphUtils INFO "lpp2" = 1
11:35:52 Py:MadGraphUtils INFO "lhaid" = 260000
11:35:52 Py:MadGraphUtils INFO "niters_FO_grid" = 4
11:35:52 Py:MadGraphUtils INFO "etaj" = -1.0
11:35:52 Py:MadGraphUtils INFO "jetradius" = 0.4
11:35:52 Py:MadGraphUtils INFO "jetalgo" = -1
11:35:52 Py:MadGraphUtils INFO "req_acc" = 0.001000
11:35:52 Py:MadGraphUtils INFO "shower_scale_factor" = 1.0
11:35:52 Py:MadGraphUtils INFO "ebeam2" = 6500
11:35:52 Py:MadGraphUtils INFO "run_tag" = tag_1
11:35:52 Py:MadGraphUtils INFO "ebeam1" = 6500
11:35:52 Py:MadGraphUtils INFO Finished checking run card - All OK!
11:35:53 [?1034hNo module named madgraph
11:35:53 INFO: ************************************************************
11:35:53 * *
11:35:53 * W E L C O M E to M A D G R A P H 5 *
11:35:53 * a M C @ N L O *
11:35:53 * *
11:35:53 * * * *
11:35:53 * * * * * *
11:35:53 * * * * * 5 * * * * *
11:35:53 * * * * * *
11:35:53 * * * *
11:35:53 * *
11:35:53 * VERSION 5.2.6.2 20xx-xx-xx *
11:35:53 * *
11:35:53 * The MadGraph5_aMC@NLO Development Team - Find us at *
11:35:53 * http://amcatnlo.cern.ch *
11:35:53 * *
11:35:53 * Type 'help' for in-line help. *
11:35:53 * *
11:35:53 ************************************************************
11:35:54 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:35:54 INFO: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
11:35:54 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:35:54 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
11:35:54 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
11:35:54 launch --name=schan
11:35:54 stty: standard input: Invalid argument
11:35:54 stty: standard input: Invalid argument
11:35:54 The following switches determine which programs are run:
11:35:54 /=============================================================================\
11:35:54 1. Type of perturbative computation  order = NLO
11:35:54 2. No MC@[N]LO matching / event generation fixed_order = OFF
11:35:54 3. Shower the generated events  shower = OFF
11:35:54 4. Decay onshell particles  madspin = ON
11:35:54 5. Add weights to events for new hypp.  reweight = OFF
11:35:54 6. Run MadAnalysis5 on the events generated madanalysis = Not Avail.
11:35:54 \=============================================================================/
11:35:54 Either type the switch number (1 to 6) to change its setting,
11:35:54 Set any switch explicitly (e.g. type 'fixed_order=ON' at the prompt)
11:35:54 Type 'help' for the list of all valid option
11:35:54 Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
11:35:54 >INFO: will run in mode: noshower
11:35:54 WARNING: You have chosen not to run a parton shower.
11:35:54 NLO events without showering are NOT physical.
11:35:54 Please, shower the LesHouches events before using them for physics analyses.
11:35:54 You have to choose NOW which parton-shower you WILL use and specify it in the run_card. 
11:35:54 Do you want to edit a card (press enter to bypass editing)?
11:35:54 /------------------------------------------------------------\
11:35:54 | 1. param : param_card.dat |
11:35:54 | 2. run : run_card.dat |
11:35:54 | 3. madspin : madspin_card.dat |
11:35:54 \------------------------------------------------------------/
11:35:54 you can also
11:35:54 - enter the path to a valid card or banner.
11:35:54 - use the 'set' command to modify a parameter directly.
11:35:54 The set option works only for param_card and run_card.
11:35:54 Type 'help set' for more information on this command.
11:35:54 - call an external program (ASperGE/MadWidth/...).
11:35:54 Type 'help' for the list of available command
11:35:54 [0, done, 1, param, 2, run, 3, madspin, enter path][90s to answer]
11:35:54 >INFO: Update the dependent parameter of the param_card.dat
11:35:55 INFO: Starting run
11:35:55 INFO: Compiling the code
11:35:55 WARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) 
11:35:55 WARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) 
11:35:56 INFO: For gauge cancellation, the width of 't' has been set to zero. 
11:35:57 INFO: Using LHAPDF v6.1.5 interface for PDFs
11:36:02 INFO: Compiling source...
11:36:19 INFO: ...done, continuing with P* directories
11:36:19 INFO: Compiling StdHEP (can take a couple of minutes) ...
11:36:46 INFO: ...done.
11:36:46 INFO: Compiling directories...
11:36:46 INFO: Compiling on 8 cores
11:36:46 INFO: Compiling P0_udx_wp_tbx...
11:36:46 INFO: Compiling P0_dxu_wp_tbx...
11:36:46 INFO: Compiling P0_dux_wm_txb...
11:36:46 INFO: Compiling P0_uxd_wm_txb...
11:37:17 INFO: P0_udx_wp_tbx done.
11:37:18 INFO: P0_dxu_wp_tbx done.
11:37:19 INFO: P0_dux_wm_txb done.
11:37:19 INFO: P0_uxd_wm_txb done.
11:37:19 INFO: Checking test output:
11:37:19 INFO: P0_udx_wp_tbx
11:37:19 INFO: Result for test_ME:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for test_MC:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for check_poles:
11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
11:37:19 INFO: P0_dxu_wp_tbx
11:37:19 INFO: Result for test_ME:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for test_MC:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for check_poles:
11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
11:37:19 INFO: P0_dux_wm_txb
11:37:19 INFO: Result for test_ME:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for test_MC:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for check_poles:
11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
11:37:19 INFO: P0_uxd_wm_txb
11:37:19 INFO: Result for test_ME:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for test_MC:
11:37:19 INFO: Passed.
11:37:19 INFO: Result for check_poles:
11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
11:37:19 INFO: Starting run
11:37:19 INFO: Using 8 cores
11:37:19 INFO: Cleaning previous results
11:37:19 INFO: Generating events without running the shower.
11:37:19 INFO: Setting up grids
11:37:20 INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 11h37 ]
11:37:30 INFO: Idle: 0, Running: 3, Completed: 1 [ 9.7s ]
11:37:30 INFO: Idle: 0, Running: 0, Completed: 4 [ 10.3s ]
11:37:30 INFO: Determining the number of unweighted events per channel
11:37:30
11:37:30 Intermediate results:
11:37:30 Random seed: 1
11:37:30 Total cross section: 1.019e+01 +- 6.4e-02 pb
11:37:30 Total abs(cross section): 1.515e+01 +- 9.1e-02 pb
11:37:30
11:37:30
11:37:30 INFO: Computing upper envelope
11:37:31 INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 11h37 ]
11:43:36 INFO: Idle: 0, Running: 3, Completed: 1 [ 6m 4s ]
11:43:53 INFO: Idle: 0, Running: 0, Completed: 4 [ 6m 21s ]
11:43:53 INFO: Updating the number of unweighted events per channel
11:43:53
11:43:53 Intermediate results:
11:43:53 Random seed: 1
11:43:53 Total cross section: 1.027e+01 +- 9.9e-03 pb
11:43:53 Total abs(cross section): 1.519e+01 +- 1.4e-02 pb
11:43:53
11:43:53
11:43:53 INFO: Generating events
11:43:54 INFO: Idle: 0, Running: 0, Completed: 1 [ current time: 11h43 ]
11:43:54 INFO: Collecting events
11:43:55 INFO:
11:43:55 --------------------------------------------------------------
11:43:55 Summary:
11:43:55 Process p p > w+ > t b~ [QCD] ; p p > w- > t~ b [QCD]
11:43:55 Run at p-p collider (6500.0 + 6500.0 GeV)
11:43:55 Number of events generated: 1
11:43:55 Total cross section: 3.163e+00 +- 6.4e-03 pb
11:43:55 --------------------------------------------------------------
11:43:55
11:43:55 INFO: The /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events.lhe.gz file has been generated.
11:43:55
11:43:55 INFO: Events generated
11:43:55 reweight -from_cards
11:43:55 decay_events -from_cards
11:43:57 INFO: Running MadSpin
11:43:57 INFO: This functionality allows for the decay of resonances
11:43:57 INFO: in a .lhe file, keeping track of the spin correlation effets.
11:43:57 INFO: BE AWARE OF THE CURRENT LIMITATIONS:
11:43:57 INFO: (1) Only a succession of 2 body decay are currently allowed
11:43:57 ************************************************************
11:43:57 * *
11:43:57 * W E L C O M E to M A D S P I N *
11:43:57 * *
11:43:57 ************************************************************
11:43:57 INFO: Extracting the banner ...
11:43:57 INFO: process: p p > w+ > t b~
11:43:57 INFO: options:
11:43:58 INFO: detected model: loop_sm-no_b_mass. Loading...
11:43:58 import Events/schan/events.lhe.gz
11:43:58 INFO: process: p p > w+ > t b~
11:43:58 INFO: options:
11:43:58 Set group_subprocesses to Auto
11:43:58 Note that you need to regenerate all processes
11:43:58 set loop optimized output to True
11:43:58 set loop color flows to False
11:43:58 Complex mass already desactivated.
11:43:58 INFO: detected model: loop_sm-no_b_mass. Loading...
11:43:58 INFO: load particles
11:43:58 INFO: load vertices
11:43:58 INFO: Restrict model loop_sm with file ../../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
11:43:58 INFO: Run "set stdout_level DEBUG" before import for more information.
11:43:58 INFO: Change particles name to pass to MG5 convention
11:43:58 Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
11:43:58 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
11:43:58 Set group_subprocesses to Auto
11:43:58 Note that you need to regenerate all processes
11:43:58 set loop optimized output to True
11:43:58 set loop color flows to False
11:43:58 Complex mass already desactivated.
11:43:58 Defined multiparticle p = g u c d s u~ c~ d~ s~
11:43:58 Defined multiparticle j = g u c d s u~ c~ d~ s~
11:43:58 Defined multiparticle l+ = e+ mu+
11:43:58 Defined multiparticle l- = e- mu-
11:43:58 Defined multiparticle vl = ve vm vt
11:43:58 Defined multiparticle vl~ = ve~ vm~ vt~
11:43:58 Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
11:43:58 Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
11:43:58 Defined multiparticle l+ = e+ mu+ ta+
11:43:58 Defined multiparticle l- = e- mu- ta-
11:43:58 set ms_dir MadSpin
11:43:58 set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
11:43:58 set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
11:43:58 decay t > w+ b, w+ > l+ vl
11:43:58 decay t~ > w- b~, w- > l- vl~
11:43:58 launch
11:43:58 INFO: Will use seed 766814452
11:43:58 INFO:
11:43:58 INFO: decay channels for t~ : ( width = 1.32 GeV )
11:43:58 INFO: BR d1 d2
11:43:58 INFO: 1.000000e+00 w- b~
11:43:58 INFO:
11:43:58 INFO:
11:43:58 INFO: decay channels for w- : ( width = 2.085 GeV )
11:43:58 INFO: BR d1 d2
11:43:58 INFO: 3.337000e-01 u~ d
11:43:58 INFO: 3.337000e-01 c~ s
11:43:58 INFO: 1.082000e-01 ve~ e-
11:43:58 INFO: 1.082000e-01 vm~ mu-
11:43:58 INFO: 1.082000e-01 vt~ ta-
11:43:58 INFO:
11:43:58 INFO:
11:43:58 INFO: decay channels for w+ : ( width = 2.085 GeV )
11:43:58 INFO: BR d1 d2
11:43:58 INFO: 3.337000e-01 u d~
11:43:58 INFO: 3.337000e-01 c s~
11:43:58 INFO: 1.082000e-01 ve e+
11:43:58 INFO: 1.082000e-01 vm mu+
11:43:58 INFO: 1.082000e-01 vt ta+
11:43:58 INFO:
11:43:58 INFO:
11:43:58 INFO: decay channels for t : ( width = 1.32 GeV )
11:43:58 INFO: BR d1 d2
11:43:58 INFO: 1.000000e+00 w+ b
11:43:58 INFO:
11:43:59 INFO: generating the production square matrix element
11:43:59 INFO: generate p p > w+ > t b~ --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w+ > t b~ pert_QCD --no_warning=duplicate;add process p p > w- > t~ b --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w- > t~ b pert_QCD --no_warning=duplicate;
11:44:03 INFO: Done 4.271
11:44:03 INFO: generating the full matrix element squared (with decay)
11:44:03 INFO: generate p p > w+ > t b~, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w+ > t b~ pert_QCD, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;add process p p > w- > t~ b, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w- > t~ b pert_QCD, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;
11:44:09 INFO: Done 6.096
11:44:09 INFO: generate matrix element for decay only (1 - > N).
11:44:09 INFO: output standalone_msF /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/MadSpin/decay_me
11:44:13 INFO: Done 3.253
11:44:13 INFO: Compiling code
11:44:41 INFO: detect independant decays
11:44:41 INFO: Done in 0.252036094666s
11:44:41 INFO:
11:44:41 INFO: Estimating the maximum weight
11:44:41 INFO: *****************************
11:44:41 INFO: Probing the first 250 events
11:44:41 INFO: with 1000 phase space points
11:44:41 INFO: For 2 decaying particle type in the final states
11:44:41 INFO:
11:44:41 INFO: Event 1/500 : 0.17s
11:44:41 INFO: Not enough events for at least one production mode.
11:44:41 INFO: This is ok as long as you don't reuse the max weight for other generations.
11:44:42 INFO:
11:44:42 INFO: Decaying the events...
11:44:42 INFO: Total number of events written: 1/1
11:44:42 INFO: Average number of trial points per production event: 1.0
11:44:42 INFO: Branching ratio to allowed decays: 0.3246
11:44:42 INFO: Number of events with weights larger than max_weight: 0
11:44:42 INFO: Number of subprocesses 1
11:44:42 INFO: Number of failures when restoring the Monte Carlo masses: 0
11:44:42 INFO: Decayed events have been written in /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events_decayed.lhe.gz
11:44:42 INFO: The decayed event file has been moved to the following location:
11:44:42 INFO: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan_decayed_1/events.lhe.gz
11:44:42 INFO: MadSpin Done
11:44:42 WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical.
11:44:42 Please, shower the Les Houches events before using them for physics analyses. 
11:44:42 INFO:
11:44:42 quit
11:44:42 INFO:
11:44:42 Py:MadGraphUtils INFO "rw_rscale" = 1.0, 2.0, 0.5
11:44:42 Py:MadGraphUtils INFO "ptl" = 0.0
11:44:42 Py:MadGraphUtils INFO "xn" = 1.0
11:44:42 Py:MadGraphUtils INFO "ptj" = 0.1
11:44:42 Py:MadGraphUtils INFO "niters_FO" = 6
11:44:42 Py:MadGraphUtils INFO "nevt_job" = -1
11:44:42 Py:MadGraphUtils INFO "drll" = 0.0
11:44:42 Py:MadGraphUtils INFO "fixed_ren_scale" = True
11:44:42 Py:MadGraphUtils INFO "npoints_FO_grid" = 5000
11:44:42 Py:MadGraphUtils INFO "iseed" = 1
11:44:42 Py:MadGraphUtils INFO "muR_ref_fixed" = 172.5
11:44:42 Py:MadGraphUtils INFO "dynamical_scale_choice" = -1
11:44:42 Py:MadGraphUtils INFO "R0gamma" = 0.4
11:44:42 Py:MadGraphUtils INFO "reweight_scale" = False
11:44:42 Py:MadGraphUtils INFO "store_rwgt_info" = False
11:44:42 Py:MadGraphUtils INFO "lpp1" = 1
11:44:42 Py:MadGraphUtils INFO "nevents" = 1
11:44:42 Py:MadGraphUtils INFO "pt_max_pdg" = {}
11:44:42 Py:MadGraphUtils INFO "ickkw" = 0
11:44:42 Py:MadGraphUtils INFO "iappl" = 0
11:44:42 Py:MadGraphUtils INFO "etal" = -1.0
11:44:42 Py:MadGraphUtils INFO "epsgamma" = 1.0
11:44:42 Py:MadGraphUtils INFO "req_acc_FO" = 0.01
11:44:42 Py:MadGraphUtils INFO "pt_min_pdg" = {}
11:44:42 Py:MadGraphUtils INFO "etagamma" = -1.0
11:44:42 Py:MadGraphUtils INFO "drll_sf" = 0.0
11:44:42 Py:MadGraphUtils INFO "muR_over_ref" = 1.0
11:44:42 Py:MadGraphUtils INFO "pdlabel" = lhapdf
11:44:42 Py:MadGraphUtils INFO "parton_shower" = HERWIGPP
11:44:42 Py:MadGraphUtils INFO "event_norm" = average
11:44:42 Py:MadGraphUtils INFO "mll_sf" = 30.0
11:44:42 Py:MadGraphUtils INFO "muF_ref_fixed" = 172.5
11:44:42 Py:MadGraphUtils INFO "bwcutoff" = 50
11:44:42 Py:MadGraphUtils INFO "ptgmin" = 20.0
11:44:42 Py:MadGraphUtils INFO "rw_fscale" = 1.0, 2.0, 0.5
11:44:42 Py:MadGraphUtils INFO "mxx_min_pdg" = {}
11:44:42 Py:MadGraphUtils INFO "npoints_FO" = 10000
11:44:42 Py:MadGraphUtils INFO "reweight_PDF" = False
11:44:42 Py:MadGraphUtils INFO "muF_over_ref" = 1.0
11:44:42 Py:MadGraphUtils INFO "isoEM" = True
11:44:42 Py:MadGraphUtils INFO "mll" = 0.0
11:44:42 Py:MadGraphUtils INFO "fixed_fac_scale" = True
11:44:42 Py:MadGraphUtils INFO "lpp2" = 1
11:44:42 Py:MadGraphUtils INFO "lhaid" = 260000
11:44:42 Py:MadGraphUtils INFO "niters_FO_grid" = 4
11:44:42 Py:MadGraphUtils INFO "etaj" = -1.0
11:44:42 Py:MadGraphUtils INFO "jetradius" = 0.4
11:44:42 Py:MadGraphUtils INFO "jetalgo" = -1
11:44:42 Py:MadGraphUtils INFO "req_acc" = 0.001000
11:44:42 Py:MadGraphUtils INFO "shower_scale_factor" = 1.0
11:44:42 Py:MadGraphUtils INFO "ebeam2" = 6500
11:44:42 Py:MadGraphUtils INFO "run_tag" = tag_1
11:44:42 Py:MadGraphUtils INFO "ebeam1" = 6500
11:44:42 Py:MadGraphUtils INFO Finished checking run card - All OK!
11:44:42 Py:MadGraphUtils INFO Modifying run card located at Cards/run_card.dat.
11:44:42 Py:MadGraphUtils INFO Keeping backup of original run card at Cards/run_card_backup.dat.
11:44:42 {'nevents': 250}
11:44:42 Py:MadGraphUtils INFO Setting nevents = 250.
11:44:42 Py:MadGraphUtils INFO Finished modification of run card.
11:44:45 [?1034hNo module named madgraph
11:44:45 INFO: ************************************************************
11:44:45 * *
11:44:45 * W E L C O M E to M A D G R A P H 5 *
11:44:45 * a M C @ N L O *
11:44:45 * *
11:44:45 * * * *
11:44:45 * * * * * *
11:44:45 * * * * * 5 * * * * *
11:44:45 * * * * * *
11:44:45 * * * *
11:44:45 * *
11:44:45 * VERSION 5.2.6.2 20xx-xx-xx *
11:44:45 * *
11:44:45 * The MadGraph5_aMC@NLO Development Team - Find us at *
11:44:45 * http://amcatnlo.cern.ch *
11:44:45 * *
11:44:45 * Type 'help' for in-line help. *
11:44:45 * *
11:44:45 ************************************************************
11:44:45 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:44:45 INFO: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
11:44:45 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
11:44:45 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
11:44:45 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
11:44:45 launch --parton --nocompile --only_generation -f --name=schan
11:44:45 WARNING: Removing old run information in
11:44:45 /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan 
11:44:45 stty: standard input: Invalid argument
11:44:45 stty: standard input: Invalid argument
11:44:45 INFO: will run in mode: noshower
11:44:45 WARNING: You have chosen not to run a parton shower.
11:44:45 NLO events without showering are NOT physical.
11:44:45 Please, shower the LesHouches events before using them for physics analyses.
11:44:45 You have to choose NOW which parton-shower you WILL use and specify it in the run_card. 
11:44:45 INFO: Starting run
11:44:45 INFO: Compiling the code
11:44:45 WARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) 
11:44:45 WARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) 
11:44:47 INFO: For gauge cancellation, the width of 't' has been set to zero. 
11:44:47 INFO: Starting run
11:44:47 INFO: Using 8 cores
11:44:47 INFO: Generating events without running the shower.
11:44:47 INFO: Updating the number of unweighted events per channel
11:44:47
11:44:47 Intermediate results:
11:44:47 Random seed: 1
11:44:47 Total cross section: 1.027e+01 +- 9.9e-03 pb
11:44:47 Total abs(cross section): 1.519e+01 +- 1.4e-02 pb
11:44:47
11:44:47
11:44:47 INFO: Generating events
11:44:48 INFO: Idle: 0, Running: 2, Completed: 2 [ current time: 11h44 ]
11:44:48 INFO: Idle: 0, Running: 1, Completed: 3 [ 0.17s ]
11:44:48 INFO: Idle: 0, Running: 0, Completed: 4 [ 0.18s ]
11:44:48 INFO: Collecting events
11:44:48 INFO:
11:44:48 --------------------------------------------------------------
11:44:48 Summary:
11:44:48 Process p p > w+ > t b~ [QCD] ; p p > w- > t~ b [QCD]
11:44:48 Run at p-p collider (6500.0 + 6500.0 GeV)
11:44:48 Number of events generated: 250
11:44:48 Total cross section: 1.027e+01 +- 9.9e-03 pb
11:44:48 --------------------------------------------------------------
11:44:48
11:44:48 INFO: The /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events.lhe.gz file has been generated.
11:44:48
11:44:48 INFO: Events generated
11:44:48 reweight -from_cards
11:44:48 decay_events -from_cards
11:44:50 INFO: Running MadSpin
11:44:50 INFO: This functionality allows for the decay of resonances
11:44:50 INFO: in a .lhe file, keeping track of the spin correlation effets.
11:44:50 INFO: BE AWARE OF THE CURRENT LIMITATIONS:
11:44:50 INFO: (1) Only a succession of 2 body decay are currently allowed
11:44:50 ************************************************************
11:44:50 * *
11:44:50 * W E L C O M E to M A D S P I N *
11:44:50 * *
11:44:50 ************************************************************
11:44:50 INFO: Extracting the banner ...
11:44:50 INFO: process: p p > w+ > t b~
11:44:50 INFO: options:
11:44:50 INFO: detected model: loop_sm-no_b_mass. Loading...
11:44:51 import Events/schan/events.lhe.gz
11:44:51 INFO: process: p p > w+ > t b~
11:44:51 INFO: options:
11:44:51 Set group_subprocesses to Auto
11:44:51 Note that you need to regenerate all processes
11:44:51 set loop optimized output to True
11:44:51 set loop color flows to False
11:44:51 Complex mass already desactivated.
11:44:51 INFO: detected model: loop_sm-no_b_mass. Loading...
11:44:51 INFO: load particles
11:44:51 INFO: load vertices
11:44:51 INFO: Restrict model loop_sm with file ../../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
11:44:51 INFO: Run "set stdout_level DEBUG" before import for more information.
11:44:51 INFO: Change particles name to pass to MG5 convention
11:44:51 Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
11:44:51 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
11:44:51 Set group_subprocesses to Auto
11:44:51 Note that you need to regenerate all processes
11:44:51 set loop optimized output to True
11:44:51 set loop color flows to False
11:44:51 Complex mass already desactivated.
11:44:51 Defined multiparticle p = g u c d s u~ c~ d~ s~
11:44:51 Defined multiparticle j = g u c d s u~ c~ d~ s~
11:44:51 Defined multiparticle l+ = e+ mu+
11:44:51 Defined multiparticle l- = e- mu-
11:44:51 Defined multiparticle vl = ve vm vt
11:44:51 Defined multiparticle vl~ = ve~ vm~ vt~
11:44:51 Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
11:44:51 Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
11:44:51 Defined multiparticle l+ = e+ mu+ ta+
11:44:51 Defined multiparticle l- = e- mu- ta-
11:44:51 set ms_dir MadSpin
11:44:51 set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
11:44:51 set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
11:44:51 decay t > w+ b, w+ > l+ vl
11:44:51 decay t~ > w- b~, w- > l- vl~
11:44:51 launch
11:44:51 INFO:
11:44:51 INFO: Decaying the events...
11:44:52 CRITICAL: over_weight: 1.23606071569 ('t_wpb_wp_veep',), occurence: 0.45871559633%, occurence_channel: 1.96078431373%
11:44:52 production_tag:P1_udx_wp_tbx [((-1, 2), (-5, 6))], decay:P1_udx_wp_tbx_t_bwp_wp_epve [('t_wpb_wp_veep',)], BW_cut: 23.0977
11:44:52
11:44:52 [decay.py at line 2297] 
11:44:52 INFO: Total number of events written: 250/250
11:44:52 INFO: Average number of trial points per production event: 7.576
11:44:52 INFO: Branching ratio to allowed decays: 0.3246
11:44:52 INFO: Number of events with weights larger than max_weight: 1
11:44:52 INFO: Number of subprocesses 18
11:44:52 INFO: Number of failures when restoring the Monte Carlo masses: 0
11:44:52 INFO: Decayed events have been written in /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events_decayed.lhe.gz
11:44:52 INFO: The decayed event file has been moved to the following location:
11:44:52 INFO: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan_decayed_2/events.lhe.gz
11:44:52 INFO: MadSpin Done
11:44:52 WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical.
11:44:52 Please, shower the Les Houches events before using them for physics analyses. 
11:44:52 INFO:
11:44:52 quit
11:44:52 INFO:
11:44:52 Py:MadGraphUtils INFO Creating gridpack (schan_gridpack.tar.gz)...
11:44:52 Py:MadGraphUtils INFO Package up proc_dir
11:45:06 Shortened traceback (most recent user call last):
11:45:06 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/jobOptions/EvgenJobTransforms/skeleton.GENtoEVGEN.py", line 225, in <module>
11:45:06 include(jo)
11:45:06 File "./MC15.412049.aMcAtNloHerwig7EG_H7UE_singletop_schan_lept.py", line 14, in <module>
11:45:06 include("./MadGraphControl_schan_NLO_Herwig7_A14.py")
11:45:06 File "./MadGraphControl_schan_NLO_Herwig7_A14.py", line 218, in <module>
11:45:06 cluster_queue=cluster_queue)
11:45:06 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/python/MadGraphControl/MadGraphUtils.py", line 584, in generate
11:45:06 raise RuntimeError('Gridpack sucessfully created, exiting the transform. IGNORE ERRORS if running gridpack generation!')
11:45:06 RuntimeError: Gridpack sucessfully created, exiting the transform. IGNORE ERRORS if running gridpack generation!
11:45:06 Py:Athena INFO leaving with code 8: "an unknown exception occurred"

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2018-08-22
Last reply:
2018-08-23

Dear Marc,

The reason of this message is that:
1) This is not bad enough to make the job to crash. (if such situation was occuring too often the job will simply crash)
2) This is worrying enough that the user need to be aware of it

In other word, you have to make the judgement call if this is good/bad enough for you.
- If this is a madspin creation run that will be reused for a large number of events, then for sure you have to discard this run and re-do it (probably with more events in the stage of the evaluation of the max-weight)
- Otherwise, this means that one event is going to have the wrong weight (1 instead of 1.24) leading to bias.
Up to you if you can leave with such error or not (low occurence of such issue are typically not a huge issue)

You can also note that your event is quite offshell and therefore the physical reason of such warning is that NWA is not valid anymore and that you are sensitive to offshell effects. So you might need to consider to use MadSpin on a more restricted range of mass (or just accept it as part of your theoretical error related to the NWA)

So as I said this is a judgement call that you have to make.

Cheers,

Olivier

> On 22 Aug 2018, at 10:52, Marc de Beurs <email address hidden> wrote:
>
> New question #672577 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/672577
>
> Dear Madspin experts,
>
> When I am running a single top sample, (t-channel or s-channel), and I want the tops to decay lepontically (t > w+ b, w+ > l+ vl) in madspin, I sometimes run into the critical message of:
>
> CRITICAL: over_weight: 1.24394983756 ('tx_bxwm_wm_tamvtx',), occurence: 0.0492610837438%, occurence_channel: 0.392156862745% production_tag:P4_us_txbuc [((1, 4), (-6, 2, 4, 5))], decay:P4_us_txbuc_tx_bxwm_wm_tamvtx [('tx_bxwm_wm_tamvtx',)], BW_cut: 13.097 [decay.py at line 2297]
>
> Is this something to worry about?
>
> I am running MG5 version 2.6.2 and below you can find the std output which also contains the cards (I am running it inside the atlas framework, so you can ignore all the messages that are unknown to you):
>
> Cheers Marc
>
> 11:34:57 Tue Aug 21 11:34:57 CEST 2018
> 11:34:57 Preloading tcmalloc_minimal.so
> 11:34:57 Preloading /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/IntelSoftware/linux/x86_64/xe2013/composer_xe_2013.3.163/compiler/lib/intel64/libintlc.so.5:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/IntelSoftware/linux/x86_64/xe2013/composer_xe_2013.3.163/compiler/lib/intel64/libimf.so
> 11:34:57 Py:Athena INFO including file "AthenaCommon/Preparation.py"
> 11:34:57 Py:Athena INFO using release [MCProd-19.2.5.33.4] [x86_64-slc6-gcc47-opt] [19.2.X.Y.Z-VAL-MCProd/rel_0] -- built on [2018-05-23 11:02]
> 11:34:57 Py:Athena INFO including file "AthenaCommon/Bootstrap.py"
> 11:34:57 Py:Athena INFO including file "AthenaCommon/Atlas.UnixStandardJob.py"
> 11:34:59 Py:Athena INFO including file "AthenaCommon/Execution.py"
> 11:34:59 Py:Athena INFO including file "runargs.generate.py"
> 11:34:59 Py:Athena INFO including file "EvgenJobTransforms/skeleton.GENtoEVGEN.py"
> 11:34:59 Py:Athena INFO including file "AthenaCommon/Atlas_Gen.UnixStandardJob.py"
> 11:34:59 Py:Athena INFO including file "PartPropSvc/PartPropSvc.py"
> 11:35:00 Py:Athena INFO including file "MC15.412049.aMcAtNloHerwig7EG_H7UE_singletop_schan_lept.py"
> 11:35:00 Py:Athena INFO including file "./MadGraphControl_schan_NLO_Herwig7_A14.py"
> 11:35:00 Py:MadGraphUtils INFO Assuming that proc_card_mg5.dat is a well-formatted process card
> 11:35:00 Py:MadGraphUtils INFO Started process generation at Tue Aug 21 11:35:00 2018
> 11:35:04 [?1034h************************************************************
> 11:35:04 * *
> 11:35:04 * W E L C O M E to *
> 11:35:04 * M A D G R A P H 5 _ a M C @ N L O *
> 11:35:04 * *
> 11:35:04 * *
> 11:35:04 * * * *
> 11:35:04 * * * * * *
> 11:35:04 * * * * * 5 * * * * *
> 11:35:04 * * * * * *
> 11:35:04 * * * *
> 11:35:04 * *
> 11:35:04 * VERSION 2.6.2 2018-04-29 *
> 11:35:04 * *
> 11:35:04 * The MadGraph5_aMC@NLO Development Team - Find us at *
> 11:35:04 * https://server06.fynu.ucl.ac.be/projects/madgraph *
> 11:35:04 * and *
> 11:35:04 * http://amcatnlo.web.cern.ch/amcatnlo/ *
> 11:35:04 * *
> 11:35:04 * Type 'help' for in-line help. *
> 11:35:04 * Type 'tutorial' to learn how MG5 works *
> 11:35:04 * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
> 11:35:04 * Type 'tutorial MadLoop' to learn how MadLoop works *
> 11:35:04 * *
> 11:35:04 ************************************************************
> 11:35:04 load MG5 configuration from ../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
> 11:35:04 set collier to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/collier/1.1/x86_64-slc6-gcc47-opt
> 11:35:04 fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
> 11:35:04 Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
> 11:35:04 MG5_aMC> set fastjet /PATH/TO/fastjet-config
> 11:35:04
> 11:35:04 set lhapdf to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.4/x86_64-slc6-gcc47-opt/bin/lhapdf-config
> 11:35:04 set ninja to /afs/.cern.ch/sw/lcg/external/MCGenerators_lcgcmt67c/gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib
> 11:35:04 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> 11:35:04 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> 11:35:04 import /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/proc_card_mg5.dat
> 11:35:04 The import format was not given, so we guess it as command
> 11:35:04 import model loop_sm-no_b_mass
> 11:35:05 INFO: reload from .py file
> 11:35:05 INFO: load particles
> 11:35:05 INFO: load vertices
> 11:35:05 INFO: Restrict model loop_sm with file ../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
> 11:35:05 INFO: Run "set stdout_level DEBUG" before import for more information.
> 11:35:05 INFO: Change particles name to pass to MG5 convention
> 11:35:05 Defined multiparticle p = g u c d s u~ c~ d~ s~
> 11:35:05 Defined multiparticle j = g u c d s u~ c~ d~ s~
> 11:35:05 Defined multiparticle l+ = e+ mu+
> 11:35:05 Defined multiparticle l- = e- mu-
> 11:35:05 Defined multiparticle vl = ve vm vt
> 11:35:05 Defined multiparticle vl~ = ve~ vm~ vt~
> 11:35:05 Pass the definition of 'j' and 'p' to 5 flavour scheme.
> 11:35:05 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
> 11:35:05 define l+ = l+ ta+
> 11:35:05 Defined multiparticle l+ = e+ mu+ ta+
> 11:35:05 define l- = l- ta-
> 11:35:05 Defined multiparticle l- = e- mu- ta-
> 11:35:05 generate p p > w+ > t b~ [QCD]
> 11:35:05 WARNING: Process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > w+> t b~ WEIGHTED=4 [ all = QCD ] can have real emission processes which are not finite.
> 11:35:05 To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
> 11:35:05 g, b, c, s, u, d, d~, u~, s~, c~, b~ 
> 11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: u d~ > w+ > t b~ [ all = QCD ] (1 / 4)
> 11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: c s~ > w+ > t b~ [ all = QCD ] (2 / 4)
> 11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: d~ u > w+ > t b~ [ all = QCD ] (3 / 4)
> 11:35:05 INFO: Generating FKS-subtracted matrix elements for born process: s~ c > w+ > t b~ [ all = QCD ] (4 / 4)
> 11:35:05 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:05 INFO: Generating virtual matrix element with MadLoop for process: u d~ > w+ > t b~ QED<=2 [ all = QCD ] (1 / 4)
> 11:35:05 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:05 INFO: Generating virtual matrix element with MadLoop for process: c s~ > w+ > t b~ QED<=2 [ all = QCD ] (2 / 4)
> 11:35:06 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:06 INFO: Generating virtual matrix element with MadLoop for process: d~ u > w+ > t b~ QED<=2 [ all = QCD ] (3 / 4)
> 11:35:06 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:06 INFO: Generating virtual matrix element with MadLoop for process: s~ c > w+ > t b~ QED<=2 [ all = QCD ] (4 / 4)
> 11:35:06 INFO: Generated 4 subprocesses with 32 real emission diagrams, 4 born diagrams and 8 virtual diagrams
> 11:35:06 add process p p > w- > t~ b [QCD]
> 11:35:06 WARNING: Process: g/u/c/d/s/u~/c~/d~/s~/b/b~ g/u/c/d/s/u~/c~/d~/s~/b/b~ > w-> t~ b WEIGHTED=4 [ all = QCD ] can have real emission processes which are not finite.
> 11:35:06 To avoid this, please use multiparticles when generating the process and be sure to include all the following particles in the multiparticle definition:
> 11:35:06 g, b, c, s, u, d, d~, u~, s~, c~, b~ 
> 11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: d u~ > w- > t~ b [ all = QCD ] (1 / 4)
> 11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: s c~ > w- > t~ b [ all = QCD ] (2 / 4)
> 11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: u~ d > w- > t~ b [ all = QCD ] (3 / 4)
> 11:35:06 INFO: Generating FKS-subtracted matrix elements for born process: c~ s > w- > t~ b [ all = QCD ] (4 / 4)
> 11:35:06 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:06 INFO: Generating virtual matrix element with MadLoop for process: d u~ > w- > t~ b QED<=2 [ all = QCD ] (1 / 4)
> 11:35:06 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:06 INFO: Generating virtual matrix element with MadLoop for process: s c~ > w- > t~ b QED<=2 [ all = QCD ] (2 / 4)
> 11:35:06 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:06 INFO: Generating virtual matrix element with MadLoop for process: u~ d > w- > t~ b QED<=2 [ all = QCD ] (3 / 4)
> 11:35:07 INFO: Generating virtual matrix elements using MadLoop:
> 11:35:07 INFO: Generating virtual matrix element with MadLoop for process: c~ s > w- > t~ b QED<=2 [ all = QCD ] (4 / 4)
> 11:35:07 INFO: Generated 4 subprocesses with 32 real emission diagrams, 4 born diagrams and 8 virtual diagrams
> 11:35:07 output -f
> 11:35:07 INFO: Writing out the aMC@NLO code, using optimized Loops
> 11:35:07 INFO: initialize a new directory: PROCNLO_loop_sm-no_b_mass_0
> 11:35:11 INFO: remove old information in PROCNLO_loop_sm-no_b_mass_0
> 11:35:18 INFO: When using the 'output_dependencies=internal' MG5_aMC option, the (optional) reduction library ninja cannot be employed because it is not distributed with the MG5_aMC code so that it cannot be copied locally.
> 11:35:18 INFO: When using the 'output_dependencies=internal' MG5_aMC option, the (optional) reduction library collier cannot be employed because it is not distributed with the MG5_aMC code so that it cannot be copied locally.
> 11:35:18 INFO: Generating real emission matrix-elements...
> 11:35:19 INFO: Generating Helas calls for FKS process: u d~ > w+ > t b~ [ all = QCD ] (1 / 8)
> 11:35:19 INFO: Processing color information for process: u d~ > w+ > t b~ [ all = QCD ]
> 11:35:19 INFO: Generating Helas calls for FKS process: c s~ > w+ > t b~ [ all = QCD ] (2 / 8)
> 11:35:19 INFO: Generating Helas calls for FKS process: d~ u > w+ > t b~ [ all = QCD ] (3 / 8)
> 11:35:20 INFO: Processing color information for process: d~ u > w+ > t b~ [ all = QCD ]
> 11:35:20 INFO: Generating Helas calls for FKS process: s~ c > w+ > t b~ [ all = QCD ] (4 / 8)
> 11:35:20 INFO: Generating Helas calls for FKS process: d u~ > w- > t~ b [ all = QCD ] (5 / 8)
> 11:35:20 INFO: Processing color information for process: d u~ > w- > t~ b [ all = QCD ]
> 11:35:20 INFO: Generating Helas calls for FKS process: s c~ > w- > t~ b [ all = QCD ] (6 / 8)
> 11:35:21 INFO: Generating Helas calls for FKS process: u~ d > w- > t~ b [ all = QCD ] (7 / 8)
> 11:35:21 INFO: Processing color information for process: u~ d > w- > t~ b [ all = QCD ]
> 11:35:21 INFO: Generating Helas calls for FKS process: c~ s > w- > t~ b [ all = QCD ] (8 / 8)
> 11:35:22 INFO: ... Done
> 11:35:22 INFO: Writing files in P0_udx_wp_tbx (1 / 4)
> 11:35:22 INFO: Creating files in directory V0_udx_wp_tbx
> 11:35:22 ALOHA: aloha creates FFV1 set of routines with options: L2,P0
> 11:35:22 INFO: Computing diagram color coefficients
> 11:35:22 INFO: Drawing loop Feynman diagrams for Process: u d~ > w+ > t b~ QED<=2 [ all = QCD ]
> 11:35:22 INFO: Generating born Feynman diagrams for Process: u d~ > w+ > t b~ QED<=2 [ all = QCD ]
> 11:35:23 INFO: Writing files in P0_dxu_wp_tbx (2 / 4)
> 11:35:23 INFO: Creating files in directory V0_dxu_wp_tbx
> 11:35:23 INFO: Computing diagram color coefficients
> 11:35:23 INFO: Drawing loop Feynman diagrams for Process: d~ u > w+ > t b~ QED<=2 [ all = QCD ]
> 11:35:23 INFO: Generating born Feynman diagrams for Process: d~ u > w+ > t b~ QED<=2 [ all = QCD ]
> 11:35:24 INFO: Writing files in P0_dux_wm_txb (3 / 4)
> 11:35:24 INFO: Creating files in directory V0_dux_wm_txb
> 11:35:24 INFO: Computing diagram color coefficients
> 11:35:24 INFO: Drawing loop Feynman diagrams for Process: d u~ > w- > t~ b QED<=2 [ all = QCD ]
> 11:35:24 INFO: Generating born Feynman diagrams for Process: d u~ > w- > t~ b QED<=2 [ all = QCD ]
> 11:35:25 INFO: Writing files in P0_uxd_wm_txb (4 / 4)
> 11:35:25 INFO: Creating files in directory V0_uxd_wm_txb
> 11:35:25 INFO: Computing diagram color coefficients
> 11:35:25 INFO: Drawing loop Feynman diagrams for Process: u~ d > w- > t~ b QED<=2 [ all = QCD ]
> 11:35:25 INFO: Generating born Feynman diagrams for Process: u~ d > w- > t~ b QED<=2 [ all = QCD ]
> 11:35:26 History written to /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/proc_card_mg5.dat
> 11:35:26 ALOHA: aloha creates FFV2 routines
> 11:35:26 ALOHA: aloha creates FFV1 routines
> 11:35:27 save configuration file to /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:35:27 INFO: Use Fortran compiler gfortran
> 11:35:27 INFO: Use c++ compiler g++
> 11:35:27 INFO: Generate jpeg diagrams
> 11:35:35 INFO: Generate web pages
> 11:35:44 Type "launch" to generate events from this process, or see
> 11:35:44 /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/README
> 11:35:44 Run "open index.html" to see more information about this process.
> 11:35:44 quit
> 11:35:44 Py:MadGraphUtils INFO Finished process generation at Tue Aug 21 11:35:44 2018
> 11:35:44 Py:MadGraphUtils WARNING Path for option fastjet not found in original config
> 11:35:44 Py:MadGraphUtils INFO Modifying config paths to avoid use of afs:
> 11:35:44 Py:MadGraphUtils INFO {'ninja': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib', 'lhapdf': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//lhapdf/6.1.5/x86_64-slc6-gcc47-opt/bin/lhapdf-config', 'syscalc_path': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//syscalc/1.1.4/x86_64-slc6-gcc47-opt', 'collier': '/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//collier/1.1/x86_64-slc6-gcc47-opt'}
> 11:35:44 Py:MadGraphUtils INFO Copying default run_card.dat from PROCNLO_loop_sm-no_b_mass_0/Cards/run_card.dat
> 11:35:44 Py:MadGraphUtils INFO Found old card at run_card.tmp.dat
> 11:35:44 Py:MadGraphUtils INFO Found old param card at PROCNLO_loop_sm-no_b_mass_0/Cards/param_card.dat
> 11:35:44 Py:MadGraphUtils INFO proc_card:
> 11:35:44
> 11:35:44 import model loop_sm-no_b_mass
> 11:35:44 define l+ = l+ ta+
> 11:35:44 define l- = l- ta-
> 11:35:44 generate p p > w+ > t b~ [QCD]
> 11:35:44 add process p p > w- > t~ b [QCD]
> 11:35:44 output -f
> 11:35:44 Py:MadGraphUtils INFO run_card:
> 11:35:44 #***********************************************************************
> 11:35:44 # MadGraph5_aMC@NLO *
> 11:35:44 # *
> 11:35:44 # run_card.dat aMC@NLO *
> 11:35:44 # *
> 11:35:44 # This file is used to set the parameters of the run. *
> 11:35:44 # *
> 11:35:44 # Some notation/conventions: *
> 11:35:44 # *
> 11:35:44 # Lines starting with a hash (#) are info or comments *
> 11:35:44 # *
> 11:35:44 # mind the format: value = variable ! comment *
> 11:35:44 # *
> 11:35:44 # Some of the values of variables can be list. These can either be *
> 11:35:44 # comma or space separated. *
> 11:35:44 # *
> 11:35:44 # To display additional parameter, you can use the command: *
> 11:35:44 # update to_full *
> 11:35:44 #***********************************************************************
> 11:35:44 #
> 11:35:44 #*******************
> 11:35:44 # Running parameters
> 11:35:44 #*******************
> 11:35:44 #
> 11:35:44 #***********************************************************************
> 11:35:44 # Tag name for the run (one word) *
> 11:35:44 #***********************************************************************
> 11:35:44 tag_1 = run_tag ! name of the run
> 11:35:44 #***********************************************************************
> 11:35:44 # Number of LHE events (and their normalization) and the required *
> 11:35:44 # (relative) accuracy on the Xsec. *
> 11:35:44 # These values are ignored for fixed order runs *
> 11:35:44 #***********************************************************************
> 11:35:44 5500 = nevents ! Number of unweighted events requested
> 11:35:44 -1.0 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 11:35:44 -1 = nevt_job! Max number of events per job in event generation.
> 11:35:44 ! (-1= no split).
> 11:35:44 #***********************************************************************
> 11:35:44 # Normalize the weights of LHE events such that they sum or average to *
> 11:35:44 # the total cross section *
> 11:35:44 #***********************************************************************
> 11:35:44 average = event_norm ! valid settings: average, sum, bias
> 11:35:44 #***********************************************************************
> 11:35:44 # Number of points per itegration channel (ignored for aMC@NLO runs) *
> 11:35:44 #***********************************************************************
> 11:35:44 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
> 11:35:44 ! number of points and iter. below)
> 11:35:44 # These numbers are ignored except if req_acc_FO is equal to -1
> 11:35:44 5000 = npoints_FO_grid ! number of points to setup grids
> 11:35:44 4 = niters_FO_grid ! number of iter. to setup grids
> 11:35:44 10000 = npoints_FO ! number of points to compute Xsec
> 11:35:44 6 = niters_FO ! number of iter. to compute Xsec
> 11:35:44 #***********************************************************************
> 11:35:44 # Random number seed *
> 11:35:44 #***********************************************************************
> 11:35:44 1 = iseed ! rnd seed (0=assigned automatically=default))
> 11:35:44 #***********************************************************************
> 11:35:44 # Collider type and energy *
> 11:35:44 #***********************************************************************
> 11:35:44 1 = lpp1 ! beam 1 type (0 = no PDF)
> 11:35:44 1 = lpp2 ! beam 2 type (0 = no PDF)
> 11:35:44 6500 = ebeam1 ! beam 1 energy in GeV
> 11:35:44 6500 = ebeam2 ! beam 2 energy in GeV
> 11:35:44 #***********************************************************************
> 11:35:44 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> 11:35:44 #***********************************************************************
> 11:35:44 lhapdf = pdlabel ! PDF set
> 11:35:44 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
> 11:35:44 ! numbers for central PDF sets are allowed. Can be a list;
> 11:35:44 ! PDF sets beyond the first are included via reweighting.
> 11:35:44 #***********************************************************************
> 11:35:44 # Include the NLO Monte Carlo subtr. terms for the following parton *
> 11:35:44 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> 11:35:44 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> 11:35:44 #***********************************************************************
> 11:35:44 HERWIGPP = parton_shower
> 11:35:44 1.0 = shower_scale_factor ! multiply default shower starting
> 11:35:44 ! scale by this factor
> 11:35:44 #***********************************************************************
> 11:35:44 # Renormalization and factorization scales *
> 11:35:44 # (Default functional form for the non-fixed scales is the sum of *
> 11:35:44 # the transverse masses divided by two of all final state particles *
> 11:35:44 # and partons. This can be changed in SubProcesses/set_scales.f or via *
> 11:35:44 # dynamical_scale_choice option) *
> 11:35:44 #***********************************************************************
> 11:35:44 True = fixed_ren_scale ! if .true. use fixed ren scale
> 11:35:44 True = fixed_fac_scale ! if .true. use fixed fac scale
> 11:35:44 172.5 = muR_ref_fixed ! fixed ren reference scale
> 11:35:44 172.5 = muF_ref_fixed ! fixed fact reference scale
> 11:35:44 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
> 11:35:44 ! dynamical choices. Can be a list; scale choices beyond the
> 11:35:44 ! first are included via reweighting
> 11:35:44 1.0 = muR_over_ref ! ratio of current muR over reference muR
> 11:35:44 1.0 = muF_over_ref ! ratio of current muF over reference muF
> 11:35:44 #***********************************************************************
> 11:35:44 # Reweight variables for scale dependence and PDF uncertainty *
> 11:35:44 #***********************************************************************
> 11:35:44 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
> 11:35:44 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
> 11:35:44 False = reweight_scale ! Reweight to get scale variation using the
> 11:35:44 ! rw_rscale and rw_fscale factors. Should be a list of
> 11:35:44 ! booleans of equal length to dynamical_scale_choice to
> 11:35:44 ! specify for which choice to include scale dependence.
> 11:35:44 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
> 11:35:44 ! list booleans of equal length to lhaid to specify for
> 11:35:44 ! which PDF set to include the uncertainties.
> 11:35:44 #***********************************************************************
> 11:35:44 # Store reweight information in the LHE file for off-line model- *
> 11:35:44 # parameter reweighting at NLO+PS accuracy *
> 11:35:44 #***********************************************************************
> 11:35:44 False = store_rwgt_info ! Store info for reweighting in LHE file
> 11:35:44 #***********************************************************************
> 11:35:44 # ickkw parameter: *
> 11:35:44 # 0: No merging *
> 11:35:44 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
> 11:35:44 # level. After showering an MLM-type merging should be applied as *
> 11:35:44 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
> 11:35:44 # 4: UNLOPS merging (with pythia8 only). No interface from within *
> 11:35:44 # MG5_aMC available, but available in Pythia8. *
> 11:35:44 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
> 11:35:44 #***********************************************************************
> 11:35:44 0 = ickkw
> 11:35:44 #***********************************************************************
> 11:35:44 #
> 11:35:44 #***********************************************************************
> 11:35:44 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
> 11:35:44 # written in the LHE event file *
> 11:35:44 #***********************************************************************
> 11:35:44 50 = bwcutoff
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts on the jets. Jet clustering is performed by FastJet. *
> 11:35:44 # - When matching to a parton shower, these generation cuts should be *
> 11:35:44 # considerably softer than the analysis cuts. *
> 11:35:44 # - More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 11:35:44 0.4 = jetradius ! The radius parameter for the jet algorithm
> 11:35:44 0.1 = ptj ! Min jet transverse momentum
> 11:35:44 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> 11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 0.0 = ptl ! Min lepton transverse momentum
> 11:35:44 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:44 0.0 = drll ! Min distance between opposite sign lepton pairs
> 11:35:44 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 11:35:44 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 11:35:44 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> 11:35:44 #***********************************************************************
> 11:35:44 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
> 11:35:44 # all the other parameters are ignored. *
> 11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 20.0 = ptgmin ! Min photon transverse momentum
> 11:35:44 -1.0 = etagamma ! Max photon abs(pseudo-rap)
> 11:35:44 0.4 = R0gamma ! Radius of isolation code
> 11:35:44 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 11:35:44 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> 11:35:44 True = isoEM ! isolate photons from EM energy (photons and leptons)
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts associated to MASSIVE particles identified by their PDG codes. *
> 11:35:44 # All cuts are applied to both particles and anti-particles, so use *
> 11:35:44 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
> 11:35:44 # {6:100, 25:200} for multiple particles *
> 11:35:44 #***********************************************************************
> 11:35:44 {} = pt_min_pdg ! Min pT for a massive particle
> 11:35:44 {} = pt_max_pdg ! Max pT for a massive particle
> 11:35:44 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
> 11:35:44 #***********************************************************************
> 11:35:44 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> 11:35:44 #***********************************************************************
> 11:35:44 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> 11:35:44 #***********************************************************************
> 11:35:44 Py:MadGraphUtils INFO param_card:
> 11:35:44 ######################################################################
> 11:35:44 ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> 11:35:44 ######################################################################
> 11:35:44 ## ##
> 11:35:44 ## Width set on Auto will be computed following the information ##
> 11:35:44 ## present in the decay.py files of the model. ##
> 11:35:44 ## See arXiv:1402.1178 for more details. ##
> 11:35:44 ## ##
> 11:35:44 ######################################################################
> 11:35:44
> 11:35:44 ###################################
> 11:35:44 ## INFORMATION FOR LOOP
> 11:35:44 ###################################
> 11:35:44 Block loop
> 11:35:44 1 9.118800e+01 # MU_R
> 11:35:44
> 11:35:44 ###################################
> 11:35:44 ## INFORMATION FOR MASS
> 11:35:44 ###################################
> 11:35:44 Block mass
> 11:35:44 6 1.725000e+02 # MT
> 11:35:44 15 1.776820e+00 # MTA
> 11:35:44 23 9.118760e+01 # MZ
> 11:35:44 25 1.250000e+02 # MH
> 11:35:44 ## Dependent parameters, given by model restrictions.
> 11:35:44 ## Those values should be edited following the
> 11:35:44 ## analytical expression. MG5 ignores those values
> 11:35:44 ## but they are important for interfacing the output of MG5
> 11:35:44 ## to external program such as Pythia.
> 11:35:44 1 0.000000 # d : 0.0
> 11:35:44 2 0.000000 # u : 0.0
> 11:35:44 3 0.000000 # s : 0.0
> 11:35:44 4 0.000000 # c : 0.0
> 11:35:44 5 0.000000 # b : 0.0
> 11:35:44 11 0.000000 # e- : 0.0
> 11:35:44 12 0.000000 # ve : 0.0
> 11:35:44 13 0.000000 # mu- : 0.0
> 11:35:44 14 0.000000 # vm : 0.0
> 11:35:44 16 0.000000 # vt : 0.0
> 11:35:44 21 0.000000 # g : 0.0
> 11:35:44 22 0.000000 # a : 0.0
> 11:35:44 24 80.399000 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
> 11:35:44
> 11:35:44 ###################################
> 11:35:44 ## INFORMATION FOR SMINPUTS
> 11:35:44 ###################################
> 11:35:44 Block sminputs
> 11:35:44 1 1.323489e+02 # aEWM1
> 11:35:44 2 1.166370e-05 # Gf
> 11:35:44 3 1.184000e-01 # aS
> 11:35:44
> 11:35:44 ###################################
> 11:35:44 ## INFORMATION FOR YUKAWA
> 11:35:44 ###################################
> 11:35:44 Block yukawa
> 11:35:44 6 1.725000e+02 # ymt
> 11:35:44 15 1.777000e+00 # ymtau
> 11:35:44
> 11:35:44 ###################################
> 11:35:44
> 11:35:44 #*************************
> 11:35:44 # Decay widths *
> 11:35:44 #*************************
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 1 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 2 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 3 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 4 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 5 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 6 1.320000e+00
> 11:35:44 # BR NDA ID1 ID2 ...
> 11:35:44 1.000000e+00 2 5 24 # 1.32
> 11:35:44 DECAY -6 1.320000e+00
> 11:35:44 # BR NDA ID1 ID2 ...
> 11:35:44 1.000000e+00 2 -5 -24 # 1.32
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 11 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 12 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 13 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 14 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 15 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 16 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 21 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 22 0.000000e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 23 2.495200e+00
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY -24 2.085000e+00
> 11:35:44 # BR NDA ID1 ID2 ...
> 11:35:44 3.3370000e-01 2 1 -2
> 11:35:44 3.337000e-01 2 3 -4
> 11:35:44 1.082000e-01 2 11 -12
> 11:35:44 1.082000e-01 2 13 -14
> 11:35:44 1.082000e-01 2 15 -16
> 11:35:44 DECAY 24 2.085000e+00
> 11:35:44 # BR NDA ID1 ID2 ...
> 11:35:44 3.337000e-01 2 -1 2
> 11:35:44 3.337000e-01 2 -3 4
> 11:35:44 1.082000e-01 2 -11 12
> 11:35:44 1.082000e-01 2 -13 14
> 11:35:44 1.082000e-01 2 -15 16
> 11:35:44 #
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 25 6.382339e-03
> 11:35:44 #
> 11:35:44 # PDG Width
> 11:35:44 DECAY 82 0.000000e+00
> 11:35:44 Py:MadGraphUtils INFO Found MadGraph/MadGraph5_aMC@NLO version 2.6.2
> 11:35:44 Py:MadGraphUtils INFO Started generating gridpack at Tue Aug 21 11:35:44 2018
> 11:35:44 Py:MadGraphUtils WARNING >>>>>> THIS KIND OF JOB SHOULD ONLY BE RUN LOCALLY - NOT IN GRID JOBS <<<<<<
> 11:35:44 Py:MadGraphUtils INFO Setting nevents = 1.
> 11:35:44 Py:MadGraphUtils INFO Setting req_acc = 0.001000
> 11:35:44 Py:MadGraphUtils INFO run_card.dat: run_card.dat
> 11:35:44 #***********************************************************************
> 11:35:44 # MadGraph5_aMC@NLO *
> 11:35:44 # *
> 11:35:44 # run_card.dat aMC@NLO *
> 11:35:44 # *
> 11:35:44 # This file is used to set the parameters of the run. *
> 11:35:44 # *
> 11:35:44 # Some notation/conventions: *
> 11:35:44 # *
> 11:35:44 # Lines starting with a hash (#) are info or comments *
> 11:35:44 # *
> 11:35:44 # mind the format: value = variable ! comment *
> 11:35:44 # *
> 11:35:44 # Some of the values of variables can be list. These can either be *
> 11:35:44 # comma or space separated. *
> 11:35:44 # *
> 11:35:44 # To display additional parameter, you can use the command: *
> 11:35:44 # update to_full *
> 11:35:44 #***********************************************************************
> 11:35:44 #
> 11:35:44 #*******************
> 11:35:44 # Running parameters
> 11:35:44 #*******************
> 11:35:44 #
> 11:35:44 #***********************************************************************
> 11:35:44 # Tag name for the run (one word) *
> 11:35:44 #***********************************************************************
> 11:35:44 tag_1 = run_tag ! name of the run
> 11:35:44 #***********************************************************************
> 11:35:44 # Number of LHE events (and their normalization) and the required *
> 11:35:44 # (relative) accuracy on the Xsec. *
> 11:35:44 # These values are ignored for fixed order runs *
> 11:35:44 #***********************************************************************
> 11:35:44 1 = nevents ! Number of unweighted events requested
> 11:35:44 0.001000 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 11:35:44 -1 = nevt_job! Max number of events per job in event generation.
> 11:35:44 ! (-1= no split).
> 11:35:44 #***********************************************************************
> 11:35:44 # Normalize the weights of LHE events such that they sum or average to *
> 11:35:44 # the total cross section *
> 11:35:44 #***********************************************************************
> 11:35:44 average = event_norm ! valid settings: average, sum, bias
> 11:35:44 #***********************************************************************
> 11:35:44 # Number of points per itegration channel (ignored for aMC@NLO runs) *
> 11:35:44 #***********************************************************************
> 11:35:44 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
> 11:35:44 ! number of points and iter. below)
> 11:35:44 # These numbers are ignored except if req_acc_FO is equal to -1
> 11:35:44 5000 = npoints_FO_grid ! number of points to setup grids
> 11:35:44 4 = niters_FO_grid ! number of iter. to setup grids
> 11:35:44 10000 = npoints_FO ! number of points to compute Xsec
> 11:35:44 6 = niters_FO ! number of iter. to compute Xsec
> 11:35:44 #***********************************************************************
> 11:35:44 # Random number seed *
> 11:35:44 #***********************************************************************
> 11:35:44 1 = iseed ! rnd seed (0=assigned automatically=default))
> 11:35:44 #***********************************************************************
> 11:35:44 # Collider type and energy *
> 11:35:44 #***********************************************************************
> 11:35:44 1 = lpp1 ! beam 1 type (0 = no PDF)
> 11:35:44 1 = lpp2 ! beam 2 type (0 = no PDF)
> 11:35:44 6500 = ebeam1 ! beam 1 energy in GeV
> 11:35:44 6500 = ebeam2 ! beam 2 energy in GeV
> 11:35:44 #***********************************************************************
> 11:35:44 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> 11:35:44 #***********************************************************************
> 11:35:44 lhapdf = pdlabel ! PDF set
> 11:35:44 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
> 11:35:44 ! numbers for central PDF sets are allowed. Can be a list;
> 11:35:44 ! PDF sets beyond the first are included via reweighting.
> 11:35:44 #***********************************************************************
> 11:35:44 # Include the NLO Monte Carlo subtr. terms for the following parton *
> 11:35:44 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> 11:35:44 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> 11:35:44 #***********************************************************************
> 11:35:44 HERWIGPP = parton_shower
> 11:35:44 1.0 = shower_scale_factor ! multiply default shower starting
> 11:35:44 ! scale by this factor
> 11:35:44 #***********************************************************************
> 11:35:44 # Renormalization and factorization scales *
> 11:35:44 # (Default functional form for the non-fixed scales is the sum of *
> 11:35:44 # the transverse masses divided by two of all final state particles *
> 11:35:44 # and partons. This can be changed in SubProcesses/set_scales.f or via *
> 11:35:44 # dynamical_scale_choice option) *
> 11:35:44 #***********************************************************************
> 11:35:44 True = fixed_ren_scale ! if .true. use fixed ren scale
> 11:35:44 True = fixed_fac_scale ! if .true. use fixed fac scale
> 11:35:44 172.5 = muR_ref_fixed ! fixed ren reference scale
> 11:35:44 172.5 = muF_ref_fixed ! fixed fact reference scale
> 11:35:44 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
> 11:35:44 ! dynamical choices. Can be a list; scale choices beyond the
> 11:35:44 ! first are included via reweighting
> 11:35:44 1.0 = muR_over_ref ! ratio of current muR over reference muR
> 11:35:44 1.0 = muF_over_ref ! ratio of current muF over reference muF
> 11:35:44 #***********************************************************************
> 11:35:44 # Reweight variables for scale dependence and PDF uncertainty *
> 11:35:44 #***********************************************************************
> 11:35:44 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
> 11:35:44 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
> 11:35:44 False = reweight_scale ! Reweight to get scale variation using the
> 11:35:44 ! rw_rscale and rw_fscale factors. Should be a list of
> 11:35:44 ! booleans of equal length to dynamical_scale_choice to
> 11:35:44 ! specify for which choice to include scale dependence.
> 11:35:44 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
> 11:35:44 ! list booleans of equal length to lhaid to specify for
> 11:35:44 ! which PDF set to include the uncertainties.
> 11:35:44 #***********************************************************************
> 11:35:44 # Store reweight information in the LHE file for off-line model- *
> 11:35:44 # parameter reweighting at NLO+PS accuracy *
> 11:35:44 #***********************************************************************
> 11:35:44 False = store_rwgt_info ! Store info for reweighting in LHE file
> 11:35:44 #***********************************************************************
> 11:35:44 # ickkw parameter: *
> 11:35:44 # 0: No merging *
> 11:35:44 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
> 11:35:44 # level. After showering an MLM-type merging should be applied as *
> 11:35:44 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
> 11:35:44 # 4: UNLOPS merging (with pythia8 only). No interface from within *
> 11:35:44 # MG5_aMC available, but available in Pythia8. *
> 11:35:44 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
> 11:35:44 #***********************************************************************
> 11:35:44 0 = ickkw
> 11:35:44 #***********************************************************************
> 11:35:44 #
> 11:35:44 #***********************************************************************
> 11:35:44 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
> 11:35:44 # written in the LHE event file *
> 11:35:44 #***********************************************************************
> 11:35:44 50 = bwcutoff
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts on the jets. Jet clustering is performed by FastJet. *
> 11:35:44 # - When matching to a parton shower, these generation cuts should be *
> 11:35:44 # considerably softer than the analysis cuts. *
> 11:35:44 # - More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 11:35:44 0.4 = jetradius ! The radius parameter for the jet algorithm
> 11:35:44 0.1 = ptj ! Min jet transverse momentum
> 11:35:44 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> 11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 0.0 = ptl ! Min lepton transverse momentum
> 11:35:44 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:44 0.0 = drll ! Min distance between opposite sign lepton pairs
> 11:35:44 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 11:35:44 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 11:35:44 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> 11:35:44 #***********************************************************************
> 11:35:44 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
> 11:35:44 # all the other parameters are ignored. *
> 11:35:44 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:44 #***********************************************************************
> 11:35:44 20.0 = ptgmin ! Min photon transverse momentum
> 11:35:44 -1.0 = etagamma ! Max photon abs(pseudo-rap)
> 11:35:44 0.4 = R0gamma ! Radius of isolation code
> 11:35:44 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 11:35:44 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> 11:35:44 True = isoEM ! isolate photons from EM energy (photons and leptons)
> 11:35:44 #***********************************************************************
> 11:35:44 # Cuts associated to MASSIVE particles identified by their PDG codes. *
> 11:35:44 # All cuts are applied to both particles and anti-particles, so use *
> 11:35:44 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
> 11:35:44 # {6:100, 25:200} for multiple particles *
> 11:35:44 #***********************************************************************
> 11:35:44 {} = pt_min_pdg ! Min pT for a massive particle
> 11:35:44 {} = pt_max_pdg ! Max pT for a massive particle
> 11:35:44 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
> 11:35:44 #***********************************************************************
> 11:35:44 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> 11:35:44 #***********************************************************************
> 11:35:44 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> 11:35:44 #***********************************************************************
> 11:35:44 Py:MadGraphUtils INFO Run schan will be performed in mode 2 with 8 jobs in parallel.
> 11:35:44 Py:MadGraphUtils INFO Decaying with MadSpin. Moving card (madspin_card.dat) into place.
> 11:35:44 Py:MadGraphUtils INFO creating local LHAPDF dir: MGC_LHAPDF/
> 11:35:44 Py:MadGraphUtils INFO Found LHAPDF ID=260000, name=NNPDF30_nlo_as_0118!
> 11:35:44 Py:MadGraphUtils INFO linking /cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current//NNPDF30_nlo_as_0118 --> MGC_LHAPDF/NNPDF30_nlo_as_0118
> 11:35:44 Py:MadGraphUtils INFO linking /cvmfs/sft.cern.ch/lcg/external/lhapdfsets/current//pdfsets.index --> MGC_LHAPDF/pdfsets.index
> 11:35:44 Py:MadGraphUtils INFO linking /cvmfs/atlas.cern.ch/repo/sw/Generators/lhapdfsets/current//lhapdf.conf --> MGC_LHAPDF/lhapdf.conf
> 11:35:44 Py:MadGraphUtils INFO Path to LHAPDF install dir:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt/share/LHAPDF
> 11:35:44 Py:MadGraphUtils INFO Path to LHAPDF data dir: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/MGC_LHAPDF
> 11:35:46 Willing to acquire file lhapdf-config from
> 11:35:46 /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/share/lhapdf-config
> 11:35:46 copy lhapdf-config
> 11:35:48 Py:MadGraphUtils INFO lhapdf-config --version: 6.1.5
> 11:35:49 Py:MadGraphUtils INFO lhapdf-config --prefix: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt
> 11:35:49 Py:MadGraphUtils INFO lhapdf-config --libdir: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/lhapdf/6.1.5/x86_64-slc6-gcc47-opt/lib
> 11:35:50 Py:MadGraphUtils INFO lhapdf-config --datadir: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/MGC_LHAPDF
> 11:35:50 Py:MadGraphUtils INFO lhapdf-config --pdfsets-path:
> 11:35:50 Py:MadGraphUtils INFO Setting lhapdf = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/lhapdf-config in PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:35:50 Py:MadGraphUtils INFO Creating links for LHAPDF
> 11:35:50 Py:MadGraphUtils INFO Available PDFs are:
> 11:35:50 Py:MadGraphUtils INFO ['NNPDF30_nlo_as_0118', 'lhapdf.conf', 'pdfsets.index']
> 11:35:50 Py:MadGraphUtils INFO For your information, the libraries available are (should include LHAPDF):
> 11:35:50 Py:MadGraphUtils INFO ['PDFsets', 'Pdfdata', 'libcts.a', 'libiregi.a', 'mpmodule.mod']
> 11:35:50 Py:MadGraphUtils INFO Path to fastjet install dir:/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc47-opt
> 11:35:52 Willing to acquire file fastjet-config from
> 11:35:52 /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/share/fastjet-config
> 11:35:52 copy fastjet-config
> 11:35:52 Py:MadGraphUtils INFO fastjet-config --version: 3.0.3
> 11:35:52 Py:MadGraphUtils INFO fastjet-config --prefix: /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/fastjet/3.0.3/x86_64-slc6-gcc47-opt
> 11:35:52 Py:MadGraphUtils INFO Setting fastjet = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/fastjet-config in PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:35:52 Py:MadGraphUtils INFO Now I will hack the make files a bit. Apologies, but there seems to be no good way around this.
> 11:35:52 Py:MadGraphUtils INFO Configuring the fancy gfortran compiler instead of g77 / f77
> 11:35:52 Py:MadGraphUtils INFO Make file hacking complete.
> 11:35:52 Py:MadGraphUtils INFO proc_card:
> 11:35:52
> 11:35:52 import model loop_sm-no_b_mass
> 11:35:52 define l+ = l+ ta+
> 11:35:52 define l- = l- ta-
> 11:35:52 generate p p > w+ > t b~ [QCD]
> 11:35:52 add process p p > w- > t~ b [QCD]
> 11:35:52 output -f
> 11:35:52 Py:MadGraphUtils INFO run_card:
> 11:35:52 #***********************************************************************
> 11:35:52 # MadGraph5_aMC@NLO *
> 11:35:52 # *
> 11:35:52 # run_card.dat aMC@NLO *
> 11:35:52 # *
> 11:35:52 # This file is used to set the parameters of the run. *
> 11:35:52 # *
> 11:35:52 # Some notation/conventions: *
> 11:35:52 # *
> 11:35:52 # Lines starting with a hash (#) are info or comments *
> 11:35:52 # *
> 11:35:52 # mind the format: value = variable ! comment *
> 11:35:52 # *
> 11:35:52 # Some of the values of variables can be list. These can either be *
> 11:35:52 # comma or space separated. *
> 11:35:52 # *
> 11:35:52 # To display additional parameter, you can use the command: *
> 11:35:52 # update to_full *
> 11:35:52 #***********************************************************************
> 11:35:52 #
> 11:35:52 #*******************
> 11:35:52 # Running parameters
> 11:35:52 #*******************
> 11:35:52 #
> 11:35:52 #***********************************************************************
> 11:35:52 # Tag name for the run (one word) *
> 11:35:52 #***********************************************************************
> 11:35:52 tag_1 = run_tag ! name of the run
> 11:35:52 #***********************************************************************
> 11:35:52 # Number of LHE events (and their normalization) and the required *
> 11:35:52 # (relative) accuracy on the Xsec. *
> 11:35:52 # These values are ignored for fixed order runs *
> 11:35:52 #***********************************************************************
> 11:35:52 1 = nevents ! Number of unweighted events requested
> 11:35:52 0.001000 = req_acc ! Required accuracy (-1=auto determined from nevents)
> 11:35:52 -1 = nevt_job! Max number of events per job in event generation.
> 11:35:52 ! (-1= no split).
> 11:35:52 #***********************************************************************
> 11:35:52 # Normalize the weights of LHE events such that they sum or average to *
> 11:35:52 # the total cross section *
> 11:35:52 #***********************************************************************
> 11:35:52 average = event_norm ! valid settings: average, sum, bias
> 11:35:52 #***********************************************************************
> 11:35:52 # Number of points per itegration channel (ignored for aMC@NLO runs) *
> 11:35:52 #***********************************************************************
> 11:35:52 0.01 = req_acc_FO ! Required accuracy (-1=ignored, and use the
> 11:35:52 ! number of points and iter. below)
> 11:35:52 # These numbers are ignored except if req_acc_FO is equal to -1
> 11:35:52 5000 = npoints_FO_grid ! number of points to setup grids
> 11:35:52 4 = niters_FO_grid ! number of iter. to setup grids
> 11:35:52 10000 = npoints_FO ! number of points to compute Xsec
> 11:35:52 6 = niters_FO ! number of iter. to compute Xsec
> 11:35:52 #***********************************************************************
> 11:35:52 # Random number seed *
> 11:35:52 #***********************************************************************
> 11:35:52 1 = iseed ! rnd seed (0=assigned automatically=default))
> 11:35:52 #***********************************************************************
> 11:35:52 # Collider type and energy *
> 11:35:52 #***********************************************************************
> 11:35:52 1 = lpp1 ! beam 1 type (0 = no PDF)
> 11:35:52 1 = lpp2 ! beam 2 type (0 = no PDF)
> 11:35:52 6500 = ebeam1 ! beam 1 energy in GeV
> 11:35:52 6500 = ebeam2 ! beam 2 energy in GeV
> 11:35:52 #***********************************************************************
> 11:35:52 # PDF choice: this automatically fixes also alpha_s(MZ) and its evol. *
> 11:35:52 #***********************************************************************
> 11:35:52 lhapdf = pdlabel ! PDF set
> 11:35:52 260000 = lhaid ! If pdlabel=lhapdf, this is the lhapdf number. Only
> 11:35:52 ! numbers for central PDF sets are allowed. Can be a list;
> 11:35:52 ! PDF sets beyond the first are included via reweighting.
> 11:35:52 #***********************************************************************
> 11:35:52 # Include the NLO Monte Carlo subtr. terms for the following parton *
> 11:35:52 # shower (HERWIG6 | HERWIGPP | PYTHIA6Q | PYTHIA6PT | PYTHIA8) *
> 11:35:52 # WARNING: PYTHIA6PT works only for processes without FSR!!!! *
> 11:35:52 #***********************************************************************
> 11:35:52 HERWIGPP = parton_shower
> 11:35:52 1.0 = shower_scale_factor ! multiply default shower starting
> 11:35:52 ! scale by this factor
> 11:35:52 #***********************************************************************
> 11:35:52 # Renormalization and factorization scales *
> 11:35:52 # (Default functional form for the non-fixed scales is the sum of *
> 11:35:52 # the transverse masses divided by two of all final state particles *
> 11:35:52 # and partons. This can be changed in SubProcesses/set_scales.f or via *
> 11:35:52 # dynamical_scale_choice option) *
> 11:35:52 #***********************************************************************
> 11:35:52 True = fixed_ren_scale ! if .true. use fixed ren scale
> 11:35:52 True = fixed_fac_scale ! if .true. use fixed fac scale
> 11:35:52 172.5 = muR_ref_fixed ! fixed ren reference scale
> 11:35:52 172.5 = muF_ref_fixed ! fixed fact reference scale
> 11:35:52 -1 = dynamical_scale_choice ! Choose one (or more) of the predefined
> 11:35:52 ! dynamical choices. Can be a list; scale choices beyond the
> 11:35:52 ! first are included via reweighting
> 11:35:52 1.0 = muR_over_ref ! ratio of current muR over reference muR
> 11:35:52 1.0 = muF_over_ref ! ratio of current muF over reference muF
> 11:35:52 #***********************************************************************
> 11:35:52 # Reweight variables for scale dependence and PDF uncertainty *
> 11:35:52 #***********************************************************************
> 11:35:52 1.0, 2.0, 0.5 = rw_rscale ! muR factors to be included by reweighting
> 11:35:52 1.0, 2.0, 0.5 = rw_fscale ! muF factors to be included by reweighting
> 11:35:52 False = reweight_scale ! Reweight to get scale variation using the
> 11:35:52 ! rw_rscale and rw_fscale factors. Should be a list of
> 11:35:52 ! booleans of equal length to dynamical_scale_choice to
> 11:35:52 ! specify for which choice to include scale dependence.
> 11:35:52 False = reweight_PDF ! Reweight to get PDF uncertainty. Should be a
> 11:35:52 ! list booleans of equal length to lhaid to specify for
> 11:35:52 ! which PDF set to include the uncertainties.
> 11:35:52 #***********************************************************************
> 11:35:52 # Store reweight information in the LHE file for off-line model- *
> 11:35:52 # parameter reweighting at NLO+PS accuracy *
> 11:35:52 #***********************************************************************
> 11:35:52 False = store_rwgt_info ! Store info for reweighting in LHE file
> 11:35:52 #***********************************************************************
> 11:35:52 # ickkw parameter: *
> 11:35:52 # 0: No merging *
> 11:35:52 # 3: FxFx Merging - WARNING! Applies merging only at the hard-event *
> 11:35:52 # level. After showering an MLM-type merging should be applied as *
> 11:35:52 # well. See http://amcatnlo.cern.ch/FxFx_merging.htm for details. *
> 11:35:52 # 4: UNLOPS merging (with pythia8 only). No interface from within *
> 11:35:52 # MG5_aMC available, but available in Pythia8. *
> 11:35:52 # -1: NNLL+NLO jet-veto computation. See arxiv:1412.8408 [hep-ph]. *
> 11:35:52 #***********************************************************************
> 11:35:52 0 = ickkw
> 11:35:52 #***********************************************************************
> 11:35:52 #
> 11:35:52 #***********************************************************************
> 11:35:52 # BW cutoff (M+/-bwcutoff*Gamma). Determines which resonances are *
> 11:35:52 # written in the LHE event file *
> 11:35:52 #***********************************************************************
> 11:35:52 50 = bwcutoff
> 11:35:52 #***********************************************************************
> 11:35:52 # Cuts on the jets. Jet clustering is performed by FastJet. *
> 11:35:52 # - When matching to a parton shower, these generation cuts should be *
> 11:35:52 # considerably softer than the analysis cuts. *
> 11:35:52 # - More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:52 #***********************************************************************
> 11:35:52 -1 = jetalgo ! FastJet jet algorithm (1=kT, 0=C/A, -1=anti-kT)
> 11:35:52 0.4 = jetradius ! The radius parameter for the jet algorithm
> 11:35:52 0.1 = ptj ! Min jet transverse momentum
> 11:35:52 -1.0 = etaj ! Max jet abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:52 #***********************************************************************
> 11:35:52 # Cuts on the charged leptons (e+, e-, mu+, mu-, tau+ and tau-) *
> 11:35:52 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:52 #***********************************************************************
> 11:35:52 0.0 = ptl ! Min lepton transverse momentum
> 11:35:52 -1.0 = etal ! Max lepton abs(pseudo-rap) (a value .lt.0 means no cut)
> 11:35:52 0.0 = drll ! Min distance between opposite sign lepton pairs
> 11:35:52 0.0 = drll_sf ! Min distance between opp. sign same-flavor lepton pairs
> 11:35:52 0.0 = mll ! Min inv. mass of all opposite sign lepton pairs
> 11:35:52 30.0 = mll_sf ! Min inv. mass of all opp. sign same-flavor lepton pairs
> 11:35:52 #***********************************************************************
> 11:35:52 # Photon-isolation cuts, according to hep-ph/9801442. When ptgmin=0, *
> 11:35:52 # all the other parameters are ignored. *
> 11:35:52 # More specific cuts can be specified in SubProcesses/cuts.f *
> 11:35:52 #***********************************************************************
> 11:35:52 20.0 = ptgmin ! Min photon transverse momentum
> 11:35:52 -1.0 = etagamma ! Max photon abs(pseudo-rap)
> 11:35:52 0.4 = R0gamma ! Radius of isolation code
> 11:35:52 1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
> 11:35:52 1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
> 11:35:52 True = isoEM ! isolate photons from EM energy (photons and leptons)
> 11:35:52 #***********************************************************************
> 11:35:52 # Cuts associated to MASSIVE particles identified by their PDG codes. *
> 11:35:52 # All cuts are applied to both particles and anti-particles, so use *
> 11:35:52 # POSITIVE PDG CODES only. Example of the syntax is {6 : 100} or *
> 11:35:52 # {6:100, 25:200} for multiple particles *
> 11:35:52 #***********************************************************************
> 11:35:52 {} = pt_min_pdg ! Min pT for a massive particle
> 11:35:52 {} = pt_max_pdg ! Max pT for a massive particle
> 11:35:52 {} = mxx_min_pdg ! inv. mass for any pair of (anti)particles
> 11:35:52 #***********************************************************************
> 11:35:52 # For aMCfast+APPLGRID use in PDF fitting (http://amcfast.hepforge.org)*
> 11:35:52 #***********************************************************************
> 11:35:52 0 = iappl ! aMCfast switch (0=OFF, 1=prepare grids, 2=fill grids)
> 11:35:52 #***********************************************************************
> 11:35:52 Py:MadGraphUtils INFO param_card:
> 11:35:52 ######################################################################
> 11:35:52 ## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
> 11:35:52 ######################################################################
> 11:35:52 ## ##
> 11:35:52 ## Width set on Auto will be computed following the information ##
> 11:35:52 ## present in the decay.py files of the model. ##
> 11:35:52 ## See arXiv:1402.1178 for more details. ##
> 11:35:52 ## ##
> 11:35:52 ######################################################################
> 11:35:52
> 11:35:52 ###################################
> 11:35:52 ## INFORMATION FOR LOOP
> 11:35:52 ###################################
> 11:35:52 Block loop
> 11:35:52 1 9.118800e+01 # MU_R
> 11:35:52
> 11:35:52 ###################################
> 11:35:52 ## INFORMATION FOR MASS
> 11:35:52 ###################################
> 11:35:52 Block mass
> 11:35:52 6 1.725000e+02 # MT
> 11:35:52 15 1.776820e+00 # MTA
> 11:35:52 23 9.118760e+01 # MZ
> 11:35:52 25 1.250000e+02 # MH
> 11:35:52 ## Dependent parameters, given by model restrictions.
> 11:35:52 ## Those values should be edited following the
> 11:35:52 ## analytical expression. MG5 ignores those values
> 11:35:52 ## but they are important for interfacing the output of MG5
> 11:35:52 ## to external program such as Pythia.
> 11:35:52 1 0.000000 # d : 0.0
> 11:35:52 2 0.000000 # u : 0.0
> 11:35:52 3 0.000000 # s : 0.0
> 11:35:52 4 0.000000 # c : 0.0
> 11:35:52 5 0.000000 # b : 0.0
> 11:35:52 11 0.000000 # e- : 0.0
> 11:35:52 12 0.000000 # ve : 0.0
> 11:35:52 13 0.000000 # mu- : 0.0
> 11:35:52 14 0.000000 # vm : 0.0
> 11:35:52 16 0.000000 # vt : 0.0
> 11:35:52 21 0.000000 # g : 0.0
> 11:35:52 22 0.000000 # a : 0.0
> 11:35:52 24 80.399000 # w+ : cmath.sqrt(MZ__exp__2/2. + cmath.sqrt(MZ__exp__4/4. - (aEW*cmath.pi*MZ__exp__2)/(Gf*sqrt__2)))
> 11:35:52
> 11:35:52 ###################################
> 11:35:52 ## INFORMATION FOR SMINPUTS
> 11:35:52 ###################################
> 11:35:52 Block sminputs
> 11:35:52 1 1.323489e+02 # aEWM1
> 11:35:52 2 1.166370e-05 # Gf
> 11:35:52 3 1.184000e-01 # aS
> 11:35:52
> 11:35:52 ###################################
> 11:35:52 ## INFORMATION FOR YUKAWA
> 11:35:52 ###################################
> 11:35:52 Block yukawa
> 11:35:52 6 1.725000e+02 # ymt
> 11:35:52 15 1.777000e+00 # ymtau
> 11:35:52
> 11:35:52 ###################################
> 11:35:52
> 11:35:52 #*************************
> 11:35:52 # Decay widths *
> 11:35:52 #*************************
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 1 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 2 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 3 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 4 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 5 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 6 1.320000e+00
> 11:35:52 # BR NDA ID1 ID2 ...
> 11:35:52 1.000000e+00 2 5 24 # 1.32
> 11:35:52 DECAY -6 1.320000e+00
> 11:35:52 # BR NDA ID1 ID2 ...
> 11:35:52 1.000000e+00 2 -5 -24 # 1.32
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 11 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 12 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 13 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 14 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 15 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 16 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 21 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 22 0.000000e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 23 2.495200e+00
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY -24 2.085000e+00
> 11:35:52 # BR NDA ID1 ID2 ...
> 11:35:52 3.3370000e-01 2 1 -2
> 11:35:52 3.337000e-01 2 3 -4
> 11:35:52 1.082000e-01 2 11 -12
> 11:35:52 1.082000e-01 2 13 -14
> 11:35:52 1.082000e-01 2 15 -16
> 11:35:52 DECAY 24 2.085000e+00
> 11:35:52 # BR NDA ID1 ID2 ...
> 11:35:52 3.337000e-01 2 -1 2
> 11:35:52 3.337000e-01 2 -3 4
> 11:35:52 1.082000e-01 2 -11 12
> 11:35:52 1.082000e-01 2 -13 14
> 11:35:52 1.082000e-01 2 -15 16
> 11:35:52 #
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 25 6.382339e-03
> 11:35:52 #
> 11:35:52 # PDG Width
> 11:35:52 DECAY 82 0.000000e+00
> 11:35:52 Py:MadGraphUtils INFO Setting automatic_html_opening = False
> 11:35:52 Py:MadGraphUtils INFO Setting run_mode = 2
> 11:35:52 Py:MadGraphUtils INFO Setting number of cores = 8 in Cards/amcatnlo_configuration.txt
> 11:35:52 Py:MadGraphUtils INFO amcatnlo_configuration.txt: Cards/amcatnlo_configuration.txt
> 11:35:52 ################################################################################
> 11:35:52 #
> 11:35:52 # Copyright (c) 2009 The MadGraph5_aMC@NLO Development team and Contributors
> 11:35:52 #
> 11:35:52 # This file is a part of the MadGraph5_aMC@NLO project, an application which
> 11:35:52 # automatically generates Feynman diagrams and matrix elements for arbitrary
> 11:35:52 # high-energy processes in the Standard Model and beyond.
> 11:35:52 #
> 11:35:52 # It is subject to the MadGraph5_aMC@NLO license which should accompany this
> 11:35:52 # distribution.
> 11:35:52 #
> 11:35:52 # For more information, visit madgraph.phys.ucl.ac.be and amcatnlo.web.cern.ch
> 11:35:52 #
> 11:35:52 ################################################################################
> 11:35:52 #
> 11:35:52 # This File contains some configuration variable for MadGraph/MadEvent
> 11:35:52 #
> 11:35:52 # Line starting by #! are comment and should remain commented
> 11:35:52 # Line starting with # should be uncommented if you want to modify the default
> 11:35:52 # value.
> 11:35:52 # Current value for all options can seen by typing "display options"
> 11:35:52 # after either ./bin/mg5_aMC or ./bin/madevent
> 11:35:52 #
> 11:35:52 # You can place this files in ~/.mg5/mg5_configuration.txt if you have more than
> 11:35:52 # one version of MG5.
> 11:35:52 #
> 11:35:52 ################################################################################
> 11:35:52 #! Prefered Fortran Compiler
> 11:35:52 #! If None: try to find g77 or gfortran on the system
> 11:35:52 #!
> 11:35:52 # fortran_compiler = None
> 11:35:52 # f2py_compiler = None
> 11:35:52
> 11:35:52
> 11:35:52 #! Prefered C++ Compiler
> 11:35:52 #! If None: try to find g++ or clang on the system
> 11:35:52 #!
> 11:35:52 # cpp_compiler = None
> 11:35:52
> 11:35:52 #! Prefered Text Editor
> 11:35:52 #! Default: use the shell default Editor
> 11:35:52 #! or try to find one available on the system
> 11:35:52 #! Be careful: Only shell based editor are allowed
> 11:35:52 # text_editor = None
> 11:35:52
> 11:35:52 #! Prefered WebBrower
> 11:35:52 #! If None: try to find one available on the system
> 11:35:52 # web_browser = None
> 11:35:52
> 11:35:52 #! Prefered PS viewer
> 11:35:52 #! If None: try to find one available on the system
> 11:35:52 # eps_viewer = None
> 11:35:52
> 11:35:52 #! Time allowed to answer question (if no answer takes default value)
> 11:35:52 #! 0: No time limit
> 11:35:52 # timeout = 60
> 11:35:52
> 11:35:52 #! Pythia8 path.
> 11:35:52 #! Defines the path to the pythia8 installation directory (i.e. the
> 11:35:52 #! on containing the lib, bin and include directories) .
> 11:35:52 #! If using a relative path, that starts from the mg5 directory
> 11:35:52 # pythia8_path = ./HEPTools/pythia8
> 11:35:52
> 11:35:52 #! MG5aMC_PY8_interface path
> 11:35:52 #! Defines the path of the C++ driver file that is used by MG5_aMC to
> 11:35:52 #! steer the Pythia8 shower.
> 11:35:52 #! Can be installed directly from within MG5_aMC with the following command:
> 11:35:52 #! MG5_aMC> install mg5amc_py8_interface
> 11:35:52 # mg5amc_py8_interface_path = ./HEPTools/MG5aMC_PY8_interface
> 11:35:52
> 11:35:52 #! Herwig++/Herwig7 paths
> 11:35:52 #! specify here the paths also to HepMC ant ThePEG
> 11:35:52 #! define the path to the herwig++, thepeg and hepmc directories.
> 11:35:52 #! paths can be absolute or relative from mg5 directory
> 11:35:52 #! WARNING: if Herwig7 has been installed with the bootstrap script,
> 11:35:52 #! then please set thepeg_path and hepmc_path to the same value as
> 11:35:52 #! hwpp_path
> 11:35:52 # hwpp_path =
> 11:35:52 # thepeg_path =
> 11:35:52 # hepmc_path =
> 11:35:52
> 11:35:52 #! Control when MG5 checks if he is up-to-date.
> 11:35:52 #! Enter the number of day between two check (0 means never)
> 11:35:52 #! A question is always asked before any update
> 11:35:52 # auto_update = 7
> 11:35:52
> 11:35:52 ################################################################################
> 11:35:52 # INFO FOR MADEVENT / aMC@NLO
> 11:35:52 ################################################################################
> 11:35:52 # If this file is in a MADEVENT Template. 'main directory' is the directory
> 11:35:52 # containing the SubProcesses directory. Otherwise this is the MadGraph5_aMC@NLO main
> 11:35:52 # directory (containing the directories madgraph and Template)
> 11:35:52
> 11:35:52 #! Allow/Forbid the automatic opening of the web browser (on the status page)
> 11:35:52 #! when launching MadEvent [True/False]
> 11:35:52 automatic_html_opening = False
> 11:35:52 #! allow notification of finished job in the notification center (Mac Only)
> 11:35:52 # notification_center = True
> 11:35:52
> 11:35:52
> 11:35:52 #! Default Running mode
> 11:35:52 #! 0: single machine/ 1: cluster / 2: multicore
> 11:35:52 run_mode = 2
> 11:35:52
> 11:35:52 #! Cluster Type [pbs|sge|condor|lsf|ge|slurm|htcaas|htcaas2] Use for cluster run only
> 11:35:52 #! And cluster queue (or partition for slurm)
> 11:35:52 #! And size of the cluster (some part of the code can adapt splitting accordingly)
> 11:35:52 # cluster_type = condor
> 11:35:52 # cluster_queue = madgraph
> 11:35:52 # cluster_size = 150
> 11:35:52
> 11:35:52 #! Path to a node directory to avoid direct writing on the central disk
> 11:35:52 #! Note that condor clusters avoid direct writing by default (therefore this
> 11:35:52 #! options does not affect condor clusters)
> 11:35:52 # cluster_temp_path = None
> 11:35:52
> 11:35:52 #! path to a node directory where local file can be found (typically pdf)
> 11:35:52 #! to avoid to send them to the node (if cluster_temp_path is on True or condor)
> 11:35:52 # cluster_local_path = None # example: /cvmfs/cp3.uclouvain.be/madgraph/
> 11:35:52
> 11:35:52 #! Cluster waiting time for status update
> 11:35:52 #! First number is when the number of waiting job is higher than the number
> 11:35:52 #! of running one (time in second). The second number is in the second case.
> 11:35:52 # cluster_status_update = 600 30
> 11:35:52
> 11:35:52 #! How to deal with failed submission (can occurs on cluster mode)
> 11:35:52 #! 0: crash, -1: print error, hangs the program up to manual instructions, N(>0) retry up to N times.
> 11:35:52 # cluster_nb_retry = 1
> 11:35:52
> 11:35:52 #! How much time to wait for the output file before resubmission/crash (filesystem can be very slow)
> 11:35:52 # cluster_retry_wait = 300
> 11:35:52
> 11:35:52 #! Nb_core to use (None = all) This is use only for multicore run
> 11:35:52 #! This correspond also to the number core used for code compilation for cluster mode
> 11:35:52 nb_core = 8
> 11:35:52
> 11:35:52 #! Pythia-PGS Package
> 11:35:52 #! relative path start from main directory
> 11:35:52 # pythia-pgs_path = ./pythia-pgs
> 11:35:52
> 11:35:52 #! Delphes Package
> 11:35:52 #! relative path start from main directory
> 11:35:52 # delphes_path = ./Delphes
> 11:35:52
> 11:35:52 #! MadAnalysis4 fortran-based package [for basic analysis]
> 11:35:52 #! relative path start from main directory
> 11:35:52 # madanalysis_path = ./MadAnalysis
> 11:35:52
> 11:35:52 #! MadAnalysis5 python-based Package [For advanced analysis]
> 11:35:52 #! relative path start from main directory
> 11:35:52 # madanalysis5_path = ./HEPTools/madanalysis5/madanalysis5
> 11:35:52
> 11:35:52 #! ExRootAnalysis Package
> 11:35:52 #! relative path start from main directory
> 11:35:52 # exrootanalysis_path = ./ExRootAnalysis
> 11:35:52
> 11:35:52 #! TOPDRAWER PATH
> 11:35:52 #! Path to the directory containing td executables
> 11:35:52 #! relative path start from main directory
> 11:35:52 # td_path = ./td
> 11:35:52
> 11:35:52 #! lhapdf-config
> 11:35:52 #! If None: try to find one available on the system
> 11:35:52 lhapdf = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/lhapdf-config
> 11:35:52
> 11:35:52 #! fastjet-config
> 11:35:52 #! If None: try to find one available on the system
> 11:35:52 fastjet = /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/fastjet-config
> 11:35:52
> 11:35:52 #! MCatNLO-utilities
> 11:35:52 #! relative path starting from main directory
> 11:35:52 # MCatNLO-utilities_path = ./MCatNLO-utilities
> 11:35:52
> 11:35:52 #! Set what OLP to use for the loop ME generation
> 11:35:52 # OLP = MadLoop
> 11:35:52
> 11:35:52 #! Set the PJFRy++ directory containing pjfry's library
> 11:35:52 #! if auto: try to find it automatically on the system (default)
> 11:35:52 #! if '' or None: disabling pjfry
> 11:35:52 #! if pjfry=/PATH/TO/pjfry/lib: use that specific installation path for PJFry++
> 11:35:52 # pjfry = auto
> 11:35:52
> 11:35:52 #! Set the Golem95 directory containing golem's library
> 11:35:52 #! It only supports version higher than 1.3.0
> 11:35:52 #! if auto: try to find it automatically on the system (default)
> 11:35:52 #! if '' or None: disabling Golem95
> 11:35:52 #! if golem=/PATH/TO/golem/lib: use that speficif installation path for Golem95
> 11:35:52 # golem = auto
> 11:35:52
> 11:35:52 #! Set the samurai directory containing samurai's library
> 11:35:52 #! It only supports version higher than 2.0.0
> 11:35:52 #! if auto: try to find it automatically on the system (default)
> 11:35:52 #! if '' or None: disabling samurai
> 11:35:52 #! if samurai=/PATH/TO/samurai/lib: use that specific installation path for samurai
> 11:35:52 # samurai = None
> 11:35:52
> 11:35:52 #! Set the Ninja directory containing ninja's library
> 11:35:52 #! if '' or None: disabling ninja
> 11:35:52 #! if ninja=/PATH/TO/ninja/lib: use that specific installation path for ninja
> 11:35:52 ninja = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//gosam_contrib/2.0/x86_64-slc6-gcc47-opt/lib
> 11:35:52
> 11:35:52 #! Set the COLLIER directory containing COLLIER's library
> 11:35:52 #! if '' or None: disabling COLLIER
> 11:35:52 #! if ninja=/PATH/TO/ninja/lib: use that specific installation path for COLLIER
> 11:35:52 # Note that it is necessary that you have generated a static library for COLLIER
> 11:35:52 collier = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//collier/1.1/x86_64-slc6-gcc47-opt
> 11:35:52
> 11:35:52 #! Set how MadLoop dependencies (such as CutTools) should be handled
> 11:35:52 #! > external : ML5 places a link to the MG5_aMC-wide libraries
> 11:35:52 #! > internal : ML5 copies all dependencies in the output so that it is independent
> 11:35:52 #! > environment_paths : ML5 searches for the dependencies in your environment path
> 11:35:52 output_dependencies = internal
> 11:35:52
> 11:35:52
> 11:35:52 #! SysCalc PATH
> 11:35:52 #! Path to the directory containing syscalc executables
> 11:35:52 #! relative path start from main directory
> 11:35:52 syscalc_path = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c//syscalc/1.1.4/x86_64-slc6-gcc47-opt
> 11:35:52
> 11:35:52 #! Absolute paths to config scripts in the bin directories for APPLgrid and aMCFast.
> 11:35:52 # applgrid = applgrid-config
> 11:35:52 # amcfast = amcfast-config
> 11:35:52
> 11:35:52
> 11:35:52
> 11:35:52 # MG5 MAIN DIRECTORY
> 11:35:52 mg5_path = /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt
> 11:35:52 Py:MadGraphUtils INFO Removing Cards/shower_card.dat to ensure we get parton level events only
> 11:35:52 Py:MadGraphUtils INFO "rw_rscale" = 1.0, 2.0, 0.5
> 11:35:52 Py:MadGraphUtils INFO "ptl" = 0.0
> 11:35:52 Py:MadGraphUtils INFO "xn" = 1.0
> 11:35:52 Py:MadGraphUtils INFO "ptj" = 0.1
> 11:35:52 Py:MadGraphUtils INFO "niters_FO" = 6
> 11:35:52 Py:MadGraphUtils INFO "nevt_job" = -1
> 11:35:52 Py:MadGraphUtils INFO "drll" = 0.0
> 11:35:52 Py:MadGraphUtils INFO "fixed_ren_scale" = True
> 11:35:52 Py:MadGraphUtils INFO "npoints_FO_grid" = 5000
> 11:35:52 Py:MadGraphUtils INFO "iseed" = 1
> 11:35:52 Py:MadGraphUtils INFO "muR_ref_fixed" = 172.5
> 11:35:52 Py:MadGraphUtils INFO "dynamical_scale_choice" = -1
> 11:35:52 Py:MadGraphUtils INFO "R0gamma" = 0.4
> 11:35:52 Py:MadGraphUtils INFO "reweight_scale" = False
> 11:35:52 Py:MadGraphUtils INFO "store_rwgt_info" = False
> 11:35:52 Py:MadGraphUtils INFO "lpp1" = 1
> 11:35:52 Py:MadGraphUtils INFO "nevents" = 1
> 11:35:52 Py:MadGraphUtils INFO "pt_max_pdg" = {}
> 11:35:52 Py:MadGraphUtils INFO "ickkw" = 0
> 11:35:52 Py:MadGraphUtils INFO "iappl" = 0
> 11:35:52 Py:MadGraphUtils INFO "etal" = -1.0
> 11:35:52 Py:MadGraphUtils INFO "epsgamma" = 1.0
> 11:35:52 Py:MadGraphUtils INFO "req_acc_FO" = 0.01
> 11:35:52 Py:MadGraphUtils INFO "pt_min_pdg" = {}
> 11:35:52 Py:MadGraphUtils INFO "etagamma" = -1.0
> 11:35:52 Py:MadGraphUtils INFO "drll_sf" = 0.0
> 11:35:52 Py:MadGraphUtils INFO "muR_over_ref" = 1.0
> 11:35:52 Py:MadGraphUtils INFO "pdlabel" = lhapdf
> 11:35:52 Py:MadGraphUtils INFO "parton_shower" = HERWIGPP
> 11:35:52 Py:MadGraphUtils INFO "event_norm" = average
> 11:35:52 Py:MadGraphUtils INFO "mll_sf" = 30.0
> 11:35:52 Py:MadGraphUtils INFO "muF_ref_fixed" = 172.5
> 11:35:52 Py:MadGraphUtils INFO "bwcutoff" = 50
> 11:35:52 Py:MadGraphUtils INFO "ptgmin" = 20.0
> 11:35:52 Py:MadGraphUtils INFO "rw_fscale" = 1.0, 2.0, 0.5
> 11:35:52 Py:MadGraphUtils INFO "mxx_min_pdg" = {}
> 11:35:52 Py:MadGraphUtils INFO "npoints_FO" = 10000
> 11:35:52 Py:MadGraphUtils INFO "reweight_PDF" = False
> 11:35:52 Py:MadGraphUtils INFO "muF_over_ref" = 1.0
> 11:35:52 Py:MadGraphUtils INFO "isoEM" = True
> 11:35:52 Py:MadGraphUtils INFO "mll" = 0.0
> 11:35:52 Py:MadGraphUtils INFO "fixed_fac_scale" = True
> 11:35:52 Py:MadGraphUtils INFO "lpp2" = 1
> 11:35:52 Py:MadGraphUtils INFO "lhaid" = 260000
> 11:35:52 Py:MadGraphUtils INFO "niters_FO_grid" = 4
> 11:35:52 Py:MadGraphUtils INFO "etaj" = -1.0
> 11:35:52 Py:MadGraphUtils INFO "jetradius" = 0.4
> 11:35:52 Py:MadGraphUtils INFO "jetalgo" = -1
> 11:35:52 Py:MadGraphUtils INFO "req_acc" = 0.001000
> 11:35:52 Py:MadGraphUtils INFO "shower_scale_factor" = 1.0
> 11:35:52 Py:MadGraphUtils INFO "ebeam2" = 6500
> 11:35:52 Py:MadGraphUtils INFO "run_tag" = tag_1
> 11:35:52 Py:MadGraphUtils INFO "ebeam1" = 6500
> 11:35:52 Py:MadGraphUtils INFO Finished checking run card - All OK!
> 11:35:53 [?1034hNo module named madgraph
> 11:35:53 INFO: ************************************************************
> 11:35:53 * *
> 11:35:53 * W E L C O M E to M A D G R A P H 5 *
> 11:35:53 * a M C @ N L O *
> 11:35:53 * *
> 11:35:53 * * * *
> 11:35:53 * * * * * *
> 11:35:53 * * * * * 5 * * * * *
> 11:35:53 * * * * * *
> 11:35:53 * * * *
> 11:35:53 * *
> 11:35:53 * VERSION 5.2.6.2 20xx-xx-xx *
> 11:35:53 * *
> 11:35:53 * The MadGraph5_aMC@NLO Development Team - Find us at *
> 11:35:53 * http://amcatnlo.cern.ch *
> 11:35:53 * *
> 11:35:53 * Type 'help' for in-line help. *
> 11:35:53 * *
> 11:35:53 ************************************************************
> 11:35:54 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:35:54 INFO: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
> 11:35:54 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:35:54 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> 11:35:54 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> 11:35:54 launch --name=schan
> 11:35:54 stty: standard input: Invalid argument
> 11:35:54 stty: standard input: Invalid argument
> 11:35:54 The following switches determine which programs are run:
> 11:35:54 /=============================================================================\
> 11:35:54 1. Type of perturbative computation  order = NLO
> 11:35:54 2. No MC@[N]LO matching / event generation fixed_order = OFF
> 11:35:54 3. Shower the generated events  shower = OFF
> 11:35:54 4. Decay onshell particles  madspin = ON
> 11:35:54 5. Add weights to events for new hypp.  reweight = OFF
> 11:35:54 6. Run MadAnalysis5 on the events generated madanalysis = Not Avail.
> 11:35:54 \=============================================================================/
> 11:35:54 Either type the switch number (1 to 6) to change its setting,
> 11:35:54 Set any switch explicitly (e.g. type 'fixed_order=ON' at the prompt)
> 11:35:54 Type 'help' for the list of all valid option
> 11:35:54 Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
> 11:35:54 >INFO: will run in mode: noshower
> 11:35:54 WARNING: You have chosen not to run a parton shower.
> 11:35:54 NLO events without showering are NOT physical.
> 11:35:54 Please, shower the LesHouches events before using them for physics analyses.
> 11:35:54 You have to choose NOW which parton-shower you WILL use and specify it in the run_card. 
> 11:35:54 Do you want to edit a card (press enter to bypass editing)?
> 11:35:54 /------------------------------------------------------------\
> 11:35:54 | 1. param : param_card.dat |
> 11:35:54 | 2. run : run_card.dat |
> 11:35:54 | 3. madspin : madspin_card.dat |
> 11:35:54 \------------------------------------------------------------/
> 11:35:54 you can also
> 11:35:54 - enter the path to a valid card or banner.
> 11:35:54 - use the 'set' command to modify a parameter directly.
> 11:35:54 The set option works only for param_card and run_card.
> 11:35:54 Type 'help set' for more information on this command.
> 11:35:54 - call an external program (ASperGE/MadWidth/...).
> 11:35:54 Type 'help' for the list of available command
> 11:35:54 [0, done, 1, param, 2, run, 3, madspin, enter path][90s to answer]
> 11:35:54 >INFO: Update the dependent parameter of the param_card.dat
> 11:35:55 INFO: Starting run
> 11:35:55 INFO: Compiling the code
> 11:35:55 WARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) 
> 11:35:55 WARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) 
> 11:35:56 INFO: For gauge cancellation, the width of 't' has been set to zero. 
> 11:35:57 INFO: Using LHAPDF v6.1.5 interface for PDFs
> 11:36:02 INFO: Compiling source...
> 11:36:19 INFO: ...done, continuing with P* directories
> 11:36:19 INFO: Compiling StdHEP (can take a couple of minutes) ...
> 11:36:46 INFO: ...done.
> 11:36:46 INFO: Compiling directories...
> 11:36:46 INFO: Compiling on 8 cores
> 11:36:46 INFO: Compiling P0_udx_wp_tbx...
> 11:36:46 INFO: Compiling P0_dxu_wp_tbx...
> 11:36:46 INFO: Compiling P0_dux_wm_txb...
> 11:36:46 INFO: Compiling P0_uxd_wm_txb...
> 11:37:17 INFO: P0_udx_wp_tbx done.
> 11:37:18 INFO: P0_dxu_wp_tbx done.
> 11:37:19 INFO: P0_dux_wm_txb done.
> 11:37:19 INFO: P0_uxd_wm_txb done.
> 11:37:19 INFO: Checking test output:
> 11:37:19 INFO: P0_udx_wp_tbx
> 11:37:19 INFO: Result for test_ME:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for test_MC:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for check_poles:
> 11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> 11:37:19 INFO: P0_dxu_wp_tbx
> 11:37:19 INFO: Result for test_ME:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for test_MC:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for check_poles:
> 11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> 11:37:19 INFO: P0_dux_wm_txb
> 11:37:19 INFO: Result for test_ME:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for test_MC:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for check_poles:
> 11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> 11:37:19 INFO: P0_uxd_wm_txb
> 11:37:19 INFO: Result for test_ME:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for test_MC:
> 11:37:19 INFO: Passed.
> 11:37:19 INFO: Result for check_poles:
> 11:37:19 INFO: Poles successfully cancel for 20 points over 20 (tolerance=1.0e-05)
> 11:37:19 INFO: Starting run
> 11:37:19 INFO: Using 8 cores
> 11:37:19 INFO: Cleaning previous results
> 11:37:19 INFO: Generating events without running the shower.
> 11:37:19 INFO: Setting up grids
> 11:37:20 INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 11h37 ]
> 11:37:30 INFO: Idle: 0, Running: 3, Completed: 1 [ 9.7s ]
> 11:37:30 INFO: Idle: 0, Running: 0, Completed: 4 [ 10.3s ]
> 11:37:30 INFO: Determining the number of unweighted events per channel
> 11:37:30
> 11:37:30 Intermediate results:
> 11:37:30 Random seed: 1
> 11:37:30 Total cross section: 1.019e+01 +- 6.4e-02 pb
> 11:37:30 Total abs(cross section): 1.515e+01 +- 9.1e-02 pb
> 11:37:30
> 11:37:30
> 11:37:30 INFO: Computing upper envelope
> 11:37:31 INFO: Idle: 0, Running: 4, Completed: 0 [ current time: 11h37 ]
> 11:43:36 INFO: Idle: 0, Running: 3, Completed: 1 [ 6m 4s ]
> 11:43:53 INFO: Idle: 0, Running: 0, Completed: 4 [ 6m 21s ]
> 11:43:53 INFO: Updating the number of unweighted events per channel
> 11:43:53
> 11:43:53 Intermediate results:
> 11:43:53 Random seed: 1
> 11:43:53 Total cross section: 1.027e+01 +- 9.9e-03 pb
> 11:43:53 Total abs(cross section): 1.519e+01 +- 1.4e-02 pb
> 11:43:53
> 11:43:53
> 11:43:53 INFO: Generating events
> 11:43:54 INFO: Idle: 0, Running: 0, Completed: 1 [ current time: 11h43 ]
> 11:43:54 INFO: Collecting events
> 11:43:55 INFO:
> 11:43:55 --------------------------------------------------------------
> 11:43:55 Summary:
> 11:43:55 Process p p > w+ > t b~ [QCD] ; p p > w- > t~ b [QCD]
> 11:43:55 Run at p-p collider (6500.0 + 6500.0 GeV)
> 11:43:55 Number of events generated: 1
> 11:43:55 Total cross section: 3.163e+00 +- 6.4e-03 pb
> 11:43:55 --------------------------------------------------------------
> 11:43:55
> 11:43:55 INFO: The /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events.lhe.gz file has been generated.
> 11:43:55
> 11:43:55 INFO: Events generated
> 11:43:55 reweight -from_cards
> 11:43:55 decay_events -from_cards
> 11:43:57 INFO: Running MadSpin
> 11:43:57 INFO: This functionality allows for the decay of resonances
> 11:43:57 INFO: in a .lhe file, keeping track of the spin correlation effets.
> 11:43:57 INFO: BE AWARE OF THE CURRENT LIMITATIONS:
> 11:43:57 INFO: (1) Only a succession of 2 body decay are currently allowed
> 11:43:57 ************************************************************
> 11:43:57 * *
> 11:43:57 * W E L C O M E to M A D S P I N *
> 11:43:57 * *
> 11:43:57 ************************************************************
> 11:43:57 INFO: Extracting the banner ...
> 11:43:57 INFO: process: p p > w+ > t b~
> 11:43:57 INFO: options:
> 11:43:58 INFO: detected model: loop_sm-no_b_mass. Loading...
> 11:43:58 import Events/schan/events.lhe.gz
> 11:43:58 INFO: process: p p > w+ > t b~
> 11:43:58 INFO: options:
> 11:43:58 Set group_subprocesses to Auto
> 11:43:58 Note that you need to regenerate all processes
> 11:43:58 set loop optimized output to True
> 11:43:58 set loop color flows to False
> 11:43:58 Complex mass already desactivated.
> 11:43:58 INFO: detected model: loop_sm-no_b_mass. Loading...
> 11:43:58 INFO: load particles
> 11:43:58 INFO: load vertices
> 11:43:58 INFO: Restrict model loop_sm with file ../../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
> 11:43:58 INFO: Run "set stdout_level DEBUG" before import for more information.
> 11:43:58 INFO: Change particles name to pass to MG5 convention
> 11:43:58 Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
> 11:43:58 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
> 11:43:58 Set group_subprocesses to Auto
> 11:43:58 Note that you need to regenerate all processes
> 11:43:58 set loop optimized output to True
> 11:43:58 set loop color flows to False
> 11:43:58 Complex mass already desactivated.
> 11:43:58 Defined multiparticle p = g u c d s u~ c~ d~ s~
> 11:43:58 Defined multiparticle j = g u c d s u~ c~ d~ s~
> 11:43:58 Defined multiparticle l+ = e+ mu+
> 11:43:58 Defined multiparticle l- = e- mu-
> 11:43:58 Defined multiparticle vl = ve vm vt
> 11:43:58 Defined multiparticle vl~ = ve~ vm~ vt~
> 11:43:58 Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
> 11:43:58 Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
> 11:43:58 Defined multiparticle l+ = e+ mu+ ta+
> 11:43:58 Defined multiparticle l- = e- mu- ta-
> 11:43:58 set ms_dir MadSpin
> 11:43:58 set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
> 11:43:58 set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
> 11:43:58 decay t > w+ b, w+ > l+ vl
> 11:43:58 decay t~ > w- b~, w- > l- vl~
> 11:43:58 launch
> 11:43:58 INFO: Will use seed 766814452
> 11:43:58 INFO:
> 11:43:58 INFO: decay channels for t~ : ( width = 1.32 GeV )
> 11:43:58 INFO: BR d1 d2
> 11:43:58 INFO: 1.000000e+00 w- b~
> 11:43:58 INFO:
> 11:43:58 INFO:
> 11:43:58 INFO: decay channels for w- : ( width = 2.085 GeV )
> 11:43:58 INFO: BR d1 d2
> 11:43:58 INFO: 3.337000e-01 u~ d
> 11:43:58 INFO: 3.337000e-01 c~ s
> 11:43:58 INFO: 1.082000e-01 ve~ e-
> 11:43:58 INFO: 1.082000e-01 vm~ mu-
> 11:43:58 INFO: 1.082000e-01 vt~ ta-
> 11:43:58 INFO:
> 11:43:58 INFO:
> 11:43:58 INFO: decay channels for w+ : ( width = 2.085 GeV )
> 11:43:58 INFO: BR d1 d2
> 11:43:58 INFO: 3.337000e-01 u d~
> 11:43:58 INFO: 3.337000e-01 c s~
> 11:43:58 INFO: 1.082000e-01 ve e+
> 11:43:58 INFO: 1.082000e-01 vm mu+
> 11:43:58 INFO: 1.082000e-01 vt ta+
> 11:43:58 INFO:
> 11:43:58 INFO:
> 11:43:58 INFO: decay channels for t : ( width = 1.32 GeV )
> 11:43:58 INFO: BR d1 d2
> 11:43:58 INFO: 1.000000e+00 w+ b
> 11:43:58 INFO:
> 11:43:59 INFO: generating the production square matrix element
> 11:43:59 INFO: generate p p > w+ > t b~ --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w+ > t b~ pert_QCD --no_warning=duplicate;add process p p > w- > t~ b --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w- > t~ b pert_QCD --no_warning=duplicate;
> 11:44:03 INFO: Done 4.271
> 11:44:03 INFO: generating the full matrix element squared (with decay)
> 11:44:03 INFO: generate p p > w+ > t b~, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w+ > t b~ pert_QCD, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;add process p p > w- > t~ b, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;define pert_QCD = -5 -4 -3 -2 -1 1 2 3 4 5 21;add process p p > w- > t~ b pert_QCD, (t~ > b~ w- , w- > l- vl~ QCD=99), (t > b w+ , w+ > l+ vl QCD=99) --no_warning=duplicate;
> 11:44:09 INFO: Done 6.096
> 11:44:09 INFO: generate matrix element for decay only (1 - > N).
> 11:44:09 INFO: output standalone_msF /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/MadSpin/decay_me
> 11:44:13 INFO: Done 3.253
> 11:44:13 INFO: Compiling code
> 11:44:41 INFO: detect independant decays
> 11:44:41 INFO: Done in 0.252036094666s
> 11:44:41 INFO:
> 11:44:41 INFO: Estimating the maximum weight
> 11:44:41 INFO: *****************************
> 11:44:41 INFO: Probing the first 250 events
> 11:44:41 INFO: with 1000 phase space points
> 11:44:41 INFO: For 2 decaying particle type in the final states
> 11:44:41 INFO:
> 11:44:41 INFO: Event 1/500 : 0.17s
> 11:44:41 INFO: Not enough events for at least one production mode.
> 11:44:41 INFO: This is ok as long as you don't reuse the max weight for other generations.
> 11:44:42 INFO:
> 11:44:42 INFO: Decaying the events...
> 11:44:42 INFO: Total number of events written: 1/1
> 11:44:42 INFO: Average number of trial points per production event: 1.0
> 11:44:42 INFO: Branching ratio to allowed decays: 0.3246
> 11:44:42 INFO: Number of events with weights larger than max_weight: 0
> 11:44:42 INFO: Number of subprocesses 1
> 11:44:42 INFO: Number of failures when restoring the Monte Carlo masses: 0
> 11:44:42 INFO: Decayed events have been written in /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events_decayed.lhe.gz
> 11:44:42 INFO: The decayed event file has been moved to the following location:
> 11:44:42 INFO: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan_decayed_1/events.lhe.gz
> 11:44:42 INFO: MadSpin Done
> 11:44:42 WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical.
> 11:44:42 Please, shower the Les Houches events before using them for physics analyses. 
> 11:44:42 INFO:
> 11:44:42 quit
> 11:44:42 INFO:
> 11:44:42 Py:MadGraphUtils INFO "rw_rscale" = 1.0, 2.0, 0.5
> 11:44:42 Py:MadGraphUtils INFO "ptl" = 0.0
> 11:44:42 Py:MadGraphUtils INFO "xn" = 1.0
> 11:44:42 Py:MadGraphUtils INFO "ptj" = 0.1
> 11:44:42 Py:MadGraphUtils INFO "niters_FO" = 6
> 11:44:42 Py:MadGraphUtils INFO "nevt_job" = -1
> 11:44:42 Py:MadGraphUtils INFO "drll" = 0.0
> 11:44:42 Py:MadGraphUtils INFO "fixed_ren_scale" = True
> 11:44:42 Py:MadGraphUtils INFO "npoints_FO_grid" = 5000
> 11:44:42 Py:MadGraphUtils INFO "iseed" = 1
> 11:44:42 Py:MadGraphUtils INFO "muR_ref_fixed" = 172.5
> 11:44:42 Py:MadGraphUtils INFO "dynamical_scale_choice" = -1
> 11:44:42 Py:MadGraphUtils INFO "R0gamma" = 0.4
> 11:44:42 Py:MadGraphUtils INFO "reweight_scale" = False
> 11:44:42 Py:MadGraphUtils INFO "store_rwgt_info" = False
> 11:44:42 Py:MadGraphUtils INFO "lpp1" = 1
> 11:44:42 Py:MadGraphUtils INFO "nevents" = 1
> 11:44:42 Py:MadGraphUtils INFO "pt_max_pdg" = {}
> 11:44:42 Py:MadGraphUtils INFO "ickkw" = 0
> 11:44:42 Py:MadGraphUtils INFO "iappl" = 0
> 11:44:42 Py:MadGraphUtils INFO "etal" = -1.0
> 11:44:42 Py:MadGraphUtils INFO "epsgamma" = 1.0
> 11:44:42 Py:MadGraphUtils INFO "req_acc_FO" = 0.01
> 11:44:42 Py:MadGraphUtils INFO "pt_min_pdg" = {}
> 11:44:42 Py:MadGraphUtils INFO "etagamma" = -1.0
> 11:44:42 Py:MadGraphUtils INFO "drll_sf" = 0.0
> 11:44:42 Py:MadGraphUtils INFO "muR_over_ref" = 1.0
> 11:44:42 Py:MadGraphUtils INFO "pdlabel" = lhapdf
> 11:44:42 Py:MadGraphUtils INFO "parton_shower" = HERWIGPP
> 11:44:42 Py:MadGraphUtils INFO "event_norm" = average
> 11:44:42 Py:MadGraphUtils INFO "mll_sf" = 30.0
> 11:44:42 Py:MadGraphUtils INFO "muF_ref_fixed" = 172.5
> 11:44:42 Py:MadGraphUtils INFO "bwcutoff" = 50
> 11:44:42 Py:MadGraphUtils INFO "ptgmin" = 20.0
> 11:44:42 Py:MadGraphUtils INFO "rw_fscale" = 1.0, 2.0, 0.5
> 11:44:42 Py:MadGraphUtils INFO "mxx_min_pdg" = {}
> 11:44:42 Py:MadGraphUtils INFO "npoints_FO" = 10000
> 11:44:42 Py:MadGraphUtils INFO "reweight_PDF" = False
> 11:44:42 Py:MadGraphUtils INFO "muF_over_ref" = 1.0
> 11:44:42 Py:MadGraphUtils INFO "isoEM" = True
> 11:44:42 Py:MadGraphUtils INFO "mll" = 0.0
> 11:44:42 Py:MadGraphUtils INFO "fixed_fac_scale" = True
> 11:44:42 Py:MadGraphUtils INFO "lpp2" = 1
> 11:44:42 Py:MadGraphUtils INFO "lhaid" = 260000
> 11:44:42 Py:MadGraphUtils INFO "niters_FO_grid" = 4
> 11:44:42 Py:MadGraphUtils INFO "etaj" = -1.0
> 11:44:42 Py:MadGraphUtils INFO "jetradius" = 0.4
> 11:44:42 Py:MadGraphUtils INFO "jetalgo" = -1
> 11:44:42 Py:MadGraphUtils INFO "req_acc" = 0.001000
> 11:44:42 Py:MadGraphUtils INFO "shower_scale_factor" = 1.0
> 11:44:42 Py:MadGraphUtils INFO "ebeam2" = 6500
> 11:44:42 Py:MadGraphUtils INFO "run_tag" = tag_1
> 11:44:42 Py:MadGraphUtils INFO "ebeam1" = 6500
> 11:44:42 Py:MadGraphUtils INFO Finished checking run card - All OK!
> 11:44:42 Py:MadGraphUtils INFO Modifying run card located at Cards/run_card.dat.
> 11:44:42 Py:MadGraphUtils INFO Keeping backup of original run card at Cards/run_card_backup.dat.
> 11:44:42 {'nevents': 250}
> 11:44:42 Py:MadGraphUtils INFO Setting nevents = 250.
> 11:44:42 Py:MadGraphUtils INFO Finished modification of run card.
> 11:44:45 [?1034hNo module named madgraph
> 11:44:45 INFO: ************************************************************
> 11:44:45 * *
> 11:44:45 * W E L C O M E to M A D G R A P H 5 *
> 11:44:45 * a M C @ N L O *
> 11:44:45 * *
> 11:44:45 * * * *
> 11:44:45 * * * * * *
> 11:44:45 * * * * * 5 * * * * *
> 11:44:45 * * * * * *
> 11:44:45 * * * *
> 11:44:45 * *
> 11:44:45 * VERSION 5.2.6.2 20xx-xx-xx *
> 11:44:45 * *
> 11:44:45 * The MadGraph5_aMC@NLO Development Team - Find us at *
> 11:44:45 * http://amcatnlo.cern.ch *
> 11:44:45 * *
> 11:44:45 * Type 'help' for in-line help. *
> 11:44:45 * *
> 11:44:45 ************************************************************
> 11:44:45 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:44:45 INFO: load configuration from /cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/input/mg5_configuration.txt
> 11:44:45 INFO: load configuration from /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Cards/amcatnlo_configuration.txt
> 11:44:45 Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> 11:44:45 Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> 11:44:45 launch --parton --nocompile --only_generation -f --name=schan
> 11:44:45 WARNING: Removing old run information in
> 11:44:45 /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan 
> 11:44:45 stty: standard input: Invalid argument
> 11:44:45 stty: standard input: Invalid argument
> 11:44:45 INFO: will run in mode: noshower
> 11:44:45 WARNING: You have chosen not to run a parton shower.
> 11:44:45 NLO events without showering are NOT physical.
> 11:44:45 Please, shower the LesHouches events before using them for physics analyses.
> 11:44:45 You have to choose NOW which parton-shower you WILL use and specify it in the run_card. 
> 11:44:45 INFO: Starting run
> 11:44:45 INFO: Compiling the code
> 11:44:45 WARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) 
> 11:44:45 WARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) 
> 11:44:47 INFO: For gauge cancellation, the width of 't' has been set to zero. 
> 11:44:47 INFO: Starting run
> 11:44:47 INFO: Using 8 cores
> 11:44:47 INFO: Generating events without running the shower.
> 11:44:47 INFO: Updating the number of unweighted events per channel
> 11:44:47
> 11:44:47 Intermediate results:
> 11:44:47 Random seed: 1
> 11:44:47 Total cross section: 1.027e+01 +- 9.9e-03 pb
> 11:44:47 Total abs(cross section): 1.519e+01 +- 1.4e-02 pb
> 11:44:47
> 11:44:47
> 11:44:47 INFO: Generating events
> 11:44:48 INFO: Idle: 0, Running: 2, Completed: 2 [ current time: 11h44 ]
> 11:44:48 INFO: Idle: 0, Running: 1, Completed: 3 [ 0.17s ]
> 11:44:48 INFO: Idle: 0, Running: 0, Completed: 4 [ 0.18s ]
> 11:44:48 INFO: Collecting events
> 11:44:48 INFO:
> 11:44:48 --------------------------------------------------------------
> 11:44:48 Summary:
> 11:44:48 Process p p > w+ > t b~ [QCD] ; p p > w- > t~ b [QCD]
> 11:44:48 Run at p-p collider (6500.0 + 6500.0 GeV)
> 11:44:48 Number of events generated: 250
> 11:44:48 Total cross section: 1.027e+01 +- 9.9e-03 pb
> 11:44:48 --------------------------------------------------------------
> 11:44:48
> 11:44:48 INFO: The /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events.lhe.gz file has been generated.
> 11:44:48
> 11:44:48 INFO: Events generated
> 11:44:48 reweight -from_cards
> 11:44:48 decay_events -from_cards
> 11:44:50 INFO: Running MadSpin
> 11:44:50 INFO: This functionality allows for the decay of resonances
> 11:44:50 INFO: in a .lhe file, keeping track of the spin correlation effets.
> 11:44:50 INFO: BE AWARE OF THE CURRENT LIMITATIONS:
> 11:44:50 INFO: (1) Only a succession of 2 body decay are currently allowed
> 11:44:50 ************************************************************
> 11:44:50 * *
> 11:44:50 * W E L C O M E to M A D S P I N *
> 11:44:50 * *
> 11:44:50 ************************************************************
> 11:44:50 INFO: Extracting the banner ...
> 11:44:50 INFO: process: p p > w+ > t b~
> 11:44:50 INFO: options:
> 11:44:50 INFO: detected model: loop_sm-no_b_mass. Loading...
> 11:44:51 import Events/schan/events.lhe.gz
> 11:44:51 INFO: process: p p > w+ > t b~
> 11:44:51 INFO: options:
> 11:44:51 Set group_subprocesses to Auto
> 11:44:51 Note that you need to regenerate all processes
> 11:44:51 set loop optimized output to True
> 11:44:51 set loop color flows to False
> 11:44:51 Complex mass already desactivated.
> 11:44:51 INFO: detected model: loop_sm-no_b_mass. Loading...
> 11:44:51 INFO: load particles
> 11:44:51 INFO: load vertices
> 11:44:51 INFO: Restrict model loop_sm with file ../../../../../../../../../../cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/sw/lcg/external/MCGenerators_lcgcmt67c/madgraph5amc/2.6.2.atlas/x86_64-slc6-gcc47-opt/models/loop_sm/restrict_no_b_mass.dat .
> 11:44:51 INFO: Run "set stdout_level DEBUG" before import for more information.
> 11:44:51 INFO: Change particles name to pass to MG5 convention
> 11:44:51 Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
> 11:44:51 Defined multiparticle all = g gh gh~ d u s c b d~ u~ s~ c~ b~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t t~ z w+ h w- ta- ta+
> 11:44:51 Set group_subprocesses to Auto
> 11:44:51 Note that you need to regenerate all processes
> 11:44:51 set loop optimized output to True
> 11:44:51 set loop color flows to False
> 11:44:51 Complex mass already desactivated.
> 11:44:51 Defined multiparticle p = g u c d s u~ c~ d~ s~
> 11:44:51 Defined multiparticle j = g u c d s u~ c~ d~ s~
> 11:44:51 Defined multiparticle l+ = e+ mu+
> 11:44:51 Defined multiparticle l- = e- mu-
> 11:44:51 Defined multiparticle vl = ve vm vt
> 11:44:51 Defined multiparticle vl~ = ve~ vm~ vt~
> 11:44:51 Defined multiparticle p = g u c d s b u~ c~ d~ s~ b~
> 11:44:51 Defined multiparticle j = g u c d s b u~ c~ d~ s~ b~
> 11:44:51 Defined multiparticle l+ = e+ mu+ ta+
> 11:44:51 Defined multiparticle l- = e- mu- ta-
> 11:44:51 set ms_dir MadSpin
> 11:44:51 set Nevents_for_max_weigth 250 # number of events for the estimate of the max. weight (default: 75)
> 11:44:51 set max_weight_ps_point 1000 # number of PS to estimate the maximum for each event (default: 400)
> 11:44:51 decay t > w+ b, w+ > l+ vl
> 11:44:51 decay t~ > w- b~, w- > l- vl~
> 11:44:51 launch
> 11:44:51 INFO:
> 11:44:51 INFO: Decaying the events...
> 11:44:52 CRITICAL: over_weight: 1.23606071569 ('t_wpb_wp_veep',), occurence: 0.45871559633%, occurence_channel: 1.96078431373%
> 11:44:52 production_tag:P1_udx_wp_tbx [((-1, 2), (-5, 6))], decay:P1_udx_wp_tbx_t_bwp_wp_epve [('t_wpb_wp_veep',)], BW_cut: 23.0977
> 11:44:52
> 11:44:52 [decay.py at line 2297] 
> 11:44:52 INFO: Total number of events written: 250/250
> 11:44:52 INFO: Average number of trial points per production event: 7.576
> 11:44:52 INFO: Branching ratio to allowed decays: 0.3246
> 11:44:52 INFO: Number of events with weights larger than max_weight: 1
> 11:44:52 INFO: Number of subprocesses 18
> 11:44:52 INFO: Number of failures when restoring the Monte Carlo masses: 0
> 11:44:52 INFO: Decayed events have been written in /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan/events_decayed.lhe.gz
> 11:44:52 INFO: The decayed event file has been moved to the following location:
> 11:44:52 INFO: /afs/cern.ch/work/m/mdebeurs/ATHENA/Gridpack_creating/s-channel/b9/PROCNLO_loop_sm-no_b_mass_0/Events/schan_decayed_2/events.lhe.gz
> 11:44:52 INFO: MadSpin Done
> 11:44:52 WARNING: You have chosen not to run a parton shower. NLO events without showering are NOT physical.
> 11:44:52 Please, shower the Les Houches events before using them for physics analyses. 
> 11:44:52 INFO:
> 11:44:52 quit
> 11:44:52 INFO:
> 11:44:52 Py:MadGraphUtils INFO Creating gridpack (schan_gridpack.tar.gz)...
> 11:44:52 Py:MadGraphUtils INFO Package up proc_dir
> 11:45:06 Shortened traceback (most recent user call last):
> 11:45:06 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/jobOptions/EvgenJobTransforms/skeleton.GENtoEVGEN.py", line 225, in <module>
> 11:45:06 include(jo)
> 11:45:06 File "./MC15.412049.aMcAtNloHerwig7EG_H7UE_singletop_schan_lept.py", line 14, in <module>
> 11:45:06 include("./MadGraphControl_schan_NLO_Herwig7_A14.py")
> 11:45:06 File "./MadGraphControl_schan_NLO_Herwig7_A14.py", line 218, in <module>
> 11:45:06 cluster_queue=cluster_queue)
> 11:45:06 File "/cvmfs/atlas.cern.ch/repo/sw/software/x86_64-slc6-gcc47-opt/19.2.5/MCProd/19.2.5.33.4/InstallArea/python/MadGraphControl/MadGraphUtils.py", line 584, in generate
> 11:45:06 raise RuntimeError('Gridpack sucessfully created, exiting the transform. IGNORE ERRORS if running gridpack generation!')
> 11:45:06 RuntimeError: Gridpack sucessfully created, exiting the transform. IGNORE ERRORS if running gridpack generation!
> 11:45:06 Py:Athena INFO leaving with code 8: "an unknown exception occurred"
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

Provide an answer of your own, or ask Marc de Beurs for more information if necessary.

To post a message you must log in.