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