I: Re: Re: [fluka-discuss]: I: 18 MeV Neutrons on Carbon

From: <francesca.belloni_at_libero.it>
Date: Thu, 16 Jul 2015 16:13:53 +0200 (CEST)

Dear Alberto, before replying I must have run the simulation with 0.1 GeV neutrons, or nowI messed up something. After lunch I decided to rmake some order: remove all the tests and restart from scratch.I downloaded the files test.inp and mgdraw.f from the link where I had uploaded them and changed only all the the LAM-BIAS values to -0.01,but now I still obtain only neutrons.My *.log file looks like: 1 6 8 0 300 1 0 0.0000000000000000 0.0000000000000000 4.1275814553889081E-002 1.0000000000000000
 2 6 8 1 300 8
 3 6 8 300 0 0
 1 21 8 0 300 1 0 0.0000000000000000 0.0000000000000000 2.0202138511377764E-002 1.0000000000000000
 2 21 8 1 300 8
 3 21 8 300 0 0
 1 63 8 0 300 1 0 0.0000000000000000 0.0000000000000000 9.0236660655996660E-003 1.0000000000000000
 2 63 8 1 300 8
 3 63 8 300 0 0
 1 71 8 0 300 1 0 0.0000000000000000 0.0000000000000000 4.5958651383185257E-002 1.0000000000000000
 2 71 8 1 300 8
 3 71 8 300 0 0
 1 89 8 0 300 1 0 0.0000000000000000 0.0000000000000000 1.7361631698946663E-002 1.0000000000000000
 2 89 8 1 300 8
 3 89 8 300 0 0
 1 89 8 0 300 1 0 -1.0320993451150469E-003 -2.3336882643914678E-002 1.6143825911649597E-002 1.0000000000000000
 2 89 8 1 300 8
 3 89 8 300 0 0
 1 95 8 0 300 1 0 0.0000000000000000 0.0000000000000000 1.0695069418869118E-002 1.0000000000000000
 2 95 8 1 300 8
 3 95 8 300 0 0

did I mess up something again,please? Sorry... regards,Francesca

----Messaggio originale----

Da: francesca.belloni_at_libero.it

Data: 16/07/2015 13.12

A: <fasso_at_slac.stanford.edu>, "fluka-discuss_at_fluka.org"<fluka-discuss_at_fluka.org>

Ogg: R: Re: Re: [fluka-discuss]: I: 18 MeV Neutrons on Carbon

Dear Alberto,

         thanks a lot for your help!
I was sincerly going in the opposite direction...I has started with a lam-bias of 10^-9 and since nothing was obtained I was
continuously increasing it.

I was not expecting antineutrons. Thought knowing it was not the case, I was trying to see if perhaps the lam-bias
for particle_id between -6 and +8 was meaning (-6,+8) ,.i.e. extremes not included, or [-6,+8], i.e. extremes included.
In the first case I thought to increase the id by one, to be sure neutrons were included.
Desperation leads me to silly thoughts :)

Thanks a lot for all the suggestions and help!


----Messaggio originale----

Da: fasso_at_slac.stanford.edu

Data: 16/07/2015 12.25

A: "fluka-discuss_at_fluka.org"<fluka-discuss_at_fluka.org>, "francesca.belloni_at_libero.it"<francesca.belloni_at_libero.it>

Ogg: Re: Re: [fluka-discuss]: I: 18 MeV Neutrons on Carbon


Dear Francesca,

there is one big issue in your input file (I confess that I have not even looked at the mgdraw routine, maybe I will have a look later).

In your LAM-BIAS command, you have set a biasing factor of 5.0E-12: that means that you have increased the cross section by a factor 2.E11!!!

If you look at the carbon cross section for 18 MeV neutrons, for instance on the Brookhaven NNDC site, you find that the total cross section

is about 1.4 b, and the inelastic cross section about 10 times lower. That means that if you set the biasing factor = 1.E-2 you have already

a 99% probability that the neutron will interact within the 5 mm of your target. That should be more than enough! Maybe you could even bias

a little less, say a biasing factor 0.05.

I did not find any other problem in your input file, but perhaps you could consider the following:

* I suggest that you could use DEFAULT PRECISION instead of NEW-DEFA

* command LOW-NEUT is not necessary. Everything is already there by default

* in the MATERIAL definitions, don't type the mass number (leave WHAT(2) blank). FLUKA has its own table of mass numbers and it is better not

   to interfere

