End code 2.0
I'm running the MadDump plugin and get the following error from MadGraph. I'd appreciate any advice you have on what to. Here is the contents of the error log file. Thank you.
generate_events run_01 -f
Traceback (most recent call last):
File "/home/
return self.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
self.
File "/home/
original_
File "/home/
self.
File "/home/
stop = Cmd.onecmd_
File "/home/
return func(arg, **opt)
File "/home/
x_improve.
File "/home/
main_
File "/home/
P_comb.
File "/home/
oneresult.
File "/home/
% (error_code, open(log).read())
Exception: Reported error: End code 2.0
Full associated log:
Process in group number 1200
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11839999999999999
*****
NNPDFDriver version 1.0.3
Grid: NNPDF23_
*****
New value of alpha_s from PDF nn23lo1: 0.13000000000000000
Define smin to 1.2327121000000
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
mdl_cabi = 0.22773599999999999
aEWM1 = 127.90000000000001
mdl_Gf = 1.1663700000000
aS = 0.11839999999999999
mdl_ymdo = 5.0400000000000
mdl_ymup = 2.5500000000000
mdl_yms = 0.10100000000000001
mdl_ymc = 1.2700000000000000
mdl_ymb = 4.7000000000000002
mdl_ymt = 172.00000000000000
mdl_yme = 5.1099999999999
mdl_ymm = 0.10566000000000000
mdl_ymtau = 1.7769999999999999
mdl_muBI = 1.0000000000000000
mdl_muChi = 0.0000000000000000
mdl_dChi = 0.0000000000000000
mdl_eps = 1.0000000000000000
mdl_MZ = 91.187600000000003
mdl_Me = 5.1099999999999
mdl_MMU = 0.10566000000000000
mdl_MTA = 1.7769999999999999
mdl_MU = 2.5500000000000
mdl_MC = 1.2700000000000000
mdl_MT = 172.00000000000000
mdl_MD = 5.0400000000000
mdl_MS = 0.10100000000000001
mdl_MB = 4.7000000000000002
mdl_MH = 125.00000000000000
mdl_mChi = 3.0000000000000
mdl_WZ = 2.4952000000000001
mdl_WW = 2.0850000000000000
mdl_WT = 1.5083359999999999
mdl_WH = 4.0699999999999
mdl_Mpi0 = 0.13497999999999999
mdl_Wpi0 = 7.8000000000000
mdl_Meta = 0.54800000000000004
mdl_Weta = 1.3100000000000
mdl_Momega = 0.78300000000000003
mdl_Womega = 8.5000000000000
mdl_Mrho = 0.77500000000000002
mdl_Wrho = 0.14699999999999999
mdl_Mphi = 1.0189999999999999
mdl_Wphi = 4.1999999999999
mdl_Mjpsi = 3.0970000000000000
mdl_Wjpsi = 9.2999999999999
Internal Params
-----
mdl_cos__cabi = 0.97418004031982097
mdl_CKM1x1 = (0.974180040319
mdl_sin__cabi = 0.22577256042856930
mdl_CKM1x2 = (0.225772560428
mdl_CKM1x3 = (0.000000000000
mdl_CKM2x1 = (-0.22577256042
mdl_CKM2x2 = (0.974180040319
mdl_CKM2x3 = (0.000000000000
mdl_CKM3x1 = (0.000000000000
mdl_CKM3x2 = (0.000000000000
mdl_CKM3x3 = (1.000000000000
mdl_MZ__exp__2 = 8315.1783937600012
mdl_MZ__exp__4 = 69142191.720053151
mdl_sqrt__2 = 1.4142135623730951
mdl_MH__exp__2 = 15625.000000000000
mdl_complexi = (0.000000000000
mdl_conjg__CKM1x1 = (0.974180040319
mdl_conjg__CKM1x2 = (0.225772560428
mdl_conjg__CKM1x3 = (0.000000000000
mdl_conjg__CKM2x1 = (-0.22577256042
At line 68 of file splev.f
Fortran runtime error: Index '4' of dimension 1 of array 't' above upper bound of 0
Error termination. Backtrace:
mdl_conjg__CKM2x2 = (0.974180040319
mdl_conjg__CKM2x3 = (0.000000000000
mdl_conjg__CKM3x1 = (0.000000000000
mdl_conjg__CKM3x2 = (0.000000000000
mdl_conjg__CKM3x3 = (1.000000000000
mdl_aEW = 7.8186082877247
mdl_MW = 79.824359746197842
mdl_sqrt__aEW = 8.8422894590285
mdl_ee = 0.31345100004952897
mdl_MW__exp__2 = 6371.9284088904105
mdl_sw2 = 0.23369913342182447
mdl_cw = 0.87538612427783857
mdl_sqrt__sw2 = 0.48342438232036300
mdl_sw = 0.48342438232036300
mdl_g1 = 0.35807170271074895
mdl_gw = 0.64839716719502682
mdl_vev = 246.22056907348590
mdl_vev__exp__2 = 60624.568634871241
mdl_lam = 0.12886689630821144
mdl_yb = 2.6995322804122
mdl_yc = 7.2944808428161
mdl_ydo = 2.8948175943144
mdl_ye = 2.9350233942354
mdl_ym = 6.0687783137949
mdl_ys = 5.8011225600348
mdl_yt = 0.98791394091683138
mdl_ytau = 1.0206529494239
mdl_yup = 1.4646398542662
mdl_muH = 88.388347648318444
mdl_ee__exp__2 = 9.8251529432049
mdl_sw__exp__2 = 0.23369913342182450
mdl_cw__exp__2 = 0.76630086657817542
Internal Params evaluated point by point
-----
mdl_sqrt__aS = 0.34409301068170506
mdl_G__exp__2 = 1.4878582807401259
Couplings of milliDipoles
-----
GC_3 = -0.00000E+00 -0.31345E+00
GC_6 = 0.00000E+00 0.31345E+00
GC_11 = 0.00000E+00 0.00000E+00
GC_12 = 0.00000E+00 0.00000E+00
Collider parameters:
------
Running at aUe machine @ 0.62330872802055315 GeV
PDF set = nn23lo1
alpha_s(Mz)= 0.1300 running at 2 loops.
alpha_s(Mz)= 0.1300 running at 2 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis
getting user params
Enter number of events and max and min iterations:
Number of events and iterations 1000 5 3
Enter desired fractional accuracy:
Desired fractional accuracy: 0.10000000000000001
Enter 0 for fixed, 2 for adjustable grid:
Suppress amplitude (0 no, 1 yes)?
Using suppressed amplitude.
Exact helicity sum (0 yes, n = number/event)?
Explicitly summing over helicities
Enter Configuration Number:
Running Configuration Number: 1
Not subdividing B.W.
Attempting mappinvarients 1 4
Completed mapping 4
about to integrate 3 1000 5 3 3 1
Using non-zero grid deformation.
3 dimensions 1000 events 3 invarients 5 iterations 1 config(s), (0.99)
Using h-tuple random number sequence.
Error opening grid
Using Uniform Grid! 16
Using uniform alpha 1.0000000000000000
Grid defined OK
Set CM energy to 0.44
Mapping Graph 1 to config 1
Setting grid 1 0.13441E-05 1
Warning: No minimum found for integration
Setting minimum to 5.1475772741132415
Transforming s_hat 1/s 3 5.1475772741132415 1.2327121000000
1 1 2 3
Masses: 0.300E-02 0.511E-03 0.300E-02 0.511E-03
Using random seed offsets 1 : 1
with seed 24
Ranmar initialization seeds 31059 9398
#0 0x7ff69d82c32a
#1 0x7ff69d82ced5
#2 0x7ff69d82d2a7
#3 0x556159e12b65
#4 0x556159dc5a8a
#5 0x556159dc5d73
#6 0x556159d9f0cb
#7 0x556159d9f9e7
#8 0x556159ddc52b
#9 0x556159d98bdc
#10 0x556159d98d4f
#11 0x7ff69ce8bb96
#12 0x556159d96a69
#13 0xffffffffffffffff
ls status:
input_app.txt
param.log
run1_app.log
automatic_
cluster_
#******
#* MadGraph5_aMC@NLO *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.7.3 2020-06-21 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https:/
#* *
#******
#* *
#* Command File for MadGraph5_aMC@NLO *
#* *
#* run as ./bin/mg5_aMC filename *
#* *
#******
set default_
set group_subprocesses Auto
set ignore_
set loop_optimized_
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_
import model milliDipoles
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_events decay /home/maroccol/
hepmc
decay phi > chi chi~
define darkmatter chi
add process interaction @electron
output output/BEBC/mQ_phi
#######
## PARAM_CARD AUTOMATICALY GENERATED BY MG5 ####
#######
#######
## INFORMATION FOR CKMBLOCK
#######
BLOCK CKMBLOCK #
1 2.277360e-01 # cabi
#######
## INFORMATION FOR FRBLOCK
#######
BLOCK FRBLOCK #
1 1.000000e+00 # mubi
2 0.000000e+00 # muchi
3 0.000000e+00 # dchi
4 1.000000e+00 # eps
#######
## INFORMATION FOR MASS
#######
BLOCK MASS #
1 5.040000e-03 # md
2 2.550000e-03 # mu
3 1.010000e-01 # ms
4 1.270000e+00 # mc
5 4.700000e+00 # mb
6 1.720000e+02 # mt
11 5.110000e-04 # me
13 1.056600e-01 # mmu
15 1.777000e+00 # mta
23 9.118760e+01 # mz
25 1.250000e+02 # mh
111 1.349800e-01 # mpi0
113 7.750000e-01 # mrho
221 5.480000e-01 # meta
223 7.830000e-01 # momega
333 1.019000e+00 # mphi
443 3.097000e+00 # mjpsi
9000005 0.003 # mchi
12 0.000000e+00 # ve : 0.0
14 0.000000e+00 # vm : 0.0
16 0.000000e+00 # vt : 0.0
21 0.000000e+00 # g : 0.0
22 0.000000e+00 # a : 0.0
24 7.982436e+01 # w+ : cmath.sqrt(
#######
## INFORMATION FOR SMINPUTS
#######
BLOCK SMINPUTS #
1 1.279000e+02 # aewm1
2 1.166370e-05 # gf
3 1.184000e-01 # as
#######
## INFORMATION FOR YUKAWA
#######
BLOCK YUKAWA #
1 5.040000e-03 # ymdo
2 2.550000e-03 # ymup
3 1.010000e-01 # yms
4 1.270000e+00 # ymc
5 4.700000e+00 # ymb
6 1.720000e+02 # ymt
11 5.110000e-04 # yme
13 1.056600e-01 # ymm
15 1.777000e+00 # ymtau
#######
## INFORMATION FOR DECAY
#######
DECAY 6 1.508336e+00 # wt
DECAY 23 2.495200e+00 # wz
DECAY 24 2.085000e+00 # ww
DECAY 25 4.070000e-03 # wh
DECAY 111 7.800000e-09 # wpi0
DECAY 113 1.470000e-01 # wrho
DECAY 221 1.310000e-06 # weta
DECAY 223 8.500000e-03 # womega
DECAY 333 4.200000e-03 # wphi
DECAY 443 9.300000e-02 # wjpsi
DECAY 1 0.000000e+00 # d : 0.0
DECAY 2 0.000000e+00 # u : 0.0
DECAY 3 0.000000e+00 # s : 0.0
DECAY 4 0.000000e+00 # c : 0.0
DECAY 5 0.000000e+00 # b : 0.0
DECAY 11 0.000000e+00 # e- : 0.0
DECAY 12 0.000000e+00 # ve : 0.0
DECAY 13 0.000000e+00 # mu- : 0.0
DECAY 14 0.000000e+00 # vm : 0.0
DECAY 15 0.000000e+00 # ta- : 0.0
DECAY 16 0.000000e+00 # vt : 0.0
DECAY 21 0.000000e+00 # g : 0.0
DECAY 22 0.000000e+00 # a : 0.0
DECAY 9000005 0.000000e+00 # chi : 0.0
#######
## INFORMATION FOR QNUMBERS 9000005
#######
BLOCK QNUMBERS 9000005 # chi
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 1 # particle/
#######
## INFORMATION FOR QNUMBERS 111
#######
BLOCK QNUMBERS 111 # pi0
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/
#######
## INFORMATION FOR QNUMBERS 221
#######
BLOCK QNUMBERS 221 # eta
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/
#######
## INFORMATION FOR QNUMBERS 223
#######
BLOCK QNUMBERS 223 # omega
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 113
#######
BLOCK QNUMBERS 113 # rho
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 333
#######
BLOCK QNUMBERS 333 # phi
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # particle/
#######
## INFORMATION FOR QNUMBERS 443
#######
BLOCK QNUMBERS 443 # jpsi
1 0 # 3 times electric charge
2 3 # number of spin states (2s+1)
3 1 # colour rep (1: singlet, 3: triplet, 8: octet)
4 0 # 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 *
# *
# To display more options, you can type the command: *
# update full_run_card *
#******
#
#******
# Tag name for the run (one word) *
#******
tag_1 = run_tag ! name of the run
#******
# Number of events and rnd seed *
# Warning: Do not generate more than 1M events in a single run *
#******
312 = 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 *
#******
9 = lpp1 ! beam 1 type
0 = lpp2 ! beam 2 type
190.07523015 = ebeam1 ! beam 1 total energy in GeV
0.000511 = ebeam2 ! beam 2 total energy in GeV
# To see polarised beam options: type "update beam_pol"
# 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
# To see heavy ion options: type "update ion_pdf"
# Renormalization and factorization scales *
#******
False = fixed_ren_scale ! if .true. use fixed ren scale
False = fixed_fac_scale ! if .true. use fixed fac scale
91.188 = scale ! fixed ren scale
91.188 = dsqrt_q2fact1 ! fixed fact scale for pdf1
91.188 = dsqrt_q2fact2 ! fixed fact scale for pdf2
-1 = dynamical_
1.0 = scalefact ! scale factor for event-by-event scales
#******
# Type and output format
#******
False = gridpack !True = setting up the grid pack
-1.0 = time_of_flight ! threshold (in mm) below which the invariant livetime is not written (-1 means not written)
average = event_norm ! average/sum. Normalization of the weight in the LHEF
# To see MLM/CKKW merging options: type "update MLM" or "update CKKW"#*
#
#******
# handling of the helicities:
# 0: sum over all helicities
# 1: importance sampling over helicities
#******
0 = nhel ! using helicities importance sampling or not.
#******
# Generation bias, check the wiki page below for more information: *
# 'cp3.irmp.
#******
None = bias_module ! Bias type of bias, [None, ptj_bias, -custom_folder-]
{} = bias_parameters ! Specifies the parameters of the module.
#
#******
# Parton level cuts definition *
#******
#
#
#******
# BW cutoff (M+/-bwcutoff*
#******
15.0 = bwcutoff ! (M+/-bwcutoff*
#IF(d)# #******
#IF(d)# # Apply pt/E/eta/
#IF(d)# # (note that etmiss/
#IF(d)# #******
False = cut_decays ! Cut decay products
#******
# Standard Cuts *
#******
# Minimum and maximum pt's (for max, -1 means no cut) *
#******
0.0 = ptj ! minimum pt for the jets
0.0 = ptb ! minimum pt for the b
0.0 = pta ! minimum pt for the photons
0.0 = ptl ! minimum pt for the charged leptons
0.0 = misset ! minimum missing Et (sum of neutrino's momenta)
-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)
{} = pt_min_pdg ! pt cut for other particles (use pdg code). Applied on particle and anti-particle
{} = pt_max_pdg ! pt cut for other particles (syntax e.g. {6: 100, 25: 50})
#
# For display option for energy cut in the partonic center of mass frame type 'update ecut'
##*****
# Maximum and minimum absolute rapidity (for max, -1 means no cut) *
#******
-1.0 = etaj ! max rap for the jets
-1.0 = etab ! max rap for the b
-1.0 = etaa ! max rap for the photons
-1.0 = etal ! max rap for the charged leptons
0.0 = etajmin ! min rap for the jets
0.0 = etabmin ! min rap for the b
0.0 = etaamin ! min rap for the photons
0.0 = etalmin ! main rap for the charged leptons
{} = eta_min_pdg ! rap cut for other particles (use pdg code). Applied on particle and anti-particle
{} = eta_max_pdg ! rap cut for other particles (syntax e.g. {6: 2.5, 23: 5})
#******
# Minimum and maximum DeltaR distance *
#******
0.0 = drjj ! min distance between jets
0.0 = drbb ! min distance between b's
0.0 = drll ! min distance between leptons
0.0 = draa ! min distance between gammas
0.0 = drbj ! min distance between b and jet
0.0 = 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.0 = 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 *
#IF(llll)# # WARNING: for four lepton final state mmll cut require to have *
#IF(llll)# # 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
{} = mxx_min_pdg ! min invariant mass of a pair of particles X/X~ (e.g. {6:250})
{'default': False} = mxx_only_
#IF(LL)# #******
#IF(LL)# # Minimum and maximum invariant mass for all letpons *
#IF(LL)# #******
0.0 = mmnl ! min invariant mass for all letpons (l+- and vl)
-1.0 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
#IF(LL)# #******
#IF(LL)# # Minimum and maximum pt for 4-momenta sum of leptons / neutrino *
#IF(LL)# # for pair of lepton includes only same flavor, opposite charge
#IF(LL)# #******
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 = ptheavy ! minimum pt for at least one heavy final state
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
#IF(jj)# #******
#IF(jj)# # Control the pt's of the jets sorted by pt *
#IF(jj)# #******
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
#IF(ll)# #******
#IF(ll)# # Control the pt's of leptons sorted by pt *
#IF(ll)# #******
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
#IF(JJ)# #******
#IF(JJ)# # Control the Ht(k)=Sum of k leading jets *
#IF(JJ)# #******
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
#IF(aj)# #******
#IF(aj)# # Photon-isolation cuts, according to hep-ph/9801442 *
#IF(aj)# # When ptgmin=0, all the other parameters are ignored *
#IF(aj)# # When ptgmin>0, pta and draj are not going to be used *
#IF(aj)# #******
#IF(aj)# %(ptgmin)s = ptgmin ! Min photon transverse momentum
#IF(aj)# %(r0gamma)s = R0gamma ! Radius of isolation code
#IF(aj)# %(xn)s = xn ! n parameter of eq.(3.4) in hep-ph/9801442
#IF(aj)# %(epsgamma)s = epsgamma ! epsilon_gamma parameter of eq.(3.4) in hep-ph/9801442
#IF(aj)# %(isoem)s = isoEM ! isolate photons from EM energy (photons and leptons)
#IF(jj)# #******
#IF(jj)# # WBF cuts *
#IF(jj)# #******
0.0 = xetamin ! minimum rapidity for two jets in the WBF case
-1.0 = deltaeta ! minimum rapidity for two jets in the WBF case
#******
# maximal pdg code for quark to be considered as a light jet *
# (otherwise b cuts are applied) *
#******
4 = maxjetflavor ! Maximum jet pdg code
#******
#
#******
# Store info for systematics studies *
# WARNING: Do not use for interference type of computation *
#******
False = use_syst ! Enable systematics studies
#
systematics = systematics_program ! none, systematics [python], SysCalc [depreceted, C++]
['--mur=0.5,1,2', '--muf=0.5,1,2', '--pdf=errorset'] = systematics_
# Syscalc is deprecated but to see the associate options type'update syscalc'
# Additional hidden parameters
#******
1.0 = xn # hidden_parameter
True = isoem # hidden_parameter
1.0 = epsgamma # hidden_parameter
0.4 = r0gamma # hidden_parameter
0.0 = ptgmin # hidden_parameter
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- Luca Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Giacomo for more information if necessary.