[fluka-discuss]: Antiproton secondaries, Energy distribution

From: Kalista, Sebastian <Sebastian.Kalista_at_oeaw.ac.at>
Date: Tue, 21 Feb 2017 09:37:11 +0000

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. 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). I am now interested in two plots:

1. The number of the produced secondaries and their energy. For example number of Pions on the y-axis and energy on the x-axis.
2. The angle between secondaries and the incoming beam direction. Number of lets say Pions again on the y-axis and on the x-axis the angle.

Here is my current input file:

TITLE
Antiproton on foil
* Set the defaults for precision simulations
DEFAULTS PRECISIO
PART-THR -1E-06 4-HELIUM _at_LASTPAR
PART-THR -1.0E-14 NEUTRON
* 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 98. 4-HELIUM -21. TARGET TARGET2 100.blabla
USRBDX 100. 1E-05 700. &
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

So for the first plot i tried to use an USRBDX detector. But i am not quite sure what the y-axis tells me and how i can convert it from fluence to total particle number.

Thank you guys in advance for your help
Sebastian


__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Tue Feb 21 2017 - 11:54:51 CET

This archive was generated by hypermail 2.3.0 : Tue Feb 21 2017 - 11:54:56 CET