Gridpack generation - Error detected in generate_events pilotrun
Dear all,
I am trying to generate a gridpack by using MG5 2.6.5, the model SMEFTsim_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set complex_mass_scheme False
define p = g u c d s b u~ c~ d~ s~ b~
define j = p
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model SMEFTsim_
generate p p > e+ e- j j j j QCD=0 NP=1 NP^2==1 SMHLOOP=0
output ZV_e_cHbox_LI
But I got the following error:
Error detected in "generate_events pilotrun"
write debug file [...]/genproduc
If you need help with this issue please contact us on https:/
MadGraph5Error : [Fail 5 times]
A compilation Error occurs when trying to compile [...]/genproduc
The compilation fails with the following output message:
gfortran -o madevent driver.o myamp.o genps.o unwgt.o setcuts.o get_color.o cuts.o cluster.o reweight.o initcluster.o addmothers.o setscales.o idenparts.o dummy_fct.o auto_dsig55.o auto_dsig128.o auto_dsig46.o auto_dsig30.o auto_dsig71.o auto_dsig3.o auto_dsig119.o auto_dsig23.o auto_dsig123.o auto_dsig87.o auto_dsig19.o auto_dsig131.o auto_dsig113.o auto_dsig103.o auto_dsig8.o auto_dsig122.o auto_dsig13.o auto_dsig.o auto_dsig21.o auto_dsig112.o auto_dsig18.o auto_dsig120.o auto_dsig102.o auto_dsig2.o auto_dsig12.o auto_dsig94.o auto_dsig79.o auto_dsig110.o auto_dsig101.o auto_dsig75.o auto_dsig11.o auto_dsig100.o auto_dsig132.o auto_dsig108.o auto_dsig118.o auto_dsig1.o auto_dsig22.o auto_dsig47.o auto_dsig130.o auto_dsig86.o auto_dsig99.o auto_dsig36.o auto_dsig109.o auto_dsig93.o auto_dsig68.o auto_dsig44.o auto_dsig121.o auto_dsig66.o auto_dsig76.o auto_dsig136.o auto_dsig77.o auto_dsig92.o auto_dsig111.o auto_dsig53.o auto_dsig91.o auto_dsig49.o auto_dsig26.o auto_dsig126.o auto_dsig85.o auto_dsig52.o auto_dsig65.o auto_dsig116.o auto_dsig51.o auto_dsig64.o auto_dsig16.o auto_dsig106.o auto_dsig88.o auto_dsig37.o auto_dsig54.o auto_dsig9.o auto_dsig58.o auto_dsig20.o auto_dsig57.o auto_dsig78.o auto_dsig84.o auto_dsig27.o auto_dsig62.o auto_dsig127.o auto_dsig90.o auto_dsig97.o auto_dsig48.o auto_dsig61.o auto_dsig56.o auto_dsig60.o auto_dsig6.o auto_dsig117.o auto_dsig35.o auto_dsig10.o auto_dsig63.o auto_dsig17.o auto_dsig39.o auto_dsig134.o auto_dsig107.o auto_dsig135.o auto_dsig139.o auto_dsig43.o auto_dsig82.o auto_dsig25.o auto_dsig125.o auto_dsig81.o auto_dsig115.o auto_dsig42.o auto_dsig80.o auto_dsig144.o auto_dsig98.o auto_dsig105.o auto_dsig50.o auto_dsig83.o auto_dsig74.o auto_dsig15.o auto_dsig59.o auto_dsig41.o auto_dsig40.o auto_dsig67.o auto_dsig73.o auto_dsig4.o auto_dsig45.o auto_dsig34.o auto_dsig38.o auto_dsig24.o auto_dsig69.o auto_dsig29.o auto_dsig129.o auto_dsig143.o auto_dsig5.o auto_dsig138.o auto_dsig95.o auto_dsig33.o auto_dsig137.o auto_dsig124.o auto_dsig142.o auto_dsig72.o auto_dsig114.o auto_dsig141.o auto_dsig140.o auto_dsig7.o auto_dsig32.o auto_dsig104.o auto_dsig96.o auto_dsig133.o auto_dsig14.o auto_dsig89.o auto_dsig70.o auto_dsig31.o auto_dsig28.o matrix133.o matrix58.o matrix30.o matrix28.o matrix3.o matrix128.o matrix94.o matrix23.o matrix119.o matrix90.o matrix132.o matrix78.o matrix19.o matrix75.o matrix113.o matrix130.o matrix22.o matrix103.o matrix122.o matrix84.o matrix67.o matrix118.o matrix54.o matrix21.o matrix112.o matrix39.o matrix20.o matrix99.o matrix18.o matrix4.o matrix13.o matrix120.o matrix102.o matrix2.o matrix93.o matrix12.o matrix110.o matrix101.o matrix47.o matrix11.o matrix100.o matrix5.o matrix109.o matrix1.o matrix92.o matrix31.o matrix53.o matrix36.o matrix91.o matrix9.o matrix126.o matrix52.o matrix63.o matrix49.o matrix74.o matrix51.o matrix136.o matrix50.o matrix45.o matrix85.o matrix62.o matrix55.o matrix61.o matrix26.o matrix96.o matrix60.o matrix98.o matrix6.o matrix73.o matrix97.o matrix116.o matrix33.o matrix37.o matrix16.o matrix106.o matrix89.o matrix131.o matrix65.o matrix69.o matrix59.o matrix68.o matrix56.o matrix70.o matrix7.o matrix137.o matrix48.o matrix71.o matrix141.o matrix14.o matrix27.o matrix66.o matrix35.o matrix79.o matrix117.o matrix83.o matrix57.o matrix17.o matrix123.o matrix29.o matrix43.o matrix82.o matrix135.o matrix72.o matrix139.o matrix81.o matrix44.o matrix42.o matrix80.o matrix25.o matrix8.o matrix125.o matrix41.o matrix40.o matrix77.o matrix115.o matrix144.o matrix105.o matrix95.o matrix15.o matrix34.o matrix76.o matrix38.o matrix108.o matrix121.o matrix143.o matrix86.o matrix107.o matrix129.o matrix64.o matrix134.o matrix111.o matrix138.o matrix24.o matrix124.o matrix142.o matrix88.o matrix32.o matrix46.o matrix114.o matrix127.o matrix140.o matrix10.o matrix87.o matrix104.o -L../../lib/ -ldhelas -ldsample -lmodel -lgeneric -lpdf -lcernlib -O3 -I/[...
/[
collect2: error: ld returned 1 exit status
make: *** [makefile:44: madevent] Error 1
In case you need it, here is the run_card.dat I used:
#******
# 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) *
#******
tag_1 = run_tag ! name of the run
#******
# Run to generate the grid pack *
#******
True = gridpack !True = setting up the grid pack
#******
# 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. *
#******
1000 = nevents ! Number of unweighted events requested
0 = 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
263000 = 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
#******
# Time of flight information. (-1 means not run)
#******
-1.0 = time_of_flight ! threshold below which info is not written
#******
# Matching - Warning! ickkw > 1 is still beta
#******
0 = ickkw ! 0 no matching, 1 MLM, 2 CKKW matching
1 = highestmult ! for ickkw=2, highest mult group
1 = ktscheme ! for ickkw=1, 1 Durham kT, 2 Pythia pTE
1.0 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
True = pdfwgt ! for ickkw=1, perform pdf reweighting
5 = asrwgtflavor ! highest quark flavor for a_s reweight
True = clusinfo ! include clustering tag in output
3.0 = lhe_version ! Change the way clustering information pass to shower.
#******
#******
#
#******
# Automatic ptj and mjj cuts if xqcut > 0
# (turn off for VBF and single top processes)
#******
True = auto_ptj_mjj ! Automatic setting of ptj and mjj
#******
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/
# (note that etmiss/
#******
False = cut_decays ! Cut decay products
#******
# Number of helicities to sum per event (0 = all helicities)
# 0 gives more stable result, but longer run time (needed for
# long decay chains e.g.).
# Use >=2 if most helicities contribute, e.g. pure QCD.
#******
0 = nhel ! Number of helicities used per event
#******
# Standard Cuts
#******
#
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
10.0 = ptj ! minimum pt for the jets
10.0 = ptb ! minimum pt for the b
0.0 = pta ! minimum pt for the photons
0.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) *
#******
6.5 = etaj ! max rap for the jets
6.5 = etab ! max rap for the b
-1.0 = etaa ! max rap for the photons
-1.0 = 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.4 = drjj ! min distance between jets
0.4 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.4 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.4 = drjl ! min distance between jet and lepton
0.0 = drab ! min distance between gamma and b
0.0 = drbl ! min distance between b and lepton
0.4 = dral ! min distance between gamma and lepton
-1.0 = drjjmax ! max distance between jets
-1.0 = drbbmax ! max distance between b's
-1.0 = drllmax ! max distance between leptons
-1.0 = draamax ! max distance between gammas
-1.0 = drbjmax ! max distance between b and jet
-1.0 = drajmax ! max distance between gamma and jet
-1.0 = drjlmax ! max distance between jet and lepton
-1.0 = drabmax ! max distance between gamma and b
-1.0 = drblmax ! max distance between b and lepton
-1.0 = 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! *
#******
100.0 = mmjj ! min invariant mass of a jet pair
100.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
#******
# KT DURHAM CUT *
#******
-1.0 = ktdurham
0.4 = dparameter
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
5 = maxjetflavor ! Maximum jet pdg code
#******
# Jet measure cuts *
#******
0.0 = xqcut ! minimum kt jet measure between partons
#******
#
#******
# Store info for systematics studies *
# WARNING: If use_syst is T, matched Pythia output is *
# meaningful ONLY if plotted taking matchscale *
# reweighting into account! *
#******
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
# 30 50 = sys_matchscale # variation of merging scale
# PDF sets and number of members (0 or none for all members).
# NNPDF30_
# MSTW2008nlo68cl
Thank you in advance for your help.
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply: