Madgraph Error detected in "generate_events run_01"
I can't generate events. What can I do? This is what I get when trying to launch:
Error detected in "generate_events run_01"
write debug file /Users/
If you need help with this issue please contact us on https:/
str : A compilation Error occurs when trying to compile /Users/
The compilation fails with the following output message:
gfortran -w -fPIC -ffixed-
gfortran -w -fPIC -ffixed-
gfortran -o gensym symmetry.o idenparts.o -L../../lib/ -lmodel -lgeneric -lpdf -lc++ -mmacosx-
AND this is the run_01_
generate_events run_01
Traceback (most recent call last):
File "/Users/
return self.onecmd_
File "/Users/
return func(arg, **opt)
File "/Users/
self.
File "/Users/
original_
File "/Users/
postcmd=False)
File "/Users/
stop = Cmd.onecmd_
File "/Users/
return func(arg, **opt)
File "/Users/
jobs, P_zero_result = ajobcreator.
File "/Users/
self.
File "/Users/
return misc.compile(
File "/Users/
raise MadGraph5Error(
MadGraph5Error: A compilation Error occurs when trying to compile /Users/
The compilation fails with the following output message:
gfortran -w -fPIC -ffixed-
gfortran -w -fPIC -ffixed-
gfortran -o gensym symmetry.o idenparts.o -L../../lib/ -lmodel -lgeneric -lpdf -lc++ -mmacosx-
ld: warning: object file (symmetry.o) was built for newer macOS version (11.1) than being linked (10.8)
ld: warning: object file (idenparts.o) was built for newer macOS version (11.1) than being linked (10.8)
ld: warning: dylib (/usr/local/
ld: warning: dylib (/usr/local/
ld: warning: ignoring file ../../lib/
ld: warning: object file (../../
ld: warning: object file (../../
ld: warning: object file (../../
ld: warning: object file (../../
ld: warning: object file (../../
ld: warning: object file (../../
Undefined symbols for architecture x86_64:
"_setpara_", referenced from:
_MAIN__ in symmetry.o
_setrun_ in libgeneric.
ld: symbol(s) not found for architecture x86_64
collect2: error: ld returned 1 exit status
make: *** [gensym] Error 1
Please try to fix this compilations issue and retry.
Help might be found at https:/
If you think that this is a bug, you can report this at https:/
automatic_
cluster_
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Revision history for this message
|
#1 |
Which model/ process are you trying to generate ?
Is it possible that you are not using a UFO model?
Cheers,
Olivier
Revision history for this message
|
#2 |
The UFO model is ok, everything runs fine in another computer. I tried with Madgraph 2.8.3 instead of 2.9.1.2 and I got this:
Command "generate_events ojpeg" interrupted with error:
OSError : [Errno 13] Permission denied
Please report this bug on https:/
More information is found in '/Users/
Please attach this file to your report.
Revision history for this message
|
#3 |
Permission denied is an issue that you do not have enough disk space (or that you have created too many files, or ...)
This is typically not something that we can fix on our side, but something that you can handle yourself (if running on a cluster, then you might need to check the cluster documentation for where you can write file/...
The question about the model/ process is mainly such that i can try to run the same generation and see if I also have the same compilation error.
Cheers,
Olivier
Revision history for this message
|
#4 |
Hi,
I already managed to run it on a cluster and got this:
Command "generate_events ojpeg" interrupted with error:
AttributeError : 'block' object has no attribute 'on_template'
Please report this bug on https:/
More information is found in '/home/
Please attach this file to your report.
This is the debug file:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.3.2 2018-06-22 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events ojpeg
Traceback (most recent call last):
File "/usr/local/
return self.onecmd_
File "/usr/local/
return func(arg, **opt)
File "/usr/local/
self.
File "/usr/local/
original_
File "/usr/local/
postcmd=False)
File "/usr/local/
stop = Cmd.onecmd_
File "/usr/local/
return func(arg, **opt)
File "/usr/local/
self.
File "/usr/local/
self.
File "/usr/local/
python_
File "/usr/local/
text += this_group.
AttributeError: 'block' object has no attribute 'on_template'
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.3.2 2018-06-22 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set low_mem_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model_v4 ssmzp_nonuniver
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > zp , zp > ta- ta+ QCD=0
output Zprime_
#******
# MadGraph/MadEvent *
# http://
# *
# param_card.dat *
# *
# Les Houches friendly file for the SM parameters of MadGraph *
# Spectrum and decay widths produced by SMCalc *
# *
# Please note the following IMPORTANT issues: *
# *
# 0. REFRAIN from editing this file by hand! Some of the parame- *
# ters are not independent *
# (such as G_Fermi, alpha_em, sin(theta_W),MZ,MW) and serious *
# problems might be encountered (such as violation of unitarity *
# or gauge invariance). Always use a calculator. *
# *
# 1. alpha_S(MZ) has been used in the calculation of the parameters *
# but, for consistency, it will be reset by madgraph to the *
# value expected IF the pdfs for collisions with hadrons are *
# used. This value is KEPT by madgraph when no pdf are used *
# lpp(i)=0 . *
# *
# 2. Values of the charm and bottom kinematic (pole) masses are *
# those used in the matrix elements and phase space UNLESS they *
# are set to ZERO from the start in the model (particles.dat) *
# This happens, for example, when using 5-flavor QCD where *
# charm and bottom are treated as partons in the initial state *
# and a zero mass might be hardwired in the model definition. *
# *
##*****
Block MODSEL # Select Model
0 1 # 0 1 = SM
Block SMINPUTS # Standard Model inputs
1 1.27900000E+02 # alpha_em(MZ)(-1) SM MSbar
2 1.16600000E-05 # G_Fermi
3 1.18000000E-01 # alpha_s(MZ) SM MSbar
4 9.15445065E+01 # Z mass (as input parameter)
Block MGSMPARAM # Standard Model parameters for MadGraph
1 2.31200000E-01 # sin(theta_W)^2
2 8.02673592E+01 # W mass (as input parameter)
Block MGYUKAWA # Yukawa masses m/v=y/sqrt(2)
# PDG YMASS
5 4.20000000E+00 # mbottom for the Yukawa y_b
4 1.25000000E+00 # mcharm for the Yukawa y_c
6 1.74300000E+02 # mtop for the Yukawa y_t
15 1.77700000E+00 # mtau for the Yukawa y_ta
Block MGCKM # CKM elements for MadGraph
1 1 9.75000000E-01 # Vud for Cabibbo matrix
Block MASS # Mass spectrum (kinematic masses)
# PDG Mass
4 1.40000000E+00 # charm pole mass
5 4.20000000E+00 # bottom pole mass
6 1.74300000E+02 # top pole mass
15 1.77700000E+00 # tau mass
23 9.15445065E+01 # Z mass
24 8.02673592E+01 # W mass
25 1.00000000E+02 # H mass
32 350.0e+00 # ZPMASS
BLOCK MGUSER
1 1 # cd11 ,first variable name
2 1 # cd12 ,first variable name
3 1 # cd21 ,first variable name
4 1 # cd22 ,first variable name
5 1 # cu11 ,first variable name
6 1 # cu12 ,first variable name
7 1 # cu21 ,first variable name
8 1 # cu22 ,first variable name
9 1 # cl11 ,first variable name
10 1 # cl12 ,first variable name
11 1 # cl21 ,first variable name
12 1 # cl22 ,first variable name
13 1 # cn11 ,first variable name
14 1 # cn12 ,first variable name
15 1 # cn21 ,first variable name
16 1 # cn22 ,first variable name
17 0 # cwwzp ,first variable name
18 0.00000000e+00 # czpzph ,first variable name
19 0.00000000e+00 # czpzhh ,first variable name
# PDG Width
DECAY 6 1.51013490E+00 # top width
DECAY 23 2.44639985E+00 # Z width
DECAY 24 2.03535570E+00 # W width
DECAY 32 10.3924 # ZPWIDTH
DECAY 25 4.27608700E-03 # H width
# BR NDA ID1 ID2
7.18385415E-02 2 4 -4 # BR( H -> c cbar )
8.03404251E-01 2 5 -5 # BR( H -> b bbar )
0.00000000E+00 2 6 -6 # BR( H -> t tbar )
4.83591503E-02 2 15 -15 # BR( H -> tau- tau+)
2.69892403E-05 2 23 23 # BR( H -> Z Z^(*))
3.67605190E-03 2 24 -24 # BR( H -> W W^(*))
2.43358656E-02 2 21 21 # BR( H -> g g )
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
Zprime_
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
100000 = nevents ! Number of unweighted events requested
17133 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
1 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
6500.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0.0 = polbeam1 ! beam polarization for beam 1
0.0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
'lhapdf' = pdlabel ! PDF set
263400 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
average = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only)
#******
1 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
0.0 = xqcut ! minimum kt jet measure between partons
#******
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15000 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
10.0 = ptj ! minimum pt for the jets
0.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
0.0 = ptheavy ! minimum pt for one heavy final state
-1.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-1.0 = ptlmax ! maximum pt for the charged leptons
-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
0.0 = etalmin ! main rap for the charged leptons
#******
# Minimum and maximum DeltaR distance *
#******
0 = drjj ! min distance between jets
0 = drbb ! min distance between b's
0 = drll ! min distance between leptons
0 = draa ! min distance between gammas
0 = drbj ! min distance between b and jet
0.1 = draj ! min distance between gamma and jet
0 = drjl ! min distance between jet and lepton
0 = drab ! min distance between gamma and b
0 = drbl ! min distance between b and lepton
0.1 = dral ! min distance between gamma and lepton
-1 = drjjmax ! max distance between jets
-1 = drbbmax ! max distance between b's
-1 = drllmax ! max distance between leptons
-1 = draamax ! max distance between gammas
-1 = drbjmax ! max distance between b and jet
-1 = drajmax ! max distance between gamma and jet
-1 = drjlmax ! max distance between jet and lepton
-1 = drabmax ! max distance between gamma and b
-1 = drblmax ! max distance between b and lepton
-1 = dralmax ! maxdistance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#******
200.0 = mmjj ! min invariant mass of a jet pair
0.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
-1.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma pair
-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#******
# Minimum and maximum invariant mass for all letpons *
#******
0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#******
# Inclusive cuts *
#******
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth lepton in pt
-1.0 = ptl1max ! maximum pt for the leading lepton in pt
-1.0 = ptl2max ! maximum pt for the second lepton in pt
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#******
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# WBF cuts *
#******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#******
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
True = use_syst ! Enable systematics studies
#
#******
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#******
#
0.5 1 2 = sys_scalefact # factorization/
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_
# MSTW2008nlo68cl
#
Revision history for this message
|
#5 |
Here my advise would be to upgrade to the latest version of the code.
This seems to be an issue with the writting of the run_card.dat which seems unlikely to be un-noticed for two years and a half.
Cheers,
Olivier
Revision history for this message
|
#6 |
Hi,
I already managed to run it on a cluster and got this:
Command "generate_events ojpeg" interrupted with error:
AttributeError : 'block' object has no attribute 'on_template'
Please report this bug on https:/
More information is found in '/home/
Please attach this file to your report.
This is the debug file:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.3.2 2018-06-22 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events ojpeg
Traceback (most recent call last):
File "/usr/local/
return self.onecmd_
File "/usr/local/
return func(arg, **opt)
File "/usr/local/
self.
File "/usr/local/
original_
File "/usr/local/
postcmd=False)
File "/usr/local/
stop = Cmd.onecmd_
File "/usr/local/
return func(arg, **opt)
File "/usr/local/
self.
File "/usr/local/
self.
File "/usr/local/
python_
File "/usr/local/
text += this_group.
AttributeError: 'block' object has no attribute 'on_template'
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.3.2 2018-06-22 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set low_mem_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model_v4 ssmzp_nonuniver
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > zp , zp > ta- ta+ QCD=0
output Zprime_
#******
# MadGraph/MadEvent *
# http://
# *
# param_card.dat *
# *
# Les Houches friendly file for the SM parameters of MadGraph *
# Spectrum and decay widths produced by SMCalc *
# *
# Please note the following IMPORTANT issues: *
# *
# 0. REFRAIN from editing this file by hand! Some of the parame- *
# ters are not independent *
# (such as G_Fermi, alpha_em, sin(theta_W),MZ,MW) and serious *
# problems might be encountered (such as violation of unitarity *
# or gauge invariance). Always use a calculator. *
# *
# 1. alpha_S(MZ) has been used in the calculation of the parameters *
# but, for consistency, it will be reset by madgraph to the *
# value expected IF the pdfs for collisions with hadrons are *
# used. This value is KEPT by madgraph when no pdf are used *
# lpp(i)=0 . *
# *
# 2. Values of the charm and bottom kinematic (pole) masses are *
# those used in the matrix elements and phase space UNLESS they *
# are set to ZERO from the start in the model (particles.dat) *
# This happens, for example, when using 5-flavor QCD where *
# charm and bottom are treated as partons in the initial state *
# and a zero mass might be hardwired in the model definition. *
# *
##*****
Block MODSEL # Select Model
0 1 # 0 1 = SM
Block SMINPUTS # Standard Model inputs
1 1.27900000E+02 # alpha_em(MZ)(-1) SM MSbar
2 1.16600000E-05 # G_Fermi
3 1.18000000E-01 # alpha_s(MZ) SM MSbar
4 9.15445065E+01 # Z mass (as input parameter)
Block MGSMPARAM # Standard Model parameters for MadGraph
1 2.31200000E-01 # sin(theta_W)^2
2 8.02673592E+01 # W mass (as input parameter)
Block MGYUKAWA # Yukawa masses m/v=y/sqrt(2)
# PDG YMASS
5 4.20000000E+00 # mbottom for the Yukawa y_b
4 1.25000000E+00 # mcharm for the Yukawa y_c
6 1.74300000E+02 # mtop for the Yukawa y_t
15 1.77700000E+00 # mtau for the Yukawa y_ta
Block MGCKM # CKM elements for MadGraph
1 1 9.75000000E-01 # Vud for Cabibbo matrix
Block MASS # Mass spectrum (kinematic masses)
# PDG Mass
4 1.40000000E+00 # charm pole mass
5 4.20000000E+00 # bottom pole mass
6 1.74300000E+02 # top pole mass
15 1.77700000E+00 # tau mass
23 9.15445065E+01 # Z mass
24 8.02673592E+01 # W mass
25 1.00000000E+02 # H mass
32 350.0e+00 # ZPMASS
BLOCK MGUSER
1 1 # cd11 ,first variable name
2 1 # cd12 ,first variable name
3 1 # cd21 ,first variable name
4 1 # cd22 ,first variable name
5 1 # cu11 ,first variable name
6 1 # cu12 ,first variable name
7 1 # cu21 ,first variable name
8 1 # cu22 ,first variable name
9 1 # cl11 ,first variable name
10 1 # cl12 ,first variable name
11 1 # cl21 ,first variable name
12 1 # cl22 ,first variable name
13 1 # cn11 ,first variable name
14 1 # cn12 ,first variable name
15 1 # cn21 ,first variable name
16 1 # cn22 ,first variable name
17 0 # cwwzp ,first variable name
18 0.00000000e+00 # czpzph ,first variable name
19 0.00000000e+00 # czpzhh ,first variable name
# PDG Width
DECAY 6 1.51013490E+00 # top width
DECAY 23 2.44639985E+00 # Z width
DECAY 24 2.03535570E+00 # W width
DECAY 32 10.3924 # ZPWIDTH
DECAY 25 4.27608700E-03 # H width
# BR NDA ID1 ID2
7.18385415E-02 2 4 -4 # BR( H -> c cbar )
8.03404251E-01 2 5 -5 # BR( H -> b bbar )
0.00000000E+00 2 6 -6 # BR( H -> t tbar )
4.83591503E-02 2 15 -15 # BR( H -> tau- tau+)
2.69892403E-05 2 23 23 # BR( H -> Z Z^(*))
3.67605190E-03 2 24 -24 # BR( H -> W W^(*))
2.43358656E-02 2 21 21 # BR( H -> g g )
#******
# MadGraph5_aMC@NLO *
# *
# run_card.dat MadEvent *
# *
# This file is used to set the parameters of the run. *
# *
# Some notation/
# *
# Lines starting with a '# ' are info or comments *
# *
# mind the format: value = variable ! comment *
#******
#
#******
# Running parameters
#******
#
#******
# Tag name for the run (one word) *
#******
Zprime_
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
# If you want to run Pythia, avoid more than 50k events in a run. *
#******
100000 = nevents ! Number of unweighted events requested
17133 = iseed ! rnd seed (0=assigned automatically=
#******
# Collider type and energy *
# lpp: 0=No PDF, 1=proton, -1=antiproton, 2=photon from proton, *
# 3=photon from electron *
#******
1 = lpp1 ! beam 1 type
1 = lpp2 ! beam 2 type
6500.0 = ebeam1 ! beam 1 total energy in GeV
6500.0 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0.0 = polbeam1 ! beam polarization for beam 1
0.0 = polbeam2 ! beam polarization for beam 2
#******
# PDF CHOICE: this automatically fixes also alpha_s and its evol. *
#******
'lhapdf' = pdlabel ! PDF set
263400 = lhaid ! if pdlabel=lhapdf, this is the lhapdf number
#******
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
3.0 = lhe_version ! Change the way clustering information pass to shower.
True = clusinfo ! include clustering tag in output
average = event_norm ! average/sum. Normalization of the weight in the LHEF
#******
# Matching parameter (MLM only)
#******
1 = ickkw ! 0 no matching, 1 MLM
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
4 = asrwgtflavor ! highest quark flavor for a_s reweight
True = auto_ptj_mjj ! Automatic setting of ptj and mjj if xqcut >0
0.0 = xqcut ! minimum kt jet measure between partons
#******
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15000 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
10.0 = ptj ! minimum pt for the jets
0.0 = ptb ! minimum pt for the b
10.0 = pta ! minimum pt for the photons
10.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
0.0 = ptheavy ! minimum pt for one heavy final state
-1.0 = ptjmax ! maximum pt for the jets
-1.0 = ptbmax ! maximum pt for the b
-1.0 = ptamax ! maximum pt for the photons
-1.0 = ptlmax ! maximum pt for the charged leptons
-1.0 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0.0 = ej ! minimum E for the jets
0.0 = eb ! minimum E for the b
0.0 = ea ! minimum E for the photons
0.0 = el ! minimum E for the charged leptons
-1.0 = ejmax ! maximum E for the jets
-1.0 = ebmax ! maximum E for the b
-1.0 = eamax ! maximum E for the photons
-1.0 = elmax ! maximum E for the charged leptons
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
2.5 = etaa ! max rap for the photons
2.5 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
0.0 = etalmin ! main rap for the charged leptons
#******
# Minimum and maximum DeltaR distance *
#******
0 = drjj ! min distance between jets
0 = drbb ! min distance between b's
0 = drll ! min distance between leptons
0 = draa ! min distance between gammas
0 = drbj ! min distance between b and jet
0.1 = draj ! min distance between gamma and jet
0 = drjl ! min distance between jet and lepton
0 = drab ! min distance between gamma and b
0 = drbl ! min distance between b and lepton
0.1 = dral ! min distance between gamma and lepton
-1 = drjjmax ! max distance between jets
-1 = drbbmax ! max distance between b's
-1 = drllmax ! max distance between leptons
-1 = draamax ! max distance between gammas
-1 = drbjmax ! max distance between b and jet
-1 = drajmax ! max distance between gamma and jet
-1 = drjlmax ! max distance between jet and lepton
-1 = drabmax ! max distance between gamma and b
-1 = drblmax ! max distance between b and lepton
-1 = dralmax ! maxdistance between gamma and lepton
#******
# Minimum and maximum invariant mass for pairs *
# WARNING: for four lepton final state mmll cut require to have *
# different lepton masses for each flavor! *
#******
200.0 = mmjj ! min invariant mass of a jet pair
0.0 = mmbb ! min invariant mass of a b pair
0.0 = mmaa ! min invariant mass of gamma gamma pair
0.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1.0 = mmjjmax ! max invariant mass of a jet pair
-1.0 = mmbbmax ! max invariant mass of a b pair
-1.0 = mmaamax ! max invariant mass of gamma gamma pair
-1.0 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#******
# Minimum and maximum invariant mass for all letpons *
#******
0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
0.0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1.0 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#******
# Inclusive cuts *
#******
0.0 = xptj ! minimum pt for at least one jet
0.0 = xptb ! minimum pt for at least one b
0.0 = xpta ! minimum pt for at least one photon
0.0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0.0 = ptj1min ! minimum pt for the leading jet in pt
0.0 = ptj2min ! minimum pt for the second jet in pt
0.0 = ptj3min ! minimum pt for the third jet in pt
0.0 = ptj4min ! minimum pt for the fourth jet in pt
-1.0 = ptj1max ! maximum pt for the leading jet in pt
-1.0 = ptj2max ! maximum pt for the second jet in pt
-1.0 = ptj3max ! maximum pt for the third jet in pt
-1.0 = ptj4max ! maximum pt for the fourth jet in pt
0 = cutuse ! reject event if fails any (0) / all (1) jet pt cuts
#******
# Control the pt's of leptons sorted by pt *
#******
0.0 = ptl1min ! minimum pt for the leading lepton in pt
0.0 = ptl2min ! minimum pt for the second lepton in pt
0.0 = ptl3min ! minimum pt for the third lepton in pt
0.0 = ptl4min ! minimum pt for the fourth lepton in pt
-1.0 = ptl1max ! maximum pt for the leading lepton in pt
-1.0 = ptl2max ! maximum pt for the second lepton in pt
-1.0 = ptl3max ! maximum pt for the third lepton in pt
-1.0 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0.0 = htjmin ! minimum jet HT=Sum(jet pt)
-1.0 = htjmax ! maximum jet HT=Sum(jet pt)
0.0 = ihtmin !inclusive Ht for all partons (including b)
-1.0 = ihtmax !inclusive Ht for all partons (including b)
0.0 = ht2min ! minimum Ht for the two leading jets
0.0 = ht3min ! minimum Ht for the three leading jets
0.0 = ht4min ! minimum Ht for the four leading jets
-1.0 = ht2max ! maximum Ht for the two leading jets
-1.0 = ht3max ! maximum Ht for the three leading jets
-1.0 = ht4max ! maximum Ht for the four leading jets
#******
# Photon-isolation cuts, according to hep-ph/9801442 *
# When ptgmin=0, all the other parameters are ignored *
# When ptgmin>0, pta and draj are not going to be used *
#******
0.0 = ptgmin ! Min photon transverse momentum
0.4 = R0gamma ! Radius of isolation code
1.0 = xn ! n parameter of eq.(3.4) in hep-ph/9801442
1.0 = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
True = isoEM ! isolate photons from EM energy (photons and leptons)
#******
# WBF cuts *
#******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
0.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# Turn on either the ktdurham or ptlund cut to activate *
# CKKW(L) merging with Pythia8 [arXiv:1410.3012, arXiv:1109.4829] *
#******
-1.0 = ktdurham
0.4 = dparameter
-1.0 = ptlund
1, 2, 3, 4, 5, 6, 21 = pdgs_for_
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
True = use_syst ! Enable systematics studies
#
#******
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#******
#
0.5 1 2 = sys_scalefact # factorization/
None = sys_alpsfact # \alpha_s emission scale factors
auto = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
NNPDF23_
# MSTW2008nlo68cl
#
Revision history for this message
|
#7 |
Could you try with 2.9.2?
Thanks,
Olivier
Can you help with this problem?
Provide an answer of your own, or ask Camila Paris for more information if necessary.