Zero events at hadron level analysis, cross section information is lost

Asked by Min He on 2017-06-15

Dear MadAnalysis team,

I am using MadGraph v2.5.5 to generate some events. The process card is shown below.

generate e+ e- > a mu+ mu-
output
launch
shower=PYTHIA8
detector=DELPHES
analysis=MADANALYSIS_5
set ebeam1 120
set ebeam2 120

The parton level analysis is OK. http://202.120.47.101:23334/PROC_sm_2/HTML/run_01/tag_1_MA5_PARTON_ANALYSIS_analysis1/HTML/index.html
But I got zero events at the hadron level analysis. http://202.120.47.101:23334/PROC_sm_2/HTML/run_01/tag_1_MA5_HADRON_ANALYSIS_analysis2_BasicReco/HTML/index.html

It seems the cross section information is lost at the hadron level. I think I have the same problem with this: https://answers.launchpad.net/madanalysis5/+question/627403
I want to know how to solve this problem.

Thanks!

Min

Question information

Language:
English Edit question
Status:
Answered
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Last query:
2017-06-19
Last reply:
2017-06-19
Benjamin Fuks (fuks) said : #1

Hi Min,

The information on the cross section is not passed through the files for the simple reason that in the hepmc format, the cross section value is updated event after event and only the last value matter. However, the xsection information in the output file is written before the reading of any event.

You can fix it (in the meantime we change this) by adding to the script:
  set.<name of the dataset>.xsection = …. [in pb]

Cheers,

Benjamin

> On 15 Jun 2017, at 11:53 , Min He <email address hidden> wrote:
>
> New question #643895 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Dear MadAnalysis team,
>
> I am using MadGraph v2.5.5 to generate some events. The process card is shown below.
>
> generate e+ e- > a mu+ mu-
> output
> launch
> shower=PYTHIA8
> detector=DELPHES
> analysis=MADANALYSIS_5
> set ebeam1 120
> set ebeam2 120
>
> The parton level analysis is OK. http://202.120.47.101:23334/PROC_sm_2/HTML/run_01/tag_1_MA5_PARTON_ANALYSIS_analysis1/HTML/index.html
> But I got zero events at the hadron level analysis. http://202.120.47.101:23334/PROC_sm_2/HTML/run_01/tag_1_MA5_HADRON_ANALYSIS_analysis2_BasicReco/HTML/index.html
>
> It seems the cross section information is lost at the hadron level. I think I have the same problem with this: https://answers.launchpad.net/madanalysis5/+question/627403
> I want to know how to solve this problem.
>
> Thanks!
>
> Min
>
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Min He (hemin-sjtu) said : #2

Hi Benjamin,

This time I added what you said to the script but I still get zero events at hadron level.
http://202.120.47.101:23334/PROC_sm_3/HTML/run_01/tag_1_MA5_HADRON_ANALYSIS_analysis2_BasicReco/HTML/index.html

generate e+ e- > a mu+ mu-
output
launch
shower=PYTHIA8
detector=DELPHES
analysis=MADANALYSIS_5
set ebeam1 120
set ebeam2 120
set tag_1_pythia8_BasicReco.xsection = 0.4334

while the dataset is named as tag_1_pythia8_BasicReco.lhe.gz

Thanks!

Min

Benjamin Fuks (fuks) said : #3

Hi,

I don't have access to the webpage…

Could you please check the hepmc flie to see whether it is full or empty, as well as the reco-lhe file? Thanks in advance.

Cheers,

Benjamin

> On 15 Jun 2017, at 15:17 , Min He <email address hidden> wrote:
>
> Question #643895 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Status: Answered => Open
>
> Min He is still having a problem:
> Hi Benjamin,
>
> This time I added what you said to the script but I still get zero events at hadron level.
> http://202.120.47.101:23334/PROC_sm_3/HTML/run_01/tag_1_MA5_HADRON_ANALYSIS_analysis2_BasicReco/HTML/index.html
>
> generate e+ e- > a mu+ mu-
> output
> launch
> shower=PYTHIA8
> detector=DELPHES
> analysis=MADANALYSIS_5
> set ebeam1 120
> set ebeam2 120
> set tag_1_pythia8_BasicReco.xsection = 0.4334
>
> while the dataset is named as tag_1_pythia8_BasicReco.lhe.gz
>
> Thanks!
>
> Min
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Min He (hemin-sjtu) said : #4

Hi Benjamin,

Both files are full of data.

Here are some lines from the hepmc file:

HepMC::Version 2.06.09
HepMC::IO_GenEvent-START_EVENT_LISTING
E 0 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 6 1 2 0 1 4.3341590000000005e-05
N 1 "Weight_MERGING=0.000"
U GEV MM
C 3.4673272000000004e-04 5.4176987500000003e-02
V -1 0 0 0 0 0 1 1 0
P 1 -11 0 0 1.1999999999891199e+02 1.1999999999999999e+02 5.1099999999999995e-04 4 0 0 -1 0
P 3 -11 0 0 1.1999999999991546e+02 1.1999999999991546e+02 0 21 0 0 -3 0

Thanks!

