[fluka-discuss]: Mgdraw.f to study processes

From: Julien Bonvalet <julien.bonvalet_at_u-bordeaux.fr>
Date: Thu, 3 Sep 2020 14:58:10 +0200 (CEST)

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

pB11.png
(image/png attachment: pB11.png)

Received on Thu Sep 03 2020 - 17:11:23 CEST

This archive was generated by hypermail 2.3.0 : Thu Sep 03 2020 - 17:11:25 CEST