pythia interface, pdg number set to 0

Hi Lisa,

Ok, I've suceed to reproduce the problem (of the missing pid). I'll look how to solve this and why this happen.

Cheers,

Olivier

On Jun 12, 2013, at 11:26 AM, Lisa <email address hidden> wrote:

> Question #230499 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/230499
>
> Description changed to:
> Hi!
> I have a question concering the MG/pythia interface. If I run pythia from madgraph and let particles decay (BR are in the parameter card, model implemented into MG via feynrules), the pythia_events.lhe file sets the pdg numbers of the mother particles to "0" (if their pdgnumber >6100011), as you can see here (3rd+4th line):
> -2 -1 0 0 0 0 0.20678670E+01 0.34615428E+00 0.10854676E+04 0.10854696E+04 0.00000000E+00 0. 0.
> 2 -1 0 0 0 0 0.45198779E+00 0.86497999E+00 -0.13690725E+04 0.13690728E+04 0.00000000E+00 0. 0.
> 0 2 1 1 0 0 0.81989152E+03 0.55163525E+03 0.36825329E+03 0.11694502E+04 0.50545001E+03 0. 0.
> 0 2 1 1 0 0 -0.81737167E+03 -0.55042412E+03 -0.65185818E+03 0.12850922E+04 0.50545001E+03 0. 0.
> 5100022 2 3 3 0 0 0.81895122E+03 0.54635823E+03 0.36595721E+03 0.11636180E+04 0.50089300E+03 0. 0.
> -13 2 3 3 0 0 0.94030035E+00 0.52770225E+01 0.22960764E+01 0.58321748E+01 0.10566000E+00 0. 0.
> ...
> The decay afterwards is fine and has the correct numbers. The pythia.log file has the corresponding pdg number:
> 11 -6100013 7 16 17 -18.78653 -76.73302 -1148.62623 1257.40321 505.45001
> 11 6100013 8 18 19 20.97614 73.18729 -111.22014 523.11185 505.45001
> 1 5100022 9 0 0 -19.49695 -71.56314 -1137.74663 1245.33634 500.89300
> 1 -13 10 0 0 0.71043 -5.16988 -10.87961 12.06687 0.10566
> This happens only for pdg numbers >61000001, for the corresponding particle with 5100001 it works and the pdg number is transfered correctly.
> Do you have any idea what might cause this problem and, besides choosing smaller pdg numbers, is there a solution?
> Cheers and thanks a lot!
> Lisa
>
> Edit: I tried to change the pdgnumbers now. It works for 4 out of ~12
> particles fine. If I change the next one (independent of the specific
> pdg number), MG gives a nice crosssection, but the pythia.lhe event file
> is empty:
>
> </LesHouchesEvents>
>
> -->
> </LesHouchesEvents>
>
> The lhe file of pythia starts with:
> <LesHouchesEvents version="1.0">
> <init>
> 0 0 0.00000E+00 0.00000E+00 0 0 0 0 3 1
> 0.00000E+00 0.00000E+00 0.10000E+01 100
> </init>
> so that even the inital particles dont get any information. I checked that this is independet of the pdg numbers I changed, only the total number of altered numbers seem to play a role. The pythia.log file is also empty.
> I forgot to mention, I use MG 5.10 and the pythia that is downloaded with it.
> Also, if I change the pdg numbers in feynrules aso, it is the same problem.
> Cheers and thanks a lot,
> Lisa
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Hi Olivier!
Could you find any solution to this problem?
Cheers,
Lisa

I am having a similar problem. My pythia_events.lhe shows for a given event:
      21 -1 0 0 0 0 -0.30595011E+02 0.77038407E+00 0.50681440E+03 0.50773761E+03 0.00000000E+00 0. 0.
       1 -1 0 0 0 0 -0.48499330E+01 0.99370577E+01 -0.53267445E+03 0.53278921E+03 0.00000000E+00 0. 0.
       1 2 1 1 0 0 -0.30821189E+03 -0.36208509E+03 -0.51083153E+02 0.47823599E+03 0.50400002E-02 0. 0.
       0 2 1 1 0 0 0.28096806E+03 0.34802600E+03 0.23458444E+02 0.45892858E+03 0.10000000E+03 0. 0.
       0 2 2 2 0 0 -0.82011117E+01 0.24766527E+02 0.17646529E+01 0.10336224E+03 0.10000000E+03 0. 0.
      21 1 0 0 0 0 0.29059903E+03 0.36290867E+03 0.23197024E+02 0.46626214E+03 0.26679143E+02 0. 0.
      21 1 0 0 0 0 -0.29153930E+03 -0.33726327E+03 -0.44092264E+02 0.45420703E+03 0.74954781E+02 0. 0.
      21 1 0 0 0 0 -0.39692442E+02 0.13113426E+03 0.34662100E+01 0.13890130E+03 0.22580411E+02 0. 0.
      21 1 0 0 0 0 0.15621148E+02 -0.47134192E+02 0.37385040E+03 0.37720860E+03 0.75202123E+01 0. 0.
      21 1 0 0 0 0 0.97622783E+01 -0.51952584E+01 0.87831879E+00 0.11966461E+02 0.44868647E+01 0. 0.
      21 1 0 0 0 0 0.52743186E+01 0.95255554E+01 0.12395879E+02 0.16916309E+02 0.37348429E+01 0. 0.
      12 1 0 0 0 0 -0.76124715E+01 -0.77500842E+01 0.00000000E+00 0.10863403E+02 0.00000000E+00 0. 0.

