MadGraph5_aMC@NLO

No Integration channels found

Asked by Jack Y. Araz

Dear Oliver

I'm having a problem with the survey part of the simulation. I'm using DMsimp model to generate mono jet events at NLO but I realised that *sometimes* I get "WARNING: No integration channels found for contribution XXXX" and it crashes by giving the following error;

File "/mt/batch/jaraz/MG5_aMC_v2_7_3/madgraph/interface/amcatnlo_run_interface.py", line 2536, in update_jobs_to_run
    jobs[i]['nevents'] += 1
IndexError: list index out of range

However, this doesn't happen all the time, if I run the exact same sample again it can pass the integration check. Thus I think it's about the allowed iteration during the survey. Is there a way to increase the number of trials? Or is there anything else I can do?

Mg5 v2.7.3

Thanks
Best regards
Jack

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Solved by:
Jack Y. Araz
Solved:
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said : 		#1

Hi Jack,
some extra information would be needed.
I suggest you to run with a fixed random seed in the run_card (iseed), and see if you still get that error.
In that case, go into the P0 dir which has no channel, and run ./gensym NLO

What process are you considering?

Cheers,

Marco

Revision history for this message
Jack Y. Araz (jackaraz) said : 		#2

Hi Marco

Thanks for the answer. I'm running DMsimp model from feynrules (tried both the combined and separate packages under the dedicated page) and running 'p p > xd xd~ j [QCD]' with different restrict cards. For instance, for a spin-1-mediated process, the gV** are on, and rest is zero etc.

I tried with a fixed seed and got the same error (currently running for spin-1-mediated one) One of the integration channels that shows up is P0_gu_xdxdxu, so I went into that folder and ran './gensym NL O' as requested. The result is attached below.

Thanks
Best regards
—
Jack Y. Araz

———————————————
 A PDF is used, so alpha_s(MZ) is going to be modified
 Old value of alpha_s from param_card: 0.11839999999999999
==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
LHAPDF 6.2.3 loading /mt/batch/jaraz/MG5_aMC_v2_7_3/HEPTools/lhapdf6/share/LHAPDF/NNPDF30_nlo_as_0118/NNPDF30_nlo_as_0118_0000.dat
NNPDF30_nlo_as_0118 PDF set, member #0, version 2; LHAPDF ID = 260000
 New value of alpha_s from PDF lhapdf : 0.11800222660863767
WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5

  External Params
  ---------------------------------

 MU_R = 91.188000000000002
 mdl_cabi = 0.22773599999999999
 mdl_gVXd = 0.25000000000000000
 aEWM1 = 127.90000000000001
 mdl_Gf = 1.1663700000000000E-005
 aS = 0.11839999999999999
 mdl_ymt = 172.00000000000000
 mdl_ymtau = 1.7769999999999999
 mdl_MZ = 91.187600000000003
 mdl_MTA = 1.7769999999999999
 mdl_MT = 172.00000000000000
 mdl_MH = 125.00000000000000
 mdl_MXr = 100.00000000000000
 mdl_MY1 = 1500.0000000000000
 mdl_WZ = 2.4952000000000001
 mdl_WW = 2.0850000000000000
 mdl_WT = 1.5083359999999999
 mdl_WH = 4.0699999999999998E-003
 mdl_WY1 = 10.000000000000000
  Internal Params
  ---------------------------------

 mdl_MXd = 100.00000000000000
 mdl_MXc = 100.00000000000000
 mdl_gVu22 = 0.25000000000000000
 mdl_gVd22 = 0.25000000000000000
 mdl_gVu11 = 0.25000000000000000
 mdl_gVd11 = 0.25000000000000000
 mdl_gAu11 = 0.0000000000000000
 mdl_gAd22 = 0.0000000000000000
 mdl_gAu22__exp__2 = 0.0000000000000000
 mdl_cos__cabi = 0.97418004031982097
 mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
 mdl_sin__cabi = 0.22577256042856930
 mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
 mdl_CKM2x1 = (-0.22577256042856930 , 0.0000000000000000 )
 mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
 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.0000000000000000 , 1.0000000000000000 )
 mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
 mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
 mdl_conjg__CKM2x1 = (-0.22577256042856930 , -0.0000000000000000 )
 mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
 mdl_MT__exp__2 = 29584.000000000000
 mdl_gVu22__exp__2 = 6.2500000000000000E-002
 mdl_gVd11__exp__2 = 6.2500000000000000E-002
 mdl_gVd22__exp__2 = 6.2500000000000000E-002
 mdl_gVu11__exp__2 = 6.2500000000000000E-002
 mdl_aEW = 7.8186082877247844E-003
 mdl_MW = 79.824359746197842
 mdl_sqrt__aEW = 8.8422894590285753E-002
 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_yt = 0.98791394091683138
 mdl_ytau = 1.0206529494239589E-002
 mdl_muH = 88.388347648318444
 mdl_I2a33 = ( 0.98791394091683138 , 0.0000000000000000 )
 mdl_I3a33 = ( 0.98791394091683138 , 0.0000000000000000 )
 mdl_ee__exp__2 = 9.8251529432049817E-002
 mdl_cw__exp__2 = 0.76630086657817542
 mdl_sw__exp__2 = 0.23369913342182450
 mdl_yt__exp__2 = 0.97597395465782455
  Internal Params evaluated point by point
  ----------------------------------------

 mdl_sqrt__aS = 0.34409301068170506
 mdl_G__exp__2 = 1.4878582807401259
 mdl_G__exp__3 = 1.8148567439626970
 mdl_G__exp__4 = 2.2137222635669636
  Couplings of DMsimp_s_spin1-vec
  ---------------------------------

        GC_10 = -0.12198E+01 0.00000E+00
        GC_11 = 0.00000E+00 0.12198E+01
   R2GC_115_5 = -0.00000E+00 -0.30647E-01
  R2GC_132_16 = -0.00000E+00 -0.62813E-02
 UVGC_161_5_1eps -0.00000E+00 -0.49802E-01
 UVGC_179_15_1ep 0.00000E+00 0.38309E-02
 UVGC_179_16_1ep -0.00000E+00 -0.27295E-01
 UVGC_179_17_1ep -0.00000E+00 -0.14366E-02
 UVGC_182_23_1ep 0.00000E+00 0.12563E-01
        GC_43 = 0.00000E+00 0.25000E+00

 Collider parameters:
 --------------------

 Running at P P machine @ 13000.000000000000 GeV
 PDF set = lhapdf
 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

Using random seed offsets: 1 , 1 , 0
  with seed 23
 Ranmar initialization seeds 29840 9397
tau_min 1 1 : 0.20000E+03 -- 0.15100E+04
tau_min 2 1 : 0.20000E+03 -- 0.15100E+04
tau_min 3 1 : 0.21000E+03 0.21000E+03 0.15100E+04
tau_min 4 1 : 0.20000E+03 -- 0.15100E+04
tau_min 5 1 : 0.20000E+03 -- 0.15100E+04
tau_min 6 1 : 0.20000E+03 -- 0.15100E+04
tau_min 7 1 : 0.20000E+03 -- 0.15100E+04
tau_min 8 1 : 0.20000E+03 -- 0.15100E+04
tau_min 9 1 : 0.20000E+03 -- 0.15100E+04
tau_min 10 1 : 0.20000E+03 -- 0.15100E+04
 ntry 1
 User-defined scale not set
STOP 1

> On Jul 30, 2020, at 9:05 AM, marco zaro <email address hidden> wrote:
>
> Your question #692086 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692086
>
> Status: Open => Needs information
>
> marco zaro requested more information:
> Hi Jack,
> some extra information would be needed.
> I suggest you to run with a fixed random seed in the run_card (iseed), and see if you still get that error.
> In that case, go into the P0 dir which has no channel, and run ./gensym NLO
>
> What process are you considering?
>
> Cheers,
>
> Marco
>
> --
> To answer this request for more information, you can either reply to
> this email or enter your reply at the following page:
> https://answers.launchpad.net/mg5amcnlo/+question/692086
>
> You received this question notification because you asked the question.

Revision history for this message
marco zaro (marco-zaro) said : 		#3

