e+ e- > j j [QCD QED] with eMELA bug with quarkphreco and lepphreco parameters
Dear mg5 team,
I tried to run NLO EW corrections of the process "e+ e- > j j" following this guide https:/
In mg5, run:
import model loop_qcd_qed_sm
generate e+ e- > j j [QCD QED] aW=10 aEW=10
output ee_jj
exit mg5, the change the following parameter in FKS_params.dat
#IRPoleChec
1.0d-5
-> #IRPoleCheckThr
-1d0
Launh mg5 again, then run:
launh ee_jj
0
set pdlabel emela
set lhaid 137002
set lpp1 -3
set lpp2 3
set ebeam 30
0
The process runs fine. By default, I have the following parameters:
quarkphreco True
Then I change the following parameter in the 'run_card.dat':
set quarkphreco False
The program runs fine, then I change the following parameter in the 'run_card.dat':
set quarkphreco True
Now, the program breaks entirely, I have tried to set everything back to default, it still broke.
The program only works when I "set quarkphreco False", but this is not the default.
And the numbers is not the same as the first time I ran the program with "quarkphreco True". I might have changed the "lepphreco" parameter too. Flipping the two parameters between "True" and "False"
Here's the error message:
INFO: Compiling source...
INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 16 cores
INFO: Compiling P0_epem_uux...
INFO: Compiling P0_epem_ddx...
INFO: Compiling P0_epem_bbx...
INFO: Compiling P0_epem_aa...
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
STOP 1
INFO: P0_epem_ddx done.
INFO: P0_epem_bbx done.
INFO: P0_epem_uux done.
INFO: P0_epem_aa done.
INFO: Checking test output:
INFO: P0_epem_uux
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: P0_epem_ddx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: P0_epem_bbx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: P0_epem_aa
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: Starting run
INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
INFO: Setting up grids
WARNING: program /afs/cern.
INFO: Idle: 0, Running: 0, Completed: 4 [ current time: 12h06 ]
Command "launch auto " interrupted with error:
Exception : program /afs/cern.
Please report this bug on https:/
More information is found in '/afs/cern.
Please attach this file to your report.
I guess, there is something to do with the "P0_epem_aa", previously, when the program was running fine, I had this message when the program tried to run "P0_epem_aa":
INFO: P0_epem_bbx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: P0_epem_aa
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for check_poles:
INFO: Poles successfully cancel for 20 points over 20 (tolerance=
INFO: Starting run
INFO: Using 16 cores
INFO: Cleaning previous results
INFO: Doing fixed order NLO
WARNING: No integration channels found for contribution P0_epem_aa
INFO: Setting up grids
INFO: Idle: 0, Running: 3, Completed: 0 [ current time: 10h54 ]
INFO: Idle: 0, Running: 2, Completed: 1 [ 3.7s ]
INFO: Idle: 0, Running: 1, Completed: 2 [ 3.9s ]
INFO: Idle: 0, Running: 0, Completed: 3 [ 4s ]
sum of cpu time of last step: 0 second
INFO:
Results after grid setup:
Total cross section: 1.267e+02 +- 1.2e+00 pb
INFO: Refining results, step 1
INFO: Idle: 0, Running: 3, Completed: 0 [ current time: 10h54 ]
INFO: Idle: 0, Running: 2, Completed: 1 [ 5.9s ]
INFO: Idle: 0, Running: 1, Completed: 2 [ 6.3s ]
INFO: Idle: 0, Running: 0, Completed: 3 [ 6.4s ]
sum of cpu time of last step: 0 second
INFO:
---
Final results and run summary:
Process e+ e- > j j [QCD QED] aS=10 aEW=10
Run at e+-e- collider (30.0 + 30.0 GeV)
Total cross section: 1.279e+02 +- 6.8e-01 pb
---
Scale variation (computed from histogram information):
Here is the content in the all_G1/log.txt . Thank you in advance.
Best, Dung.
LHAPDF 6.5.4 loading ../../.
LL_ALPMZ PDF set, member #0, version 1
Reading INFO from grid
2508966
===
INFO: MadFKS read these parameters from FKS_params.dat
===
> IRPoleCheckThre
> PrecisionVirtua
> SelectedContrib
> VetoedContribut
> QCD_squared_
> QED_squared_
> SelectedCouplin
> NHelForMCoverHels = 4
> VirtualFraction = 1.0000000000000000
> MinVirtualFraction = 5.0000000000000
> SeparateFlavour
> UsePolyVirtual = F
===
SPLIT TYPE USED: F T
A PDF is used, so alpha_s(MZ) is going to be modified
Old value of alpha_s from param_card: 0.11900000000000002
======= Evolution parameters =======
Perturbative order: LO
Factorisation scheme: Delta
Renormalisation scheme: alpha(MZ)
method: magnus
Evolution scheme: VFNS
NLMAX = 3, NUMAX = 2, NDMAX = 3
Mass thresholds:
me = 0.000510999 GeV
mu = 0.002160000 GeV
md = 0.004670000 GeV
ms = 0.093000000 GeV
mm = 0.105660000 GeV
mc = 1.270000000 GeV
mt = 1.776860000 GeV
mb = 4.180000000 GeV
MW = 80.379000000 GeV
MZ = 91.187600000 GeV
Alpha FIXED with reference value:
Alpha(Qref) = 0.007755544
W effects included in evolution
New value of alpha_s from PDF emela : 0.11900000000000002
WARNING: the value of maxjetflavorspe
Hence it will be set to: 5
***
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
***
External Params
---
MU_R = 1000.0000000000000
aEWM1 = 128.93000000000001
mdl_MW = 80.418999999999997
aS = 0.11900000000000002
mdl_ymt = 173.30000000000001
mdl_MZ = 91.188000000000002
mdl_MT = 173.30000000000001
mdl_MH = 125.00000000000000
mdl_WZ = 2.4888300000000001
mdl_WW = 2.0846499999999999
mdl_WT = 1.3672800000000001
mdl_WH = 4.1871800000000
Internal Params
---
mdl_lhv = 1.0000000000000000
mdl_SCKM22 = 1.0000000000000000
mdl_epsIR = 1.0000000000000000
mdl_epsUV = 1.0000000000000000
mdl_CKM33 = 1.0000000000000000
mdl_MC__exp__2 = 0.0000000000000000
mdl_MS__exp__2 = 0.0000000000000000
mdl_LogSwitch = 1.0000000000000000
mdl_MB__exp__4 = 0.0000000000000000
mdl_Ql__exp__4 = 1.0000000000000000
mdl_vep = 0.0000000000000000
mdl_CMSParam = 0.0000000000000000
mdl_Me__exp__2 = 0.0000000000000000
mdl_
mdl_Ql__exp__2 = 1.0000000000000000
mdl_CKM22 = 1.0000000000000000
mdl_yb = 0.0000000000000000
mdl_MU__exp__2 = 0.0000000000000000
mdl_MB__exp__2 = 0.0000000000000000
mdl_
mdl_
mdl_
mdl_
mdl_SCKM11 = 1.0000000000000000
mdl_SCKM33 = 1.0000000000000000
mdl_ep = 1.0000000000000
mdl_Ncol = 3.0000000000000000
mdl_CA = 3.0000000000000000
mdl_TF = 0.50000000000000000
mdl_CF = 1.3333333333333333
mdl_complexi = (0.000000000000
mdl_MZ__exp__2 = 8315.2513440000002
mdl_MW__exp__2 = 6467.2155609999991
mdl_sqrt__2 = 1.4142135623730951
mdl_sw2 = 0.22224653309289089
mdl_cw = 0.88190332061236121
mdl_sqrt__sw2 = 0.47143030565810135
mdl_sw = 0.47143030565810135
mdl_MH__exp__2 = 15625.000000000000
mdl_I3d = -0.500000000000
mdl_I3u = 0.50000000000000000
mdl_I3l = -0.500000000000
mdl_I3v = 0.50000000000000000
mdl_Qd = -0.333333333333
mdl_Qu = 0.66666666666666663
mdl_Ql = -1.0000000000000000
mdl_
mdl_sw__exp__2 = 0.22224653309289086
mdl_MT__exp__2 = 30032.890000000003
mdl_cw__exp__2 = 0.77775346690710911
mdl_MT__exp__4 = 901974481.75210023
mdl_sw__exp__4 = 4.9393521471809
mdl_MW__exp__4 = 41824877.112440534
mdl_MT__exp__3 = 5204699.8370000012
mdl_MH__exp__4 = 244140625.00000000
mdl_MZ__exp__4 = 69143404.913893804
mdl_MT__exp__6 = 27088900393267.840
mdl_MW__exp__6 = 270490496098.48813
mdl_MT__exp__8 = 8.1355796573196
mdl_MW__exp__8 = 1749320345470752.0
mdl_MT__exp__10 = 2.4433496893452
mdl_MW__exp__10 = 1.1313231759402
mdl_MZ__exp__6 = 574944790638.99182
mdl_sw__exp__3 = 0.10477375102743487
mdl_MZ__exp__3 = 758251.13955667207
mdl_cw__exp__4 = 0.60490045528602765
mdl_cw__exp__6 = 0.47046342623239684
mdl_MH__exp__6 = 3814697265625.0000
mdl_cw__exp__3 = 0.68590336508315575
mdl_Qd__exp__3 = -3.703703703703
mdl_Qd__exp__2 = 0.11111111111111110
mdl_I3d__exp__2 = 0.25000000000000000
mdl_I3d__exp__3 = -0.125000000000
mdl_Qu__exp__3 = 0.29629629629629622
mdl_Qu__exp__2 = 0.44444444444444442
mdl_I3u__exp__2 = 0.25000000000000000
mdl_I3u__exp__3 = 0.12500000000000000
mdl_Ql__exp__3 = -1.0000000000000000
mdl_I3l__exp__2 = 0.25000000000000000
mdl_I3l__exp__3 = -0.125000000000
mdl_I3v__exp__3 = 0.12500000000000000
mdl_MW__exp__3 = 520087.00820005889
mdl_sw__exp__6 = 1.0977538904358
mdl_sw__exp__8 = 2.4397199633860
mdl_Qu__exp__4 = 0.19753086419753083
mdl_Qd__exp__4 = 1.2345679012345
mdl_I3l__exp__4 = 6.2500000000000
mdl_I3u__exp__4 = 6.2500000000000
mdl_I3d__exp__4 = 6.2500000000000
mdl_I3v__exp__4 = 6.2500000000000
mdl_I3v__exp__2 = 0.25000000000000000
mdl_aEW = 7.7561467462964
mdl_Gf = (1.198749835046
mdl_sqrt__aEW = 8.8068988561788
mdl_ee = 0.31219643584339185
mdl_g1 = 0.35400301659666522
mdl_gw = 0.66223242777652103
mdl_vev = 242.87243157213206
mdl_vev__exp__2 = 58987.018017759969
mdl_lam = 0.13244439645428069
mdl_yt = 1.0091026337275697
mdl_muH = 88.388347648318430
mdl_I233 = (1.009102633727
mdl_I333 = (1.009102633727
mdl_AxialZUp = -0.187728181847
mdl_AxialZDown = 0.18772818184783885
mdl_VectorZUp = 7.6469681635801
mdl_VectorZDown = -0.132098931741
mdl_AxialG0Up3 = -0.713543315221
mdl_VectorHUp3 = (-0.00000000000
mdl_VectorAUp = 0.20813095722892791
mdl_VectorADown = -0.104065478614
mdl_VectorWmDxU = 0.23413452020120429
mdl_AxialWmDxU = -0.234134520201
mdl_VectorWpUxD = 0.23413452020120429
mdl_AxialWpUxD = -0.234134520201
mdl_
mdl_
mdl_
mdl_
mdl_ee__exp__2 = 9.7466614553317
mdl_R2VV_FIN_ = 9.8754327521535
mdl_ee__exp__3 = 3.0428729677267
mdl_R2SFF_FIN_ = 3.0830749076334
mdl_ee__exp__4 = 9.4997409524848
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
Internal Params evaluated point by point
---
mdl_sqrt__aS = 0.34496376621320685
mdl_G__exp__4 = 2.2362154867812229
mdl_G__exp__2 = 1.4953981031087418
mdl_
mdl_G__exp__3 = 1.8286695702272273
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
mdl_
Couplings of loop_qcd_qed_sm
---
GC_11 = 0.00000E+00 0.12229E+01
R2_DDA = 0.00000E+00 0.26279E-02
R2_UUA = -0.00000E+00 -0.52558E-02
R2_DDZ_V2 = 0.00000E+00 0.73741E-02
R2_DDZ_V3 = 0.00000E+00 0.70239E-03
c_UVZMass1EW = -0.00000E+00 -0.10195E-01
c_UVZMass2EW = -0.00000E+00 -0.46990E+03
c_UVAMass1EW = 0.00000E+00 0.58992E-02
c_UVAZMass1EW = 0.00000E+00 0.66190E-02
c_UVAZMass2EW = 0.00000E+00 -0.96797E+02
c_UVemepMass1EW -0.00000E+00 0.74872E-02
c_UVemepMass2EW 0.00000E+00 -0.75658E-03
c_UVepemA1EW = 0.00000E+00 0.11789E-02
c_UVepemA2EW = 0.00000E+00 -0.52931E-02
c_UVcxcA1EW = -0.00000E+00 -0.23915E-02
c_UVcxcA2EW = 0.00000E+00 0.36162E-02
c_UVdxdA1EW = 0.00000E+00 0.34729E-02
c_UVbxbA1EW = 0.00000E+00 0.47466E-02
c_UVbxbA2EW = 0.00000E+00 -0.18343E-02
c_UVepemZ1EW = 0.00000E+00 0.18389E-01
c_UVepemZ2EW = 0.00000E+00 0.50798E-02
c_UVcxcZ1EW = 0.00000E+00 -0.17259E-01
c_UVcxcZ2EW = 0.00000E+00 -0.34333E-02
c_UVdxdZ1EW = 0.00000E+00 0.16262E-01
c_UVbxbZ1EW = 0.00000E+00 0.20176E-01
c_UVbxbZ2EW = 0.00000E+00 0.17307E-02
GC_1 = -0.00000E+00 -0.10407E+00
GC_2 = 0.00000E+00 0.20813E+00
GC_3 = -0.00000E+00 -0.31220E+00
GC_5 = 0.00000E+00 0.97467E-01
GC_6 = 0.00000E+00 0.19493E+00
GC_30 = -0.10091E+01 -0.00000E+00
GC_54 = 0.00000E+00 0.34109E+00
GC_68 = -0.00000E+00 -0.29201E+00
GC_70 = -0.00000E+00 -0.58402E+00
GC_75 = -0.00000E+00 -0.36466E+00
GC_76 = -0.00000E+00 -0.27815E-01
GC_77 = 0.00000E+00 0.83444E-01
GC_78 = -0.37546E+00 0.00000E+00
GC_79 = 0.00000E+00 -0.20857E+00
GC_80 = 0.00000E+00 0.37546E+00
GC_81 = 0.00000E+00 0.13023E+00
GC_82 = 0.00000E+00 0.87002E-01
GC_83 = 0.00000E+00 0.28193E+00
GC_84 = -0.13421E+02 0.00000E+00
GC_92 = -0.25107E+02 0.00000E+00
GC_99 = 0.00000E+00 0.68474E+02
GC_124 = 0.00000E+00 0.46827E+00
R2_AAboson1 = 0.00000E+00 0.61721E-03
R2_AAboson2 = -0.00000E+00 -0.41148E-03
R2_AAboson3 = -0.00000E+00 -0.79833E+01
R2_AAU = 0.00000E+00 0.54864E-03
R2_AAD = 0.00000E+00 0.13716E-03
R2_AAt3 = -0.00000E+00 -0.98863E+02
R2_AZboson1 = 0.00000E+00 0.11546E-02
R2_AZboson2 = -0.00000E+00 -0.76975E-03
R2_AZboson3 = -0.00000E+00 -0.14934E+02
R2_AZl = 0.00000E+00 0.27468E-04
R2_AZU = 0.00000E+00 0.20157E-03
R2_AZD = 0.00000E+00 0.17411E-03
R2_AZt3 = -0.00000E+00 -0.36323E+02
R2_ZZboson1 = 0.00000E+00 0.21599E-02
R2_ZZboson2 = -0.00000E+00 -0.14400E-02
R2_ZZboson3 = -0.00000E+00 -0.27938E+02
R2_ZZl = 0.00000E+00 0.15061E-03
R2_ZZU = 0.00000E+00 0.52040E-03
R2_ZZD = 0.00000E+00 0.66735E-03
R2_ZZt3 = -0.00000E+00 -0.93775E+02
R2_ZZv = 0.00000E+00 0.29756E-03
R2_LLCm = 0.00000E+00 0.79359E-03
R2_LLCp0 = 0.00000E+00 0.89269E-03
R2_QQCpcs = 0.00000E+00 0.13886E-02
R2_bbA2Cp = 0.00000E+00 0.14909E-03
R2_bbA2Cp_t = 0.00000E+00 -0.38205E-03
R2_bbA2Cm = 0.00000E+00 0.18352E-04
R2_ccA2Cp = 0.00000E+00 -0.29818E-03
R2_ccA2Cp_s = 0.00000E+00 -0.57801E-03
R2_ccA2Cm = 0.00000E+00 -0.14682E-03
R2_eeACp = 0.00000E+00 0.55739E-03
R2_eeACm = 0.00000E+00 0.49551E-03
R2_ddZ2Cp_u = -0.00000E+00 0.88821E-03
R2_bbZ2Cp = 0.00000E+00 0.45821E-03
R2_bbZ2Cp_t = -0.00000E+00 0.12469E-02
R2_bbZ2Cm = 0.00000E+00 -0.98104E-05
R2_ccZ2Cp = 0.00000E+00 -0.37851E-03
R2_ccZ2Cp_s = -0.00000E+00 -0.73372E-03
R2_ccZ2Cm = 0.00000E+00 0.78483E-04
R2_eeZCp = 0.00000E+00 0.37237E-03
R2_llZCp = -0.00000E+00 0.57923E-03
R2_eeZCm = 0.00000E+00 -0.26488E-03
R2_AWWlv = 0.00000E+00 -0.28901E-03
R2_AWWcs = 0.00000E+00 -0.86702E-03
c_UVZMass1EW_1e -0.00000E+00 -0.30052E-02
c_UVZMass2EW_1e -0.00000E+00 -0.15682E+03
c_UVAMass1EW_1e 0.00000E+00 0.75437E-03
c_UVAZMass1EW_1 0.00000E+00 0.12924E-02
c_UVAZMass2EW_1 0.00000E+00 -0.19202E+02
c_UVemepMass1EW 0.00000E+00 0.16641E-02
c_UVemepMass2EW -0.00000E+00 -0.17637E-03
c_UVepemA1EW_1e 0.00000E+00 0.53012E-03
c_UVepemA2EW_1e 0.00000E+00 -0.80135E-03
c_UVcxcA1EW_1ep -0.00000E+00 -0.67708E-03
c_UVcxcA2EW_1ep 0.00000E+00 0.55462E-03
c_UVdxdA1EW_1ep 0.00000E+00 0.79346E-03
c_UVbxbA1EW_1ep 0.00000E+00 0.11290E-02
c_UVbxbA2EW_1ep 0.00000E+00 -0.28343E-03
c_UVepemZ1EW_1e 0.00000E+00 0.19690E-02
c_UVepemZ2EW_1e -0.00000E+00 -0.29434E-04
c_UVcxcZ1EW_1ep 0.00000E+00 -0.21056E-02
c_UVcxcZ2EW_1ep 0.00000E+00 0.87213E-05
c_UVdxdZ1EW_1ep 0.00000E+00 0.22732E-02
c_UVbxbZ1EW_1ep 0.00000E+00 0.33044E-02
c_UVbxbZ2EW_1ep -0.00000E+00 -0.10902E-05
Collider parameters:
---
Running at e+e- machine @ 60.000000000000000 GeV
PDF set = emela
alpha_s(Mz)= 0.1190 running at 2 loops.
alpha_s(Mz)= 0.1190 running at 2 loops.
Renormalization scale set on event-by-event basis
Factorization scale set on event-by-event basis
INFO: orders_tag_plot is computed as: + QCD * 1 + QED * 100
orders_
AMP_SPLIT: 1 correspond to S.O. 0 4
orders_
AMP_SPLIT: 2 correspond to S.O. 0 6
getting user params
Number of phase-space points per iteration: -1
Maximum number of iterations is: 6
Desired accuracy is: 5.0000000000000
Using adaptive grids: 2
Using Multi-channel integration
Do MC over helicities for the virtuals
Number of channels to integrate together: 4
Running Configuration Number(s): 1 1 2 2
initial-
Splitting channel: 0
Weight multiplier: 1.0000000000000000
doing the all of this channel
Normal integration (Sfunction != 1)
RESTART: Fresh run
about to integrate 7 -1 6
imode is 0
channel 1 : 1 T 0 0 0.1000E+01 0.0000E+00 0.1000E+01
channel 2 : 1 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
channel 3 : 2 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
channel 4 : 2 T 0 0 0.0000E+00 0.0000E+00 0.1000E+01
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG
STOP 1
Thanks for using LHAPDF 6.5.4. Please make sure to cite the paper:
Eur.Phys.J. C75 (2015) 3, 132 (http://
------- iteration 1
Update # PS points (even_rn): 280 --> 256
Using random seed offsets: 0 , 4 , 0
with seed 42
Ranmar initialization seeds 22919 9420
initial-final FKS maps:
0 : 8 1 2 3 4 5 6 7 8
1 : 4 5 6 7 8 0 0 0 0
2 : 4 1 2 3 4 0 0 0 0
Total number of FKS directories is 8
For the Born we use nFKSprocesses:
1 2 1 2 2 2 2 2
tau_min 1 1 : 0.00000E+00 -- 0.20000E+02
tau_min 2 1 : 0.00000E+00 -- 0.20000E+02
tau_min 3 1 : 0.00000E+00 -- 0.20000E+02
tau_min 4 1 : 0.00000E+00 -- 0.20000E+02
tau_min 5 1 : 0.20000E+02 0.20000E+02 0.20000E+02
tau_min 6 1 : 0.20000E+02 0.20000E+02 0.20000E+02
tau_min 7 1 : 0.20000E+02 0.20000E+02 0.20000E+02
tau_min 8 1 : 0.20000E+02 0.20000E+02 0.20000E+02
ERROR, too many photons 2
Time in seconds: 0
Question information
- Language:
- English Edit question
- Status:
- Open
- Assignee:
- marco zaro Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Le Van Dung for more information if necessary.