error:invalid particle format while hadronisation in Pythia

Asked by Antara Sarkar on 2020-05-30

Dear experts,

I was trying to generate vbf process and shower and hadronise with Pythia using my model file as follows:

define v = w+ w- z
define uq = u c d s t b
generate p p > h j j $$v /a g h uq NP=0 SMHLOOP=0, h > a a NP=0, SMHLOOP=1

The SMHLOOP interation is an effective one between higgs and photons. As showering and hadronisation is occuring, I get the following error:

INFO: Combining Events
  === Results Summary for run: run_01 tag: tag_1 ===

     Cross-section : 0.003275 +- 1.184e-06 pb
     Nb of events : 1000000

INFO: Creating Plots for parton level
INFO: Plots for parton level generated, see /home/dell/my_computation/MG5_aMC_v2_6_5/sm_hjj_500000_deta15/HTML/run_01/plots_parton.html
INFO: End Plots for parton level
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
decay_events -from_cards
INFO: Running Pythia
INFO: Creating shower LHE File (for plot)
INFO: Warning! Never use this lhe file for detector studies!
INFO: Creating Pythia LHE Root File
Warning in <TTree::Bronch>: Using split mode on a class: TRootWeight with a custom Streamer
** Reading pythia_events.lhe
** ERROR: invalid particle format#####------------------------------] (53.50%)
** [################################################################] (100.00%)
** Exiting...
INFO: Creating Plots for Pythia level
INFO: Plots for pythia level generated, see /home/dell/my_computation/MG5_aMC_v2_6_5/sm_hjj_500000_deta15/HTML/run_01/plots_pythia_tag_1.html
INFO: End Plots for pythia level
INFO: finish
INFO: No pgs_card detected, so not run pgs
INFO: No delphes_card detected, so not run Delphes
  === Results Summary for run: run_01 tag: tag_1 ===

     Cross-section : 0.003275 +- 1.184e-06 pb
     Nb of events : 1000000

INFO: storing files of previous run
INFO: Storing Pythia files of previous run

The process is completed. But I am not sure if it completely hadronised the events or not. Could you please explain this error and tell if the process is hadronised or not?

Regards,
Antara

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2020-05-31
Last reply:
2020-06-02

This question was reopened

The error seems related to ExRootAnalysis.
Not at all related to pythia.
The log of pythia is actually not displayed. So I can not comment on it (anyway I'm not the expert in pythia log)

The issue with ExRootAnalysis is not critical since you typically do not want to use such weird file
(this is a root conversion of the lhe file created by connverting the hepmc file to a lhef file when applying a very basic clustering). So I will not worry about it.

Cheers,

Olivier

> On 30 May 2020, at 12:50, Antara Sarkar <email address hidden> wrote:
>
> New question #691050 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/691050
>
> Dear experts,
>
> I was trying to generate vbf process and shower and hadronise with Pythia using my model file as follows:
>
> define v = w+ w- z
> define uq = u c d s t b
> generate p p > h j j $$v /a g h uq NP=0 SMHLOOP=0, h > a a NP=0, SMHLOOP=1
>
> The SMHLOOP interation is an effective one between higgs and photons. As showering and hadronisation is occuring, I get the following error:
>
> INFO: Combining Events
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 0.003275 +- 1.184e-06 pb
> Nb of events : 1000000
>
> INFO: Creating Plots for parton level
> INFO: Plots for parton level generated, see /home/dell/my_computation/MG5_aMC_v2_6_5/sm_hjj_500000_deta15/HTML/run_01/plots_parton.html
> INFO: End Plots for parton level
> store_events
> INFO: Storing parton level results
> INFO: End Parton
> reweight -from_cards
> decay_events -from_cards
> INFO: Running Pythia
> INFO: Creating shower LHE File (for plot)
> INFO: Warning! Never use this lhe file for detector studies!
> INFO: Creating Pythia LHE Root File
> Warning in <TTree::Bronch>: Using split mode on a class: TRootWeight with a custom Streamer
> ** Reading pythia_events.lhe
> ** ERROR: invalid particle format#####------------------------------] (53.50%)
> ** [################################################################] (100.00%)
> ** Exiting...
> INFO: Creating Plots for Pythia level
> INFO: Plots for pythia level generated, see /home/dell/my_computation/MG5_aMC_v2_6_5/sm_hjj_500000_deta15/HTML/run_01/plots_pythia_tag_1.html
> INFO: End Plots for pythia level
> INFO: finish
> INFO: No pgs_card detected, so not run pgs
> INFO: No delphes_card detected, so not run Delphes
> === Results Summary for run: run_01 tag: tag_1 ===
>
> Cross-section : 0.003275 +- 1.184e-06 pb
> Nb of events : 1000000
>
> INFO: storing files of previous run
> INFO: Storing Pythia files of previous run
>
> The process is completed. But I am not sure if it completely hadronised the events or not. Could you please explain this error and tell if the process is hadronised or not?
>
> Regards,
> Antara
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Antara Sarkar (antara-s) said : #2

Thanks Olivier Mattelaer, that solved my question.

Antara Sarkar (antara-s) said : #3

Hi Olivier,

Could you please suggest me in one regard? I do not have a cluster in working condition right now in my University. However, I need to generate and study very large number of events (500M) with carying the parameters which is taking a lot of time to run on my computer. Could you please suggest if some help can be provided by the Madgraph in this regard?

Thanks and regards,
Antara

Hi,

For that number of events, I would suggest to use the gridpack strategy.
For computing resources, our web-service is unlikely to support such type of generation.

Cheers,

Olivier

Can you help with this problem?

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

To post a message you must log in.