RE: [fluka-discuss]: Flux/fluence

From: Andrea Fontana <andrea.fontana_at_pv.infn.it>
Date: Tue, 14 Aug 2018 17:18:22 +0200 (CEST)

Dear Mihaela,
     to have realistic results, please make sure that you enable also
the PHOTONUC and MUPHOTON cards. Neutrons will also be produced
through photonuclear interactions by photons whenever the muons
produce electromagnetic showers.

Perhaps this paper can be useful:
https://journals.aps.org/prd/abstract/10.1103/PhysRevD.64.013012

Andrea

> Dear Andrea,
>
> Now the program is working. I have one more question. When I use USRTRACK
> for muons and neutrons what I obtain is { d Phi } over { dE }. What seems
> weird is the fact that I get three orders of magnitude more neutrons than
> muons. Maybe I didn't understand correcty the output of this card...? Or
> these are the stopped particles?
> ________________________________________
> From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on
> behalf of Mihaela Parvu [mihaela.parvu_at_cern.ch]
> Sent: 13 August 2018 15:48
> To: andrea.fontana_at_pv.infn.it
> Cc: andrea.fontana_at_pv.infn.it; fluka-discuss_at_fluka.org
> Subject: RE: [fluka-discuss]: Flux/fluence
>
> I do not think that this is the problem because I made a new folder with a
> new .flair file.
> ________________________________________
> From: Andrea Fontana [andrea.fontana_at_pv.infn.it]
> Sent: 13 August 2018 14:43
> To: Mihaela Parvu
> Cc: andrea.fontana_at_pv.infn.it; fluka-discuss_at_fluka.org
> Subject: RE: [fluka-discuss]: Flux/fluence
>
> Hi Mihaela,
> the error is simply due to the fact that you are using the
> renamed .inp file (MIN_v2.inp), with the old MIN.flair file
> (which still try to load the old file MIN.inp).
>
> Simply remove the "_v2" from the filename and use MIN.inp
> or create a new .flair file by reading the .inp from an empty
> Flair project.
>
> Andrea
>
>> Dear Andrea,
>>
>> Thank you so much for your help. Unfortunately, now I get this plot
>> error.
>> It's like I do not have neutrons, but it is impossible. Also, When I try
>> to plot the geometry I get another error:
>>
>> plot 'MIN_v2_geometry.dat' ind 0 us 1:2 w l ls 1 notitle
>> <e< line 0: warning: Skipping data file with no valid points
>>
>> NEUTRONS USRBDX
>>
>> Started 2018.08.13 14:19:31
>> reset
>> set terminal wxt 0 enhanced raise title 'MIN_try1.flair-0: USR-1D
>> ../MIN
>> v2 31'
>>
>> set title 'USR-1D /MIN v2 31'
>> unset grid
>> set xlabel ''
>> set xtics
>> set ylabel ''
>> set ytics
>> set logscale x
>> set logscale y
>> unset logscale z
>> unset logscale cb
>> unset logscale x2
>> unset logscale y2
>> set key default
>> plot '../modified_CERN/MIN_v2_31_tab.lis' ind 0 us 1:3 w steps lw 2 lc
>> rgb
>> 'red' t 'Detector 1','' ind 0 us (sqrt($1*$2)):3:($3*$4/100.) w
>> errorbars
>> lw 2 lc rgb 'red' pt 0 ps 1 notitle
>> <e< gnuplot> plot '../modified_CERN/MIN_v2_31_tab.lis' ind 0 us 1:3 w
>> steps lw 2 lc rgb 'red' t 'Detector 1','' ind 0 us
>> (sqrt($1*$2)):3:($3*$4/100.) w errorbars lw 2 lc rgb 'red' pt 0 ps 1
>> notitle
>>
>> <e<
>>
>> ^
>>
>> <e< line 0: all points y value undefined!
>>
>> Gnuplot error
>>
>> ________________________________________
>> From: Andrea Fontana [andrea.fontana_at_pv.infn.it]
>> Sent: 13 August 2018 13:46
>> To: Mihaela Parvu
>> Cc: fluka-discuss_at_fluka.org
>> Subject: Re: [fluka-discuss]: Flux/fluence
>>
>> Dear Mihaela,
>> I have run your example and I see that you have forgotten the
>> normalization to the area of the first USRBDX card (output 25).
>>
>> Moreover this normalization should read 360000, while the
>> normalization of your second USRBDX card should read 360720.36
>> (and not 360000), to take into account the 6mm thickness of the wall.
>>
>> Even with this normalization correctly done, you will still
>> see differences because of the muons interactions in the steel of
>> the wall and in the argon of the target. If you put both materials
>> to vacuum, you will get the same fluxes.
>>
>> In attachment a corrected version of your input (I have added two
>> USRBIN cards to better visualize your results).
>>
>> Perhaps this discussion is useful:
>>
>> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/9211.html
>>
>> Hope this helps,
>> Andrea
>>
>>> Dear Fluka experts,
>>>
>>> I am trying to simulate muon induced neutrons. For this I am using my
>>> own
>>> predefined source of muons with flux as a function of energy. When I
>>> use
>>> USRBDX, the fluence I am getting from VOID to wall does not coincide
>>> with
>>> the incident one.
>>> Can you please tell me help me? I have attached the input files
>>>
>>
>>
>> --
>> ========================================================================
>> Dr. Andrea Fontana tel: +39 0382 987991
>> Istituto Nazionale fax: +39 0382 423241
>> di Fisica Nucleare
>> Sezione di Pavia e-mail: andrea.fontana_at_pv.infn.it
>> Via Bassi 6 web : www.pv.infn.it/~fontana
>> 27100 PAVIA, Italy
>> ========================================================================
>>
>>
>>
>
>
> --
> ========================================================================
> Dr. Andrea Fontana tel: +39 0382 987991
> Istituto Nazionale fax: +39 0382 423241
> di Fisica Nucleare
> Sezione di Pavia e-mail: andrea.fontana_at_pv.infn.it
> Via Bassi 6 web : www.pv.infn.it/~fontana
> 27100 PAVIA, Italy
> ========================================================================
>
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id?c_info
>
>
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at
> https://www.fluka.org/fluka.php?id=acc_info
>
>


-- 
========================================================================
Dr. Andrea Fontana                    tel: +39 0382 987991
Istituto Nazionale                    fax: +39 0382 423241
di Fisica Nucleare
Sezione di Pavia                      e-mail: andrea.fontana_at_pv.infn.it
Via Bassi 6                           web   : www.pv.infn.it/~fontana
27100 PAVIA, Italy
========================================================================
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Aug 14 2018 - 18:38:20 CEST

This archive was generated by hypermail 2.3.0 : Tue Aug 14 2018 - 18:38:22 CEST