* just a curiosity: why have you biased also the antineutrons? Do you expect any?

* you must be an old "customer", since in most of your input you use numbers instead of names. That was the way we used FLUKA many years

   ago, but now I have lost the habit and have had myself some difficulty in understanding your input. Try the names! They will make your life

   much easier :-)

Kind regards,


From: francesca.belloni_at_libero.it <francesca.belloni_at_libero.it>

Sent: Thursday, July 16, 2015 1:45 AM

To: Fasso, Alberto; fluka-discuss_at_fluka.org

Subject: R: Re: [fluka-discuss]: I: 18 MeV Neutrons on Carbon

Dear Alberto,

         I apologize for sending files this way. The browser I use at work is preventing me to attach files, and I

am not allowed to make any change on the computer I am using.

I had interfaced the FLUKA output with ROOT because I am terrible in fortran...in one

command line I insert 5 mistakes :)

A TFile is practically a folder and a TTree is a way to store data, like ntuple.

But you are right that the TTRee is not necesary (in my case it is even empty...).

I cleaned the input and put test.inp and mgdraw.f at the follwoing link


If i use instead than 18 MeV neutron, 0.1 GeV neutrons I obtain also residual nuclei....

Is the fact that I obtain practically nothing with 18 MeV neutrons linked to the low energy neutron treatment, please?

Or do i make some other basic mistake, please?

Thanks a lot,



----Messaggio originale----

Da: fasso_at_slac.stanford.edu

Data: 16/07/2015 10.19

A: "fluka-discuss_at_fluka.org"<fluka-discuss_at_fluka.org>, "francesca.belloni_at_libero.it"<francesca.belloni_at_libero.it>

Ogg: Re: [fluka-discuss]: I: 18 MeV Neutrons on Carbon

Dear Francesca,

please forgive a Paleolithic person...

1) What is a TTree? I looked it up on Google and it looks that it is something to do with ROOT. But please consider that not everybody is using ROOT! In particular,

it is not necessary to know ROOT for using FLUKA

2) I had a look at the link you sent, and I was able to download the following:

NeutCarb.tar.gz_files/ (a directory)


The stuff seems to gzipped, but the .gz is an extension in the middle of the file (or directory) name. How does one gunzip them?

The bottom line is: can you please provide input file and user routines in a less sophisticated way, if you want to be helped by Paleolithic people? :-):-):-)

Just two simple ASCII files would be ok...


From: owner-fluka-discuss_at_mi.infn.it <owner-fluka-discuss_at_mi.infn.it> on behalf of francesca.belloni_at_libero.it <francesca.belloni_at_libero.it>

Sent: Wednesday, July 15, 2015 11:29 AM

To: fluka-discuss_at_fluka.org

Subject: [fluka-discuss]: I: 18 MeV Neutrons on Carbon

Dear all,

     I made a mistake. The link is





----Messaggio originale----

Da: francesca.belloni_at_libero.it

Data: 15/07/2015 19.35

A: <fluka-discuss_at_fluka.org>

Ogg: 18 MeV Neutrons on Carbon

Dear all,

           I am trying to simulate the response of a diamond detector to a monochromatic pencil-like neutron beam.

I was expecting different reactions to take place and, as a consequence, to be able to score with mgdraw the

energy deposited by charged particles, for example alphas, in the carbon.

But I must be doing something completely wrong, since the energy deposited in the diamond is always 0

(not to get a large TTree for nothing, in the mgdraw routine I asked to dump the energy deposited only if it is higher than 0.

As a cosnequence nothing is written on the TTree).

To understand what'sgoing on, I edited also usdraw by following


and dump the results on the *.log file.

But I get always only neutrons.

I have problems in attaching files using the browser at work, therefore I uploaded the simulation at the link


(this time everyone should have direct access, at least for the download).

Could anyone give a look please?






To run the simulation:

-expand the file in a folder

-add to your .bashrc

        export LD_LIBRARY_PATH=$LD_LIBRARY_PATH:/folder_where_you_put_the_simulation

-source .bashrc

-digit make

-digit $FLUPRO/flutil/rfluka -e rootfluka -N0 -M5 test &amp;

You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Jul 16 2015 - 17:36:02 CEST

This archive was generated by hypermail 2.3.0 : Thu Jul 16 2015 - 17:36:05 CEST