Sorry there was a typo in my previous Email. In order to set hadron/muon thresholds to 1 keV, the PART-THR cards must be:
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
PART-THR -0.000001 -6.0 _at_LASTPAR
PART-THR -1.0E-14 NEUTRON
By mistake I had set it to 100 keV in my previous Email and not to 1 keV (thanks to Francesco for noticing).
Cheers, Anton
________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Anton Lechner [Anton.Lechner_at_cern.ch]
Sent: 02 February 2017 12:07
To: Kalista, Sebastian; fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: RE: USRBDX Flair
Dear Sebastian,
Thanks for the input file. Please be so kind and always reply with fluka-discuss in cc such that other users can also profit from the answers.
In order to plot USRBDX spectra with Flair, you need to do the following:
1) You need to request unformatted data and not formatted, i.e. you need to specify "-21." on WHAT(3) instead of "21.".
2) Then you need to convert the binary file to ascii files in the data frame of Flair (in the Run tab). Alternatively, you can also do this conversion on the command line using the $FLUPRO/flutil/usxsuw command (see Section 2.5.1. of the FLUKA manual for details). This conversion will give you 3 files (*.bnx, *_tab.lis, *_sum.lis), the two latter being in human readable format. Note, that the format is different from the one you get if you directly select formatted data on the USRBDX card.
3) Under the PLOT tab in Flair, you can use the *_tab.lis file to plot the spectrum under USR-1D.
For illustration, I attach a screenshot of the plot I got with your input file by following the above procedure. You might want to adjust the binning of the spectrum (I would rather go for a logarithmic energy binning on the USRBDX card).
Another important remark: you have a beam energy of 10 keV, but PRECISIO imposes by a default a cut-off energy of 100 keV for hadrons. The reason you still see pions crossing the boundary is because you have a very thin material layer and there is some sub-cut ranging out implemented for hadrons in FLUKA. For a correct simulation you need to reduce the threshold with the PART-THR card:
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
PART-THR -0.0001 -6.0 _at_LASTPAR
PART-THR -1.0E-14 NEUTRON
I use 1 keV, which is more or less the lower advisable energy cut. Since the first line overwrites the the threshold for neutrons, I reset it in the second line to the neutron default threshold. Note that the PART-THR only acts on hadrons and muons, while electron/positron and photon production/transport thresholds can be set via the EMFCUT card.
Cheers, Anton
________________________________
From: Kalista, Sebastian [Sebastian.Kalista_at_oeaw.ac.at]
Sent: 01 February 2017 08:10
To: Anton Lechner
Subject: AW: USRBDX Flair
Dear Anton,
Thank you for your quick answer.
TITLE
Antiproton on foil
* Set the defaults for precision simulations
DEFAULTS PRECISIO
* Antiprotonen
BEAM -1E-05 APROTON
* Define the beam position
BEAMPOS -10.
GEOBEGIN COMBNAME
0 0
* Black body
SPH blkbody 0.0 0.0 0.0 100000.0
* Void sphere
SPH void 0.0 0.0 0.0 10000.0
* Cylindrical target
RPP target -10.0 10.0 -10.0 10.0 -0.5E-03 0.5E-03
* Cylindrical target
RPP target1 -10.0 10.0 -10.0 10.0 0.5E-03 1E-03
END
* Black hole
BLKBODY 5 +blkbody -void
* Void around
VOID 5 +void -target-target1
* Target
TARGET 5 +target
TARGET2 5 +target1
END
GEOEND
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA BLCKHOLE BLKBODY
ASSIGNMA VACUUM VOID
ASSIGNMA ALUMINUM TARGET
ASSIGNMA ALUMINUM TARGET2
USRBDX 101. PIONS+- 21. TARGET TARGET2 100.blabla
USRBDX 5. 0. &
USRBIN 10. APROTON -50. 20.0 20.0 20.0Enedepap
USRBIN -20.0 -20.0 -20.0 100. 100. 100. &
USRBIN 10. PIONS+- -50. 20.0 20.0 20.0Enedeppi
USRBIN -20.0 -20.0 -20.0 100. 100. 100. &
USRBIN 10. PROTON -50. 20.0 20.0 20.0Enedepp
USRBIN -20.0 -20.0 -20.0 100. 100. 100. &
USRBIN 10. DEUTERON -50. 20.0 20.0 20.0Enedepd
USRBIN -20.0 -20.0 -20.0 100. 100. 100. &
USRBIN 10. 4-HELIUM -50. 20.0 20.0 20.0Enedephe
USRBIN -20.0 -20.0 -20.0 100. 100. 100. &
* Set the random number seed
RANDOMIZ 1.0
* Set the number of primary histories to be simulated in the run
START 10000.
STOP
Sebastian
________________________________
Von: Anton Lechner [Anton.Lechner_at_cern.ch]
Gesendet: Mittwoch, 1. Februar 2017 07:39
An: Kalista, Sebastian; fluka-discuss_at_fluka.org
Betreff: RE: USRBDX Flair
Dear Sebastian,
Could you be so kind and provide us the whole input file.
Thanks!
Anton
________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Kalista, Sebastian [Sebastian.Kalista_at_oeaw.ac.at]
Sent: 31 January 2017 21:03
To: fluka-discuss_at_fluka.org
Subject: [fluka-discuss]: USRBDX Flair
Dear FLUKA users and experts,
I am currently working on an antiproton annihilation experiment and I am trying to simulate the experimental setup. For my first few steps I chose an easy set-up and started implementing it in FLAIR. I am interested in the fluence of different particles through a surface, so i chose the USRBDX detector. To be exact antiprotons impinging on an aluminum foil.
I defined two regions TARGET and TARGET2 (both thin aluminum foils touching each other in the x-y-plane) and implemented the detector as follows:
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....
USRBDX 101. PIONS+- 21. TARGET TARGET2 100.blabla
USRBDX 5. 0. &
I am not sure where the problem is, but it is impossible for me to process the outcome and to plot it. Does anyone know the solution or has similar problems?
Best Regards
Sebastian Kalista
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Feb 02 2017 - 15:51:19 CET