More on heavy ions

From: Joel DeWitt <joel.m.dewitt@gmail.com>
Date: Tue Jul 01 2008 - 21:21:03 CEST
fluka-discuss:

The suggestions given by Francesco have been very helpful (http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1622.html), of which I am extremely grateful.  However, I have encountered a fundamental problem.  As an example, I have attached my input for 1 GeV/n 56-Fe on aluminum.  The resulting spectra obtained with USRYIELD do not integrate to 1.  Our experimental data is binned LET, normalized per unit ion per unit area.  Multiplying a fluence contribution by the area and summing all the contributions of course gives 1.  Perhaps I am handling the double differential yield inappropriately.

Suggestions will be most appreciated.

Joel DeWitt

TITLE
BB4 (NSRL-5) 956 MeV/n 56-Fe, Al
* Blackhole
ASSIGNMA BLCKHOLE BLKHOLE
* Air
MATERIAL 0.0012048 AIR
* Air
COMPOUND -9.256D-4 NITROGEN -2.837D-4 OXYGEN -1.572D-5 ARGONAIR
* Air
STERNHEI 10.5961 1.7418 4.2759 0.10914 3.3994 AIR
* Air
ASSIGNMA AIR VAC
DEFAULTS NEW-DEFA
* 956 MeV/n 56-Fe
BEAM -0.956 -0.082425 -1.7 10.0 10.0 1.0HEAVYION
HI-PROPE 26.0 56.0
* The user requests transport of heavy ions with nuclear interactions
EVENTYPE 2.0 DPMJET
* 1.815 g/cm2 air (E_i = 1 GeV/n, E_f = 956 MeV/n) from Dreamweaver: 1.815/1.2048E-3 = 1506.47 cm
BEAMPOS 0.0 0.0 -1506.47 0.0 0.0
GEOBEGIN COMBNAME
    0 0 BB4 (NSRL-5) 956 MeV/n 56-Fe, Al
SPH BLK 0.0 0.0 0.0 10000.0
RPP VOI -2000.0 2000.0 -2000.0 2000.0 -2000.0 2000.0
* CR-39, front (1.40 g/cm3, 0.05 cm)
RPP CR3900 -2.0 2.0 -2.0 2.0 0.0 0.05
* Al, 5 g/cm2 (2.6989 g/cm3, 1.85 cm)
RPP Al05 -5.0 5.0 -5.0 5.0 0.06 1.91
* CR-39, 5 g/cm2
RPP CR3905 -2.0 2.0 -2.0 2.0 1.92 1.97
* Al, 10 g/cm2
RPP Al10 -5.0 5.0 -5.0 5.0 1.98 3.83
* CR-39, 10 g/cm2
RPP CR3910 -2.0 2.0 -2.0 2.0 3.84 3.89
* Al, 15 g/cm2
RPP Al15 -5.0 5.0 -5.0 5.0 3.9 5.75
* CR-39, 15 g/cm2
RPP CR3915 -2.0 2.0 -2.0 2.0 5.76 5.81
* Al, 20 g/cm2
RPP Al20 -5.0 5.0 -5.0 5.0 5.82 7.67
* CR-39, 20 g/cm2
RPP CR3920 -2.0 2.0 -2.0 2.0 7.68 7.73
* Al, 30 g/cm2 (2.6989 g/cm3, +3.71 cm)
RPP Al30 -5.0 5.0 -5.0 5.0 7.74 11.45
* CR-39, 30 g/cm2
RPP CR3930 -2.0 2.0 -2.0 2.0 11.46 11.51
END
BLKHOLE 5 +BLK -VOI
* Front detector
PNTD00 5 +VOI +CR3900
* 5 g/cm2 aluminum
Alum05 5 +VOI +Al05
* 5 g/cm2 detector
PNTD05 5 +VOI +CR3905
* 10 g/cm2 aluminum
Alum10 5 +VOI +Al10
* 10 g/cm2 detector
PNTD10 5 +VOI +CR3910
* 15 g/cm2 aluminum
Alum15 5 +VOI +Al15
* 15 g/cm2 detector
PNTD15 5 +VOI +CR3915
* 20 g/cm2 aluminum
Alum20 5 +VOI +Al20
* 20 g/cm2 detector
PNTD20 5 +VOI +CR3920
* 30 g/cm2 aluminum
Alum30 5 +VOI +Al30
* 30 g/cm2 detector
PNTD30 5 +VOI +CR3930
* System
VAC 5 +VOI -CR3900 -CR3905 -CR3910 -CR3915 -CR3920 -CR3930 -Al05 -Al10
               -Al15 -Al20 -Al30