And my pythia.log shows for the same event:
    1 !p+! 21 2212 0 0 0 0.00000 0.00000 3999.99989 4000.00000 0.93827
    2 !p+! 21 2212 0 0 0 0.00000 0.00000 -3999.99989 4000.00000 0.93827
 ==================================================================================================================================
    3 !d! 21 1 1 0 0 4.20584 -0.37510 1167.31128 1167.31892 0.00000
    4 !d! 21 1 2 0 0 1.89846 1.51675 -710.18881 710.19297 0.00000
    5 !g! 21 21 3 0 0 -30.59501 0.77038 506.81440 507.73761 0.00000
    6 !d! 21 1 4 0 0 -4.84993 9.93706 -532.67445 532.78921 0.00000
    7 !d! 21 1 0 0 0 -308.21189 -362.08509 -51.08315 478.23599 0.00504
    8 !! 21 51 0 0 0 280.96806 348.02600 23.45844 458.92858 100.00000
    9 !! 21 -51 0 0 0 -8.20111 24.76653 1.76465 103.36224 100.00000
 ==================================================================================================================================
   10 1 51 8 0 0 280.96806 348.02600 23.45844 458.92858 100.00000
   11 1 -51 9 0 0 -8.20111 24.76653 1.76465 103.36224 100.00000
...

Note that the PDG code +/-51 should not decay, but in the lhe file it is listed as an intermediate, rather than final state, particle. Also, interestingly the particle should have a name, which is not showing up in the log file. Any idea how to fix this? Just change the PDG number?

Hi,

That's a question for the Pythia author.

Cheers,

Olivier

On Jul 25, 2013, at 4:21 PM, Nick <email address hidden> wrote:

> Question #230499 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/230499
>
> Nick requested more information:
> I am having a similar problem. My pythia_events.lhe shows for a given event:
> 21 -1 0 0 0 0 -0.30595011E+02 0.77038407E+00 0.50681440E+03 0.50773761E+03 0.00000000E+00 0. 0.
> 1 -1 0 0 0 0 -0.48499330E+01 0.99370577E+01 -0.53267445E+03 0.53278921E+03 0.00000000E+00 0. 0.
> 1 2 1 1 0 0 -0.30821189E+03 -0.36208509E+03 -0.51083153E+02 0.47823599E+03 0.50400002E-02 0. 0.
> 0 2 1 1 0 0 0.28096806E+03 0.34802600E+03 0.23458444E+02 0.45892858E+03 0.10000000E+03 0. 0.
> 0 2 2 2 0 0 -0.82011117E+01 0.24766527E+02 0.17646529E+01 0.10336224E+03 0.10000000E+03 0. 0.
> 21 1 0 0 0 0 0.29059903E+03 0.36290867E+03 0.23197024E+02 0.46626214E+03 0.26679143E+02 0. 0.
> 21 1 0 0 0 0 -0.29153930E+03 -0.33726327E+03 -0.44092264E+02 0.45420703E+03 0.74954781E+02 0. 0.
> 21 1 0 0 0 0 -0.39692442E+02 0.13113426E+03 0.34662100E+01 0.13890130E+03 0.22580411E+02 0. 0.
> 21 1 0 0 0 0 0.15621148E+02 -0.47134192E+02 0.37385040E+03 0.37720860E+03 0.75202123E+01 0. 0.
> 21 1 0 0 0 0 0.97622783E+01 -0.51952584E+01 0.87831879E+00 0.11966461E+02 0.44868647E+01 0. 0.
> 21 1 0 0 0 0 0.52743186E+01 0.95255554E+01 0.12395879E+02 0.16916309E+02 0.37348429E+01 0. 0.
> 12 1 0 0 0 0 -0.76124715E+01 -0.77500842E+01 0.00000000E+00 0.10863403E+02 0.00000000E+00 0. 0.
>
>
> And my pythia.log shows for the same event:
> 1 !p+! 21 2212 0 0 0 0.00000 0.00000 3999.99989 4000.00000 0.93827
> 2 !p+! 21 2212 0 0 0 0.00000 0.00000 -3999.99989 4000.00000 0.93827
> ==================================================================================================================================
> 3 !d! 21 1 1 0 0 4.20584 -0.37510 1167.31128 1167.31892 0.00000
> 4 !d! 21 1 2 0 0 1.89846 1.51675 -710.18881 710.19297 0.00000
> 5 !g! 21 21 3 0 0 -30.59501 0.77038 506.81440 507.73761 0.00000
> 6 !d! 21 1 4 0 0 -4.84993 9.93706 -532.67445 532.78921 0.00000
> 7 !d! 21 1 0 0 0 -308.21189 -362.08509 -51.08315 478.23599 0.00504
> 8 !! 21 51 0 0 0 280.96806 348.02600 23.45844 458.92858 100.00000
> 9 !! 21 -51 0 0 0 -8.20111 24.76653 1.76465 103.36224 100.00000
> ==================================================================================================================================
> 10 1 51 8 0 0 280.96806 348.02600 23.45844 458.92858 100.00000
> 11 1 -51 9 0 0 -8.20111 24.76653 1.76465 103.36224 100.00000
> ...
>
>
> Note that the PDG code +/-51 should not decay, but in the lhe file it is listed as an intermediate, rather than final state, particle. Also, interestingly the particle should have a name, which is not showing up in the log file. Any idea how to fix this? Just change the PDG number?
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

The discussion at the following link, particularly comment #6, may be of some help. I am not sure that it fixed my problem but it may fix yours, Lisa:

https://answers.launchpad.net/madgraph5/+question/222285