ee collision with ISR in v3.3.0

Asked by lingfeng-li

Dear MG5 guru,

I am trying to generate some Higgs factory diboson process with ISR effects included (thys leading to lower average Ecm). As mentioned in arXiv:2108.10261, I am supposed to use lpp= 3/-3 for the process with ISR. I also adopted the "set group_subprocesses False" commant from arXiv:2111.02442. The test VBF run is fine. However, when I generate e+ e- > w+ w- process using pdflavel=eva, I got the error that asks to use lpp=0 only.

Here is the full message:
Error when reading /home/lingfeng/MG5_aMC_v3_3_0/eeww/SubProcesses/P1_epem_wpwm/G1/results.dat
Command "generate_events run_01" interrupted with error:
Exception : Reported error: End code 1.0
  Full associated log:
  Process in group number 1
  A PDF is used, so alpha_s(MZ) is going to be modified
  Old value of alpha_s from param_card: 0.11799999999999999
  New value of alpha_s from PDF eva : 0.11799999999999999
  Define smin to 25868.863817482143
  *****************************************************
  * MadGraph/MadEvent *
  * -------------------------------- *
  * http://madgraph.hep.uiuc.edu *
  * http://madgraph.phys.ucl.ac.be *
  * http://madgraph.roma2.infn.it *
  * -------------------------------- *
  * *
  * PARAMETER AND COUPLING VALUES *
  * *
  *****************************************************

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

  mdl_MB = 4.7000000000000002
  mdl_MT = 173.00000000000000
  mdl_MTA = 1.7769999999999999
  mdl_MZ = 91.188000000000002
  mdl_MH = 125.00000000000000
  aEWM1 = 132.50700000000001
  mdl_Gf = 1.1663900000000000E-005
  aS = 0.11799999999999999
  mdl_ymb = 4.7000000000000002
  mdl_ymt = 173.00000000000000
  mdl_ymtau = 1.7769999999999999
  mdl_WT = 1.4915000000000000
  mdl_WZ = 2.4414039999999999
  mdl_WW = 2.0476000000000001
  mdl_WH = 6.3823389999999999E-003
   Internal Params
   ---------------------------------

  mdl_conjg__CKM3x3 = 1.0000000000000000
  mdl_CKM3x3 = 1.0000000000000000
  mdl_conjg__CKM1x1 = 1.0000000000000000
  mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
  mdl_MZ__exp__2 = 8315.2513440000002
  mdl_MZ__exp__4 = 69143404.913893804
  mdl_sqrt__2 = 1.4142135623730951
  mdl_MH__exp__2 = 15625.000000000000
  mdl_aEW = 7.5467711139788835E-003
  mdl_MW = 80.419002445756163
  mdl_sqrt__aEW = 8.6872153846781555E-002
  mdl_ee = 0.30795376724436879
  mdl_MW__exp__2 = 6467.2159543705357
  mdl_sw2 = 0.22224648578577766
  mdl_cw = 0.88190334743339216
  mdl_sqrt__sw2 = 0.47143025548407230
  mdl_sw = 0.47143025548407230
  mdl_g1 = 0.34919219678733299
  mdl_gw = 0.65323293034757990
  mdl_vev = 246.21845810181637
  mdl_vev__exp__2 = 60623.529110035903
  mdl_lam = 0.12886910601690263
  mdl_yb = 2.6995554250465490E-002
  mdl_yt = 0.99366614581500623
  mdl_ytau = 1.0206617000654717E-002
  mdl_muH = 88.388347648318430
  mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
  mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
  mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
  mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
  mdl_ee__exp__2 = 9.4835522759998875E-002
  mdl_sw__exp__2 = 0.22224648578577769
  mdl_cw__exp__2 = 0.77775351421422245
   Internal Params evaluated point by point
   ----------------------------------------

  mdl_sqrt__aS = 0.34351128074635334
  mdl_G__exp__2 = 1.4828317324943823
   Couplings of sm
   ---------------------------------

          GC_3 = -0.00000E+00 -0.30795E+00
         GC_50 = -0.00000E+00 -0.28804E+00
         GC_53 = 0.00000E+00 0.57609E+00
         GC_59 = 0.00000E+00 0.82310E-01
        GC_100 = 0.00000E+00 0.46191E+00

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

  Running at e+e- machine @ 240.00000000000000 GeV
  PDF set = eva
  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 2000 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)?
  Monte-Carlo over helicities

  ********************************************
  * You are using the DiscreteSampler module *
  * part of the MG5_aMC framework *
  * Author: Valentin Hirschi *
  ********************************************

 Enter Configuration Number:
 Running Configuration Number: 1
  Not subdividing B.W.
  Attempting mappinvarients 1 4
  Determine nb_t
  T-channel found: 0
  Completed mapping 4
  about to integrate 4 2000 5 3 4 1
  Using non-zero grid deformation.
   4 dimensions 2000 events 4 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 240.00
  Mapping Graph 1 to config 1
  Determine nb_t
  T-channel found: 0
 Setting grid 1 0.44911E+00 1
  Transforming s_hat 1/s 3 0.44911221905350945 25868.863817482143 57599.999999999993
    1 1 2 3 4
  Masses: 0.000E+00 0.000E+00 0.804E+02 0.804E+02
 Using random seed offsets 1 : 1
   with seed 24
  Ranmar initialization seeds 31059 9398
   Particle 3 4
       Et > 0.0 0.0
        E > 0.0 0.0
      Eta < -1.0 -1.0
    xqcut: 0.0 0.0
 d R # 3 > -0.0 0.0
 s min # 3> 0.0 0.0
 xqcutij # 3> 0.0 0.0
  alpha_s for scale 95.290809227980077 is 0.11722214083519961
  PDF not supported for pdg -11
  For lepton colliders, please set the lpp* variables to 0 in the run_card current is -3
 STOP 1

 ls status:
 input_app.txt
 run1_app.log

Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/home/lingfeng/MG5_aMC_v3_3_0/eeww/run_01_tag_1_debug.log'.
Please attach this file to your report.

