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 [1;34mPass the definition of 'j' and 'p' to 5 flavour scheme.[0m 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 [1;34mWARNING: 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~ [0m 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 [1;34mWARNING: 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~ [0m 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 [1m order[0m = [32mNLO[0m 11:35:54 2. No MC@[N]LO matching / event generation [1mfixed_order[0m = [32mOFF[0m 11:35:54 3. Shower the generated events [1m shower[0m = [32mOFF[0m 11:35:54 4. Decay onshell particles [1m madspin[0m = [32mON[0m 11:35:54 5. Add weights to events for new hypp. [1m reweight[0m = [31mOFF[0m 11:35:54 6. Run MadAnalysis5 on the events generated [1mmadanalysis[0m = [01mNot Avail.[0m 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 [1;34mWARNING: 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. [0m 11:35:54 Do you want to edit a card (press enter to bypass editing)? 11:35:54 /------------------------------------------------------------\ 11:35:54 | [31m1[0m. param : [32mparam_card.dat[0m | 11:35:54 | [31m2[0m. run : [32mrun_card.dat[0m | 11:35:54 | [31m3[0m. madspin : [32mmadspin_card.dat[0m | 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 [[4m0[0m, 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 [1;34mWARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) [0m 11:35:55 [1;34mWARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) [0m 11:35:56 [1mINFO: For gauge cancellation, the width of 't' has been set to zero. [0m 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 [1;34mWARNING: 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. [0m 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 [1;34mWARNING: 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 [0m 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 [1;34mWARNING: 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. [0m 11:44:45 INFO: Starting run 11:44:45 INFO: Compiling the code 11:44:45 [1;34mWARNING: Missing mass in the lhef file (-6) . Please fix this (use the "update missing" command if needed) [0m 11:44:45 [1;34mWARNING: Missing mass in the lhef file (-24) . Please fix this (use the "update missing" command if needed) [0m 11:44:47 [1mINFO: For gauge cancellation, the width of 't' has been set to zero. [0m 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 [1;31mCRITICAL: 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 [1;30m[decay.py at line 2297][0m [0m 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 [1;34mWARNING: 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. [0m 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 11:45:06 include(jo) 11:45:06 File "./MC15.412049.aMcAtNloHerwig7EG_H7UE_singletop_schan_lept.py", line 14, in 11:45:06 include("./MadGraphControl_schan_NLO_Herwig7_A14.py") 11:45:06 File "./MadGraphControl_schan_NLO_Herwig7_A14.py", line 218, in 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"

MadGraph5_aMC@NLO

History of question #672577