<init> NAN values when using tW with MadSTRPlugin with istr==5

Asked by Javier

Dear All,

So in other question in this website we set MG to photoproduce tW out of pp collisions at 13 TeV. The photon is associated with one of the protons. We run along the MadSTRPlugin and are interested in the following istr versions:

 istr==1
 istr==5

 that we call 'DR' and 'DS' respectively. So with DR version things go fine but with DS version we have a <init> section that looks like:

 <init>
2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
 NAN NAN NAN 0
</init>

in the .LHE file. Also there are other section with 'nan' entries as shown below. So these issues are not present in the DR one, and the two samples only differ in the istr value.

We are interested in knowing what could be the reason for these 'nan' entries in the DS case, as we are processing these LHE files to add detector effects.

Thanks
Javier

<event>
 10 0 +nan 1.36274300e+02 7.54677110e-03 8.75997870e-02
       22 -1 0 0 0 0 +0.00000000000e+00 +0.00000000000e+00 +6.56586360030e+02 6.56586360030e+02 0.00000000000e+00 0.0000e+00 1.0000e+00
       -5 -1 0 0 0 501 +0.00000000000e+00 +0.00000000000e+00 -3.53228680303e+01 3.56680179697e+01 4.95000000000e+00 0.0000e+00 1.0000e+00
       24 2 1 2 0 0 +1.36383218254e+01 -4.47047477119e+01 +8.18144869133e+01 1.23705324620e+02 8.01553415573e+01 0.0000e+00 0.0000e+00
       12 1 3 3 0 0 +1.89779193914e+01 +1.91552660486e+01 +1.39244965106e+01 3.03476069055e+01 0.00000000000e+00 0.0000e+00 -1.0000e+00
      -11 1 3 3 0 0 -5.33959756597e+00 -6.38600137604e+01 +6.78899904027e+01 9.33577177141e+01 5.10999000000e-04 0.0000e+00 1.0000e+00
       -6 2 1 2 0 501 -1.36383218254e+01 +4.47047477119e+01 +5.39449005087e+02 5.68549053380e+02 1.73373235306e+02 0.0000e+00 0.0000e+00
       -5 1 6 6 0 501 -2.43028769382e+00 +7.65557538783e+01 +1.40161287262e+02 1.59801059749e+02 4.95000000000e+00 0.0000e+00 1.0000e+00
      -24 2 6 6 0 0 -1.12080341316e+01 -3.18510061664e+01 +3.99287717825e+02 4.08747993631e+02 8.06482118129e+01 0.0000e+00 0.0000e+00
       11 1 8 8 0 0 -4.55325824030e+01 -2.08362364371e+01 +1.96410796115e+02 2.02693279711e+02 5.10999000000e-04 0.0000e+00 -1.0000e+00
      -12 1 8 8 0 0 +3.43245482714e+01 -1.10147697293e+01 +2.02876921710e+02 2.06054713920e+02 0.00000000000e+00 0.0000e+00 1.0000e+00
#aMCatNLO 1 0 0 2 4 0.95788887E+03 0.00000000E+00 9 5 0 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00
  <rwgt>
</rwgt>
</event>

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
marco zaro Edit question
Last query:
Last reply:
Revision history for this message
marco zaro (marco-zaro) said :
#1

Dear Javier,
thanks for reporting your issue.
Can I ask what version of MG5_aMC are you using, and the commands for the generation?
Thanks

Marco

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

Hi,
just to let you know that I was not able to reproduce your error.
I am using MG5_aMC v2_8_3_2, for the process

import model loop_sm-no_b_mass
define p = 21 2 4 1 3 -2 -4 -1 -3 5 -5 # pass to 5 flavors
define j = p
generate p p > t w- [QCD]
output tw-str

and I get sensible numbers in the event files for istr 1-->5 (did not try 6).

Did you change any parameter in the cards? I used the default cards...

Thanks

Cheers,

Marco