Min

Benjamin Fuks (fuks) said : #5

Hi,

Do you mind generating a file with maybe 100 events and share it? I cannot reproduce your issue.

Cheers,

Benj

> On 15 Jun 2017, at 19:03 , Min He <email address hidden> wrote:
>
> Question #643895 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Status: Answered => Open
>
> Min He is still having a problem:
> Hi Benjamin,
>
> Both files are full of data.
>
> Here are some lines from the hepmc file:
>
> HepMC::Version 2.06.09
> HepMC::IO_GenEvent-START_EVENT_LISTING
> E 0 -1 -1.0000000000000000e+00 -1.0000000000000000e+00 -1.0000000000000000e+00 0 0 6 1 2 0 1 4.3341590000000005e-05
> N 1 "Weight_MERGING=0.000"
> U GEV MM
> C 3.4673272000000004e-04 5.4176987500000003e-02
> V -1 0 0 0 0 0 1 1 0
> P 1 -11 0 0 1.1999999999891199e+02 1.1999999999999999e+02 5.1099999999999995e-04 4 0 0 -1 0
> P 3 -11 0 0 1.1999999999991546e+02 1.1999999999991546e+02 0 21 0 0 -3 0
>
> Thanks!
>
> Min
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Min He (hemin-sjtu) said : #6

Hi,

Here are these two files.
https://www.dropbox.com/sh/dhjfnvsuasg1xy1/AAB0Mt_rBPP3WQqlVibgmQvPa?dl=0

I generated 10000 events using script shown below.

generate e+ e- > a mu+ mu-
output
launch
shower=PYTHIA8
detector=DELPHES
analysis=MADANALYSIS_5
set ebeam1 120
set ebeam2 120
set tag_1_pythia8_BasicReco.xsection = 0.4334

Min

Benjamin Fuks (fuks) said : #7

Hi Min,

I cannot reproduce your bug. MA5 reads your events and display them correctly. Do you mind sharing the report?

Cheers,

Benjamin

> On 16 Jun 2017, at 10:44 , Min He <email address hidden> wrote:
>
> Question #643895 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Status: Answered => Open
>
> Min He is still having a problem:
> Hi,
>
> Here are these two files.
> https://www.dropbox.com/sh/dhjfnvsuasg1xy1/AAB0Mt_rBPP3WQqlVibgmQvPa?dl=0
>
> I generated 10000 events using script shown below.
>
> generate e+ e- > a mu+ mu-
> output
> launch
> shower=PYTHIA8
> detector=DELPHES
> analysis=MADANALYSIS_5
> set ebeam1 120
> set ebeam2 120
> set tag_1_pythia8_BasicReco.xsection = 0.4334
>
> Min
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Min He (hemin-sjtu) said : #8

Hi Benjamin,

I added the report to this dropbox link.
https://www.dropbox.com/sh/dhjfnvsuasg1xy1/AAB0Mt_rBPP3WQqlVibgmQvPa?dl=0

Thanks!

Min

Benjamin Fuks (fuks) said : #9

But those reports are fine… Histograms are not empty...

> On 16 Jun 2017, at 11:35 , Min He <email address hidden> wrote:
>
> Question #643895 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Status: Answered => Open
>
> Min He is still having a problem:
> Hi Benjamin,
>
> I added the report to this dropbox link.
> https://www.dropbox.com/sh/dhjfnvsuasg1xy1/AAB0Mt_rBPP3WQqlVibgmQvPa?dl=0
>
> Thanks!
>
> Min
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Min He (hemin-sjtu) said : #10

Hi Benjamin,

For the parton level analysis report(Parton.pdf), it seems OK. But for the hadron level analysis report(Hadron.pdf), in the section Datasets, we can see

2.1
 tag_1_pythia8_basicreco
• Sample consisting of: signal events.
• Generated events: 3202 events.
• Normalization to the luminosity: 0+/- 0 events.
• Ratio (event weight): 0.0 .

The events number is zero, and the histograms are empty if there is a photon or a muon. I do not understand why this happened.

Thanks!

Min

Benjamin Fuks (fuks) said : #11

Hi Min,

No, the event number is not zero:
  > • Generated events: 3202 events.
The cross section is probably set to zero, which means that each of the 3202 events has a zero weight.

Cheers,

Benjamin

> On 19 Jun 2017, at 17:38 , Min He <email address hidden> wrote:
>
> Question #643895 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/643895
>
> Status: Answered => Open
>
> Min He is still having a problem:
> Hi Benjamin,
>
> For the parton level analysis report(Parton.pdf), it seems OK. But for
> the hadron level analysis report(Hadron.pdf), in the section Datasets,
> we can see
>
> 2.1
> tag_1_pythia8_basicreco
> • Sample consisting of: signal events.
> • Generated events: 3202 events.
> • Normalization to the luminosity: 0+/- 0 events.
> • Ratio (event weight): 0.0 .
>
> The events number is zero, and the histograms are empty if there is a
> photon or a muon. I do not understand why this happened.
>
> Thanks!
>
> Min
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Can you help with this problem?

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

To post a message you must log in.