Error when generating e+ e- > h2 in a 2HDM
Hello,
I am trying to generate e+ e- > h2 in a 2HDM but I am running into the error
"Command "generate_events run_01" interrupted with error:
Exception : Reported error: End code 2.0 "
In the run_01_
"At line 199 of file ranmar.f
Fortran runtime error: Index '0' of dimension 1 of array 'ranu' below lower bound of 1"
I'm not sure what the issue here is. Do you have any insight/suggestions to fix this?
The rest of the debug log is below.
Thanks!
Exception: Reported error: End code 2.0
Full associated log:
Process in group number 1
No PDF is used, alpha_s(MZ) from param_card is used
Value of alpha_s from param_card: 0.11799999999999999
The default order of alpha_s running is fixed to 2
Define smin to 16900.000000000000
******
* MadGraph/MadEvent *
* -------
* http://
* http://
* http://
* -------
* *
* PARAMETER AND COUPLING VALUES *
* *
******
External Params
-----
mdl_cabi = 0.22773599999999999
mdl_l2 = 0.50000000000000000
mdl_l3 = 1.0000000000000000
mdl_lR7 = 0.10000000000000001
mdl_lI7 = 0.20000000000000001
mdl_mixh = 0.29999999999999999
mdl_mixh2 = 0.10000000000000001
mdl_mixh3 = 0.20000000000000001
aEWM1 = 127.90000000000001
mdl_Gf = 1.1663900000000
aS = 0.11799999999999999
mdl_ymb = 4.7000000000000002
mdl_ymt = 172.00000000000000
mdl_ymtau = 1.7769999999999999
mdl_GDI1x1 = 0.0000000000000000
mdl_GDI1x2 = 0.0000000000000000
mdl_GDI1x3 = 0.0000000000000000
mdl_GDI2x1 = 0.0000000000000000
mdl_GDI2x2 = 0.0000000000000000
mdl_GDI2x3 = 0.0000000000000000
mdl_GDI3x1 = 0.0000000000000000
mdl_GDI3x2 = 0.0000000000000000
mdl_GDI3x3 = 0.0000000000000000
mdl_GDR1x1 = 1.0000000000000000
mdl_GDR1x2 = 1.0000000000000
mdl_GDR1x3 = 1.0000000000000
mdl_GDR2x1 = 1.0000000000000
mdl_GDR2x2 = 1.0000000000000000
mdl_GDR2x3 = 1.0000000000000
mdl_GDR3x1 = 1.0000000000000
mdl_GDR3x2 = 1.0000000000000
mdl_GDR3x3 = 1.0000000000000000
mdl_GLI1x1 = 0.0000000000000000
mdl_GLI1x2 = 0.0000000000000000
mdl_GLI1x3 = 0.0000000000000000
mdl_GLI2x1 = 0.0000000000000000
mdl_GLI2x2 = 0.0000000000000000
mdl_GLI2x3 = 0.0000000000000000
mdl_GLI3x1 = 0.0000000000000000
mdl_GLI3x2 = 0.0000000000000000
mdl_GLI3x3 = 0.0000000000000000
mdl_GLR1x1 = 1.0000000000000000
mdl_GLR1x2 = 1.0000000000000
mdl_GLR1x3 = 1.0000000000000
mdl_GLR2x1 = 1.0000000000000
mdl_GLR2x2 = 1.0000000000000000
mdl_GLR2x3 = 1.0000000000000
mdl_GLR3x1 = 1.0000000000000
mdl_GLR3x2 = 1.0000000000000
mdl_GLR3x3 = 1.0000000000000000
mdl_GUI1x1 = 0.0000000000000000
mdl_GUI1x2 = 0.0000000000000000
mdl_GUI1x3 = 0.0000000000000000
mdl_GUI2x1 = 0.0000000000000000
mdl_GUI2x2 = 0.0000000000000000
mdl_GUI2x3 = 0.0000000000000000
mdl_GUI3x1 = 0.0000000000000000
mdl_GUI3x2 = 0.0000000000000000
mdl_GUI3x3 = 0.0000000000000000
mdl_GUR1x1 = 1.0000000000000000
mdl_GUR1x2 = 1.0000000000000
mdl_GUR1x3 = 1.0000000000000
mdl_GUR2x1 = 1.0000000000000
mdl_GUR2x2 = 1.0000000000000000
mdl_GUR2x3 = 1.0000000000000
mdl_GUR3x1 = 1.0000000000000
mdl_GUR3x2 = 1.0000000000000
mdl_GUR3x3 = 1.0000000000000000
mdl_MZ = 91.187600000000003
mdl_MT = 172.00000000000000
mdl_MB = 4.7000000000000002
mdl_mhc = 150.00000000000000
mdl_mh1 = 120.00000000000000
mdl_mh2 = 130.00000000000000
mdl_mh3 = 140.00000000000000
mdl_WZ = 2.4952000000000001
mdl_WW = 2.0850000000000000
mdl_WT = 1.5083359999999999
mdl_whc = 1.0000000000000000
mdl_Wh1 = 1.0000000000000000
mdl_Wh2 = 1.0000000000000000
mdl_Wh3 = 1.0000000000000000
Internal Params
-----
mdl_cos__cabi = 0.97418004031982097
mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
mdl_sin__cabi = 0.22577256042856930
mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
mdl_CKM1x3 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_CKM2x1 = (-0.22577256042
mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
mdl_CKM2x3 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_CKM3x1 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_CKM3x2 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_CKM3x3 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_cos__mixh = 0.95533648912560598
mdl_cos__mixh2 = 0.99500416527802571
mdl_TH1x1 = 0.95056378592206325
mdl_sin__mixh = 0.29552020666133955
mdl_TH1x2 = 0.29404383655185579
mdl_sin__mixh2 = 9.9833416646828
mdl_TH1x3 = -9.983341664682
mdl_cos__mixh3 = 0.98006657784124163
mdl_sin__mixh3 = 0.19866933079506122
mdl_TH2x1 = -0.270681488391
mdl_TH2x2 = 0.94215466351136834
mdl_TH2x3 = 0.19767681165408385
mdl_TH3x1 = 0.15218416716418803
mdl_TH3x2 = -0.160881360665
mdl_TH3x3 = 0.97517032720181585
mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
mdl_GD1x1 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GD1x2 = ( 1.0000000000000
mdl_GD1x3 = ( 1.0000000000000
mdl_GD2x1 = ( 1.0000000000000
mdl_GD2x2 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GD2x3 = ( 1.0000000000000
mdl_GD3x1 = ( 1.0000000000000
mdl_GD3x2 = ( 1.0000000000000
mdl_GD3x3 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GL1x1 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GL1x2 = ( 1.0000000000000
mdl_GL1x3 = ( 1.0000000000000
mdl_GL2x1 = ( 1.0000000000000
mdl_GL2x2 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GL2x3 = ( 1.0000000000000
mdl_GL3x1 = ( 1.0000000000000
mdl_GL3x2 = ( 1.0000000000000
mdl_GL3x3 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GU1x1 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GU1x2 = ( 1.0000000000000
mdl_GU1x3 = ( 1.0000000000000
mdl_GU2x1 = ( 1.0000000000000
mdl_GU2x2 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_GU2x3 = ( 1.0000000000000
mdl_GU3x1 = ( 1.0000000000000
mdl_GU3x2 = ( 1.0000000000000
mdl_GU3x3 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_l7 = ( 0.10000000000000001 , 0.20000000000000001 )
mdl_MZ__exp__2 = 8315.1783937600012
mdl_MZ__exp__4 = 69142191.720053151
mdl_sqrt__2 = 1.4142135623730951
mdl_mhc__exp__2 = 22500.000000000000
mdl_mh1__exp__2 = 14400.000000000000
mdl_mh2__exp__2 = 16900.000000000000
mdl_mh3__exp__2 = 19600.000000000000
mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
mdl_conjg__GD1x1 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
mdl_conjg__GD2x1 = ( 1.0000000000000
mdl_conjg__CKM1x3 = ( 0.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GD3x1 = ( 1.0000000000000
mdl_I1a11 = ( 0.97643776592410669 , -0.0000000000000000 )
mdl_conjg__CKM2x1 = (-0.22577256042
mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
mdl_conjg__CKM2x3 = ( 0.0000000000000000 , -0.0000000000000000 )
mdl_I1a12 = (-0.21603076002
mdl_conjg__CKM3x1 = ( 0.0000000000000000 , -0.0000000000000000 )
mdl_conjg__CKM3x2 = ( 0.0000000000000000 , -0.0000000000000000 )
mdl_conjg__CKM3x3 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_I1a13 = ( 1.0000000000000
mdl_conjg__GD1x2 = ( 1.0000000000000
mdl_conjg__GD2x2 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GD3x2 = ( 1.0000000000000
mdl_I1a21 = ( 0.23551436083176752 , -0.0000000000000000 )
mdl_I1a22 = ( 0.97192231471553525 , 0.0000000000000000 )
mdl_I1a23 = ( 1.0000000000000
mdl_conjg__GD1x3 = ( 1.0000000000000
mdl_conjg__GD2x3 = ( 1.0000000000000
mdl_conjg__GD3x3 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_I1a31 = ( 1.1999526007483
mdl_I1a32 = ( 7.4840747989125
mdl_I1a33 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_I2a11 = ( 0.97192231471553525 , 0.0000000000000000 )
mdl_I2a12 = (-0.21603076002
mdl_I2a13 = ( 7.4840747989125
mdl_I2a21 = ( 0.23551436083176752 , 0.0000000000000000 )
mdl_I2a22 = ( 0.97643776592410669 , 0.0000000000000000 )
mdl_I2a23 = ( 1.1999526007483
mdl_I2a31 = ( 1.0000000000000
mdl_I2a32 = ( 1.0000000000000
mdl_I2a33 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_conjg__GU1x1 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GU2x1 = ( 1.0000000000000
mdl_conjg__GU3x1 = ( 1.0000000000000
mdl_I3a11 = ( 0.97192231471553525 , 0.0000000000000000 )
mdl_I3a12 = ( 0.23551436083176752 , 0.0000000000000000 )
mdl_I3a13 = ( 1.0000000000000
mdl_conjg__GU1x2 = ( 1.0000000000000
mdl_conjg__GU2x2 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GU3x2 = ( 1.0000000000000
mdl_I3a21 = (-0.21603076002
mdl_I3a22 = ( 0.97643776592410669 , 0.0000000000000000 )
mdl_I3a23 = ( 1.0000000000000
mdl_conjg__GU1x3 = ( 1.0000000000000
mdl_conjg__GU2x3 = ( 1.0000000000000
mdl_conjg__GU3x3 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_I3a31 = ( 7.4840747989125
mdl_I3a32 = ( 1.1999526007483
mdl_I3a33 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_I4a11 = ( 0.97643776592410669 , 0.0000000000000000 )
mdl_I4a12 = ( 0.23551436083176752 , 0.0000000000000000 )
mdl_I4a13 = ( 1.1999526007483
mdl_I4a21 = (-0.21603076002
mdl_I4a22 = ( 0.97192231471553525 , 0.0000000000000000 )
mdl_I4a23 = ( 7.4840747989125
mdl_I4a31 = ( 1.0000000000000
mdl_I4a32 = ( 1.0000000000000
mdl_I4a33 = ( 1.0000000000000000 , 0.0000000000000000 )
mdl_TH1x1__exp__2 = 0.90357151110648615
mdl_TH2x1__exp__2 = 7.3268468158056
mdl_TH3x1__exp__2 = 2.3160020735457
mdl_TH1x2__exp__2 = 8.6461777814134
mdl_TH2x2__exp__2 = 0.88765540997621972
mdl_TH3x2__exp__2 = 2.5882812209645
mdl_TH1x3__exp__2 = 9.9667110793791
mdl_TH2x3__exp__2 = 3.9076121865724
mdl_TH3x3__exp__2 = 0.95095716705489652
mdl_conjg__GL1x1 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GL1x2 = ( 1.0000000000000
mdl_conjg__GL1x3 = ( 1.0000000000000
mdl_conjg__GL2x1 = ( 1.0000000000000
mdl_conjg__GL2x2 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__GL2x3 = ( 1.0000000000000
mdl_conjg__GL3x1 = ( 1.0000000000000
mdl_conjg__GL3x2 = ( 1.0000000000000
mdl_conjg__GL3x3 = ( 1.0000000000000000 , -0.0000000000000000 )
mdl_conjg__l7 = ( 0.10000000000000001 ,-0.20000000000
mdl_TH2x1__exp__3 = -1.983241801321
mdl_TH3x1__exp__3 = 3.5245884671309
mdl_TH1x1__exp__4 = 0.81644147568325887
mdl_TH1x1__exp__3 = 0.85890235644870105
mdl_TH2x1__exp__4 = 5.3682684262280
mdl_TH3x1__exp__4 = 5.3638656046682
mdl_TH2x2__exp__3 = 0.83630868410019099
mdl_TH3x2__exp__3 = -4.164062046142
mdl_TH1x2__exp__3 = 2.5423552863562
mdl_TH1x2__exp__4 = 7.4756390227807
mdl_TH2x2__exp__4 = 0.78793212686005076
mdl_TH3x2__exp__4 = 6.6991996787978
mdl_TH2x3__exp__3 = 7.7244431822227
mdl_TH3x3__exp__3 = 0.92734521175183526
mdl_TH1x3__exp__3 = -9.950108197862
mdl_TH1x3__exp__4 = 9.9335329739819
mdl_TH2x3__exp__4 = 1.5269433000649
mdl_TH3x3__exp__4 = 0.90431953357307437
mdl_aEW = 7.8186082877247
mdl_MW = 79.824660036055974
mdl_sqrt__aEW = 8.8422894590285
mdl_ee = 0.31345100004952897
mdl_MW__exp__2 = 6371.9763498719121
mdl_sw2 = 0.23369336794341478
mdl_cw = 0.87538941737753784
mdl_sqrt__sw2 = 0.48341841911889827
mdl_sw = 0.48341841911889827
mdl_g1 = 0.35807035569216145
mdl_gw = 0.64840516548963911
mdl_vev = 246.21845810181625
mdl_vev__exp__2 = 60623.529110035844
mdl_mu2 = -7811.7645550179222
mdl_l1 = 0.12097597712865511
mdl_l4 = -0.144630005393
mdl_lI5 = -1.028410447757
mdl_lI6 = 1.0301433873426
mdl_lR5 = -2.133948632280
mdl_lR6 = 9.7316377085732
mdl_yb = 2.6995554250465
mdl_yt = 0.98792241086809907
mdl_ytau = 1.0206617000654
mdl_mu1 = -7333.9906710740534
mdl_l5 = ( -2.133948632280
mdl_l6 = ( 9.7316377085732
mdl_mu3 = ( -294.98311095700723 , -312.25463815040297 )
mdl_I5a31 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I5a32 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I5a33 = ( 2.6995554250465
mdl_I6a13 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I6a23 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I6a33 = ( 0.98792241086809907 , 0.0000000000000000 )
mdl_I7a31 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I7a32 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I7a33 = ( 0.98792241086809907 , 0.0000000000000000 )
mdl_I8a13 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I8a23 = ( 0.0000000000000000 , 0.0000000000000000 )
mdl_I8a33 = ( 2.6995554250465
mdl_gw__exp__2 = 0.42042925863364627
mdl_cw__exp__2 = 0.76630663205658511
mdl_sw__exp__2 = 0.23369336794341478
At line 199 of file ranmar.f
Fortran runtime error: Index '0' of dimension 1 of array 'ranu' below lower bound of 1
Error termination. Backtrace:
#0 0x10d7c3729
#1 0x10d7c43f5
#2 0x10d7c479a
#3 0x108856ab1
#4 0x10883e427
#5 0x108851d74
#6 0x10881faed
#7 0x108887523
mdl_g1__exp__2 = 0.12821437962551102
mdl_conjg__l5 = ( -2.133948632280
mdl_conjg__l6 = ( 9.7316377085732
Internal Params evaluated point by point
-----
mdl_sqrt__aS = 0.34351128074635334
mdl_G__exp__2 = 1.4828317324943823
Couplings of 2HDM
-----
GC_179 = -0.11376E+00 -0.66620E+00
GC_359 = 0.11376E+00 -0.66620E+00
Collider parameters:
------
Running at e e machine @ 1000.0000000000000 GeV
PDF set = none
alpha_s(Mz)= 0.1180 running at 2 loops.
alpha_s(Mz)= 0.1180 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 3
Completed mapping 3
about to integrate -1 1000 5 3 -1 1
Using non-zero grid deformation.
-1 dimensions 1000 events 2 invarients 5 iterations 1 config(s), (0.99)
Using h-tuple random number sequence.
Error opening grid
Using Uniform Grid! 12
Using uniform alpha 1.0000000000000000
Grid defined OK
Set CM energy to 1000.00
Mapping Graph 1 to config 1
1
Masses: 0.000E+00 0.000E+00 0.130E+03
******
* You are using the DiscreteSampler module *
* part of the MG5_aMC framework *
* Author: Valentin Hirschi *
******
Particle 3
Et > 0.0
E > 0.0
Eta < -1.0
xqcut: 0.0
alpha_s for scale 1000.0000000000000 is 8.6843671986092
ls status:
input_app.txt
run1_app.log
automatic_
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Douglas Tuckler for more information if necessary.