Re: [fluka-discuss]: some questions about deuteron inpfile

From: Francesc Salvat-Pujol <francesc.salvat.pujol_at_cern.ch>
Date: Fri, 9 Nov 2018 11:46:10 +0100

Dear Lucy,

Deuteron interactions at low energies (say below ~150 MeV/n) are an area
of current code development and are not yet included in the production
version.

At higher energies there is nothing deuteron-specific in the models, so
it is rather discouraged to employ them.

With kind regards,

Cesc

On Mon, Apr 30 2018, at 11:41 +0800, zhanghaiqing_at_sinap.ac.cn wrote:
>
>Dear all,
>I want to simulate the nuclear interaction of O-16(D, p) O-17, a Deuteron related discussions in Oct 17 2006 - 13:22:14 CEST (topic: Re: Help with deuteron beam) said that I have to activate nuclear interactions by insert an EVENTYPE card with SDUM=DPMJET and WHAT(3)=2.0, and use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries. But when I insert the EVENTYPE card with SDUM=DPMJET and WHAT(3)=2.0, then the flair program gave a warning: Obsolete card: ENENTYPE. And I don’t know how to use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries. My question is:
>1、 now that the ENENTYPE has been obsolete, when I want to simulate nuclear interactions like O-16(D, p) O-17, do I need to use other instead CARDS to activate the nuclear interaction?
>2、 Do I need to use the executable (named flukadpm3) produced by ldpm3qmd (in flutil), which links the needed libraries? How to do it?
>3、 My input file and the results are below, . I want so simulate a deuteron beam bombarding on different materials containing Oxygen, and then use two Au-Si surface barrier detectors to scoring proton. The two Au-Si surface barrier detectors are placed at 0 degrees and 165 degrees respectively to the beam direction, the distance from detectors to target are 5.5cm. The D beam energy is 6GeV.
> The scoring results are too small and the errors are too large. Is there anything wrong with my input file?
>
>TITLE
>Charged proton fluence in
>* Set the defaults for precision simulations
>DEFAULTS PRECISIO
>* Define the beam characteristics
>* WARNING: Beam Energy/momentum per nucleon
>BEAM -6. 0.0 0.1 0.1 D+
>* Define the beam position
>BEAMPOS 0.0 0.0 -50.
>*DPMJET
>EVENTYPE 2. DPMJET
>*DPMJET
>GEOBEGIN COMBNAME
> 0 0 A simple Be target inside vacuum
>SPH blakhole 0.0 0.0 0.0 200.
>SPH vacumbox 0.0 0.0 0.0 100.
>#if 0
>SPH vacumbo2 0.0 0.0 0.0 45.
>#endif
>#if 0
>* Cylindrical target
>RPP UO2TARG -1. 1. -1. 1. 0.0 0.3
>#endif
>RCC UO2TARG 0.0 0.0 0.0 0.0 0.0 0.3 1.27
>RCC Dec 1.423504748 0.0 5.312592045 2.58819E-5 0.0 9.65926E-5 0.375
>PLA cutplane 1.423504748 0.0 5.312592045 1.423504751 0.0 5.312592054
>RCC Dec1 0.0 0.0 -5.5 0.0 0.0 -1E-4 0.375
>PLA cutplan1 0.0 0.0 -5.5 0.0 0.0 -5.500001
>END
>* Void around
>BLCKHOLE 5 +blakhole -vacumbox
>vacarun2 5 +vacumbox -UO2TARG -Dec -Dec1
>* Target
>UO2TARG 5 +UO2TARG
>* Dec 1st harf
>updec 5 +Dec +cutplane
>* Dec 2st harf
>downdec 5 +Dec -cutplane
>* Dec 1st harf
>updec1 5 +Dec1 +cutplan1
>* Dec 2st harf
>downdec1 5 +Dec1 -cutplan1
>END
>GEOEND
>MATERIAL 8. 16. 0.00143 8. OXYGEN
>MATERIAL 13. 27. 2.6989 10. ALUMINUM
>MATERIAL 79. 196.96655 19.320 15. GOLD
>MATERIAL 14. 28.0885 2.329 14. SILICON
>MATERIAL 3.7 26. AL2O3
>COMPOUND 2. ALUMINUM 3. OXYGEN AL2O3
>LOW-MAT OXYGEN 8. 16. 296. OXYGEN
>LOW-MAT ALUMINUM 13. 27. 296. ALUMINUM
>* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
>ASSIGNMA AL2O3 UO2TARG
>ASSIGNMA BLCKHOLE BLCKHOLE
>ASSIGNMA VACUUM vacarun2
>ASSIGNMA SILICON downdec
>ASSIGNMA GOLD updec
>ASSIGNMA SILICON downdec1
>ASSIGNMA GOLD updec1
>SCORE PROTON
>USRBDX 9. PROTON -47. updec downdec 0.4417886prcurrair
>USRBDX 4E-3 50. 0.0 36. &
>USRYIELD 124. PROTON -56. updec downdec 1.protyield
>USRYIELD 180. 0.0 18. 4E-3 2603. &
>USRBDX 9. PROTON -48. updec1 downdec1 0.4417886prcurrair
>USRBDX 4E-3 50. 0.0 36. &
>USRYIELD 124. PROTON -57. updec1 downdec1 1.protyield
>USRYIELD 180. 0.0 18. 4E-3 2603. &
>RANDOMIZ 1. 86898.
>* Set the number of primary histories to be simulated in the run
>START 4.5E8
>STOP
>
>some results:(1)
> **** protony35_56.yie ****
>
>
> Total primaries run: 2250000000
>
> Total weight of the primaries run: 2.24999987E+09
>
>
> Detector n: 1 protyield
> (User norm: 1.
> sigma: 1. mb
> distr. scored: 1 ,
> from reg. 4 to reg. 5 )
> linear 1st variable (x1) binning from 0.0000E+00 to 3.1416E+00 18 bins ( 1.7453E-01 wide)
> 2nd variable (x2) ranges from 0.0000E+00 to 4.0000E-03
> 1st variable (x1) is: Laboratory Angle (radians)
> 2nd variable (x2) is: Laboratory Kinetic Energy
> The scored double differential yield is (normalized per primary particle):
> plain d2 N / d x1 d x2 where x1, x2 are
> the first and second variables
>
>
> 1st variable (x1) differential is a solid angle one
>
> Tot. response (integrated over x1) 1.5555557E-06 +/- 17.49635 %
> ( --> Tot. response without normaliz. 1.5555557E-06 +/- 17.49635 % )
>
> **** Double differ. yield as a function of x1 ****
> **** (not integrated) ****
>
>
> x1 boundaries:
>
> 3.141593 2.967060 2.792527 2.617994 2.443461
> 2.268928 2.094395 1.919862 1.745329 1.570796
> 1.396263 1.221730 1.047198 0.8726646 0.6981317
> 0.5235988 0.3490658 0.1745329
> Lowest x1 boundary : 0.000000
>
> Yield (double differ.):
>
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 1.9204074E-06 +/- 15.30931 %
> 1.9598615E-06 +/- 44.72136 % 1.1640070E-06 +/- 99.00000 %
> ……
>some results:(2)
> **** protony35_47.bnx ****
>
>
> Charged proton fluence in
>
> Total primaries run: 2250000000
> Total weight of the primaries run: 2.24999987E+09
>
>
> Detector n: 1( 1) prcurrair
> (Area: 0.441788614 cmq,
> distr. scored: 1 ,
> from reg. 4 to 5,
> two ways scoring,
> current scoring)
>
> Tot. resp. (Part/cmq/pr) 1.7102215E-08 +/- 14.99710 %
> ( --> (Part/pr) 7.5555642E-09 +/- 14.99710 % )
>
> **** Different. Fluxes as a function of energy ****
> **** (integrated over solid angle) ****
>
>
> Energy boundaries (GeV):
>
> 4.0000030E-03 3.8906229E-03 3.7842335E-03 3.6807533E-03 3.5801029E-03
> 3.4822049E-03 3.3869836E-03 3.2943664E-03 3.2042817E-03 3.1166603E-03
> 3.0314352E-03 2.9485403E-03 2.8679122E-03 2.7894890E-03 2.7132102E-03
> 2.6390173E-03 2.5668533E-03 2.4966625E-03 2.4283910E-03 2.3619866E-03
> 2.2973979E-03 2.2345753E-03 2.1734708E-03 2.1140371E-03 2.0562285E-03
> 2.0000008E-03 1.9453107E-03 1.8921160E-03 1.8403759E-03 1.7900508E-03
> 1.7411017E-03 1.6934912E-03 1.6471825E-03 1.6021403E-03 1.5583296E-03
> 1.5157170E-03 1.4742695E-03 1.4339556E-03 1.3947440E-03 1.3566046E-03
> 1.3195082E-03 1.2834262E-03 1.2483308E-03 1.2141950E-03 1.1809928E-03
> 1.1486985E-03 1.1172872E-03 1.0867349E-03 1.0570182E-03 1.0281139E-03
> Lowest boundary (GeV): 1.0000000E-03
>
> Flux (Part/GeV/cmq/pr):
>
> 9.1973816E-06 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 9.7217981E-06 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 1.1167414E-05 +/- 99.00000 %
> 0.000000 +/- 0.000000 % 1.1804158E-05 +/- 99.00000 %
> 1.2136017E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 2.5655980E-05 +/- 61.23724 % 1.3188634E-05 +/- 99.00000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 1.4332549E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 1.6926635E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 1.9443596E-05 +/- 99.00000 %
> 1.9990230E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 2.2334827E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 2.4272034E-05 +/- 99.00000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 2.9470983E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 3.1151358E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> 0.000000 +/- 0.000000 % 0.000000 +/- 0.000000 %
> 3.4805016E-05 +/- 99.00000 % 0.000000 +/- 0.000000 %
> ……
>
>best regards
>
>Lucy Zhang
>











--
Francesc Salvat Pujol
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
Tel: +41 22 76 64011
Fax: +41 22 76 69474
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
Received on Fri Nov 09 2018 - 13:22:33 CET

This archive was generated by hypermail 2.3.0 : Fri Nov 09 2018 - 13:22:37 CET