Dear fluka experts,
I am studying the proton Boron 11 fusion and I would like to discriminate ways to produce alphas particles between:
1) p + 11B --> 12C* --> 8Be* + alpha0 --> alpha1 + alpha2 +alpha0
2) p + 11B --> 12C* --> 8Be + alpha00 --> alpha01 + alpha02 +alpha00
3) p + 11B --> 12C* --> alpha11 + alpha22 +alpha20
To achieve that, i use mgdraw.f routine:
ENTRY SODRAW
RETURN
*
ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
IF ( UDFIRST ) THEN
UDFIRST = .FALSE.
ENDIF
do 20 ip = 1, NP
IF (KPART(ip) .EQ. -6 .OR. KPART(ip) .EQ. 7 .OR.
& KPART(ip) .EQ. 8) THEN
WRITE(IODRAW,*)'2NDY ',NCASE,JTRACK,ZSCO,
& KPART(ip),ZPART,APART,TKI(ip)
ENDIF
20 continue
do 30 ip=1,NPHEAV
WRITE (IODRAW,*) 'Heav', NCASE, JTRACK,ZSCO,
& ICODE,KHEAVY(ip),ICHEAV(KHEAVY(ip)),IBHEAV(KHEAVY(ip)),
& TKHEAV(ip)
30 continue
IF (ICODE.NE.100) THEN
WRITE(IODRAW,*) 'nuc', NCASE,JTRACK,ICODE,ICHTAR,
& IBTAR, ICESTR, IBESTR,ICRES,IBRES
ENDIF
RETURN
I obtain these kind of informations:
2NDY 1882 1 0.13088856811273053 -6 0.0000000000000000 0.0000000000000000 2.2024593427923289E-002
Heav 1882 1 0.13088856811273053 101 6 2 4 5.8379209697727672E-003
Heav 1882 1 0.13088856811273053 101 6 2 4 7.9018346964389828E-003
nuc 1882 1 101 5 11 4 8 0 0
Here for the Ncase=1882, it seems we have 1 alpha as secondary and 2 alphas as heavy secondaries from 8Be evaporation (produced at the same place in the simulation). OK. But how can i make difference between reactions 1) and 2) : Be8 is in an excited state before evaporation or not?
For reaction 3):
Heav 5523 1 8.4310820700689224E-002 101 6 2 4 6.5847353447588833E-003
Heav 5523 1 8.4310820700689224E-002 101 6 2 4 6.6214411305522847E-003
Heav 5523 1 8.4310820700689224E-002 101 6 2 4 9.2596922984330380E-004
nuc 5523 1 101 5 11 6 12 0 0
Here it seems that correspond to equation 3. 12C directly break up into 3 alphas (3 heavy secondaries).
Whhen I count the cases for each type of reaction, I obtain the proportions in the attached picture pB11.png
It seems that reaction 3 contributes in reality for less than 5% with protons from 0 to 5MeV (it is this kind of spectrum I use in FLUKA). So this reaction with direct break up of C12 in 3 alphas particles seems to be overestimated in FLUKA. Is it possible or did I make a mistake in my routine mgdraw?
Thanks in advance
Julien
CELIA FRANCE
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Sep 03 2020 - 17:11:23 CEST