Now I wonder if it is only the bosons are available in the current version. If no, what should I do to simulate correctly?

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

pdlabel=eva

is for W/Z/A pdf so indeed asking for
e+ e- > w+ w- is not supported with that PDF.
Since that mode does not have an electron pdf

For ISR, you have the following PDF: isronlyll
or one combined with beamsthralung:
cepc240ll clic3000ll fcce240ll fcce365ll ilc500ll

Cheers,

Olivier

> On 13 Nov 2021, at 08:01, lingfeng-li <email address hidden> wrote:
>
> New question #699475 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/699475
>
> Dear MG5 guru,
>
> I am trying to generate some Higgs factory diboson process with ISR effects included (thys leading to lower average Ecm). As mentioned in arXiv:2108.10261, I am supposed to use lpp= 3/-3 for the process with ISR. I also adopted the "set group_subprocesses False" commant from arXiv:2111.02442. The test VBF run is fine. However, when I generate e+ e- > w+ w- process using pdflavel=eva, I got the error that asks to use lpp=0 only.
>
> Here is the full message:
> Error when reading /home/lingfeng/MG5_aMC_v3_3_0/eeww/SubProcesses/P1_epem_wpwm/G1/results.dat
> Command "generate_events run_01" interrupted with error:
> Exception : Reported error: End code 1.0
> Full associated log:
> Process in group number 1
> A PDF is used, so alpha_s(MZ) is going to be modified
> Old value of alpha_s from param_card: 0.11799999999999999
> New value of alpha_s from PDF eva : 0.11799999999999999
> Define smin to 25868.863817482143
> *****************************************************
> * MadGraph/MadEvent *
> * -------------------------------- *
> * http://madgraph.hep.uiuc.edu *
> * http://madgraph.phys.ucl.ac.be *
> * http://madgraph.roma2.infn.it *
> * -------------------------------- *
> * *
> * PARAMETER AND COUPLING VALUES *
> * *
> *****************************************************
>
> External Params
> ---------------------------------
>
> mdl_MB = 4.7000000000000002
> mdl_MT = 173.00000000000000
> mdl_MTA = 1.7769999999999999
> mdl_MZ = 91.188000000000002
> mdl_MH = 125.00000000000000
> aEWM1 = 132.50700000000001
> mdl_Gf = 1.1663900000000000E-005
> aS = 0.11799999999999999
> mdl_ymb = 4.7000000000000002
> mdl_ymt = 173.00000000000000
> mdl_ymtau = 1.7769999999999999
> mdl_WT = 1.4915000000000000
> mdl_WZ = 2.4414039999999999
> mdl_WW = 2.0476000000000001
> mdl_WH = 6.3823389999999999E-003
> Internal Params
> ---------------------------------
>
> mdl_conjg__CKM3x3 = 1.0000000000000000
> mdl_CKM3x3 = 1.0000000000000000
> mdl_conjg__CKM1x1 = 1.0000000000000000
> mdl_complexi = ( 0.0000000000000000 , 1.0000000000000000 )
> mdl_MZ__exp__2 = 8315.2513440000002
> mdl_MZ__exp__4 = 69143404.913893804
> mdl_sqrt__2 = 1.4142135623730951
> mdl_MH__exp__2 = 15625.000000000000
> mdl_aEW = 7.5467711139788835E-003
> mdl_MW = 80.419002445756163
> mdl_sqrt__aEW = 8.6872153846781555E-002
> mdl_ee = 0.30795376724436879
> mdl_MW__exp__2 = 6467.2159543705357
> mdl_sw2 = 0.22224648578577766
> mdl_cw = 0.88190334743339216
> mdl_sqrt__sw2 = 0.47143025548407230
> mdl_sw = 0.47143025548407230
> mdl_g1 = 0.34919219678733299
> mdl_gw = 0.65323293034757990
> mdl_vev = 246.21845810181637
> mdl_vev__exp__2 = 60623.529110035903
> mdl_lam = 0.12886910601690263
> mdl_yb = 2.6995554250465490E-002
> mdl_yt = 0.99366614581500623
> mdl_ytau = 1.0206617000654717E-002
> mdl_muH = 88.388347648318430
> mdl_I1x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_I2x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I3x33 = ( 0.99366614581500623 , 0.0000000000000000 )
> mdl_I4x33 = ( 2.6995554250465490E-002, 0.0000000000000000 )
> mdl_ee__exp__2 = 9.4835522759998875E-002
> mdl_sw__exp__2 = 0.22224648578577769
> mdl_cw__exp__2 = 0.77775351421422245
> Internal Params evaluated point by point
> ----------------------------------------
>
> mdl_sqrt__aS = 0.34351128074635334
> mdl_G__exp__2 = 1.4828317324943823
> Couplings of sm
> ---------------------------------
>
> GC_3 = -0.00000E+00 -0.30795E+00
> GC_50 = -0.00000E+00 -0.28804E+00
> GC_53 = 0.00000E+00 0.57609E+00
> GC_59 = 0.00000E+00 0.82310E-01
> GC_100 = 0.00000E+00 0.46191E+00
>
> Collider parameters:
> --------------------
>
> Running at e+e- machine @ 240.00000000000000 GeV
> PDF set = eva
> 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 2000 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)?
> Monte-Carlo over helicities
>
> ********************************************
> * You are using the DiscreteSampler module *
> * part of the MG5_aMC framework *
> * Author: Valentin Hirschi *
> ********************************************
>
> Enter Configuration Number:
> Running Configuration Number: 1
> Not subdividing B.W.
> Attempting mappinvarients 1 4
> Determine nb_t
> T-channel found: 0
> Completed mapping 4
> about to integrate 4 2000 5 3 4 1
> Using non-zero grid deformation.
> 4 dimensions 2000 events 4 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 240.00
> Mapping Graph 1 to config 1
> Determine nb_t
> T-channel found: 0
> Setting grid 1 0.44911E+00 1
> Transforming s_hat 1/s 3 0.44911221905350945 25868.863817482143 57599.999999999993
> 1 1 2 3 4
> Masses: 0.000E+00 0.000E+00 0.804E+02 0.804E+02
> Using random seed offsets 1 : 1
> with seed 24
> Ranmar initialization seeds 31059 9398
> Particle 3 4
> Et > 0.0 0.0
> E > 0.0 0.0
> Eta < -1.0 -1.0
> xqcut: 0.0 0.0
> d R # 3 > -0.0 0.0
> s min # 3> 0.0 0.0
> xqcutij # 3> 0.0 0.0
> alpha_s for scale 95.290809227980077 is 0.11722214083519961
> PDF not supported for pdg -11
> For lepton colliders, please set the lpp* variables to 0 in the run_card current is -3
> STOP 1
>
> ls status:
> input_app.txt
> run1_app.log
>
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in '/home/lingfeng/MG5_aMC_v3_3_0/eeww/run_01_tag_1_debug.log'.
> Please attach this file to your report.
>
> Now I wonder if it is only the bosons are available in the current version. If no, what should I do to simulate correctly?
>
>
> --
> 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 lingfeng-li for more information if necessary.

To post a message you must log in.