GSI Forum
GSI Helmholtzzentrum für Schwerionenforschung

Home » R3BRoot » Simulation Issues » Deviation in multi-neutron simulations (Deviation in multi-neutron simulations)
Deviation in multi-neutron simulations [message #19268] Tue, 24 May 2016 14:56 Go to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
Dear Mr. Kresan,

I am trying to reproduce your results on multi-neutron simulations from your page https://www.r3broot.gsi.de/multi-neutron-simulation

I follow the steps indicated by r3blandsim.C and precalibr.C: First I run Monte Carlo simulations for 1-5 neutrons on 600 MeV per particle (I used r3ball.C)
I fire the particles from the origin with a box generator with Ptot = 600 MeV + rest energy (consistent with relativity), theta = 0 - 80 mrad, phi = 0-360 deg.
I only have NeuLAND at 13 m from the origin and no other detectors in the cave. Then I run the digitizer, cluster finder and R3B pre-calibrator as indicated by precalibr.C
(I use calibr_mini.C to compute the cuts).

However, my histograms do not look at all like your histograms. I do not mean the efficiencies and the cuts, I mean the shape and color of the blobs themselves.
Initially I thought the difference was due to the use of land/neuland digitizer, so then I decided to test both of them (50000 events), but the histograms are
nearly identical (land digitizer is just a lot slower).

What I find even more strange is that the picture of your histograms on https://www.r3broot.gsi.de/multi-neutron-simulation
also does not look like the histograms on page 56 of the NeuLAND technical design report. And the pictures on Jan Mayer his talk of March 1st (slide 7) look again
different. I do not mean the cuts, efficiencies or even the amount of red in the blobs (this just depends on how much you zoom in the z-cirection). I mean the x and y location of the blob center and the relative pattern of the colors and of the dispersion of the hits. This is somewhat different in all distinct cases.

Can you explain the nature and reasons for these differences to me? And maybe help me to improve my histograms?

Yours sincerely,
Christiaan Douma.

PS: The LAND_Digitizer.png picture I will upload in a response since I can only upload one picture at a time.
Re: Deviation in multi-neutron simulations [message #19269 is a reply to message #19268] Tue, 24 May 2016 15:00 Go to previous messageGo to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
second picture.

Christiaan.
Re: Deviation in multi-neutron simulations [message #19270 is a reply to message #19269] Tue, 24 May 2016 15:11 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
Dear Christiaan,

can you please repeat your simulations with Geant3 and input ASCII files (link on the https://www.r3broot.gsi.de/multi-neutron-simulation). Then we can validate your results. Treatment of neutrons in Geant3 has better agreement with experimental results.

Best regards,
Dima
Re: Deviation in multi-neutron simulations [message #19271 is a reply to message #19270] Tue, 24 May 2016 15:45 Go to previous messageGo to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
I already used Geant3 for my simulations. I can repeat my simulations with your ascii input file (And I will),
but this still does not explain the differences between the results of your web page and the TDR and Jan Mayer his results.

Christiaan.
Re: Deviation in multi-neutron simulations [message #19272 is a reply to message #19271] Tue, 24 May 2016 15:53 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
Status at the time of TDR: in order to have an agreement with old LAND experiments, we introduced 0.4 factor scaling to the deposited energy/light per MC particle.

Afterwards, in the simulations of Jan Mayer, we removed this factor, and introduced PMT saturation model together with final QDC resolution of 5%. Maybe Jan can comment on status of agreement of this stage with experimental data.

This is where the differences in deposited energy are coming from.

And picture on the web is just an example, do not take it as quantitative recommendation. If I remember correctly, when I produced it - there was no PMT saturation and no scaling factor.

Best regards,
Dima
Re: Deviation in multi-neutron simulations [message #19273 is a reply to message #19272] Tue, 24 May 2016 16:25 Go to previous messageGo to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
Thank you for this explanation. It clarifies much to know which differences there are
in the software between the pictures. In order to get a good agreement with your picture
on the web, I am now running the simulations with your ascii-input files. From initial
view I can already see that the angular spread of the neutrons in your ascii-files is
much smaller than with me. So for the complete picture for me: can you explain precisely which
distributions in energy, momentum, etc. you used for the neutrons to generate your pictures
on the web?

PS: I am using Jan Mayer his NeuLAND digitizer for this run, so that would mean that I have the
PMT saturation model together with final QDC resolution of 5% in the simulation.

Christiaan.
Re: Deviation in multi-neutron simulations [message #19279 is a reply to message #19273] Wed, 25 May 2016 09:45 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
ASCII files for multi-neutron simulation contain realistic distribution from interaction of 132Sn beam on a hydrogen target. The angular distribution of neutrons and a fragment is calculated assuming 500 keV relative energy between this particles - in transversal direction. This relative energy is delta-function like - to test reconstruction resolution. Longitudinal momentum comes from beam energy of 600 MeV.

Yes, digitiser of Jan contains saturation and 5% gaussian smearing.

Best regards,
Dima
Re: Deviation in multi-neutron simulations [message #19282 is a reply to message #19279] Wed, 25 May 2016 11:03 Go to previous messageGo to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
With the ASCII-files (600 MeV) from the web page I still do not get the same as (or close to) your pictures. Do you have any Idea what I'm doing wrong?
I just use r3ball.C (2500 Ampere ALADIN field) and then I use Jan Mayer his digitizer, then the clusterfinder and then NeutronCalibr2D.

Also about the distribution of the neutrons: I understand the concept of a 'realistic distribution', but do you know what this distribution approximately is?
Is it uniform in angular spread or Gaussian, etc?

Christiaan.
Re: Deviation in multi-neutron simulations [message #19284 is a reply to message #19282] Wed, 25 May 2016 11:46 Go to previous messageGo to next message
Dmytro Kresan is currently offline  Dmytro Kresan
Messages: 166
Registered: June 2004
first-grade participant
From: *gsi.de
I did not generate the ASCII files myself. I wrote already everything I know about this input.

How do your pictures compare with Jan's results?

If you want to find what goes wrong, and compare quantitatively your plots with ones from r3broot web page, please use also the same digitizer, like in the step-by-step instructions on the web. Cluster finder algorithm depends on digis input it gets.

Best regards,
Dima
Re: Deviation in multi-neutron simulations [message #19286 is a reply to message #19284] Wed, 25 May 2016 11:48 Go to previous messageGo to next message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
The pictures do not reproduce Janhis results. I already tried to use the land digitizer instead of the neuland digitizer.
(not for the Sn-ASCII-files, but for my initial try with box generator). There is almost no difference in the histograms (I uploaded them).

Christiaan.
Re: Deviation in multi-neutron simulations [message #19288 is a reply to message #19268] Wed, 25 May 2016 14:31 Go to previous messageGo to next message
Jan Mayer is currently offline  Jan Mayer
Messages: 20
Registered: August 2015
occasional visitor
From: *ikp.uni-koeln.de
First let me say that the last set of Plots (NeutronCalibration_0p.png ) does look reasonable - positioning of the cuts is debatable, we had that topic before.

It is quite difficult to find the cause of deviations from afar. It is not the "fault" of the Digitizer.

I propose the following:
- Please go to R3BRoot/macros/r3b/neuland/matrices/
- In the *.sh files for Step 0, 1 and 3 in
for NDOUBLEPLANES in $(seq xx xx); do
replace the
seq(xx xx)
with
seq(30 30)
so you only do the simulations for 30dps
- Run all the scripts in order
- Post the result.
Re: Deviation in multi-neutron simulations [message #19314 is a reply to message #19284] Thu, 02 June 2016 11:58 Go to previous message
C. A. Douma is currently offline  C. A. Douma
Messages: 88
Registered: September 2015
Location: Groningen
continuous participant
From: *kvi-cart.rug.nl
Dear all,

Thanks for your help and suggestions. The problem is now solved.
In my simulations I did not use the standard orientation of NeuLAND. The digitizer was not equipped to
handle this. After adding a piece of code for this, the histograms look fine.

Christiaan.
Previous Topic: beta in the gamma generator in r3ball.C
Next Topic: Problem with TGeant4?
Goto Forum:
  


Current Time: Thu Dec 05 04:47:58 CET 2024

Total time taken to generate the page: 0.00719 seconds