Hi Jack,
ok, thanks.
It looks like you are using
   10 = dynamical_scale_choice ! Choose one (or more) of the predefined
in the run_card, but you have not defined the corresponding scale inside setscales.f (around line 564)

Let me know if this helps

Marco

> On 30 Jul 2020, at 10:55, Jack Y. Araz <email address hidden> wrote:
>
> Question #692086 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692086
>
> Status: Needs information => Open
>
> Jack Y. Araz gave more information on the question:
> Hi Marco
>
> Thanks for the answer. I'm running DMsimp model from feynrules (tried
> both the combined and separate packages under the dedicated page) and
> running 'p p > xd xd~ j [QCD]' with different restrict cards. For
> instance, for a spin-1-mediated process, the gV** are on, and rest is
> zero etc.
>
> I tried with a fixed seed and got the same error (currently running for
> spin-1-mediated one) One of the integration channels that shows up is
> P0_gu_xdxdxu, so I went into that folder and ran './gensym NL O' as
> requested. The result is attached below.
>
> Thanks
> Best regards
> —
> Jack Y. Araz
>
> ———————————————
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11839999999999999
> ==== LHAPDF6 USING DEFAULT-TYPE LHAGLUE INTERFACE ====
> LHAPDF 6.2.3 loading /mt/batch/jaraz/MG5_aMC_v2_7_3/HEPTools/lhapdf6/share/LHAPDF/NNPDF30_nlo_as_0118/NNPDF30_nlo_as_0118_0000.dat
> NNPDF30_nlo_as_0118 PDF set, member #0, version 2; LHAPDF ID = 260000
> New value of alpha_s from PDF lhapdf : 0.11800222660863767
> WARNING: the value of maxjetflavorspecified in the run_card ( 4) is inconsistent with the number of light flavours inthe model. Hence it will be set to: 5
>
>
> External Params
> ---------------------------------
>
> MU_R = 91.188000000000002
> mdl_cabi = 0.22773599999999999
> mdl_gVXd = 0.25000000000000000
> aEWM1 = 127.90000000000001
> mdl_Gf = 1.1663700000000000E-005
> aS = 0.11839999999999999
> mdl_ymt = 172.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_MZ = 91.187600000000003
> mdl_MTA = 1.7769999999999999
> mdl_MT = 172.00000000000000
> mdl_MH = 125.00000000000000
> mdl_MXr = 100.00000000000000
> mdl_MY1 = 1500.0000000000000
> mdl_WZ = 2.4952000000000001
> mdl_WW = 2.0850000000000000
> mdl_WT = 1.5083359999999999
> mdl_WH = 4.0699999999999998E-003
> mdl_WY1 = 10.000000000000000
> Internal Params
> ---------------------------------
>
> mdl_MXd = 100.00000000000000
> mdl_MXc = 100.00000000000000
> mdl_gVu22 = 0.25000000000000000
> mdl_gVd22 = 0.25000000000000000
> mdl_gVu11 = 0.25000000000000000
> mdl_gVd11 = 0.25000000000000000
> mdl_gAu11 = 0.0000000000000000
> mdl_gAd22 = 0.0000000000000000
> mdl_gAu22__exp__2 = 0.0000000000000000
> mdl_cos__cabi = 0.97418004031982097
> mdl_CKM1x1 = ( 0.97418004031982097 , 0.0000000000000000 )
> mdl_sin__cabi = 0.22577256042856930
> mdl_CKM1x2 = ( 0.22577256042856930 , 0.0000000000000000 )
> mdl_CKM2x1 = (-0.22577256042856930 , 0.0000000000000000 )
> mdl_CKM2x2 = ( 0.97418004031982097 , 0.0000000000000000 )
> 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.0000000000000000 , 1.0000000000000000 )
> mdl_conjg__CKM1x1 = ( 0.97418004031982097 , -0.0000000000000000 )
> mdl_conjg__CKM1x2 = ( 0.22577256042856930 , -0.0000000000000000 )
> mdl_conjg__CKM2x1 = (-0.22577256042856930 , -0.0000000000000000 )
> mdl_conjg__CKM2x2 = ( 0.97418004031982097 , -0.0000000000000000 )
> mdl_MT__exp__2 = 29584.000000000000
> mdl_gVu22__exp__2 = 6.2500000000000000E-002
> mdl_gVd11__exp__2 = 6.2500000000000000E-002
> mdl_gVd22__exp__2 = 6.2500000000000000E-002
> mdl_gVu11__exp__2 = 6.2500000000000000E-002
> mdl_aEW = 7.8186082877247844E-003
> mdl_MW = 79.824359746197842
> mdl_sqrt__aEW = 8.8422894590285753E-002
> 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_yt = 0.98791394091683138
> mdl_ytau = 1.0206529494239589E-002
> mdl_muH = 88.388347648318444
> mdl_I2a33 = ( 0.98791394091683138 , 0.0000000000000000 )
> mdl_I3a33 = ( 0.98791394091683138 , 0.0000000000000000 )
> mdl_ee__exp__2 = 9.8251529432049817E-002
> mdl_cw__exp__2 = 0.76630086657817542
> mdl_sw__exp__2 = 0.23369913342182450
> mdl_yt__exp__2 = 0.97597395465782455
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34409301068170506
> mdl_G__exp__2 = 1.4878582807401259
> mdl_G__exp__3 = 1.8148567439626970
> mdl_G__exp__4 = 2.2137222635669636
> Couplings of DMsimp_s_spin1-vec
> ---------------------------------
>
> GC_10 = -0.12198E+01 0.00000E+00
> GC_11 = 0.00000E+00 0.12198E+01
> R2GC_115_5 = -0.00000E+00 -0.30647E-01
> R2GC_132_16 = -0.00000E+00 -0.62813E-02
> UVGC_161_5_1eps -0.00000E+00 -0.49802E-01
> UVGC_179_15_1ep 0.00000E+00 0.38309E-02
> UVGC_179_16_1ep -0.00000E+00 -0.27295E-01
> UVGC_179_17_1ep -0.00000E+00 -0.14366E-02
> UVGC_182_23_1ep 0.00000E+00 0.12563E-01
> GC_43 = 0.00000E+00 0.25000E+00
>
> Collider parameters:
> --------------------
>
> Running at P P machine @ 13000.000000000000 GeV
> PDF set = lhapdf
> 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
>
>
> Using random seed offsets: 1 , 1 , 0
> with seed 23
> Ranmar initialization seeds 29840 9397
> tau_min 1 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 2 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 3 1 : 0.21000E+03 0.21000E+03 0.15100E+04
> tau_min 4 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 5 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 6 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 7 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 8 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 9 1 : 0.20000E+03 -- 0.15100E+04
> tau_min 10 1 : 0.20000E+03 -- 0.15100E+04
> ntry 1
> User-defined scale not set
> STOP 1
>
>
>> On Jul 30, 2020, at 9:05 AM, marco zaro <email address hidden> wrote:
>>
>> Your question #692086 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/692086
>>
>> Status: Open => Needs information
>>
>> marco zaro requested more information:
>> Hi Jack,
>> some extra information would be needed.
>> I suggest you to run with a fixed random seed in the run_card (iseed), and see if you still get that error.
>> In that case, go into the P0 dir which has no channel, and run ./gensym NLO
>>
>> What process are you considering?
>>
>> Cheers,
>>
>> Marco
>>
>> --
>> To answer this request for more information, you can either reply to
>> this email or enter your reply at the following page:
>> https://answers.launchpad.net/mg5amcnlo/+question/692086
>>
>> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are assigned to this
> question.

Revision history for this message
Jack Y. Araz (jackaraz) said : 		#4

Hi Marco

My setscales.f between L564-575 (without the comments) is as follows

      elseif(dynamical_scale_choice.eq.10) then
  write(*,*) "User-defined scale not set"
         stop 1
         temp_scale_id='User-defined dynamical scale' ! use a meaningful string
  tmp = 1500d0

Which is the mediator mass. As you suggested I used -1 for the dynamical scale choice and it works without a problem. I realised that I forgot to remove

  write(*,*) "User-defined scale not set"
         stop 1

from the setscales.f without those, it worked nicely.

Thank you very much
Best regards
Jack Y. Araz