MadGraph5_aMC@NLO

Error when generating process Z + jets

Asked by Gilberto

Hi Madgraph team

I am trying to generate the process

generate p p > Z jp, ( Z > mu+ mu- )

where

define jp = g u c d s b u~ c~ d~ s~ b~

To do this I am considering the mass of the top to 170 GeV , so I allow madgraph to determine the width automatically. These are the relevant modifications inside the run card:

    50000.0 = ebeam1 ! beam 1 total energy in GeV
    50000.0 = ebeam2 ! beam 2 total energy in GeV
.....
......

#*********************************************************************
# BW cutoff (M+/-bwcutoff*Gamma) ! Define on/off-shell for "$" and decay
#*********************************************************************
  10.0 = bwcutoff ! (M+/-bwcutoff*Gamma)

.....
....

#*********************************************************************
   True = cut_decays ! Cut decay products

Then no cuts afterwards but:

70.0 = mmll ! min invariant mass of l+l- (same flavour) lepton pair

-----------------------------------------------------------------------------------------------------------------------------------------------------
*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/

When I run my session I got:

Current estimate of cross-section: 202235.23 +- 21790.4797481
    P1_gq_zq_z_ll
    P1_qq_zg_z_ll
INFO: Idle: 1, Running: 12, Completed: 0 [ current time: 11h05 ]
INFO: Idle: 0, Running: 13, Completed: 0 [ 5m 0s ]
INFO: Idle: 0, Running: 12, Completed: 1 [ 6m 38s ]
At line 174 of file driver.f (unit = 27, file = '/tmp/gfortrantmp9dWDTL')
Fortran runtime error: Sequential READ or WRITE not allowed after EOF marker, possibly use REWIND or BACKSPACE
ERROR DETECTED
INFO: Idle: 0, Running: 11, Completed: 2 [ 6m 42s ]
At line 174 of file driver.f (unit = 27, file = '/tmp/gfortrantmpyJ4ckI')
Fortran runtime error: Sequential READ or WRITE not allowed after EOF marker, possibly use REWIND or BACKSPACE
ERROR DETECTED
INFO: Idle: 0, Running: 10, Completed: 3 [ 6m 48s ]
At line 378 of file unwgt.f
Fortran runtime error: Array reference out of bounds for array 'store_event', upper bound of dimension 1 exceeded (538976302 > 200000)

