Message: Ntuple fills per-event instead of per-hit Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question Ntuple fills per-event instead of per-hit 

Forum: Analysis
Date: 09 Apr, 2013
From: Matt Taggart <Matt Taggart>


I've been working on producing a simulation of spontaneous fission based on the Livermore fission model built on top of the B4d example structure.

The simulation works as expected, however, when filling the ntuple to analyse the output using root I find that I am getting only a sample of the total data. It appears that the ntuple is being filled only once per fission event, rather than per particle/gamma. For my simulation of 10,000 fissions I should expect (and indeed see when reading out the neutron/gamma emission) ~37k neutrons and ~82k gammas, but each leaf has exactly 10,000 entries.

I believe the problem is caused by my implementation of the AddNtupleRow() method, as the B4d example includes this in the EventAction class file (thus explaining the "once per event" issue). The issue arises where my primary neutrons and gammas are built within separate loops determined by n/g multiplicity generated by the fission. Including AddNtupleRow at this point allows the extraction of particle energy etc. as desired, but introduces a multiple count when including the method in each loop. I also wish to extract hit data, which is called in my EventAction class (again, only occurs once per event at present).

I wonder whether the analysis from B4d is not sophisticated enough to fill different sized leaves? In which case is there a more complete/developed example that could be cannibalised? Any help/suggestions would be greatly appreciated.

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: Ntuple fills per-event instead of per-hit   (Guy Barrand - 09 Apr, 2013)
2 None: Re: Ntuple fills per-event instead of per-hit   (Ivana Hrivnacova - 09 Apr, 2013)
3 None: Re: Ntuple fills per-event instead of per-hit   (Matt Taggart - 09 Apr, 2013)
1 None: Re: Ntuple fills per-event instead of per-hit   (Ivana Hrivnacova - 11 Apr, 2013)
 Add Message Add Message
to: "Ntuple fills per-event instead of per-hit"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews

[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]