Particle.PID lost

Asked by Maosen Zhou

Dear experts,
I'm generating pp>H>hh>WWWW samples.
First, I produced pp>H>hh with MadGraph5;
Second, i passed lhe files produced in last step to standalone Pythia8. In this step, actually i just modified main43.cc in examples/ by adding four lines in Initialization: pythia.readString("25:onMode = off");
           pythia.readString("25:onIfMatch = 24 -24");
           pythia.readString("24:onMode = 3");
           pythia.readString("24:onIfAny = 11 12 13 14 15 16")
which should be hh > w+w-w+w-, w+ > l+ vl, w- > x x(x means anything). And i checked the production log, the event list proved this.
Third, I pass HepMC out to Delphes-3.0.10 for fast simulation.
Then I checked the output root file, like:
root -l delphes_out.root
Delphes->Show(0)
In Particle Container, Particle.PID = 21, 21, 21, 1560, 21, 21, 21, 1560, 21, 21, 21, 1560, 21, 21, 2, 2, 21, 2, 21, 21
I did not see 25, 24 and 11~16. I don't know what's going on here.
Maybe this is due to Pythia8 or Delphes3, but it's highly appreciated if you could drop some sentences. Thanks!

Cheers,
Maosen

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

Hi,

Please contact Delphi’s author.
> I did not see 25, 24
This sounds normal to me since this is detector level.

Cheers,

Olivier

> On Oct 26, 2015, at 06:13, Maosen Zhou <email address hidden> wrote:
>
> New question #273057 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/273057
>
> Dear experts,
> I'm generating pp>H>hh>WWWW samples.
> First, I produced pp>H>hh with MadGraph5;
> Second, i passed lhe files produced in last step to standalone Pythia8. In this step, actually i just modified main43.cc in examples/ by adding four lines in Initialization: pythia.readString("25:onMode = off");
> pythia.readString("25:onIfMatch = 24 -24");
> pythia.readString("24:onMode = 3");
> pythia.readString("24:onIfAny = 11 12 13 14 15 16")
> which should be hh > w+w-w+w-, w+ > l+ vl, w- > x x(x means anything). And i checked the production log, the event list proved this.
> Third, I pass HepMC out to Delphes-3.0.10 for fast simulation.
> Then I checked the output root file, like:
> root -l delphes_out.root
> Delphes->Show(0)
> In Particle Container, Particle.PID = 21, 21, 21, 1560, 21, 21, 21, 1560, 21, 21, 21, 1560, 21, 21, 2, 2, 21, 2, 21, 21
> I did not see 25, 24 and 11~16. I don't know what's going on here.
> Maybe this is due to Pythia8 or Delphes3, but it's highly appreciated if you could drop some sentences. Thanks!
>
> Cheers,
> Maosen
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Maosen Zhou (maosen-zhou) said :
#2

Hi Olivier,
Thanks. I already contact Delphi's author.
But i am still confused about one issue. I tried to generate p p > h h > b b~ y y at pythia level with MadGraph5, then pass hep output to Delphes3. Then I checked the delphes out like what i did above, I did see 6, -6 and 25. Do you have any ideas? Thanks!

Cheers,
Maosen

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

No idea of what this function does inside Delphes so I can not help you.

Cheers,

Olivier
> On Oct 26, 2015, at 08:08, Maosen Zhou <email address hidden> wrote:
>
> Question #273057 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/273057
>
> Maosen Zhou posted a new comment:
> Hi Olivier,
> Thanks. I already contact Delphi's author.
> But i am still confused about one issue. I tried to generate p p > h h > b b~ y y at pythia level with MadGraph5, then pass hep output to Delphes3. Then I checked the delphes out like what i did above, I did see 6, -6 and 25. Do you have any ideas? Thanks!
>
> Cheers,
> Maosen
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Maosen Zhou (maosen-zhou) said :
#4

Hi Olivier,
Thanks anyway!

Cheers,
Maosen

Can you help with this problem?

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

To post a message you must log in.