-----------------------------------------------------------------------------------------------------------------------------------------------------
*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/
*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/*/

-This happens several times until

Error when reading /data/theorie/gtx/WorktbWZ/samples/Z_jets_mup_mum_v2/SubProcesses/P1_gq_zq_z_ll/G2a0/results.dat
Command "generate_events run_01" interrupted with error:
ValueError : could not convert string to float:
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in '/data/theorie/gtx/WorktbWZ/samples/Z_jets_mup_mum_v2/run_01_tag_1_debug.log'.
Please attach this file to your report.

Any idea of what is the problem?.

Thanks
Gilberto

Question information

Language:
English Edit question
Status:
Solved
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

Hi Gilberto,

Do I understand correctly that what you try to do is equivalent to the following script?

define jp = g u c d s b u~ c~ d~ s~ b~
generate p p > Z jp, ( Z > mu+ mu- )
output
launch
set LHC 100
set bwcutoff 10
set cut_decays T
set mmll 70
set mt 170
set wt Auto

I have run this and it run perfectly fine on my laptop (in less than a minute)
#***************************************************************************
#
# original cross-section: 1863.12111955
# scale variation: +31.2% -25.5%
# central scheme variation: +3.72% -6.53%
# PDF variation: +3.27% -3.27%
#
# dynamical scheme # 1 : 1886.05 +30% -24.6% # \sum ET
# dynamical scheme # 2 : 1886.05 +30% -24.6% # \sum\sqrt{m^2+pt^2}
# dynamical scheme # 3 : 1741.47 +36.6% -28.9% # 0.5 \sum\sqrt{m^2+pt^2}
# dynamical scheme # 4 : 1932.37 +27.6% -23% # \sqrt{\hat s}
#***************************************************************************

I guess that this is not the case, since your problem seems related to a divergency of your matrix-element which is not handle by the cut that you provide.
The absence of those cut make the phase-space integration to fail to converge. The code then increase the number of phase-space points to force the convergence (This will actually not help if you miss a cut but ok)
Then at some point, you hit a pure local problem that your scratch disk space is too small to handle that number of PS point.

Cheers,

Olivier

On 20 Feb 2017, at 11:43, Gilberto <<email address hidden><mailto:<email address hidden>>> wrote:

  True = cut_decays

Revision history for this message
Gilberto (gilbtx) said : 		#2

Hi Oliver

Thanks a lot for your quick answer! I was making some experiments. And the problem seems to be not imposing a lower bound on the pt of the jets. Once this is done things seem to be ok. During your run I assume some lower bound on the Pt was used.

Cheers
Gilberto

On Monday, February 20, 2017 13:38 CET, Olivier Mattelaer <email address hidden> wrote:

> Your question #466518 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi Gilberto,
>
> Do I understand correctly that what you try to do is equivalent to the
> following script?
>
> define jp = g u c d s b u~ c~ d~ s~ b~
> generate p p > Z jp, ( Z > mu+ mu- )
> output
> launch
> set LHC 100
> set bwcutoff 10
> set cut_decays T
> set mmll 70
> set mt 170
> set wt Auto
>
> I have run this and it run perfectly fine on my laptop (in less than a minute)
> #***************************************************************************
> #
> # original cross-section: 1863.12111955
> # scale variation: +31.2% -25.5%
> # central scheme variation: +3.72% -6.53%
> # PDF variation: +3.27% -3.27%
> #
> # dynamical scheme # 1 : 1886.05 +30% -24.6% # \sum ET
> # dynamical scheme # 2 : 1886.05 +30% -24.6% # \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 3 : 1741.47 +36.6% -28.9% # 0.5 \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 4 : 1932.37 +27.6% -23% # \sqrt{\hat s}
> #***************************************************************************
>
>
> I guess that this is not the case, since your problem seems related to a divergency of your matrix-element which is not handle by the cut that you provide.
> The absence of those cut make the phase-space integration to fail to converge. The code then increase the number of phase-space points to force the convergence (This will actually not help if you miss a cut but ok)
> Then at some point, you hit a pure local problem that your scratch disk space is too small to handle that number of PS point.
>
>
> Cheers,
>
> Olivier
>
>
>
> On 20 Feb 2017, at 11:43, Gilberto <<email address hidden><mailto:<email address hidden>>> wrote:
>
> True = cut_decays
>
> -- If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/mg5amcnlo/+question/466518/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> You received this question notification because you asked the question.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#3

Hi Gilberto,

You indeed need a cut on the pt of the jet.
If you are interested in the low pt regime as well.
Then you need to use some matching/merging method (like MLM).

Cheers,

Olivier

On 20 Feb 2017, at 15:27, Gilberto <<email address hidden><mailto:<email address hidden>>> wrote:

Question #466518 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/466518

   Status: Answered => Open

Gilberto is still having a problem:
Hi Oliver

Thanks a lot for your quick answer! I was making some experiments. And
the problem seems to be not imposing a lower bound on the pt of the
jets. Once this is done things seem to be ok. During your run I assume
some lower bound on the Pt was used.

Cheers
Gilberto

On Monday, February 20, 2017 13:38 CET, Olivier Mattelaer <<email address hidden><mailto:<email address hidden>>> wrote:

Your question #466518 on MadGraph5_aMC@NLO changed:
https://answers.launchpad.net/mg5amcnlo/+question/466518

   Status: Open => Answered

Olivier Mattelaer proposed the following answer:
Hi Gilberto,

Do I understand correctly that what you try to do is equivalent to the
following script?

define jp = g u c d s b u~ c~ d~ s~ b~
generate p p > Z jp, ( Z > mu+ mu- )
output
launch
set LHC 100
set bwcutoff 10
set cut_decays T
set mmll 70
set mt 170
set wt Auto

I have run this and it run perfectly fine on my laptop (in less than a minute)
#***************************************************************************
#
# original cross-section: 1863.12111955
# scale variation: +31.2% -25.5%
# central scheme variation: +3.72% -6.53%
# PDF variation: +3.27% -3.27%
#
# dynamical scheme # 1 : 1886.05 +30% -24.6% # \sum ET
# dynamical scheme # 2 : 1886.05 +30% -24.6% # \sum\sqrt{m^2+pt^2}
# dynamical scheme # 3 : 1741.47 +36.6% -28.9% # 0.5 \sum\sqrt{m^2+pt^2}
# dynamical scheme # 4 : 1932.37 +27.6% -23% # \sqrt{\hat s}
#***************************************************************************

I guess that this is not the case, since your problem seems related to a divergency of your matrix-element which is not handle by the cut that you provide.
The absence of those cut make the phase-space integration to fail to converge. The code then increase the number of phase-space points to force the convergence (This will actually not help if you miss a cut but ok)
Then at some point, you hit a pure local problem that your scratch disk space is too small to handle that number of PS point.

Cheers,

Olivier

On 20 Feb 2017, at 11:43, Gilberto <<email address hidden><mailto:<email address hidden>><mailto:<email address hidden>>> wrote:

 True = cut_decays

-- If this answers your question, please go to the following page to let us
know that it is solved:
https://answers.launchpad.net/mg5amcnlo/+question/466518/+confirm?answer_id=0

If you still need help, you can reply to this email or go to the
following page to enter your feedback:
https://answers.launchpad.net/mg5amcnlo/+question/466518

You received this question notification because you asked the question.

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

Revision history for this message
Gilberto (gilbtx) said : 		#4

Hi Oliver

Thanks a lot for the information!

Cheers
Gilberto

On Monday, February 20, 2017 20:17 CET, Olivier Mattelaer <email address hidden> wrote:

> Your question #466518 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi Gilberto,
>
> You indeed need a cut on the pt of the jet.
> If you are interested in the low pt regime as well.
> Then you need to use some matching/merging method (like MLM).
>
> Cheers,
>
> Olivier
>
> On 20 Feb 2017, at 15:27, Gilberto
> <<email address hidden><mailto:<email address hidden>>>
> wrote:
>
> Question #466518 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> Status: Answered => Open
>
> Gilberto is still having a problem:
> Hi Oliver
>
> Thanks a lot for your quick answer! I was making some experiments. And
> the problem seems to be not imposing a lower bound on the pt of the
> jets. Once this is done things seem to be ok. During your run I assume
> some lower bound on the Pt was used.
>
> Cheers
> Gilberto
>
> On Monday, February 20, 2017 13:38 CET, Olivier Mattelaer
> <<email address hidden><mailto:<email address hidden>>>
> wrote:
>
> Your question #466518 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi Gilberto,
>
> Do I understand correctly that what you try to do is equivalent to the
> following script?
>
> define jp = g u c d s b u~ c~ d~ s~ b~
> generate p p > Z jp, ( Z > mu+ mu- )
> output
> launch
> set LHC 100
> set bwcutoff 10
> set cut_decays T
> set mmll 70
> set mt 170
> set wt Auto
>
> I have run this and it run perfectly fine on my laptop (in less than a minute)
> #***************************************************************************
> #
> # original cross-section: 1863.12111955
> # scale variation: +31.2% -25.5%
> # central scheme variation: +3.72% -6.53%
> # PDF variation: +3.27% -3.27%
> #
> # dynamical scheme # 1 : 1886.05 +30% -24.6% # \sum ET
> # dynamical scheme # 2 : 1886.05 +30% -24.6% # \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 3 : 1741.47 +36.6% -28.9% # 0.5 \sum\sqrt{m^2+pt^2}
> # dynamical scheme # 4 : 1932.37 +27.6% -23% # \sqrt{\hat s}
> #***************************************************************************
>
>
> I guess that this is not the case, since your problem seems related to a divergency of your matrix-element which is not handle by the cut that you provide.
> The absence of those cut make the phase-space integration to fail to converge. The code then increase the number of phase-space points to force the convergence (This will actually not help if you miss a cut but ok)
> Then at some point, you hit a pure local problem that your scratch disk space is too small to handle that number of PS point.
>
>
> Cheers,
>
> Olivier
>
>
> On 20 Feb 2017, at 11:43, Gilberto
> <<email address hidden><mailto:<email address hidden>><mailto:<email address hidden>>>
> wrote:
>
> True = cut_decays
>
> -- If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/mg5amcnlo/+question/466518/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> You received this question notification because you asked the question.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.
>
> -- If this answers your question, please go to the following page to let us
> know that it is solved:
> https://answers.launchpad.net/mg5amcnlo/+question/466518/+confirm?answer_id=2
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/mg5amcnlo/+question/466518
>
> You received this question notification because you asked the question.

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#5