janmayer
commented
10 months ago Hi,
writing here is fine. I have not worked on the topic for quite a while now.
I used only Neuland in the simulation (without 4mm steel window as mentioned in the Fig. description, btw. where did you put it and why?)
This geometry is quite important. The 4mm steel represents the exit window of the reaction chamber (its the only part of the reaction chamber the neutrons should see, so as long as you are not studying beam induced background, it should suffice to only have this part of the chamber there. Also, the air between the reaction chamber and the detector is slightly important.
and QGSP_INCLXX_EMV physics list.
The physics list should not be that important, but you probably want to test that. Note that cutoff energies (in the simulation, different from PMT thresholds) are quite important.
You should be able to reproduce everything using this repository - run the notebooks one by one. (There still might be changes to to newer software versions.)
Best, Jan
Hello Jan,
I have been trying out calorimetric method for neutron reconstruction and for start I wanted to reproduce one of your results (one in the center; 600 MeV) mult12-5.pdf . Below I copy my results which turned out quite different. I used only Neuland in the simulation (without 4mm steel window as mentioned in the Fig. description, btw. where did you put it and why?) and QGSP_INCLXX_EMV physics list. Other things were as in the Fig. description. I tried to change (increase) PMT threshold since my efficiencies seem to be too high but it didn't change much. Also, I tried to change initial parameters for the calorimetric cuts but it didn't change much (on a level of 1%). Do you have an idea where could the differences come from? Did you use some additional conditions which are not explicitly described here?
My results:
Number of entries in chain: 50000 Number of generated events with 1n is 10 0000. For 1n generated, 0 detected = 21.47 For 1n generated, 1 detected = 65.08 For 1n generated, 2 detected = 13.34 ................................ Number of generated events with 2n is 10 0000. For 2n generated, 0 detected = 4.46 For 2n generated, 1 detected = 44.52 For 2n generated, 2 detected = 30.92 For 2n generated, 3 detected = 16.66 ................................ Number of generated events with 3n is 10 0000. For 3n generated, 0 detected = 0.84 For 3n generated, 1 detected = 21.3 For 3n generated, 2 detected = 28.07 For 3n generated, 3 detected = 27.72 For 3n generated, 4 detected = 15.36 ................................ Number of generated events with 4n is 10 0000. For 4n generated, 0 detected = 0.18 For 4n generated, 1 detected = 9.34 For 4n generated, 2 detected = 17.06 For 4n generated, 3 detected = 24.56 For 4n generated, 4 detected = 23.96 For 4n generated, 5 detected = 24.9 ................................ Number of generated events with 5n is 10 0000. For 5n generated, 0 detected = 0.03 For 5n generated, 1 detected = 3.59 For 5n generated, 2 detected = 9.42 For 5n generated, 3 detected = 17.8 For 5n generated, 4 detected = 22.17 For 5n generated, 5 detected = 46.99
Thank you, Ivana
P.S. I hope it is fine to write here, if you would prefer to continue discussion in email or elsewhere, please let me know.