Events/run_01/events.lhe.gz-259-# istr = 6 -> DS with reshuffling on all FS particles, running BW *
Events/run_01/events.lhe.gz-260-#***********************************************************************
Events/run_01/events.lhe.gz:261: 1 = istr ! strategy to be used to remove resonances
Events/run_01/events.lhe.gz-262- ! appearing in real emissions
Events/run_01/events.lhe.gz-263- True = str_include_pdf ! compensate for PDFs when doing reshuffling
Events/run_01/events.lhe.gz-361- 2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
Events/run_01/events.lhe.gz-362- 0.34620229E+02 0.15352445E+00 0.53064232E+02 0
Events/run_01/events.lhe.gz:363: </init>
Events/run_01/events.lhe.gz-364- <event>
Events/run_01/events.lhe.gz-365- 4 0 0.53064232E+02 0.13522646E+03 0.75467711E-02 0.10316577E+00
--
--
Events/run_02/events.lhe.gz-259-# istr = 6 -> DS with reshuffling on all FS particles, running BW *
Events/run_02/events.lhe.gz-260-#***********************************************************************
Events/run_02/events.lhe.gz:261: 2 = istr ! strategy to be used to remove resonances
Events/run_02/events.lhe.gz-262- ! appearing in real emissions
Events/run_02/events.lhe.gz-263- True = str_include_pdf ! compensate for PDFs when doing reshuffling
Events/run_02/events.lhe.gz-361- 2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
Events/run_02/events.lhe.gz-362- 0.31279916E+02 0.31132120E+00 0.62802172E+02 0
Events/run_02/events.lhe.gz:363: </init>
Events/run_02/events.lhe.gz-364- <event>
Events/run_02/events.lhe.gz-365- 4 0 -.62802172E+02 0.12683800E+03 0.75467711E-02 0.11331579E+00
--
--
Events/run_03/events.lhe.gz-259-# istr = 6 -> DS with reshuffling on all FS particles, running BW *
Events/run_03/events.lhe.gz-260-#***********************************************************************
Events/run_03/events.lhe.gz:261: 3 = istr ! strategy to be used to remove resonances
Events/run_03/events.lhe.gz-262- ! appearing in real emissions
Events/run_03/events.lhe.gz-263- True = str_include_pdf ! compensate for PDFs when doing reshuffling
Events/run_03/events.lhe.gz-361- 2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
Events/run_03/events.lhe.gz-362- 0.33381083E+02 0.29314340E+00 0.59806139E+02 0
Events/run_03/events.lhe.gz:363: </init>
Events/run_03/events.lhe.gz-364- <event>
Events/run_03/events.lhe.gz-365- 4 0 0.59806139E+02 0.14370389E+03 0.75467711E-02 0.11099743E+00
--
--
Events/run_04/events.lhe.gz-259-# istr = 6 -> DS with reshuffling on all FS particles, running BW *
Events/run_04/events.lhe.gz-260-#***********************************************************************
Events/run_04/events.lhe.gz:261: 4 = istr ! strategy to be used to remove resonances
Events/run_04/events.lhe.gz-262- ! appearing in real emissions
Events/run_04/events.lhe.gz-263- True = str_include_pdf ! compensate for PDFs when doing reshuffling
Events/run_04/events.lhe.gz-361- 2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
Events/run_04/events.lhe.gz-362- 0.28805000E+02 0.33320549E+00 0.66276736E+02 0
Events/run_04/events.lhe.gz:363: </init>
Events/run_04/events.lhe.gz-364- <event>
Events/run_04/events.lhe.gz-365- 5 0 0.66276736E+02 0.34282886E+03 0.75467711E-02 0.98435175E-01
--
--
Events/run_05/events.lhe.gz-259-# istr = 6 -> DS with reshuffling on all FS particles, running BW *
Events/run_05/events.lhe.gz-260-#***********************************************************************
Events/run_05/events.lhe.gz:261: 5 = istr ! strategy to be used to remove resonances
Events/run_05/events.lhe.gz-262- ! appearing in real emissions
Events/run_05/events.lhe.gz-263- True = str_include_pdf ! compensate for PDFs when doing reshuffling
Events/run_05/events.lhe.gz-361- 2212 2212 0.65000000E+04 0.65000000E+04 -1 -1 244800 244800 -4 1
Events/run_05/events.lhe.gz-362- 0.34268113E+02 0.29328843E+00 0.60412768E+02 0
Events/run_05/events.lhe.gz:363: </init>
Events/run_05/events.lhe.gz-364- <event>
Events/run_05/events.lhe.gz-365- 5 0 -.60412768E+02 0.13724408E+03 0.75467711E-02 0.11151760E+00

Revision history for this message
Javier (jamkoons) said :
#3

Hello Marco,

Thanks for replying.

so the only thing I change is istr from 1 to 5 and also set lpp1=2 by hand in the runcard afterwards, so tW can be photoproduced and and the photon can be associated with one of the p's.

Thanks
Javier

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

Ok, which syntax are you using for generating the process?

On Thu, 7 Oct 2021, 22:30 Javier, <email address hidden>
wrote:

> Question #698976 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/698976
>
> Status: Needs information => Open
>
> Javier gave more information on the question:
> Hello Marco,
>
> Thanks for replying.
>
> so the only thing I change is istr from 1 to 5 and also set lpp1=2 by
> hand in the runcard afterwards, so tW can be photoproduced and and the
> photon can be associated with one of the p's.
>
> Thanks
> Javier
>
> --
> You received this question notification because you are assigned to this
> question.
>

Can you help with this problem?

Provide an answer of your own, or ask Javier for more information if necessary.

To post a message you must log in.