Processing less no. of events while generating root file

Asked by shibasipu on 2017-08-14

Hi Benjamin,
                         I have generated 1 Lakh events using MG5+PY8 chain. But while generating the root file in MA5, it is processing less no. of events (i.e. 2*46874=93748). It is not showing any error message. For your reference, I am coping the screen output here.

==========================================
      with the MA5-tuned card: ../../Inputma5_cms.tcl
    * Running over files ...
    * 1/2 /home/shibananda/pp2WW4b_PY8ON_14_2Lakh_06August2017_iseed7/Events/run_01_0/tag_1_pythia8_events.hepmc.gz
        => file size: 1.95 Go
        => sample format: HEPMC file produced by an unknown generator (cross section assumed in pb).
WARNING: --------------------------------------------------------------------------------
WARNING: Msg | Number of mothers greather than 2
WARNING: Details |
WARNING: Where | Function = FinalizeEvent ; File=Reader/HEPMCReader.cpp ; Line=170
WARNING: --------------------------------------------------------------------------------
        => progress: [=============================> ]Warning in <TVector3::PseudoRapidity>: transvers momentum = 0! return +/- 10e10
Warning in <TVector3::PseudoRapidity>: transvers momentum = 0! return +/- 10e10
        => progress: [===================================]
        => total number of events: 46874 ( analyzed: 46874 ; skipped: 0 )
    * 2/2 /home/shibananda/pp2WW4b_PY8ON_14_2Lakh_06August2017_iseed7/Events/run_01_1/tag_1_pythia8_events.hepmc.gz
        => file size: 1.95 Go
        => sample format: HEPMC file produced by an unknown generator (cross section assumed in pb).
WARNING: --------------------------------------------------------------------------------
WARNING: Msg | Number of mothers greather than 2
WARNING: Details |
WARNING: Where | Function = FinalizeEvent ; File=Reader/HEPMCReader.cpp ; Line=170
WARNING: --------------------------------------------------------------------------------
        => progress: [===================================]
        => total number of events: 46874 ( analyzed: 46874 ; skipped: 0 )
    * Finalizing all components ...
    * Total number of processed events: 93748.
+----------------------------------------------------------------------------------------------------------------------+
| LogReport-Warning |
+----------------------------------------------------------------------------------------------------------------------+
| Message NIterations @ File Line |
|----------------------------------------------------------------------------------------------------------------------|
| Number of mothers greather than 2 2 Reader/HEPMCReader.cpp 170 |
+----------------------------------------------------------------------------------------------------------------------+
    * Goodbye.
MA5: *******************************************************
MA5: Checking SampleAnalyzer output...
MA5: Extracting data from the output files...
MA5: Preparing data for the reports ...
MA5: Generating all plots ...
MA5: Generating the HMTL report ...
MA5: -> To open this HTML report, please type 'open'.
MA5: Generating the PDF report ...
MA5: -> To open this PDF report, please type 'open ANALYSIS_5/PDF'.
MA5: Generating the DVI report ...
MA5: Well done! Elapsed time = 2608 seconds
===============================================================
Please let me know where could be the problems.

Regards,
Shiba

Question information

Language:
English Edit question
Status:
Answered
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Last query:
2017-08-16
Last reply:
2017-08-16
shibasipu (shibasipu) said : #1

Just to add the above thread,
 the file size of each hepmc file is 2.1 GB.

But MA5 showing it is 1.95 GB.
=====================
* 1/2 /home/shibananda/pp2WW4b_PY8ON_14_2Lakh_06August2017_iseed7/Events/run_01_0/tag_1_pythia8_events.hepmc.gz
        => file size: 1.95 Go
====================================

Regards,
Shiba

Benjamin Fuks (fuks) said : #2

Hi Shiba,

Have you checked the number oif events included in the hepmc file? (check the last line starting with "E xxxxx" where xxxxx is a number). I do not think ma5 would be the problem here.

Cheers,

Benjamin

> On 14 Aug 2017, at 12:03 , shibasipu <email address hidden> wrote:
>
> Question #655867 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/655867
>
> shibasipu gave more information on the question:
> Just to add the above thread,
> the file size of each hepmc file is 2.1 GB.
>
> But MA5 showing it is 1.95 GB.
> =====================
> * 1/2 /home/shibananda/pp2WW4b_PY8ON_14_2Lakh_06August2017_iseed7/Events/run_01_0/tag_1_pythia8_events.hepmc.gz
> => file size: 1.95 Go
> ====================================
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

shibasipu (shibasipu) said : #3

Hi Benj,
                 Actually after extracting the hepmc file, the file size is too large around 6 GB. So, when I am opening any editor, the system is getting hanged. Is there any other way to check the above ?

Regards,
Shiba

Benjamin Fuks (fuks) said : #4

Bash commands?

Cheers,

Benjamin

> On 16 Aug 2017, at 08:23 , shibasipu <email address hidden> wrote:
>
> Question #655867 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/655867
>
> Status: Answered => Open
>
> shibasipu is still having a problem:
> Hi Benj,
> Actually after extracting the hepmc file, the file size is too large around 6 GB. So, when I am opening any editor, the system is getting hanged. Is there any other way to check the above ?
>
>
>
>
>
> Regards,
> Shiba
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

shibasipu (shibasipu) said : #5

thanks. sorry. I was checking in VI editor i.e. vi filename.

Regards,
Shiba

Manfred Hampl (m-hampl) said : #6

Remarks:

1. The size difference may be caused by different multipliers (base 1000 or base 2^10 = 1024)
1.95 * 1.024 * 1.024 * 1.024 = 2.09
1.95 GiB = approx. 2.1 GB

2. You could use grep to display "... the last line starting with "E xxxxx" where xxxxx is a number ...":
Something like
grep '^E ' tag_1_pythia8_events.hepmc | tail

Can you help with this problem?

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

To post a message you must log in.