Some of LandFirstHits are empty - Neutrons simulation. [message #18523] |
Mon, 21 September 2015 11:18 |
Dmytro Kresan
Messages: 166 Registered: June 2004
|
first-grade participant |
From: *gsi.de
|
|
Dear Christiaan,
The answer to your problem with fraction of NeuLAND First Hits being not filled, when performing a neutron simulation.
I have also observed this issue.
The first thing I did - place only NeuLAND geometry inside of cave filled with vacuum. The fraction of empty data was reduced. Obviously, there are scattered neutrons on the way to NeuLAND, which are not considered as primaries.
Second - I made the wrapping volumes also active. This has shown that all FirstHits are recorded.
Explanation: stepping function is called only for active volumes (scintillators), and not for passive ones.
Best regards,
Dima
|
|
|
|
Re: Some of LandFirstHits are empty - Neutrons simulation. [message #18525 is a reply to message #18524] |
Mon, 21 September 2015 13:09 |
Dmytro Kresan
Messages: 166 Registered: June 2004
|
first-grade participant |
From: *gsi.de
|
|
It is defined in function CheckIfSensitive() - file R3BLand.cxx, line number 451. You have to implement the check whether string volName contains "wrapping", and return kTRUE in this case.
NOTE: this will produce LandPoints also for wrapping volumes, which does not really make sense. To avoid it: implement at line 151 (after call to StepHistory()) check whether volume name contains "wrapping". If yes - make return kTRUE. This will store FirstHits but not the points.
Best regards,
Dima
|
|
|