pilotrun error generating gridpacks
Hi
I am trying to create a gridpack for a process that works fine without the gridpack. My process is this:
import model RPVMSSM_UFO -modelname
define j = g u c d s u~ c~ d~ s~ b b~
define excl = su1 su2 su3 su4 su5 su6 sd1 sd2 sd3 sd4 sd5 sd6 su1~ su2~ su3~ su4~ su5~ su6~ sd1~ sd2~ sd3~ sd4~ sd5~ sd6~ x1+ x2+ x1- x2-
generate p p > su1 su1~ /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @0
add process p p > su1 su1~ j /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @1
add process p p > su1 su1~ j j /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @2
output RPVSt100_
and the error is:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.2.2 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events pilotrun
Traceback (most recent call last):
File "/afs/cern.
return self.onecmd_
File "/afs/cern.
return func(arg, **opt)
File "/afs/cern.
postcmd=False)
File "/afs/cern.
stop = Cmd.onecmd_
File "/afs/cern.
return func(arg, **opt)
File "/afs/cern.
self.
File "/afs/cern.
self.
File "/afs/cern.
return f(self, *args, **opt)
File "/afs/cern.
idle, run, finish, fail = self.control(
File "/afs/cern.
raise error.__class__, '[Fail %i times] \n %s ' % (i+1, error)
ClusterManagmen
Job exhausted all retries.
automatic_
cluster_
any idea how to fix it?
thank you for your help.
Alejandro
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- Olivier Mattelaer Edit question
- Last query:
- Last reply:
- Whiteboard:
- This bug is not related to gridpack running. Even without gridpack there are plenty of errors like these: Error: Color combination 1 2 -> 0 not implemented in addmothers.f Confirmed in version 2.2.2 (the version used by the original reporter of the problem) as well as in the most recent version (2.3.0beta). Best, Rikkert
|
Revision history for this message
|
#1 |
Dear Alejandro,
I'll try to reproduce your bug.
Can you please uploads your input parameters here? That is, param_card.dat, run_card.dat and other non-standard settings. The model you are using is this one: http://
Best,
Rikkert
|
Revision history for this message
|
#2 |
Hi,
yes that is the model.
This is the run card, and the param_card:
#******
# 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 *
#******
.false. = 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. *
#******
1500 = 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 = ebeam1 ! beam 1 total energy in GeV
6500 = ebeam2 ! beam 2 total energy in GeV
#******
# Beam polarization from -100 (left-handed) to 100 (right-handed) *
#******
0 = polbeam1 ! beam polarization for beam 1
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 *
#******
F = fixed_ren_scale ! if .true. use fixed ren scale
F = fixed_fac_scale ! if .true. use fixed fac scale
91.1880 = scale ! fixed ren scale
91.1880 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.1880 = dsqrt_q2fact2 ! fixed fact scale for pdf2
1 = scalefact ! scale factor for event-by-event scales
#******
# Matching - Warning! ickkw > 1 is still beta
#******
1 = 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 = alpsfact ! scale factor for QCD emission vx
F = chcluster ! cluster only according to channel diag
F = pdfwgt ! for ickkw=1, perform pdf reweighting
5 = asrwgtflavor ! highest quark flavor for a_s reweight
T = 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)
#******
T = auto_ptj_mjj ! Automatic setting of ptj and mjj
#******
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15 = bwcutoff ! (M+/-bwcutoff*
#******
# Apply pt/E/eta/dr/mij cuts on decay products or not
# (note that etmiss/
#******
F = 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) *
#******
0 = ptj ! minimum pt for the jets
0 = ptb ! minimum pt for the b
10 = pta ! minimum pt for the photons
0 = ptl ! minimum pt for the charged leptons
0 = misset ! minimum missing Et (sum of neutrino's momenta)
0 = ptheavy ! minimum pt for one heavy final state
0 = ptonium ! minimum pt for the quarkonium states
-1 = ptjmax ! maximum pt for the jets
-1 = ptbmax ! maximum pt for the b
-1 = ptamax ! maximum pt for the photons
-1 = ptlmax ! maximum pt for the charged leptons
-1 = missetmax ! maximum missing Et (sum of neutrino's momenta)
#******
# Minimum and maximum E's (in the center of mass frame) *
#******
0 = ej ! minimum E for the jets
0 = eb ! minimum E for the b
0 = ea ! minimum E for the photons
0 = el ! minimum E for the charged leptons
-1 = ejmax ! maximum E for the jets
-1 = ebmax ! maximum E for the b
-1 = eamax ! maximum E for the photons
-1 = elmax ! maximum E for the charged leptons
#******
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
5 = etaj ! max rap for the jets
-1 = etab ! max rap for the b
-1 = etaa ! max rap for the photons
-1 = etal ! max rap for the charged leptons
-1 = etaonium ! max rap for the quarkonium states
0 = etajmin ! min rap for the jets
0 = etabmin ! min rap for the b
0 = etaamin ! min rap for the photons
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! *
#******
0 = mmjj ! min invariant mass of a jet pair
0 = mmbb ! min invariant mass of a b pair
0 = mmaa ! min invariant mass of gamma gamma pair
50 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
-1 = mmjjmax ! max invariant mass of a jet pair
-1 = mmbbmax ! max invariant mass of a b pair
-1 = mmaamax ! max invariant mass of gamma gamma pair
-1 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#******
# Minimum and maximum invariant mass for all letpons *
#******
0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#******
# Minimum and maximum pt for 4-momenta sum of leptons *
#******
0 = ptllmin ! Minimum pt for 4-momenta sum of leptons(l and vl)
-1 = ptllmax ! Maximum pt for 4-momenta sum of leptons(l and vl)
#******
# Inclusive cuts *
#******
0 = xptj ! minimum pt for at least one jet
0 = xptb ! minimum pt for at least one b
0 = xpta ! minimum pt for at least one photon
0 = xptl ! minimum pt for at least one charged lepton
#******
# Control the pt's of the jets sorted by pt *
#******
0 = ptj1min ! minimum pt for the leading jet in pt
0 = ptj2min ! minimum pt for the second jet in pt
0 = ptj3min ! minimum pt for the third jet in pt
0 = ptj4min ! minimum pt for the fourth jet in pt
-1 = ptj1max ! maximum pt for the leading jet in pt
-1 = ptj2max ! maximum pt for the second jet in pt
-1 = ptj3max ! maximum pt for the third jet in pt
-1 = 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 = ptl1min ! minimum pt for the leading lepton in pt
0 = ptl2min ! minimum pt for the second lepton in pt
0 = ptl3min ! minimum pt for the third lepton in pt
0 = ptl4min ! minimum pt for the fourth lepton in pt
-1 = ptl1max ! maximum pt for the leading lepton in pt
-1 = ptl2max ! maximum pt for the second lepton in pt
-1 = ptl3max ! maximum pt for the third lepton in pt
-1 = ptl4max ! maximum pt for the fourth lepton in pt
#******
# Control the Ht(k)=Sum of k leading jets *
#******
0 = htjmin ! minimum jet HT=Sum(jet pt)
-1 = htjmax ! maximum jet HT=Sum(jet pt)
0 = ihtmin !inclusive Ht for all partons (including b)
-1 = ihtmax !inclusive Ht for all partons (including b)
0 = ht2min ! minimum Ht for the two leading jets
0 = ht3min ! minimum Ht for the three leading jets
0 = ht4min ! minimum Ht for the four leading jets
-1 = ht2max ! maximum Ht for the two leading jets
-1 = ht3max ! maximum Ht for the three leading jets
-1 = 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 = 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 = xetamin ! minimum rapidity for two jets in the WBF case
0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# KT DURHAM CUT *
#******
-1 = 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 *
#******
20 = 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! *
#******
T = use_syst ! Enable systematics studies
#
#******
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
#######
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
#######
#######
## INFORMATION FOR DSQMIX
#######
Block dsqmix
1 3 9.862151e-01 # RRd1x3
1 6 1.654681e-01 # RRd1x6
2 3 -1.654681e-01 # RRd2x3
2 6 9.862151e-01 # RRd2x6
3 5 1.000000e+00 # RRd3x5
4 4 1.000000e+00 # RRd4x4
5 2 1.000000e+00 # RRd5x2
6 1 1.000000e+00 # RRd6x1
#######
## INFORMATION FOR FRALPHA
#######
Block fralpha
1 -1.075727e-01 # alp
#######
## INFORMATION FOR HMIX
#######
Block hmix
1 4.642435e+02 # RMUH
2 9.749107e+00 # tb
4 3.251216e+05 # MA2
#######
## INFORMATION FOR IMVCKM
#######
Block imvckm
1 3 -3.349042e-03 # ICKM1x3
2 1 -1.377152e-04 # ICKM2x1
2 2 -3.212913e-05 # ICKM2x2
3 1 -3.258550e-03 # ICKM3x1
3 2 -7.602238e-04 # ICKM3x2
#######
## INFORMATION FOR MASS
#######
Block mass
5 4.832388e+00 # MB
6 1.731000e+02 # MT
15 1.776990e+00 # Mta
23 9.118760e+01 # MZ
24 7.982901e+01 # MW
25 1.139668e+02 # MH01
35 5.643097e+02 # MH02
36 5.646897e+02 # MA0
37 5.704586e+02 # MH
1000001 5.378740e+02 # Msd1
1000002 100 # Msu1
1000003 6.105113e+02 # Msd2
1000004 6.133992e+02 # Msu2
1000005 6.155144e+02 # Msd3
1000006 6.157041e+02 # Msu3
1000011 3.037530e+02 # Msl1
1000012 3.336805e+02 # Msn1
1000013 3.168669e+02 # Msl2
1000014 3.370975e+02 # Msn2
1000015 3.169185e+02 # Msl3
1000016 3.371096e+02 # Msn3
1000021 6.161578e+02 # Mgo
1000022 9.948601e+01 # Mneu1
1000023 1.895047e+02 # Mneu2
1000024 1.894680e+02 # Mch1
1000025 -4.695446e+02 # Mneu3
1000035 4.773438e+02 # Mneu4
1000037 4.805668e+02 # Mch2
2000001 6.155223e+02 # Msd4
2000002 6.157276e+02 # Msu4
2000003 6.349333e+02 # Msd5
2000004 6.302198e+02 # Msu5
2000005 6.349335e+02 # Msd6
2000006 6.302357e+02 # Msu6
2000011 3.464667e+02 # Msl4
2000013 3.464780e+02 # Msl5
2000015 3.481667e+02 # Msl6
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
4 0.000000 # c : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
#######
## INFORMATION FOR MSD2
#######
Block msd2
1 1 3.543979e+05 # RmD21x1
2 2 3.543929e+05 # RmD22x2
3 3 3.461351e+05 # RmD23x3
#######
## INFORMATION FOR MSE2
#######
Block mse2
1 1 9.819512e+04 # RmE21x1
2 2 9.817883e+04 # RmE22x2
3 3 9.361293e+04 # RmE23x3
#######
## INFORMATION FOR MSL2
#######
Block msl2
1 1 1.174128e+05 # RmL21x1
2 2 1.174047e+05 # RmL22x2
3 3 1.151411e+05 # RmL23x3
#######
## INFORMATION FOR MSOFT
#######
Block msoft
1 1.021024e+02 # RMx1
2 1.926534e+02 # RMx2
3 5.884243e+02 # RMx3
21 1.021411e+05 # mHu2
22 -2.154519e+05 # mHd2
#######
## INFORMATION FOR MSQ2
#######
Block msq2
1 1 3.752466e+05 # RmQ21x1
2 2 3.752418e+05 # RmQ22x2
3 3 2.699498e+05 # RmQ23x3
#######
## INFORMATION FOR MSU2
#######
Block msu2
1 1 3.564881e+05 # RmU21x1
2 2 3.564835e+05 # RmU22x2
3 3 1.494114e+05 # RmU23x3
#######
## INFORMATION FOR NMIX
#######
Block nmix
1 1 -9.927936e-01 # RNN1x1
1 2 3.701124e-02 # RNN1x2
1 3 -1.087277e-01 # RNN1x3
1 4 3.419427e-02 # RNN1x4
2 1 -6.145626e-02 # RNN2x1
2 2 -9.731574e-01 # RNN2x2
2 3 1.992759e-01 # RNN2x3
2 4 -9.734984e-02 # RNN2x4
3 1 4.928328e-02 # RNN3x1
3 2 -7.521552e-02 # RNN3x2
3 3 -6.987909e-01 # RNN3x3
3 4 -7.096514e-01 # RNN3x4
4 1 9.030553e-02 # RNN4x1
4 2 -2.143305e-01 # RNN4x2
4 3 -6.783499e-01 # RNN4x3
4 4 6.969568e-01 # RNN4x4
#######
## INFORMATION FOR RVLAMLLE
#######
Block rvlamlle
1 2 1 2.000000e-01 # RLLE1x2x1
1 2 2 2.000000e-01 # RLLE1x2x2
1 2 3 2.000000e-01 # RLLE1x2x3
1 3 1 2.000000e-01 # RLLE1x3x1
1 3 2 2.000000e-01 # RLLE1x3x2
1 3 3 2.000000e-01 # RLLE1x3x3
2 1 1 -2.000000e-01 # RLLE2x1x1
2 1 2 -2.000000e-01 # RLLE2x1x2
2 1 3 -2.000000e-01 # RLLE2x1x3
2 3 1 2.000000e-01 # RLLE2x3x1
2 3 2 2.000000e-01 # RLLE2x3x2
2 3 3 2.000000e-01 # RLLE2x3x3
3 1 1 -2.000000e-01 # RLLE3x1x1
3 1 2 -2.000000e-01 # RLLE3x1x2
3 1 3 -2.000000e-01 # RLLE3x1x3
3 2 1 -2.000000e-01 # RLLE3x2x1
3 2 2 -2.000000e-01 # RLLE3x2x2
3 2 3 -2.000000e-01 # RLLE3x2x3
#######
## INFORMATION FOR RVLAMLQD
#######
Block rvlamlqd
1 1 1 1.000000e-01 # RLQD1x1x1
1 1 2 1.000000e-01 # RLQD1x1x2
1 1 3 1.000000e-01 # RLQD1x1x3
1 2 1 1.000000e-01 # RLQD1x2x1
1 2 2 1.000000e-01 # RLQD1x2x2
1 2 3 1.000000e-01 # RLQD1x2x3
1 3 1 1.000000e-01 # RLQD1x3x1
1 3 2 1.000000e-01 # RLQD1x3x2
1 3 3 1.000000e-01 # RLQD1x3x3
2 1 1 1.000000e-01 # RLQD2x1x1
2 1 2 1.000000e-01 # RLQD2x1x2
2 1 3 1.000000e-01 # RLQD2x1x3
2 2 1 1.000000e-01 # RLQD2x2x1
2 2 2 1.000000e-01 # RLQD2x2x2
2 2 3 1.000000e-01 # RLQD2x2x3
2 3 1 1.000000e-01 # RLQD2x3x1
2 3 2 1.000000e-01 # RLQD2x3x2
2 3 3 1.000000e-01 # RLQD2x3x3
3 1 1 1.000000e-01 # RLQD3x1x1
3 1 2 1.000000e-01 # RLQD3x1x2
3 1 3 1.000000e-01 # RLQD3x1x3
3 2 1 1.000000e-01 # RLQD3x2x1
3 2 2 1.000000e-01 # RLQD3x2x2
3 2 3 1.000000e-01 # RLQD3x2x3
3 3 1 1.000000e-01 # RLQD3x3x1
3 3 2 1.000000e-01 # RLQD3x3x2
3 3 3 1.000000e-01 # RLQD3x3x3
#######
## INFORMATION FOR RVLAMUDD
#######
Block rvlamudd
1 1 2 2.000000e-01 # RLDD1x1x2
1 1 3 2.000000e-01 # RLDD1x1x3
1 2 1 -2.000000e-01 # RLDD1x2x1
1 2 3 2.000000e-01 # RLDD1x2x3
1 3 1 -2.000000e-01 # RLDD1x3x1
1 3 2 -2.000000e-01 # RLDD1x3x2
2 1 2 2.000000e-01 # RLDD2x1x2
2 1 3 2.000000e-01 # RLDD2x1x3
2 2 1 -2.000000e-01 # RLDD2x2x1
2 2 3 2.000000e-01 # RLDD2x2x3
2 3 1 -2.000000e-01 # RLDD2x3x1
2 3 2 -2.000000e-01 # RLDD2x3x2
3 1 2 2.000000e-01 # RLDD3x1x2
3 1 3 2.000000e-01 # RLDD3x1x3
3 2 1 -2.000000e-01 # RLDD3x2x1
3 2 3 2.000000e-01 # RLDD3x2x3
3 3 1 -2.000000e-01 # RLDD3x3x1
3 3 2 -2.000000e-01 # RLDD3x3x2
#######
## INFORMATION FOR RVTLLE
#######
Block rvtlle
1 2 1 2.000000e-01 # RTLE1x2x1
1 2 2 2.000000e-01 # RTLE1x2x2
1 2 3 2.000000e-01 # RTLE1x2x3
1 3 1 2.000000e-01 # RTLE1x3x1
1 3 2 2.000000e-01 # RTLE1x3x2
1 3 3 2.000000e-01 # RTLE1x3x3
2 1 1 -2.000000e-01 # RTLE2x1x1
2 1 2 -2.000000e-01 # RTLE2x1x2
2 1 3 -2.000000e-01 # RTLE2x1x3
2 3 1 2.000000e-01 # RTLE2x3x1
2 3 2 2.000000e-01 # RTLE2x3x2
2 3 3 2.000000e-01 # RTLE2x3x3
3 1 1 -2.000000e-01 # RTLE3x1x1
3 1 2 -2.000000e-01 # RTLE3x1x2
3 1 3 -2.000000e-01 # RTLE3x1x3
3 2 1 -2.000000e-01 # RTLE3x2x1
3 2 2 -2.000000e-01 # RTLE3x2x2
3 2 3 -2.000000e-01 # RTLE3x2x3
#######
## INFORMATION FOR RVTLQD
#######
Block rvtlqd
1 1 1 1.000000e-01 # RTQD1x1x1
1 1 2 1.000000e-01 # RTQD1x1x2
1 1 3 1.000000e-01 # RTQD1x1x3
1 2 1 1.000000e-01 # RTQD1x2x1
1 2 2 1.000000e-01 # RTQD1x2x2
1 2 3 1.000000e-01 # RTQD1x2x3
1 3 1 1.000000e-01 # RTQD1x3x1
1 3 2 1.000000e-01 # RTQD1x3x2
1 3 3 1.000000e-01 # RTQD1x3x3
2 1 1 1.000000e-01 # RTQD2x1x1
2 1 2 1.000000e-01 # RTQD2x1x2
2 1 3 1.000000e-01 # RTQD2x1x3
2 2 1 1.000000e-01 # RTQD2x2x1
2 2 2 1.000000e-01 # RTQD2x2x2
2 2 3 1.000000e-01 # RTQD2x2x3
2 3 1 1.000000e-01 # RTQD2x3x1
2 3 2 1.000000e-01 # RTQD2x3x2
2 3 3 1.000000e-01 # RTQD2x3x3
3 1 1 1.000000e-01 # RTQD3x1x1
3 1 2 1.000000e-01 # RTQD3x1x2
3 1 3 1.000000e-01 # RTQD3x1x3
3 2 1 1.000000e-01 # RTQD3x2x1
3 2 2 1.000000e-01 # RTQD3x2x2
3 2 3 1.000000e-01 # RTQD3x2x3
3 3 1 1.000000e-01 # RTQD3x3x1
3 3 2 1.000000e-01 # RTQD3x3x2
3 3 3 1.000000e-01 # RTQD3x3x3
#######
## INFORMATION FOR RVTUDD
#######
Block rvtudd
1 1 2 2.000000e-01 # RTDD1x1x2
1 1 3 2.000000e-01 # RTDD1x1x3
1 2 1 -2.000000e-01 # RTDD1x2x1
1 2 3 2.000000e-01 # RTDD1x2x3
1 3 1 -2.000000e-01 # RTDD1x3x1
1 3 2 -2.000000e-01 # RTDD1x3x2
2 1 2 2.000000e-01 # RTDD2x1x2
2 1 3 2.000000e-01 # RTDD2x1x3
2 2 1 -2.000000e-01 # RTDD2x2x1
2 2 3 2.000000e-01 # RTDD2x2x3
2 3 1 -2.000000e-01 # RTDD2x3x1
2 3 2 -2.000000e-01 # RTDD2x3x2
3 1 2 2.000000e-01 # RTDD3x1x2
3 1 3 2.000000e-01 # RTDD3x1x3
3 2 1 -2.000000e-01 # RTDD3x2x1
3 2 3 2.000000e-01 # RTDD3x2x3
3 3 1 -2.000000e-01 # RTDD3x3x1
3 3 2 -2.000000e-01 # RTDD3x3x2
#######
## INFORMATION FOR SELMIX
#######
Block selmix
1 3 3.474325e-01 # RRl1x3
1 6 9.377050e-01 # RRl1x6
2 5 1.000000e+00 # RRl2x5
3 4 1.000000e+00 # RRl3x4
4 1 1.000000e+00 # RRl4x1
5 2 1.000000e+00 # RRl5x2
6 3 -9.377050e-01 # RRl6x3
6 6 3.474325e-01 # RRl6x6
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.279340e+02 # aEWM1
3 1.176000e-01 # aS
#######
## INFORMATION FOR SNUMIX
#######
Block snumix
1 3 1.000000e+00 # RRn1x3
2 2 1.000000e+00 # RRn2x2
3 1 1.000000e+00 # RRn3x1
#######
## INFORMATION FOR TD
#######
Block td
1 1 -2.551888e-01 # Rtd1x1
2 2 -4.374642e+00 # Rtd2x2
3 3 -1.677180e+02 # Rtd3x3
#######
## INFORMATION FOR TE
#######
Block te
1 1 -2.155743e-02 # Rte1x1
2 2 -4.457265e+00 # Rte2x2
3 3 -7.436274e+01 # Rte3x3
#######
## INFORMATION FOR TU
#######
Block tu
1 1 -8.900702e-03 # Rtu1x1
2 2 -3.560256e+00 # Rtu2x2
3 3 -5.887617e+02 # Rtu3x3
#######
## INFORMATION FOR UMIX
#######
Block umix
1 1 -9.576948e-01 # RUU1x1
1 2 2.877859e-01 # RUU1x2
2 1 2.877859e-01 # RUU2x1
2 2 9.576948e-01 # RUU2x2
#######
## INFORMATION FOR UPMNS
#######
Block upmns
1 1 1.000000e+00 # RMNS1x1
2 2 1.000000e+00 # RMNS2x2
3 3 1.000000e+00 # RMNS3x3
#######
## INFORMATION FOR USQMIX
#######
Block usqmix
1 3 5.111614e-01 # RRu1x3
1 6 8.594848e-01 # RRu1x6
2 3 -8.594848e-01 # RRu2x3
2 6 5.111614e-01 # RRu2x6
3 5 1.000000e+00 # RRu3x5
4 4 1.000000e+00 # RRu4x4
5 1 1.000000e+00 # RRu5x1
6 2 1.000000e+00 # RRu6x2
#######
## INFORMATION FOR VCKM
#######
Block vckm
1 1 9.738404e-01 # RCKM1x1
1 2 2.271982e-01 # RCKM1x2
1 3 2.173989e-03 # RCKM1x3
2 1 -2.270868e-01 # RCKM2x1
2 2 9.729587e-01 # RCKM2x2
2 3 4.222469e-02 # RCKM2x3
3 1 7.478279e-03 # RCKM3x1
3 2 -4.161426e-02 # RCKM3x2
3 3 9.991002e-01 # RCKM3x3
#######
## INFORMATION FOR VMIX
#######
Block vmix
1 1 -9.899065e-01 # RVV1x1
1 2 1.417221e-01 # RVV1x2
2 1 1.417221e-01 # RVV2x1
2 2 9.899065e-01 # RVV2x2
#######
## INFORMATION FOR YD
#######
Block yd
1 1 1.884671e-04 # Ryd1x1
2 2 3.230863e-03 # Ryd2x2
3 3 1.360783e-01 # Ryd3x3
#######
## INFORMATION FOR YE
#######
Block ye
1 1 2.882278e-05 # Rye1x1
2 2 5.959638e-03 # Rye2x2
3 3 1.002358e-01 # Rye3x3
#######
## INFORMATION FOR YU
#######
Block yu
1 1 8.789391e-06 # Ryu1x1
2 2 3.515756e-03 # Ryu2x2
3 3 8.829349e-01 # Ryu3x3
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.463305e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.141000e+00 # WW
DECAY 25 3.518961e-03 # WH01
DECAY 35 1.389974e+00 # WH02
DECAY 36 1.742413e+00 # WA0
DECAY 37 1.176094e+00 # WH
DECAY 1000001 7.060640e+00 # Wsd1
DECAY 1000002 0.01 # Wsu1
DECAY 1000003 6.043630e-01 # Wsd2
DECAY 1000004 8.637400e+00 # Wsu2
DECAY 1000005 3.208630e-01 # Wsd3
DECAY 1000006 1.283740e+00 # Wsu3
DECAY 1000011 1.333692e+00 # Wsl1
DECAY 1000012 2.205383e+00 # Wsn1
DECAY 1000013 1.315543e+00 # Wsl2
DECAY 1000014 2.263504e+00 # Wsn2
DECAY 1000015 1.315543e+00 # Wsl3
DECAY 1000016 2.263504e+00 # Wsn3
DECAY 1000021 2.524968e+00 # Wgo
DECAY 1000023 4.107365e-05 # Wneu2
DECAY 1000024 1.056560e-03 # Wch1
DECAY 1000025 3.091347e+00 # Wneu3
DECAY 1000035 3.193022e+00 # Wneu4
DECAY 1000037 5.260644e+00 # Wch2
DECAY 2000001 3.208630e-01 # Wsd4
DECAY 2000002 1.283740e+00 # Wsu4
DECAY 2000003 6.304480e+00 # Wsd5
DECAY 2000004 6.376850e+00 # Wsu5
DECAY 2000005 6.304480e+00 # Wsd6
DECAY 2000006 6.376850e+00 # Wsu6
DECAY 2000011 2.332189e+00 # Wsl4
DECAY 2000013 2.332189e+00 # Wsl5
DECAY 2000015 2.187202e+00 # Wsl6
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # tau- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
DECAY 1000022 0.000000 # n1 : 0.0
#======
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#======
Block QNUMBERS 1000022 # n1
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000023 # n2
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000025 # n3
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000035 # n4
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000024 # x1+
1 3 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000037 # x2+
1 3 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000021 # go
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 35 # h02
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 36 # a0
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 37 # h+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000012 # sv1
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000014 # sv2
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000016 # sv3
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000011 # sl1-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000013 # sl2-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000015 # sl3-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000011 # sl4-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000013 # sl5-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000015 # sl6-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000002 # su1
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000004 # su2
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000006 # su3
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000002 # su4
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000004 # su5
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000006 # su6
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000001 # sd1
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000003 # sd2
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000005 # sd3
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000001 # sd4
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000003 # sd5
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000005 # sd6
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
|
Revision history for this message
|
#3 |
Rikkert point me to the fact that this is back to the problem of Question #265473
So looks like you did not solve your problem at all. I will investigate this.
Cheers,
Olivier
|
Revision history for this message
|
#4 |
Hi
actually I think is not. Because in 265473 I was producing the process
import model RPVMSSM_UFO -modelname
define j = g u c d s u~ c~ d~ s~ b b~
generate p p > su1 su1~ /su2 QED=0, su1 > d~ s~, su1~ > d s @0
add process p p > su1 su1~ j /su2 QED=0, su1 > d~ s~, su1~ > d s @1
add process p p > su1 su1~ j j /su2 QED=0, su1 > d~ s~, su1~ > d s @2
output RPVSt100_
using only MG. I couldn't produce at all the lhe file. But now, I reduce the number of diagrams by doing this:
import model RPVMSSM_UFO -modelname
define j = g u c d s u~ c~ d~ s~ b b~
define excl = su1 su2 su3 su4 su5 su6 sd1 sd2 sd3 sd4 sd5 sd6 su1~ su2~ su3~ su4~ su5~ su6~ sd1~ sd2~ sd3~ sd4~ sd5~ sd6~ x1+ x2+ x1- x2-
generate p p > su1 su1~ /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @0
add process p p > su1 su1~ j /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @1
add process p p > su1 su1~ j j /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @2
output RPVSt100_
and now, using only MG and no gridpack included I have my lhe file.
Now I am trying to produce a gridpack using this card and I got that error.
I hope this info helps.
cheers,
|
Revision history for this message
|
#5 |
Hi,
I tried the exact same syntax and clearly the color error occur and should therefore be discarded.
So let focus on that first, and when it will be fix. We can check if this problem still occur or not.
Cheers,
Olivier
|
Revision history for this message
|
#6 |
Great Olivier, thank you for your time and help.
cheers,
Alejandro
|
Revision history for this message
|
#7 |
Hi Olivier
unfortunately the fix that you suggested here: https:/
it is not working for me.
Now, if I download the latest version of MG5_aMC_v2_3_0_beta and run the gridpack generation of my process WITHOUT adding the addmothers.f file , I am getting this warnings:
CRITICAL: Fail to run correctly job 22170.
with option: {'log': None, 'stdout': None, 'argument': ['0', '1', '2'], 'nb_submit': 1, 'stderr': None, 'prog': '/cms/gomez/
file missing: /cms/gomez/
Fails 1 times
No resubmition.
CRITICAL: Fail to run correctly job 22171.
with option: {'log': None, 'stdout': None, 'argument': ['0', '1'], 'nb_submit': 1, 'stderr': None, 'prog': '/cms/gomez/
file missing: /cms/gomez/
Fails 1 times
No resubmition.
INFO: All jobs finished
The job finishes without crushing and I got a tar file at the end. However, when I try to run the ./run script I got this error:
gomez@hexcms:
Updating grid_card.dat...
Now generating 10000 events with random seed 4738294 and granularity -1
P2_gg_su1su1xgg
P2_gg_su1su1xqq
P2_gg_su1su1xbb
P2_gq_su1su1xgq
P2_qq_su1su1xqq
P2_qq_su1su1xgg
P2_qq_su1su1xbb
P1_gg_su1su1xg_
P1_gq_su1su1xq_
P1_qq_su1su1xg_
P0_gg_su1su1x_
P0_qq_su1su1x_
Traceback (most recent call last):
File "./bin/gridrun", line 84, in <module>
cmd_line = cmd_interface.
File "/cms/gomez/
self.
File "/cms/gomez/
self.
File "/cms/gomez/
cwd=Pdir)
File "/cms/gomez/
% arg
Exception: ['./../
Error: event file not found !
On the other hand, if I download the latest MG version and I include the addmothers.f file as you suggested it, the gridpack generation is interrupted with this error:
INFO: Idle: 0, Running: 0, Completed: 312 [ 11m 28s ]
Command "generate_events " interrupted with error:
IOError : [Errno 2] No such file or directory: '/cms/gomez/
Please report this bug on https:/
More information is found in '/cms/gomez/
Please attach this file to your report.
and here is the file:
gomez@hexcms:
#******
#* MadGraph5_
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.2.3.0.beta *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#******
generate_events
Traceback (most recent call last):
File "/cms/gomez/
return self.onecmd_
File "/cms/gomez/
return func(arg, **opt)
File "/cms/gomez/
postcmd=False)
File "/cms/gomez/
stop = Cmd.onecmd_
File "/cms/gomez/
return func(arg, **opt)
File "/cms/gomez/
cross, error = sum_html.
File "/cms/gomez/
Presults = collect_result(cmd, folder_
File "/cms/gomez/
P_comb.
File "/cms/gomez/
oneresult.
File "/cms/gomez/
finput = open(filepath)
IOError: [Errno 2] No such file or directory: '/cms/gomez/
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.3.0.beta 2015-04-10 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set group_subprocesses Auto
set ignore_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model RPVMSSM_UFO -modelname
define j = g u c d s u~ c~ d~ s~ b b~
define excl = su1 su2 su3 su4 su5 su6 sd1 sd2 sd3 sd4 sd5 sd6 su1~ su2\
~ su3~ su4~ su5~ su6~ sd1~ sd2~ sd3~ sd4~ sd5~ sd6~ x1+ x2+ x1- x2-
generate p p > su1 su1~ /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s @0
add process p p > su1 su1~ j /su2 $ excl QED=0, su1 > d~ s~, su1~ > d \
s @1
add process p p > su1 su1~ j j /su2 $ excl QED=0, su1 > d~ s~, su1~ > \
d s @2
output RPVSt100_
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 FOLLOWING UFO MODEL ####
#######
## ##
## Width set on Auto will be computed following the information ##
## present in the decay.py files of the model. ##
## See arXiv:1402.1178 for more details. ##
## ##
#######
#######
## INFORMATION FOR DSQMIX
#######
Block dsqmix
1 3 9.862151e-01 # RRd1x3
1 6 1.654681e-01 # RRd1x6
2 3 -1.654681e-01 # RRd2x3
2 6 9.862151e-01 # RRd2x6
3 5 1.000000e+00 # RRd3x5
4 4 1.000000e+00 # RRd4x4
5 2 1.000000e+00 # RRd5x2
6 1 1.000000e+00 # RRd6x1
#######
## INFORMATION FOR FRALPHA
#######
Block fralpha
1 -1.075727e-01 # alp
#######
## INFORMATION FOR HMIX
#######
Block hmix
1 4.642435e+02 # RMUH
2 9.749107e+00 # tb
4 3.251216e+05 # MA2
#######
## INFORMATION FOR IMVCKM
#######
Block imvckm
1 3 -3.349042e-03 # ICKM1x3
2 1 -1.377152e-04 # ICKM2x1
2 2 -3.212913e-05 # ICKM2x2
3 1 -3.258550e-03 # ICKM3x1
3 2 -7.602238e-04 # ICKM3x2
#######
## INFORMATION FOR MASS
#######
Block mass
5 4.832388e+00 # MB
6 1.731000e+02 # MT
15 1.776990e+00 # Mta
23 9.118760e+01 # MZ
24 7.982901e+01 # MW
25 1.139668e+02 # MH01
35 5.643097e+02 # MH02
36 5.646897e+02 # MA0
37 5.704586e+02 # MH
1000001 5.378740e+02 # Msd1
1000002 100 # Msu1
1000003 6.105113e+02 # Msd2
1000004 6.133992e+02 # Msu2
1000005 6.155144e+02 # Msd3
1000006 6.157041e+02 # Msu3
1000011 3.037530e+02 # Msl1
1000012 3.336805e+02 # Msn1
1000013 3.168669e+02 # Msl2
1000014 3.370975e+02 # Msn2
1000015 3.169185e+02 # Msl3
1000016 3.371096e+02 # Msn3
1000021 6.161578e+02 # Mgo
1000022 9.948601e+01 # Mneu1
1000023 1.895047e+02 # Mneu2
1000024 1.894680e+02 # Mch1
1000025 -4.695446e+02 # Mneu3
1000035 4.773438e+02 # Mneu4
1000037 4.805668e+02 # Mch2
2000001 6.155223e+02 # Msd4
2000002 6.157276e+02 # Msu4
2000003 6.349333e+02 # Msd5
2000004 6.302198e+02 # Msu5
2000005 6.349335e+02 # Msd6
2000006 6.302357e+02 # Msu6
2000011 3.464667e+02 # Msl4
2000013 3.464780e+02 # Msl5
2000015 3.481667e+02 # Msl6
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
1 0.000000 # d : 0.0
2 0.000000 # u : 0.0
3 0.000000 # s : 0.0
4 0.000000 # c : 0.0
11 0.000000 # e- : 0.0
12 0.000000 # ve : 0.0
13 0.000000 # mu- : 0.0
14 0.000000 # vm : 0.0
16 0.000000 # vt : 0.0
21 0.000000 # g : 0.0
22 0.000000 # a : 0.0
#######
## INFORMATION FOR MSD2
#######
Block msd2
1 1 3.543979e+05 # RmD21x1
2 2 3.543929e+05 # RmD22x2
3 3 3.461351e+05 # RmD23x3
#######
## INFORMATION FOR MSE2
#######
Block mse2
1 1 9.819512e+04 # RmE21x1
2 2 9.817883e+04 # RmE22x2
3 3 9.361293e+04 # RmE23x3
#######
## INFORMATION FOR MSL2
#######
Block msl2
1 1 1.174128e+05 # RmL21x1
2 2 1.174047e+05 # RmL22x2
3 3 1.151411e+05 # RmL23x3
#######
## INFORMATION FOR MSOFT
#######
Block msoft
1 1.021024e+02 # RMx1
2 1.926534e+02 # RMx2
3 5.884243e+02 # RMx3
21 1.021411e+05 # mHu2
22 -2.154519e+05 # mHd2
#######
## INFORMATION FOR MSQ2
#######
Block msq2
1 1 3.752466e+05 # RmQ21x1
2 2 3.752418e+05 # RmQ22x2
3 3 2.699498e+05 # RmQ23x3
#######
## INFORMATION FOR MSU2
#######
Block msu2
1 1 3.564881e+05 # RmU21x1
2 2 3.564835e+05 # RmU22x2
3 3 1.494114e+05 # RmU23x3
#######
## INFORMATION FOR NMIX
#######
Block nmix
1 1 -9.927936e-01 # RNN1x1
1 2 3.701124e-02 # RNN1x2
1 3 -1.087277e-01 # RNN1x3
1 4 3.419427e-02 # RNN1x4
2 1 -6.145626e-02 # RNN2x1
2 2 -9.731574e-01 # RNN2x2
2 3 1.992759e-01 # RNN2x3
2 4 -9.734984e-02 # RNN2x4
3 1 4.928328e-02 # RNN3x1
3 2 -7.521552e-02 # RNN3x2
3 3 -6.987909e-01 # RNN3x3
3 4 -7.096514e-01 # RNN3x4
4 1 9.030553e-02 # RNN4x1
4 2 -2.143305e-01 # RNN4x2
4 3 -6.783499e-01 # RNN4x3
4 4 6.969568e-01 # RNN4x4
#######
## INFORMATION FOR RVLAMLLE
#######
Block rvlamlle
1 2 1 2.000000e-01 # RLLE1x2x1
1 2 2 2.000000e-01 # RLLE1x2x2
1 2 3 2.000000e-01 # RLLE1x2x3
1 3 1 2.000000e-01 # RLLE1x3x1
1 3 2 2.000000e-01 # RLLE1x3x2
1 3 3 2.000000e-01 # RLLE1x3x3
2 1 1 -2.000000e-01 # RLLE2x1x1
2 1 2 -2.000000e-01 # RLLE2x1x2
2 1 3 -2.000000e-01 # RLLE2x1x3
2 3 1 2.000000e-01 # RLLE2x3x1
2 3 2 2.000000e-01 # RLLE2x3x2
2 3 3 2.000000e-01 # RLLE2x3x3
3 1 1 -2.000000e-01 # RLLE3x1x1
3 1 2 -2.000000e-01 # RLLE3x1x2
3 1 3 -2.000000e-01 # RLLE3x1x3
3 2 1 -2.000000e-01 # RLLE3x2x1
3 2 2 -2.000000e-01 # RLLE3x2x2
3 2 3 -2.000000e-01 # RLLE3x2x3
#######
## INFORMATION FOR RVLAMLQD
#######
Block rvlamlqd
1 1 1 1.000000e-01 # RLQD1x1x1
1 1 2 1.000000e-01 # RLQD1x1x2
1 1 3 1.000000e-01 # RLQD1x1x3
1 2 1 1.000000e-01 # RLQD1x2x1
1 2 2 1.000000e-01 # RLQD1x2x2
1 2 3 1.000000e-01 # RLQD1x2x3
1 3 1 1.000000e-01 # RLQD1x3x1
1 3 2 1.000000e-01 # RLQD1x3x2
1 3 3 1.000000e-01 # RLQD1x3x3
2 1 1 1.000000e-01 # RLQD2x1x1
2 1 2 1.000000e-01 # RLQD2x1x2
2 1 3 1.000000e-01 # RLQD2x1x3
2 2 1 1.000000e-01 # RLQD2x2x1
2 2 2 1.000000e-01 # RLQD2x2x2
2 2 3 1.000000e-01 # RLQD2x2x3
2 3 1 1.000000e-01 # RLQD2x3x1
2 3 2 1.000000e-01 # RLQD2x3x2
2 3 3 1.000000e-01 # RLQD2x3x3
3 1 1 1.000000e-01 # RLQD3x1x1
3 1 2 1.000000e-01 # RLQD3x1x2
3 1 3 1.000000e-01 # RLQD3x1x3
3 2 1 1.000000e-01 # RLQD3x2x1
3 2 2 1.000000e-01 # RLQD3x2x2
3 2 3 1.000000e-01 # RLQD3x2x3
3 3 1 1.000000e-01 # RLQD3x3x1
3 3 2 1.000000e-01 # RLQD3x3x2
3 3 3 1.000000e-01 # RLQD3x3x3
#######
## INFORMATION FOR RVLAMUDD
#######
Block rvlamudd
1 1 2 2.000000e-01 # RLDD1x1x2
1 1 3 2.000000e-01 # RLDD1x1x3
1 2 1 -2.000000e-01 # RLDD1x2x1
1 2 3 2.000000e-01 # RLDD1x2x3
1 3 1 -2.000000e-01 # RLDD1x3x1
1 3 2 -2.000000e-01 # RLDD1x3x2
2 1 2 2.000000e-01 # RLDD2x1x2
2 1 3 2.000000e-01 # RLDD2x1x3
2 2 1 -2.000000e-01 # RLDD2x2x1
2 2 3 2.000000e-01 # RLDD2x2x3
2 3 1 -2.000000e-01 # RLDD2x3x1
2 3 2 -2.000000e-01 # RLDD2x3x2
3 1 2 2.000000e-01 # RLDD3x1x2
3 1 3 2.000000e-01 # RLDD3x1x3
3 2 1 -2.000000e-01 # RLDD3x2x1
3 2 3 2.000000e-01 # RLDD3x2x3
3 3 1 -2.000000e-01 # RLDD3x3x1
3 3 2 -2.000000e-01 # RLDD3x3x2
#######
## INFORMATION FOR RVTLLE
#######
Block rvtlle
1 2 1 2.000000e-01 # RTLE1x2x1
1 2 2 2.000000e-01 # RTLE1x2x2
1 2 3 2.000000e-01 # RTLE1x2x3
1 3 1 2.000000e-01 # RTLE1x3x1
1 3 2 2.000000e-01 # RTLE1x3x2
1 3 3 2.000000e-01 # RTLE1x3x3
2 1 1 -2.000000e-01 # RTLE2x1x1
2 1 2 -2.000000e-01 # RTLE2x1x2
2 1 3 -2.000000e-01 # RTLE2x1x3
2 3 1 2.000000e-01 # RTLE2x3x1
2 3 2 2.000000e-01 # RTLE2x3x2
2 3 3 2.000000e-01 # RTLE2x3x3
3 1 1 -2.000000e-01 # RTLE3x1x1
3 1 2 -2.000000e-01 # RTLE3x1x2
3 1 3 -2.000000e-01 # RTLE3x1x3
3 2 1 -2.000000e-01 # RTLE3x2x1
3 2 2 -2.000000e-01 # RTLE3x2x2
3 2 3 -2.000000e-01 # RTLE3x2x3
#######
## INFORMATION FOR RVTLQD
#######
Block rvtlqd
1 1 1 1.000000e-01 # RTQD1x1x1
1 1 2 1.000000e-01 # RTQD1x1x2
1 1 3 1.000000e-01 # RTQD1x1x3
1 2 1 1.000000e-01 # RTQD1x2x1
1 2 2 1.000000e-01 # RTQD1x2x2
1 2 3 1.000000e-01 # RTQD1x2x3
1 3 1 1.000000e-01 # RTQD1x3x1
1 3 2 1.000000e-01 # RTQD1x3x2
1 3 3 1.000000e-01 # RTQD1x3x3
2 1 1 1.000000e-01 # RTQD2x1x1
2 1 2 1.000000e-01 # RTQD2x1x2
2 1 3 1.000000e-01 # RTQD2x1x3
2 2 1 1.000000e-01 # RTQD2x2x1
2 2 2 1.000000e-01 # RTQD2x2x2
2 2 3 1.000000e-01 # RTQD2x2x3
2 3 1 1.000000e-01 # RTQD2x3x1
2 3 2 1.000000e-01 # RTQD2x3x2
2 3 3 1.000000e-01 # RTQD2x3x3
3 1 1 1.000000e-01 # RTQD3x1x1
3 1 2 1.000000e-01 # RTQD3x1x2
3 1 3 1.000000e-01 # RTQD3x1x3
3 2 1 1.000000e-01 # RTQD3x2x1
3 2 2 1.000000e-01 # RTQD3x2x2
3 2 3 1.000000e-01 # RTQD3x2x3
3 3 1 1.000000e-01 # RTQD3x3x1
3 3 2 1.000000e-01 # RTQD3x3x2
3 3 3 1.000000e-01 # RTQD3x3x3
#######
## INFORMATION FOR RVTUDD
#######
Block rvtudd
1 1 2 2.000000e-01 # RTDD1x1x2
1 1 3 2.000000e-01 # RTDD1x1x3
1 2 1 -2.000000e-01 # RTDD1x2x1
1 2 3 2.000000e-01 # RTDD1x2x3
1 3 1 -2.000000e-01 # RTDD1x3x1
1 3 2 -2.000000e-01 # RTDD1x3x2
2 1 2 2.000000e-01 # RTDD2x1x2
2 1 3 2.000000e-01 # RTDD2x1x3
2 2 1 -2.000000e-01 # RTDD2x2x1
2 2 3 2.000000e-01 # RTDD2x2x3
2 3 1 -2.000000e-01 # RTDD2x3x1
2 3 2 -2.000000e-01 # RTDD2x3x2
3 1 2 2.000000e-01 # RTDD3x1x2
3 1 3 2.000000e-01 # RTDD3x1x3
3 2 1 -2.000000e-01 # RTDD3x2x1
3 2 3 2.000000e-01 # RTDD3x2x3
3 3 1 -2.000000e-01 # RTDD3x3x1
3 3 2 -2.000000e-01 # RTDD3x3x2
#######
## INFORMATION FOR SELMIX
#######
Block selmix
1 3 3.474325e-01 # RRl1x3
1 6 9.377050e-01 # RRl1x6
2 5 1.000000e+00 # RRl2x5
3 4 1.000000e+00 # RRl3x4
4 1 1.000000e+00 # RRl4x1
5 2 1.000000e+00 # RRl5x2
6 3 -9.377050e-01 # RRl6x3
6 6 3.474325e-01 # RRl6x6
#######
## INFORMATION FOR SMINPUTS
#######
Block sminputs
1 1.279340e+02 # aEWM1
3 1.176000e-01 # aS
#######
## INFORMATION FOR SNUMIX
#######
Block snumix
1 3 1.000000e+00 # RRn1x3
2 2 1.000000e+00 # RRn2x2
3 1 1.000000e+00 # RRn3x1
#######
## INFORMATION FOR TD
#######
Block td
1 1 -2.551888e-01 # Rtd1x1
2 2 -4.374642e+00 # Rtd2x2
3 3 -1.677180e+02 # Rtd3x3
#######
## INFORMATION FOR TE
#######
Block te
1 1 -2.155743e-02 # Rte1x1
2 2 -4.457265e+00 # Rte2x2
3 3 -7.436274e+01 # Rte3x3
#######
## INFORMATION FOR TU
#######
Block tu
1 1 -8.900702e-03 # Rtu1x1
2 2 -3.560256e+00 # Rtu2x2
3 3 -5.887617e+02 # Rtu3x3
#######
## INFORMATION FOR UMIX
#######
Block umix
1 1 -9.576948e-01 # RUU1x1
1 2 2.877859e-01 # RUU1x2
2 1 2.877859e-01 # RUU2x1
2 2 9.576948e-01 # RUU2x2
#######
## INFORMATION FOR UPMNS
#######
Block upmns
1 1 1.000000e+00 # RMNS1x1
2 2 1.000000e+00 # RMNS2x2
3 3 1.000000e+00 # RMNS3x3
#######
## INFORMATION FOR USQMIX
#######
Block usqmix
1 3 5.111614e-01 # RRu1x3
1 6 8.594848e-01 # RRu1x6
2 3 -8.594848e-01 # RRu2x3
2 6 5.111614e-01 # RRu2x6
3 5 1.000000e+00 # RRu3x5
4 4 1.000000e+00 # RRu4x4
5 1 1.000000e+00 # RRu5x1
6 2 1.000000e+00 # RRu6x2
#######
## INFORMATION FOR VCKM
#######
Block vckm
1 1 9.738404e-01 # RCKM1x1
1 2 2.271982e-01 # RCKM1x2
1 3 2.173989e-03 # RCKM1x3
2 1 -2.270868e-01 # RCKM2x1
2 2 9.729587e-01 # RCKM2x2
2 3 4.222469e-02 # RCKM2x3
3 1 7.478279e-03 # RCKM3x1
3 2 -4.161426e-02 # RCKM3x2
3 3 9.991002e-01 # RCKM3x3
#######
## INFORMATION FOR VMIX
#######
Block vmix
1 1 -9.899065e-01 # RVV1x1
1 2 1.417221e-01 # RVV1x2
2 1 1.417221e-01 # RVV2x1
2 2 9.899065e-01 # RVV2x2
#######
## INFORMATION FOR YD
#######
Block yd
1 1 1.884671e-04 # Ryd1x1
2 2 3.230863e-03 # Ryd2x2
3 3 1.360783e-01 # Ryd3x3
#######
## INFORMATION FOR YE
#######
Block ye
1 1 2.882278e-05 # Rye1x1
2 2 5.959638e-03 # Rye2x2
3 3 1.002358e-01 # Rye3x3
#######
## INFORMATION FOR YU
#######
Block yu
1 1 8.789391e-06 # Ryu1x1
2 2 3.515756e-03 # Ryu2x2
3 3 8.829349e-01 # Ryu3x3
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.463305e+00 # WT
DECAY 23 2.495200e+00 # WZ
DECAY 24 2.141000e+00 # WW
DECAY 25 3.518961e-03 # WH01
DECAY 35 1.389974e+00 # WH02
DECAY 36 1.742413e+00 # WA0
DECAY 37 1.176094e+00 # WH
DECAY 1000001 7.060640e+00 # Wsd1
DECAY 1000002 0.001 # Wsu1
DECAY 1000003 6.043630e-01 # Wsd2
DECAY 1000004 8.637400e+00 # Wsu2
DECAY 1000005 3.208630e-01 # Wsd3
DECAY 1000006 1.283740e+00 # Wsu3
DECAY 1000011 1.333692e+00 # Wsl1
DECAY 1000012 2.205383e+00 # Wsn1
DECAY 1000013 1.315543e+00 # Wsl2
DECAY 1000014 2.263504e+00 # Wsn2
DECAY 1000015 1.315543e+00 # Wsl3
DECAY 1000016 2.263504e+00 # Wsn3
DECAY 1000021 2.524968e+00 # Wgo
DECAY 1000023 4.107365e-05 # Wneu2
DECAY 1000024 1.056560e-03 # Wch1
DECAY 1000025 3.091347e+00 # Wneu3
DECAY 1000035 3.193022e+00 # Wneu4
DECAY 1000037 5.260644e+00 # Wch2
DECAY 2000001 3.208630e-01 # Wsd4
DECAY 2000002 1.283740e+00 # Wsu4
DECAY 2000003 6.304480e+00 # Wsd5
DECAY 2000004 6.376850e+00 # Wsu5
DECAY 2000005 6.304480e+00 # Wsd6
DECAY 2000006 6.376850e+00 # Wsu6
DECAY 2000011 2.332189e+00 # Wsl4
DECAY 2000013 2.332189e+00 # Wsl5
DECAY 2000015 2.187202e+00 # Wsl6
## Dependent parameters, given by model restrictions.
## Those values should be edited following the
## analytical expression. MG5 ignores those values
## but they are important for interfacing the output of MG5
## to external program such as Pythia.
DECAY 1 0.000000 # d : 0.0
DECAY 2 0.000000 # u : 0.0
DECAY 3 0.000000 # s : 0.0
DECAY 4 0.000000 # c : 0.0
DECAY 5 0.000000 # b : 0.0
DECAY 11 0.000000 # e- : 0.0
DECAY 12 0.000000 # ve : 0.0
DECAY 13 0.000000 # mu- : 0.0
DECAY 14 0.000000 # vm : 0.0
DECAY 15 0.000000 # tau- : 0.0
DECAY 16 0.000000 # vt : 0.0
DECAY 21 0.000000 # g : 0.0
DECAY 22 0.000000 # a : 0.0
DECAY 1000022 0.000000 # n1 : 0.0
#======
# QUANTUM NUMBERS OF NEW STATE(S) (NON SM PDG CODE)
#======
Block QNUMBERS 1000022 # n1
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000023 # n2
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000025 # n3
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000035 # n4
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 1000024 # x1+
1 3 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000037 # x2+
1 3 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000021 # go
1 0 # 3 times electric charge
2 2 # number of spin states (2S+1)
3 8 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 35 # h02
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 36 # a0
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # Particle/
Block QNUMBERS 37 # h+
1 3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000012 # sv1
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000014 # sv2
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000016 # sv3
1 0 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000011 # sl1-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000013 # sl2-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000015 # sl3-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000011 # sl4-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000013 # sl5-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000015 # sl6-
1 -3 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000002 # su1
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000004 # su2
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000006 # su3
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000002 # su4
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000004 # su5
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000006 # su6
1 2 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000001 # sd1
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000003 # sd2
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 1000005 # sd3
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000001 # sd4
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000003 # sd5
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
Block QNUMBERS 2000005 # sd6
1 -1 # 3 times electric charge
2 1 # number of spin states (2S+1)
3 3 # colour rep (1: singlet, 3: triplet, 8: octet)
4 1 # Particle/
#******
# 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. *
#******
10000 = 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. *
#******
nn23lo1 = pdlabel ! PDF set
230000 = 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
#******
# Matching - Warning! ickkw > 1 is still beta
#******
1 = 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 = alpsfact ! scale factor for QCD emission vx
False = chcluster ! cluster only according to channel diag
True = pdfwgt ! for ickkw=1, perform pdf reweighting
4 = 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/dr/mij cuts on decay products or not
# (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) *
#******
20.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.0 = drjj ! min distance between jets
0.0 = drbb ! min distance between b's
0.4 = drll ! min distance between leptons
0.4 = draa ! min distance between gammas
0.0 = drbj ! min distance between b and jet
0.4 = draj ! min distance between gamma and jet
0.0 = 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! *
#******
0.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
#******
# 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) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
# Jet measure cuts *
#******
30.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! *
#******
False = use_syst ! Enable systematics studies
#
#******
# Parameter of the systematics study
# will be used by SysCalc (if installed)
#******
#
0.5 1 2 = sys_scalefact # factorization/
0.5 1 2 = 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).
Ct10nlo.LHgrid = sys_pdf # matching scales
# MSTW2008nlo68cl
I hope that this info helps.
cheers,
|
Revision history for this message
|
#8 |
Could you try with:
lp:~maddevelopers/mg5amcnlo/2.3.1_fix_addmothers
It works perfectly fine for me in that version.
Cheers,
Olivier
|
Revision history for this message
|
#9 |
Hi Olivier
unfortunately I still have issues with this. First, I dont know if I should open another threat or can I added it here.
I discover that, if I send the jobs to my cluster (and also to CERN cluster) I cannot generate the gridpack, even with your fix.
BUT, if I run it locally, I have the gridpack at the end, without errors in the process.
(Well actually just some warnings like:
INFO: Idle: 67, Running: 8, Completed: 68 [ 14m 41s ]
STOP 1
STOP 1
STOP 1
STOP 1
INFO: Idle: 64, Running: 8, Completed: 71 [ 14m 46s ]
but I think there is no problem with them right?
Anyway, the issue continues even after I got the gridpack. When I try to run the run.sh script I am getting this error:
Updating grid_card.dat...
Now generating 100 events with random seed 789 and granularity -1
P2_gg_su1su1xgg
P2_gg_su1su1xqq
P2_gq_su1su1xgq
P2_qq_su1su1xqq
P2_qq_su1su1xgg
P1_gg_su1su1xg_
P1_gq_su1su1xq_
P1_qq_su1su1xg_
P0_gg_su1su1x_
P0_qq_su1su1x_
Traceback (most recent call last):
File "./bin/gridrun", line 84, in <module>
cmd_line = cmd_interface.
File "/cms/gomez/
self.
File "/cms/gomez/
self.
File "/cms/gomez/
cwd=Pdir)
File "/cms/gomez/
% arg
Exception: ['./../
Error: event file not found !
any suggestion about it?
thanks again for all your time.
|
Revision history for this message
|
#10 |
Hi,
> STOP 1
> STOP 1
> STOP 1
> INFO: Idle: 64, Running: 8, Completed: 71 [ 14m 46s ]
> but I think there is no problem with them right?
No clearly you can not trust those results.
In the latest version, I have try to crash the code properly instead of having it continuing to run.
So Which version of the code are you using?
Are you using a version from launchpad?
I guess not since in the development version that I use, it goes trough.
Cheers,
Olivier
On 09 Jun 2015, at 10:51, alefisico <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> alefisico is still having a problem:
> Hi Olivier
>
> unfortunately I still have issues with this. First, I dont know if I should open another threat or can I added it here.
> I discover that, if I send the jobs to my cluster (and also to CERN cluster) I cannot generate the gridpack, even with your fix.
> BUT, if I run it locally, I have the gridpack at the end, without errors in the process.
> (Well actually just some warnings like:
> INFO: Idle: 67, Running: 8, Completed: 68 [ 14m 41s ]
> STOP 1
> STOP 1
> STOP 1
> STOP 1
> INFO: Idle: 64, Running: 8, Completed: 71 [ 14m 46s ]
> but I think there is no problem with them right?
>
> Anyway, the issue continues even after I got the gridpack. When I try to run the run.sh script I am getting this error:
> Updating grid_card.dat...
> Now generating 100 events with random seed 789 and granularity -1
> P2_gg_su1su1xgg
> P2_gg_su1su1xqq
> P2_gq_su1su1xgq
> P2_qq_su1su1xqq
> P2_qq_su1su1xgg
> P1_gg_su1su1xg_
> P1_gq_su1su1xq_
> P1_qq_su1su1xg_
> P0_gg_su1su1x_
> P0_qq_su1su1x_
> Traceback (most recent call last):
> File "./bin/gridrun", line 84, in <module>
> cmd_line = cmd_interface.
> File "/cms/gomez/
> self.launch(
> File "/cms/gomez/
> self.refine4gri
> File "/cms/gomez/
> cwd=Pdir)
> File "/cms/gomez/
> % arg
> Exception: ['./../
> Error: event file not found !
>
> any suggestion about it?
>
> thanks again for all your time.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#11 |
Hi Olivier
I have tried with several versions of MG (2.2.2, 2.2.3, 2.3.0.beta 2.3.1) including the one from launchpad. Just to make sure I will explain all what I am doing
##### For the version in the launchpad:
I am downloading the version of the launchpad like this:
bzr branch lp:~maddevelopers/mg5amcnlo/2.3.1_fix_addmothers
there I am running the usual cards:
https:/
https:/
https:/
and then MG if I use the cluster to generate the gridpack I have this error:
https:/
and the gridpack is not created. However, if I run it without the cluster, with multiple cores it runs for ~8hours, I got this STOP1 warnings.
https:/
but I have a gridpack at the end.
Then, to fix the issue in the message #9, once I untar the gridpack, I just copied a gen_ximprove.f file from another version of madgraph and put the file inside madevent/Source/, then I just go to the Source directory and
make ../bin/
This solves the problem and I have an lhe file at the end. But also I noticed that once I create an lhe file from this gridpack I am having events that look like this:
11 2 0.9308000E+01 0.2172364E+03 0.7816531E-02 0.1151592E+00
21 -1 0 0 507 506 0.00000000000E+00 0.00000000000E+00 0.19779054729E+04 0.19779054729E+04 0.00000000000E+00 0. -1.
3 -1 0 0 508 0 0.00000000000E+00 0.00000000000E+00 -0.45643327256E+03 0.45643327256E+03 0.00000000000E+00 0. -1.
1000002 2 1 2 508 0 -0.84219622959E+02 0.31204399495E+02 -0.91462372374E+02 0.22651703404E+03 0.18675636994E+03 0. 0.
1000002 2 3 3 508 0 -0.85463213764E+01 0.56497649818E+02 -0.90908042272E+02 0.14672665639E+03 0.99997077464E+02 0. 0.
-1000002 2 1 2 0 506 0.11087088756E+03 0.12118608604E+02 0.19052421139E+04 0.19111217424E+04 0.99997937121E+02 0. 0.
-1 1 4 4 0 502 -0.53954020986E+02 0.38621826254E+02 -0.51063232565E+02 0.83726552322E+02 0.00000000000E+00 0. -1.
-3 1 4 4 0 501 0.45407699609E+02 0.17875823564E+02 -0.39844809707E+02 0.63000104067E+02 0.00000000000E+00 0. -1.
1 1 5 5 505 0 0.58264518820E+02 -0.43511230204E+02 0.84926181828E+03 0.85236941363E+03 0.00000000000E+00 0. 1.
3 1 5 5 504 0 0.52606368736E+02 0.55629838808E+02 0.10559802956E+04 0.10587523288E+04 0.00000000000E+00 0. 1.
21 1 3 3 508 508 -0.75673301582E+02 -0.25293250324E+02 -0.55433010133E+00 0.79790377654E+02 0.00000000000E+00 0. -1.
3 1 1 2 507 0 -0.26651264597E+02 -0.43323008099E+02 -0.29230754121E+03 0.29669996896E+03 0.00000000000E+00 0. -1.
that pythia8 cannot handle them.
#### Run with 2.2.2
Something weird happened here. I notice that I can create an the lhe file with the cluster, using the same cards as above in this version of MG. I got a lot of this messages:
CRITICAL: Fail to run correctly job 36147.
with option: {'log': None, 'stdout': None, 'argument': [], 'nb_submit': 1, 'stderr': None, 'prog': 'ajob1', 'output_files': ['G1', 'G3'], 'time_check': 1434535204.7950001, 'cwd': '/cms/gomez/
file missing: /cms/gomez/
Fails 1 times
No resubmition
but at the end I got an lhe file, with events that by eye look fine and I can run this using pythia8 without an issue.
However, when I try to run the same for a gridpack, it doesn't work again.
Hopefully this helps.
cheers,
Alejandro
|
Revision history for this message
|
#12 |
Hi,
> I have tried with several versions of MG (2.2.2, 2.2.3, 2.3.0.beta
> 2.3.1) including the one from launchpad. Just to make sure I will
> explain all what I am doing
You can really forget all version but this one now:
2.3.1
(2.3.1_
When I look at your log:
> https:/
I see:
https:/
So it seems to me that you are using the 2.3.0 version and not 2.3.1. Is this correct?
> but at the end I got an lhe file, with events that by eye look fine and I can run this using pythia8 without an issue.
> However, when I try to run the same for a gridpack, it doesn’t work again.
You really should not use that lhe file. Since he is clearly wrong.
Looks like the gridpack is more secure, it is really stupid to generate a wrong LHE file.
I have just push a new version of 2.3.1, where the header is different of the current 2.3.0.beta.
so simply by the header we will be able to be sure which version you are using.
You should have the following line in the header of mg5::
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.3.1.beta 2015-06-17 *
* BZR 2.3.1 400 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* and *
* http://
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
*******
The important line are:
* VERSION 2.3.1.beta 2015-06-17 *
* BZR 2.3.1 400 *
If you do not have those lines, it means that you use a wrong version.
Cheers,
Olivier
On 17 Jun 2015, at 11:26, alefisico <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> alefisico is still having a problem:
> Hi Olivier
>
> I have tried with several versions of MG (2.2.2, 2.2.3, 2.3.0.beta
> 2.3.1) including the one from launchpad. Just to make sure I will
> explain all what I am doing
>
> ##### For the version in the launchpad:
>
> I am downloading the version of the launchpad like this:
> bzr branch lp:~maddevelopers/mg5amcnlo/2.3.1_fix_addmothers
>
> there I am running the usual cards:
> https:/
> https:/
> https:/
>
> and then MG if I use the cluster to generate the gridpack I have this error:
> https:/
>
> and the gridpack is not created. However, if I run it without the cluster, with multiple cores it runs for ~8hours, I got this STOP1 warnings.
> https:/
> but I have a gridpack at the end.
> Then, to fix the issue in the message #9, once I untar the gridpack, I just copied a gen_ximprove.f file from another version of madgraph and put the file inside madevent/Source/, then I just go to the Source directory and
> make ../bin/
>
> This solves the problem and I have an lhe file at the end. But also I
> noticed that once I create an lhe file from this gridpack I am having
> events that look like this:
>
> 11 2 0.9308000E+01 0.2172364E+03 0.7816531E-02 0.1151592E+00
> 21 -1 0 0 507 506 0.00000000000E+00 0.00000000000E+00 0.19779054729E+04 0.19779054729E+04 0.00000000000E+00 0. -1.
> 3 -1 0 0 508 0 0.00000000000E+00 0.00000000000E+00 -0.45643327256E+03 0.45643327256E+03 0.00000000000E+00 0. -1.
> 1000002 2 1 2 508 0 -0.84219622959E+02 0.31204399495E+02 -0.91462372374E+02 0.22651703404E+03 0.18675636994E+03 0. 0.
> 1000002 2 3 3 508 0 -0.85463213764E+01 0.56497649818E+02 -0.90908042272E+02 0.14672665639E+03 0.99997077464E+02 0. 0.
> -1000002 2 1 2 0 506 0.11087088756E+03 0.12118608604E+02 0.19052421139E+04 0.19111217424E+04 0.99997937121E+02 0. 0.
> -1 1 4 4 0 502 -0.53954020986E+02 0.38621826254E+02 -0.51063232565E+02 0.83726552322E+02 0.00000000000E+00 0. -1.
> -3 1 4 4 0 501 0.45407699609E+02 0.17875823564E+02 -0.39844809707E+02 0.63000104067E+02 0.00000000000E+00 0. -1.
> 1 1 5 5 505 0 0.58264518820E+02 -0.43511230204E+02 0.84926181828E+03 0.85236941363E+03 0.00000000000E+00 0. 1.
> 3 1 5 5 504 0 0.52606368736E+02 0.55629838808E+02 0.10559802956E+04 0.10587523288E+04 0.00000000000E+00 0. 1.
> 21 1 3 3 508 508 -0.75673301582E+02 -0.25293250324E+02 -0.55433010133E+00 0.79790377654E+02 0.00000000000E+00 0. -1.
> 3 1 1 2 507 0 -0.26651264597E+02 -0.43323008099E+02 -0.29230754121E+03 0.29669996896E+03 0.00000000000E+00 0. -1.
>
> that pythia8 cannot handle them.
>
>
> #### Run with 2.2.2
>
> Something weird happened here. I notice that I can create an the lhe
> file with the cluster, using the same cards as above in this version of
> MG. I got a lot of this messages:
>
> CRITICAL: Fail to run correctly job 36147.
> with option: {'log': None, 'stdout': None, 'argument': [], 'nb_submit': 1, 'stderr': None, 'prog': 'ajob1', 'output_files': ['G1', 'G3'], 'time_check': 1434535204.7950001, 'cwd': '/cms/gomez/
> file missing: /cms/gomez/
> Fails 1 times
> No resubmition
>
> but at the end I got an lhe file, with events that by eye look fine and I can run this using pythia8 without an issue.
> However, when I try to run the same for a gridpack, it doesn't work again.
>
> Hopefully this helps.
>
> cheers,
> Alejandro
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#13 |
Hi Olivier
>When I look at your log:
>> https:/
>I see:
https:/
>So it seems to me that you are using the 2.3.0 version and not 2.3.1. Is this correct?
Yes, but this was only one of my several test. For the others I am sure that I was using the one with 2.3.1.fix_
Then, I tried to download this: https:/
but everytime that I opened ./bin/mg5_aMC
I got this:
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 2.3.0.beta 2015-04-10 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* and *
* http://
* *
* Type 'help' for in-line help. *
* Type 'tutorial' to learn how MG5 works *
* Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
* Type 'tutorial MadLoop' to learn how MadLoop works *
* *
that as you mentioned is a previous version. What am I doing wrong?
thanks for all your patience.
cheers,
Alejandro
|
Revision history for this message
|
#14 |
Hi,
> Then, I tried to download this:
> https:/
Looks like my modification failed to be loaded to launchpad since someone was pushing at the same time.
This is solve now.
Can you retry?
> thanks for all your patience.
Thanks for yours.
Olivier
On 17 Jun 2015, at 14:16, alefisico <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> alefisico posted a new comment:
> Hi Olivier
>
>> When I look at your log:
>>> https:/
>
>> I see:
> https:/
>
>> So it seems to me that you are using the 2.3.0 version and not 2.3.1.
> Is this correct?
>
> Yes, but this was only one of my several test. For the others I am sure
> that I was using the one with 2.3.1.fix_
>
> Then, I tried to download this:
> https:/
>
> but everytime that I opened ./bin/mg5_aMC
> I got this:
> *******
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 2.3.0.beta 2015-04-10 *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https:/
> * and *
> * http://
> * *
> * Type 'help' for in-line help. *
> * Type 'tutorial' to learn how MG5 works *
> * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
> * Type 'tutorial MadLoop' to learn how MadLoop works *
> * *
>
>
> that as you mentioned is a previous version. What am I doing wrong?
>
> thanks for all your patience.
> cheers,
> Alejandro
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#15 |
HI
ok yes, downloading this version seems to work for this error, however now I am getting a different error. This is the output:
*******
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 5.2.3.1.beta *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https:/
* *
* Type 'help' for in-line help. *
* *
*******
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
INFO: load configuration from /afs/cern.
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
generate_events pilotrun
The following switches determine which programs are run:
1 Run the pythia shower/
2 Run PGS as detector simulator: pgs=NOT INSTALLED
3 Run Delphes as detector simulator: delphes=NOT INSTALLED
4 Decay particles with the MadSpin module: madspin=OFF
5 Add weight to events based on coupling parameters: reweight=OFF
Either type the switch number (1 to 5) to change its default setting,
or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
Type '0', 'auto', 'done' or just press enter when you are done.
[0, 4, 5, auto, done, madspin=ON, madspin=OFF, madspin, reweight=ON, ... ][60s to answer]
>Do you want to edit a card (press enter to bypass editing)?
1 / param : param_card.dat
2 / run : run_card.dat
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path][60s to answer]
>INFO: modify parameter gridpack of the run_card.dat to .true.
Do you want to edit a card (press enter to bypass editing)?
1 / param : param_card.dat
2 / run : run_card.dat
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path]
>INFO: modify param_card information BLOCK mass with id (1000002,) set to 100.0
Do you want to edit a card (press enter to bypass editing)?
1 / param : param_card.dat
2 / run : run_card.dat
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path]
>INFO: modify param_card information BLOCK decay with id (1000002,) set to 0.01
Do you want to edit a card (press enter to bypass editing)?
1 / param : param_card.dat
2 / run : run_card.dat
you can also
- enter the path to a valid card or banner.
- use the 'set' command to modify a parameter directly.
The set option works only for param_card and run_card.
Type 'help set' for more information on this command.
- call an external program (ASperGE/
Type 'help' for the list of available command
[0, done, 1, param, 2, run, enter path]
>Generating gridpack with run name pilotrun
survey pilotrun --accuracy=0.01 --points=2000 --iterations=8 --gridpack=.true.
INFO: compile directory
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P2_gg_su1su1xgg
INFO: P2_gg_su1su1xqq
INFO: P2_gg_su1su1xbb
INFO: P2_gq_su1su1xgq
INFO: P2_qq_su1su1xqq
INFO: P2_qq_su1su1xgg
INFO: P2_qq_su1su1xbb
INFO: P1_gg_su1su1xg_
INFO: P1_gq_su1su1xq_
INFO: P1_qq_su1su1xg_
INFO: P0_gg_su1su1x_
INFO: P0_qq_su1su1x_
INFO: Idle: 139, Running: 17, Completed: 0 [ 0.11s ]
INFO: Idle: 139, Running: 17, Completed: 0 [ 0.22s ]
INFO: Idle: 125, Running: 26, Completed: 5 [ 30.8s ]
INFO: Idle: 110, Running: 29, Completed: 17 [ 1m 1s ]
INFO: Idle: 102, Running: 24, Completed: 30 [ 1m 31s ]
INFO: Idle: 65, Running: 42, Completed: 49 [ 2m 1s ]
INFO: Start to wait 60s between checking status.
Note that you can change this time in the configuration file.
Press ctrl-C to force the update.
INFO: Idle: 0, Running: 58, Completed: 98 [ 3m 1s ]
INFO: Idle: 0, Running: 47, Completed: 109 [ 3m 32s ]
INFO: Idle: 0, Running: 9, Completed: 147 [ 4m 2s ]
INFO: Idle: 0, Running: 7, Completed: 149 [ 4m 32s ]
INFO: Idle: 0, Running: 5, Completed: 151 [ 5m 2s ]
INFO: Idle: 0, Running: 3, Completed: 153 [ 5m 32s ]
INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 2s ]
INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 33s ]
INFO: Idle: 0, Running: 1, Completed: 155 [ 7m 3s ]
INFO: All jobs finished
INFO: Idle: 0, Running: 0, Completed: 156 [ 7m 33s ]
Error when reading /afs/cern.
Command "generate_events pilotrun" interrupted with error:
ValueError : need more than 5 values to unpack
Please report this bug on https:/
More information is found in '/afs/cern.
Please attach this file to your report.
quit
INFO:
INFO:
and this is the pilotrun_tag.log:
https:/
cheers,
|
Revision history for this message
|
#16 |
Thanks,
I have run g g > su1 su1~ g g (with decay) since looks like that’s the responsible of the problem.
without gridpack it works (at least up to the survey, I stop the job after that).
I’m rerunning it now in gridpack mode to see if it changes anything.
(it should not since the normal mode achieve already a very good precision in the channel crashing)
if it goes trough, I will force the matching/merging to be on True.
Thanks,
Olivier
On 17 Jun 2015, at 22:41, alefisico <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Answered => Open
>
> alefisico is still having a problem:
> HI
>
> ok yes, downloading this version seems to work for this error, however
> now I am getting a different error. This is the output:
>
> *******
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * M A D E V E N T *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 5.2.3.1.beta *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https:/
> * *
> * Type 'help' for in-line help. *
> * *
> *******
> INFO: load configuration from /afs/cern.
> INFO: load configuration from /afs/cern.
> INFO: load configuration from /afs/cern.
> Using default eps viewer "evince". Set another one in ./input/
> Using default web browser "firefox". Set another one in ./input/
> generate_events pilotrun
> The following switches determine which programs are run:
> 1 Run the pythia shower/
> 2 Run PGS as detector simulator: pgs=NOT INSTALLED
> 3 Run Delphes as detector simulator: delphes=NOT INSTALLED
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weight to events based on coupling parameters: reweight=OFF
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 4, 5, auto, done, madspin=ON, madspin=OFF, madspin, reweight=ON, ... ][60s to answer]
>> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, enter path][60s to answer]
>> INFO: modify parameter gridpack of the run_card.dat to .true.
> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, enter path]
>> INFO: modify param_card information BLOCK mass with id (1000002,) set to 100.0
> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, enter path]
>> INFO: modify param_card information BLOCK decay with id (1000002,) set to 0.01
> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> you can also
> - enter the path to a valid card or banner.
> - use the 'set' command to modify a parameter directly.
> The set option works only for param_card and run_card.
> Type 'help set' for more information on this command.
> - call an external program (ASperGE/
> Type 'help' for the list of available command
> [0, done, 1, param, 2, run, enter path]
>> Generating gridpack with run name pilotrun
> survey pilotrun --accuracy=0.01 --points=2000 --iterations=8 --gridpack=.true.
> INFO: compile directory
> compile Source Directory
> Using random number seed offset = 21
> INFO: Running Survey
> Creating Jobs
> Working on SubProcesses
> INFO: P2_gg_su1su1xgg
> INFO: P2_gg_su1su1xqq
> INFO: P2_gg_su1su1xbb
> INFO: P2_gq_su1su1xgq
> INFO: P2_qq_su1su1xqq
> INFO: P2_qq_su1su1xgg
> INFO: P2_qq_su1su1xbb
> INFO: P1_gg_su1su1xg_
> INFO: P1_gq_su1su1xq_
> INFO: P1_qq_su1su1xg_
> INFO: P0_gg_su1su1x_
> INFO: P0_qq_su1su1x_
> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.11s ]
> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.22s ]
> INFO: Idle: 125, Running: 26, Completed: 5 [ 30.8s ]
> INFO: Idle: 110, Running: 29, Completed: 17 [ 1m 1s ]
> INFO: Idle: 102, Running: 24, Completed: 30 [ 1m 31s ]
> INFO: Idle: 65, Running: 42, Completed: 49 [ 2m 1s ]
> INFO: Start to wait 60s between checking status.
> Note that you can change this time in the configuration file.
> Press ctrl-C to force the update.
> INFO: Idle: 0, Running: 58, Completed: 98 [ 3m 1s ]
> INFO: Idle: 0, Running: 47, Completed: 109 [ 3m 32s ]
> INFO: Idle: 0, Running: 9, Completed: 147 [ 4m 2s ]
> INFO: Idle: 0, Running: 7, Completed: 149 [ 4m 32s ]
> INFO: Idle: 0, Running: 5, Completed: 151 [ 5m 2s ]
> INFO: Idle: 0, Running: 3, Completed: 153 [ 5m 32s ]
> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 2s ]
> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 33s ]
> INFO: Idle: 0, Running: 1, Completed: 155 [ 7m 3s ]
> INFO: All jobs finished
> INFO: Idle: 0, Running: 0, Completed: 156 [ 7m 33s ]
> Error when reading /afs/cern.
> Command "generate_events pilotrun" interrupted with error:
> ValueError : need more than 5 values to unpack
> Please report this bug on https:/
> More information is found in '/afs/cern.
> Please attach this file to your report.
> quit
> INFO:
>
> INFO:
>
>
> and this is the pilotrun_tag.log:
> https:/
>
> cheers,
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#17 |
Hi,
Ok i have run the gridpack generation with both matching/merging on/off
and both succeed…
I’m going to run it on a linux cluster as a last hope…
Cheers,
Olivier
On 18 Jun 2015, at 00:16, Olivier Mattelaer <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Thanks,
>
> I have run g g > su1 su1~ g g (with decay) since looks like that’s the responsible of the problem.
> without gridpack it works (at least up to the survey, I stop the job after that).
>
> I’m rerunning it now in gridpack mode to see if it changes anything.
> (it should not since the normal mode achieve already a very good precision in the channel crashing)
> if it goes trough, I will force the matching/merging to be on True.
>
> Thanks,
>
> Olivier
>
> On 17 Jun 2015, at 22:41, alefisico
> <email address hidden> wrote:
>
>> Question #265751 on MadGraph5_aMC@NLO changed:
>> https:/
>>
>> Status: Answered => Open
>>
>> alefisico is still having a problem:
>> HI
>>
>> ok yes, downloading this version seems to work for this error, however
>> now I am getting a different error. This is the output:
>>
>> *******
>> * *
>> * W E L C O M E to *
>> * M A D G R A P H 5 _ a M C @ N L O *
>> * M A D E V E N T *
>> * *
>> * * * *
>> * * * * * *
>> * * * * * 5 * * * * *
>> * * * * * *
>> * * * *
>> * *
>> * VERSION 5.2.3.1.beta *
>> * *
>> * The MadGraph5_aMC@NLO Development Team - Find us at *
>> * https:/
>> * *
>> * Type 'help' for in-line help. *
>> * *
>> *******
>> INFO: load configuration from /afs/cern.
>> INFO: load configuration from /afs/cern.
>> INFO: load configuration from /afs/cern.
>> Using default eps viewer "evince". Set another one in ./input/
>> Using default web browser "firefox". Set another one in ./input/
>> generate_events pilotrun
>> The following switches determine which programs are run:
>> 1 Run the pythia shower/
>> 2 Run PGS as detector simulator: pgs=NOT INSTALLED
>> 3 Run Delphes as detector simulator: delphes=NOT INSTALLED
>> 4 Decay particles with the MadSpin module: madspin=OFF
>> 5 Add weight to events based on coupling parameters: reweight=OFF
>> Either type the switch number (1 to 5) to change its default setting,
>> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
>> Type '0', 'auto', 'done' or just press enter when you are done.
>> [0, 4, 5, auto, done, madspin=ON, madspin=OFF, madspin, reweight=ON, ... ][60s to answer]
>>> Do you want to edit a card (press enter to bypass editing)?
>> 1 / param : param_card.dat
>> 2 / run : run_card.dat
>> you can also
>> - enter the path to a valid card or banner.
>> - use the 'set' command to modify a parameter directly.
>> The set option works only for param_card and run_card.
>> Type 'help set' for more information on this command.
>> - call an external program (ASperGE/
>> Type 'help' for the list of available command
>> [0, done, 1, param, 2, run, enter path][60s to answer]
>>> INFO: modify parameter gridpack of the run_card.dat to .true.
>> Do you want to edit a card (press enter to bypass editing)?
>> 1 / param : param_card.dat
>> 2 / run : run_card.dat
>> you can also
>> - enter the path to a valid card or banner.
>> - use the 'set' command to modify a parameter directly.
>> The set option works only for param_card and run_card.
>> Type 'help set' for more information on this command.
>> - call an external program (ASperGE/
>> Type 'help' for the list of available command
>> [0, done, 1, param, 2, run, enter path]
>>> INFO: modify param_card information BLOCK mass with id (1000002,) set to 100.0
>> Do you want to edit a card (press enter to bypass editing)?
>> 1 / param : param_card.dat
>> 2 / run : run_card.dat
>> you can also
>> - enter the path to a valid card or banner.
>> - use the 'set' command to modify a parameter directly.
>> The set option works only for param_card and run_card.
>> Type 'help set' for more information on this command.
>> - call an external program (ASperGE/
>> Type 'help' for the list of available command
>> [0, done, 1, param, 2, run, enter path]
>>> INFO: modify param_card information BLOCK decay with id (1000002,) set to 0.01
>> Do you want to edit a card (press enter to bypass editing)?
>> 1 / param : param_card.dat
>> 2 / run : run_card.dat
>> you can also
>> - enter the path to a valid card or banner.
>> - use the 'set' command to modify a parameter directly.
>> The set option works only for param_card and run_card.
>> Type 'help set' for more information on this command.
>> - call an external program (ASperGE/
>> Type 'help' for the list of available command
>> [0, done, 1, param, 2, run, enter path]
>>> Generating gridpack with run name pilotrun
>> survey pilotrun --accuracy=0.01 --points=2000 --iterations=8 --gridpack=.true.
>> INFO: compile directory
>> compile Source Directory
>> Using random number seed offset = 21
>> INFO: Running Survey
>> Creating Jobs
>> Working on SubProcesses
>> INFO: P2_gg_su1su1xgg
>> INFO: P2_gg_su1su1xqq
>> INFO: P2_gg_su1su1xbb
>> INFO: P2_gq_su1su1xgq
>> INFO: P2_qq_su1su1xqq
>> INFO: P2_qq_su1su1xgg
>> INFO: P2_qq_su1su1xbb
>> INFO: P1_gg_su1su1xg_
>> INFO: P1_gq_su1su1xq_
>> INFO: P1_qq_su1su1xg_
>> INFO: P0_gg_su1su1x_
>> INFO: P0_qq_su1su1x_
>> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.11s ]
>> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.22s ]
>> INFO: Idle: 125, Running: 26, Completed: 5 [ 30.8s ]
>> INFO: Idle: 110, Running: 29, Completed: 17 [ 1m 1s ]
>> INFO: Idle: 102, Running: 24, Completed: 30 [ 1m 31s ]
>> INFO: Idle: 65, Running: 42, Completed: 49 [ 2m 1s ]
>> INFO: Start to wait 60s between checking status.
>> Note that you can change this time in the configuration file.
>> Press ctrl-C to force the update.
>> INFO: Idle: 0, Running: 58, Completed: 98 [ 3m 1s ]
>> INFO: Idle: 0, Running: 47, Completed: 109 [ 3m 32s ]
>> INFO: Idle: 0, Running: 9, Completed: 147 [ 4m 2s ]
>> INFO: Idle: 0, Running: 7, Completed: 149 [ 4m 32s ]
>> INFO: Idle: 0, Running: 5, Completed: 151 [ 5m 2s ]
>> INFO: Idle: 0, Running: 3, Completed: 153 [ 5m 32s ]
>> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 2s ]
>> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 33s ]
>> INFO: Idle: 0, Running: 1, Completed: 155 [ 7m 3s ]
>> INFO: All jobs finished
>> INFO: Idle: 0, Running: 0, Completed: 156 [ 7m 33s ]
>> Error when reading /afs/cern.
>> Command "generate_events pilotrun" interrupted with error:
>> ValueError : need more than 5 values to unpack
>> Please report this bug on https:/
>> More information is found in '/afs/cern.
>> Please attach this file to your report.
>> quit
>> INFO:
>>
>> INFO:
>>
>>
>> and this is the pilotrun_tag.log:
>> https:/
>>
>> cheers,
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#18 |
Works also on the linux system on which I have access.
So I can not reproduce this …
I’m sorry but I do not know how I can help you with that.
Cheers,
Olivier
On 18 Jun 2015, at 16:06, Olivier Mattelaer <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> Ok i have run the gridpack generation with both matching/merging on/off
> and both succeed…
>
> I’m going to run it on a linux cluster as a last hope…
>
> Cheers,
>
> Olivier
>
>
> On 18 Jun 2015, at 00:16, Olivier Mattelaer
> <email address hidden> wrote:
>
>> Question #265751 on MadGraph5_aMC@NLO changed:
>> https:/
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Thanks,
>>
>> I have run g g > su1 su1~ g g (with decay) since looks like that’s the responsible of the problem.
>> without gridpack it works (at least up to the survey, I stop the job after that).
>>
>> I’m rerunning it now in gridpack mode to see if it changes anything.
>> (it should not since the normal mode achieve already a very good precision in the channel crashing)
>> if it goes trough, I will force the matching/merging to be on True.
>>
>> Thanks,
>>
>> Olivier
>>
>> On 17 Jun 2015, at 22:41, alefisico
>> <email address hidden> wrote:
>>
>>> Question #265751 on MadGraph5_aMC@NLO changed:
>>> https:/
>>>
>>> Status: Answered => Open
>>>
>>> alefisico is still having a problem:
>>> HI
>>>
>>> ok yes, downloading this version seems to work for this error, however
>>> now I am getting a different error. This is the output:
>>>
>>> *******
>>> * *
>>> * W E L C O M E to *
>>> * M A D G R A P H 5 _ a M C @ N L O *
>>> * M A D E V E N T *
>>> * *
>>> * * * *
>>> * * * * * *
>>> * * * * * 5 * * * * *
>>> * * * * * *
>>> * * * *
>>> * *
>>> * VERSION 5.2.3.1.beta *
>>> * *
>>> * The MadGraph5_aMC@NLO Development Team - Find us at *
>>> * https:/
>>> * *
>>> * Type 'help' for in-line help. *
>>> * *
>>> *******
>>> INFO: load configuration from /afs/cern.
>>> INFO: load configuration from /afs/cern.
>>> INFO: load configuration from /afs/cern.
>>> Using default eps viewer "evince". Set another one in ./input/
>>> Using default web browser "firefox". Set another one in ./input/
>>> generate_events pilotrun
>>> The following switches determine which programs are run:
>>> 1 Run the pythia shower/
>>> 2 Run PGS as detector simulator: pgs=NOT INSTALLED
>>> 3 Run Delphes as detector simulator: delphes=NOT INSTALLED
>>> 4 Decay particles with the MadSpin module: madspin=OFF
>>> 5 Add weight to events based on coupling parameters: reweight=OFF
>>> Either type the switch number (1 to 5) to change its default setting,
>>> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
>>> Type '0', 'auto', 'done' or just press enter when you are done.
>>> [0, 4, 5, auto, done, madspin=ON, madspin=OFF, madspin, reweight=ON, ... ][60s to answer]
>>>> Do you want to edit a card (press enter to bypass editing)?
>>> 1 / param : param_card.dat
>>> 2 / run : run_card.dat
>>> you can also
>>> - enter the path to a valid card or banner.
>>> - use the 'set' command to modify a parameter directly.
>>> The set option works only for param_card and run_card.
>>> Type 'help set' for more information on this command.
>>> - call an external program (ASperGE/
>>> Type 'help' for the list of available command
>>> [0, done, 1, param, 2, run, enter path][60s to answer]
>>>> INFO: modify parameter gridpack of the run_card.dat to .true.
>>> Do you want to edit a card (press enter to bypass editing)?
>>> 1 / param : param_card.dat
>>> 2 / run : run_card.dat
>>> you can also
>>> - enter the path to a valid card or banner.
>>> - use the 'set' command to modify a parameter directly.
>>> The set option works only for param_card and run_card.
>>> Type 'help set' for more information on this command.
>>> - call an external program (ASperGE/
>>> Type 'help' for the list of available command
>>> [0, done, 1, param, 2, run, enter path]
>>>> INFO: modify param_card information BLOCK mass with id (1000002,) set to 100.0
>>> Do you want to edit a card (press enter to bypass editing)?
>>> 1 / param : param_card.dat
>>> 2 / run : run_card.dat
>>> you can also
>>> - enter the path to a valid card or banner.
>>> - use the 'set' command to modify a parameter directly.
>>> The set option works only for param_card and run_card.
>>> Type 'help set' for more information on this command.
>>> - call an external program (ASperGE/
>>> Type 'help' for the list of available command
>>> [0, done, 1, param, 2, run, enter path]
>>>> INFO: modify param_card information BLOCK decay with id (1000002,) set to 0.01
>>> Do you want to edit a card (press enter to bypass editing)?
>>> 1 / param : param_card.dat
>>> 2 / run : run_card.dat
>>> you can also
>>> - enter the path to a valid card or banner.
>>> - use the 'set' command to modify a parameter directly.
>>> The set option works only for param_card and run_card.
>>> Type 'help set' for more information on this command.
>>> - call an external program (ASperGE/
>>> Type 'help' for the list of available command
>>> [0, done, 1, param, 2, run, enter path]
>>>> Generating gridpack with run name pilotrun
>>> survey pilotrun --accuracy=0.01 --points=2000 --iterations=8 --gridpack=.true.
>>> INFO: compile directory
>>> compile Source Directory
>>> Using random number seed offset = 21
>>> INFO: Running Survey
>>> Creating Jobs
>>> Working on SubProcesses
>>> INFO: P2_gg_su1su1xgg
>>> INFO: P2_gg_su1su1xqq
>>> INFO: P2_gg_su1su1xbb
>>> INFO: P2_gq_su1su1xgq
>>> INFO: P2_qq_su1su1xqq
>>> INFO: P2_qq_su1su1xgg
>>> INFO: P2_qq_su1su1xbb
>>> INFO: P1_gg_su1su1xg_
>>> INFO: P1_gq_su1su1xq_
>>> INFO: P1_qq_su1su1xg_
>>> INFO: P0_gg_su1su1x_
>>> INFO: P0_qq_su1su1x_
>>> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.11s ]
>>> INFO: Idle: 139, Running: 17, Completed: 0 [ 0.22s ]
>>> INFO: Idle: 125, Running: 26, Completed: 5 [ 30.8s ]
>>> INFO: Idle: 110, Running: 29, Completed: 17 [ 1m 1s ]
>>> INFO: Idle: 102, Running: 24, Completed: 30 [ 1m 31s ]
>>> INFO: Idle: 65, Running: 42, Completed: 49 [ 2m 1s ]
>>> INFO: Start to wait 60s between checking status.
>>> Note that you can change this time in the configuration file.
>>> Press ctrl-C to force the update.
>>> INFO: Idle: 0, Running: 58, Completed: 98 [ 3m 1s ]
>>> INFO: Idle: 0, Running: 47, Completed: 109 [ 3m 32s ]
>>> INFO: Idle: 0, Running: 9, Completed: 147 [ 4m 2s ]
>>> INFO: Idle: 0, Running: 7, Completed: 149 [ 4m 32s ]
>>> INFO: Idle: 0, Running: 5, Completed: 151 [ 5m 2s ]
>>> INFO: Idle: 0, Running: 3, Completed: 153 [ 5m 32s ]
>>> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 2s ]
>>> INFO: Idle: 0, Running: 2, Completed: 154 [ 6m 33s ]
>>> INFO: Idle: 0, Running: 1, Completed: 155 [ 7m 3s ]
>>> INFO: All jobs finished
>>> INFO: Idle: 0, Running: 0, Completed: 156 [ 7m 33s ]
>>> Error when reading /afs/cern.
>>> Command "generate_events pilotrun" interrupted with error:
>>> ValueError : need more than 5 values to unpack
>>> Please report this bug on https:/
>>> More information is found in '/afs/cern.
>>> Please attach this file to your report.
>>> quit
>>> INFO:
>>>
>>> INFO:
>>>
>>>
>>> and this is the pilotrun_tag.log:
>>> https:/
>>>
>>> cheers,
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> You received this question notification because you are an answer
>> contact for MadGraph5_aMC@NLO.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
|
Revision history for this message
|
#19 |
Hi Olivier
It is very weird. Just as last check, can you send me the cards that you are using?
thanks
|
Revision history for this message
|
#20 |
Sure here is my script:
import model ../RPVMSSM_UFO --modelname
define excl = su1 su2 su3 su4 su5 su6 sd1 sd2 sd3 sd4 sd5 sd6 su1~ su2\
~ su3~ su4~ su5~ su6~ sd1~ sd2~ sd3~ sd4~ sd5~ sd6~ x1+ x2+ x1- x2-
generate g g > su1 su1~ g g /su2 $ excl QED=0, su1 > d~ s~, su1~ > d s
output
launch
set gridpack T
set mass 1000002 100
set width 1000002 0.01
set ickkw 1
set xqcut 30
Cheers,
Olivier
On 19 Jun 2015, at 11:37, alefisico <email address hidden> wrote:
> Question #265751 on MadGraph5_aMC@NLO changed:
> https:/
>
> alefisico posted a new comment:
> Hi Olivier
>
> It is very weird. Just as last check, can you send me the cards that you
> are using?
>
> thanks
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
Can you help with this problem?
Provide an answer of your own, or ask alefisico for more information if necessary.