END
GEOEND
* Aluminum (Al)
STERNHEI 4.2395 0.1708 3.0127 0.08024 3.6345 0.12ALUMINUM
* Al, 5 g/cm2
ASSIGNMA ALUMINUM Alum05
* Al, 10 g/cm2
ASSIGNMA ALUMINUM Alum10
* Al, 15 g/cm2
ASSIGNMA ALUMINUM Alum15
* Al, 20 g/cm2
ASSIGNMA ALUMINUM Alum20
* Al, 30 g/cm2
ASSIGNMA ALUMINUM Alum30
* CR-39
MATERIAL 1.4 CR39
* CR-39
COMPOUND 12.0 CARBON 18.0 HYDROGEN 7.0 OXYGENCR39
* CR-39
* Used -C = Cbar for polycarbonate (Sternheimer)
STERNHEI 3.3201 0.1562 2.6507 0.12679 3.3076 CR39
* CR-39, front
ASSIGNMA CR39 PNTD00
* CR-39, 5 g/cm2
ASSIGNMA CR39 PNTD05
* CR-39, 10 g/cm2
ASSIGNMA CR39 PNTD10
* CR-39, 15 g/cm2
ASSIGNMA CR39 PNTD15
* CR-39, 20 g/cm2
ASSIGNMA CR39 PNTD20
* CR-39, 30 g/cm2
ASSIGNMA CR39 PNTD30
* System
USRBIN 10.0 ENERGY -21. 8.0 8.0 40.0Energy
USRBIN -8.0 -8.0 -5.0 100.0 100.0 100.0&
* CR-39, front, 2.32 cm2 scanned area; (1 LET bin width x 100 charge bin width / 2.32 cm2)^-1 = 0.0232
* Half-integer limits on Min2 and Max2 ensures what is included and what is not
USRYIELD 2223.0 HEAVYION -22. VAC PNTD00 0.0232FLUKA00
USRYIELD 180.0 5.0 175.0 100.5 0.5 2703.&
* CR-39, 5 g/cm2
USRYIELD 2223.0 HEAVYION -23. VAC PNTD05 0.0232FLUKA05
USRYIELD 180.0 5.0 175.0 100.5 0.5 2703.&
* CR-39, 10 g/cm2
USRYIELD 2223.0 HEAVYION -24. VAC PNTD10 0.0232FLUKA10
USRYIELD 230.0 5.0 225.0 100.5 0.5 2703.&
* CR-39, 15 g/cm2
USRYIELD 2223.0 HEAVYION -25. VAC PNTD15 0.0232FLUKA15
USRYIELD 255.0 5.0 0250.0 100.5 0.5 2703.&
* CR-39, 20 g/cm2
USRYIELD 2223.0 HEAVYION -26. VAC PNTD20 0.0232FLUKA20
USRYIELD 280.0 5.0 275.0 100.5 0.5 2703.&
* CR-39, 30 g/cm2
USRYIELD 2223.0 HEAVYION -27. VAC PNTD30 0.0232FLUKA30
USRYIELD 555.0 5.0 550.0 100.5 0.5 2703.&
* CR-39, front, 2.32 cm2 scanned area
USRBIN 10.0 ENERGY 28. 0.76 0.76 0.05CR-3900
USRBIN -0.76 -0.76 0.0 1.0 1.0 1.0&
* CR-39, 5 g/cm2
USRBIN 10.0 ENERGY 29. 0.76 0.76 1.97CR-3905
USRBIN -0.76 -0.76 1.92 1.0 1.0 1.0&
* CR-39, 10 g/cm2
USRBIN 10.0 ENERGY 30. 0.76 0.76 3.89CR-3910
USRBIN -0.76 -0.76 3.84 1.0 1.0 1.0&
* CR-39, 15 g/cm2
USRBIN 10.0 ENERGY 31. 0.76 0.76 5.81CR-3915
USRBIN -0.76 -0.76 5.76 1.0 1.0 1.0&
* CR-39, 20 g/cm2
USRBIN 10.0 ENERGY 32. 0.76 0.76 7.73CR-3920
USRBIN -0.76 -0.76 7.68 1.0 1.0 1.0&
* CR-39, 30 g/cm2
USRBIN 10.0 ENERGY 33. 0.76 0.76 11.51CR-3930
USRBIN -0.76 -0.76 11.46 1.0 1.0 1.0&
RANDOMIZ 1.0
START 82803.0
STOP
Received on Thu Jul 3 09:11:18 2008

This archive was generated by hypermail 2.1.8 : Thu Jul 03 2008 - 09:11:20 CEST