Dear Shirley,
Electronuclear interactions are not accessible in the current production
version (2011.2x.5).
SDUM=ELECTNUC in the PHOTONUC card will have no effect until the next
major release.
With kind regards,
Cesc
On Thu, Feb 07 2019, at 14:50 -0800, Shirley Li wrote:
>
>Dear FLUKA experts,
>
>I'm trying to simulate electron scattering in FLUKA, i.e., e+C -> e +
>nucleons + mesons. I inject 2 GeV electrons on a thin carbon target. I
>turn on electronuclear interaction with PHOTONUC card, and bias the
>inelastic interaction by 1e-6 (also tried 5e-4, didn't work either). I
>see 0 electron inelastic interaction in the output, neither in USRYIELD
>card nor in mgdraw (JTRACK=3 and ICODE=101). Are there any problems in my
>input to turn on this interaction? Inputs and results attached.
>
>Thank you very much!
>
>--
>Shirley Li
>Postdoctoral Research Associate
>Theory Group, MS 81
>SLAC National Accelerator Laboratory
>2575 Sand Hill Rd, Menlo Park, CA 94025
>*$ CREATE MGDRAW.FOR
>*COPY MGDRAW
>* *
>*=== mgdraw ===========================================================*
>* *
> SUBROUTINE MGDRAW ( ICODE, MREG )
>
> INCLUDE '(DBLPRC)'
> INCLUDE '(DIMPAR)'
> INCLUDE '(IOUNIT)'
>*
>*----------------------------------------------------------------------*
>* *
>* Copyright (C) 1990-2006 by Alfredo Ferrari *
>* All Rights Reserved. *
>* *
>* *
>* MaGnetic field trajectory DRAWing: actually this entry manages *
>* all trajectory dumping for *
>* drawing *
>* *
>* Created on 01 march 1990 by Alfredo Ferrari *
>* INFN - Milan *
>* Last change 05-may-06 by Alfredo Ferrari *
>* INFN - Milan *
>* *
>*----------------------------------------------------------------------*
>*
> INCLUDE '(CASLIM)'
> INCLUDE '(COMPUT)'
> INCLUDE '(SOURCM)'
> INCLUDE '(FHEAVY)'
> INCLUDE '(FLKSTK)'
> INCLUDE '(GENSTK)'
> INCLUDE '(MGDDCM)'
> INCLUDE '(PAPROP)'
> INCLUDE '(QUEMGD)'
> INCLUDE '(SUMCOU)'
> INCLUDE '(TRACKR)'
> INCLUDE '(EMFSTK)'
> INCLUDE '(RESNUC)'
>
>*
> DIMENSION DTQUEN ( MXTRCK, MAXQMG )
>*
> CHARACTER*20 FILNAM
> LOGICAL LFCOPE
> SAVE LFCOPE
> DATA LFCOPE / .FALSE. /
>
>
>*
>*----------------------------------------------------------------------*
>* *
>* Icode = 1: call from Kaskad *
>* Icode = 2: call from Emfsco *
>* Icode = 3: call from Kasneu *
>* Icode = 4: call from Kashea *
>* Icode = 5: call from Kasoph *
>* *
>*----------------------------------------------------------------------*
>* *
>c IF ( .NOT. LFCOPE ) THEN
>c LFCOPE = .TRUE.
>c IF ( KOMPUT .EQ. 2 ) THEN
>c FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
>c ELSE
>c FILNAM = CFDRAW
>c END IF
>c OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
>c & 'UNFORMATTED' )
>c END IF
>c
>c WRITE (IODRAW,*) NTRACK, MTRACK, JTRACK, SNGL (ETRACK),
>c & SNGL (WTRACK)
>c WRITE (IODRAW,*) ( SNGL (XTRACK (I)), SNGL (YTRACK (I)),
>c & SNGL (ZTRACK (I)), I = 0, NTRACK ),
>c & ( SNGL (DTRACK (I)), I = 1, MTRACK ),
>c & SNGL (CTRACK)
>
>* +-------------------------------------------------------------------*
>* | Quenching is activated
> IF ( LQEMGD ) THEN
> IF ( MTRACK .GT. 0 ) THEN
> RULLL = ZERZER
> CALL QUENMG ( ICODE, MREG, RULLL, DTQUEN )
> WRITE (IODRAW,*) ( ( SNGL (DTQUEN (I,JBK)), I = 1, MTRACK ),
> & JBK = 1, NQEMGD )
> END IF
> END IF
>* | End of quenching
>* +-------------------------------------------------------------------*
> RETURN
>*
>*======================================================================*
>* *
>* Boundary-(X)crossing DRAWing: *
>* *
>* Icode = 1x: call from Kaskad *
>* 19: boundary crossing *
>* Icode = 2x: call from Emfsco *
>* 29: boundary crossing *
>* Icode = 3x: call from Kasneu *
>* 39: boundary crossing *
>* Icode = 4x: call from Kashea *
>* 49: boundary crossing *
>* Icode = 5x: call from Kasoph *
>* 59: boundary crossing *
>* *
>*======================================================================*
>* *
> ENTRY BXDRAW ( ICODE, MREG, NEWREG, XSCO, YSCO, ZSCO )
> WRITE (41, 800) NCASE, JTRACK, ETRACK, MREG, NEWREG,
> & ZSCO, XSCO, YSCO
> 800 FORMAT (I7, I4, F15.5, I3, I3, F15.3, F15.3, F15.3)
> RETURN
>*
>*======================================================================*
>* *
>* Event End DRAWing: *
>* *
>*======================================================================*
>* *
> ENTRY EEDRAW ( ICODE )
> RETURN
>*
>*======================================================================*
>* *
>* ENergy deposition DRAWing: *
>* *
>* Icode = 1x: call from Kaskad *
>* 10: elastic interaction recoil *
>* 11: inelastic interaction recoil *
>* 12: stopping particle *
>* 13: pseudo-neutron deposition *
>* 14: escape *
>* 15: time kill *
>* Icode = 2x: call from Emfsco *
>* 20: local energy deposition (i.e. photoelectric) *
>* 21: below threshold, iarg=1 *
>* 22: below threshold, iarg=2 *
>* 23: escape *
>* 24: time kill *
>* Icode = 3x: call from Kasneu *
>* 30: target recoil *
>* 31: below threshold *
>* 32: escape *
>* 33: time kill *
>* Icode = 4x: call from Kashea *
>* 40: escape *
>* 41: time kill *
>* 42: delta ray stack overflow *
>* Icode = 5x: call from Kasoph *
>* 50: optical photon absorption *
>* 51: escape *
>* 52: time kill *
>* *
>*======================================================================*
>* *
> ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
> IF ( .NOT. LFCOPE ) THEN
> LFCOPE = .TRUE.
> IF ( KOMPUT .EQ. 2 ) THEN
> FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
> ELSE
> FILNAM = CFDRAW
> END IF
> OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
> & 'FORMATTED' )
> END IF
>
>
>* +-------------------------------------------------------------------*
>* | Quenching is activated : calculate quenching factor
>* | and store quenched energy in DTQUEN(1, jbk)
> IF ( LQEMGD ) THEN
> RULLL = RULL
> CALL QUENMG ( ICODE, MREG, RULLL, DTQUEN )
> WRITE (IODRAW,*) ( SNGL (DTQUEN(1, JBK)), JBK = 1, NQEMGD )
> END IF
>* | end quenching
>* +-------------------------------------------------------------------*
> RETURN
>*
>*======================================================================*
>* *
>* SOurce particle DRAWing: *
>* *
>*======================================================================*
>*
> ENTRY SODRAW
> IF ( .NOT. LFCOPE ) THEN
> LFCOPE = .TRUE.
> IF ( KOMPUT .EQ. 2 ) THEN
> FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
> ELSE
> FILNAM = CFDRAW
> END IF
> OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
> & 'FORMATTED' )
> END IF
> WRITE (60, *) -NCASE, NPFLKA, ILOFLK(NPFLKA), TKEFLK(NPFLKA),
> & TXFLK(NPFLKA), TYFLK(NPFLKA),
> & NSTMAX, SNGL (TKESUM), SNGL (WEIPRI)
>* +-------------------------------------------------------------------*
>* | (Radioactive) isotope: it works only for 1 source particle on
>* | the stack for the time being
> IF ( ILOFLK (NPFLKA) .GE. 100000 .AND. LRADDC (NPFLKA) ) THEN
> IARES = MOD ( ILOFLK (NPFLKA), 100000 ) / 100
> IZRES = MOD ( ILOFLK (NPFLKA), 10000000 ) / 100000
> IISRES = ILOFLK (NPFLKA) / 10000000
> IONID = ILOFLK (NPFLKA)
> WRITE (IODRAW,*) ( IONID,SNGL(-TKEFLK(I)),
> & SNGL (WTFLK(I)), SNGL (XFLK (I)),
> & SNGL (YFLK (I)), SNGL (ZFLK (I)),
> & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
> & SNGL (TZFLK(I)), I = 1, NPFLKA )
>* |
>* +-------------------------------------------------------------------*
>* | Patch for heavy ions: it works only for 1 source particle on
>* | the stack for the time being
> ELSE IF ( ABS (ILOFLK (NPFLKA)) .GE. 10000 ) THEN
> IONID = ILOFLK (NPFLKA)
> CALL DCDION ( IONID )
> WRITE (IODRAW,*) ( IONID,SNGL(TKEFLK(I)+AMNHEA(-IONID)),
> & SNGL (WTFLK(I)), SNGL (XFLK (I)),
> & SNGL (YFLK (I)), SNGL (ZFLK (I)),
> & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
> & SNGL (TZFLK(I)), I = 1, NPFLKA )
>* |
>* +-------------------------------------------------------------------*
>* | Patch for heavy ions: ???
> ELSE IF ( ILOFLK (NPFLKA) .LT. -6 ) THEN
> WRITE(IODRAW,*)(IONID,SNGL(TKEFLK(I)+AMNHEA(-ILOFLK(NPFLKA))),
> & SNGL (WTFLK(I)), SNGL (XFLK (I)),
> & SNGL (YFLK (I)), SNGL (ZFLK (I)),
> & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
> & SNGL (TZFLK(I)), I = 1, NPFLKA )
>* |
>* +-------------------------------------------------------------------*
>* |
> ELSE
>c WRITE (IODRAW,*) ( ILOFLK(I), SNGL (TKEFLK(I)+AM(ILOFLK(I))),
>c & SNGL (WTFLK(I)), SNGL (XFLK (I)),
>c & SNGL (YFLK (I)), SNGL (ZFLK (I)),
>c & SNGL (TXFLK(I)), SNGL (TYFLK(I)),
>c & SNGL (TZFLK(I)), I = 1, NPFLKA )
> END IF
>* |
>* +-------------------------------------------------------------------*
> RETURN
>*
>*======================================================================*
>* *
>* USer dependent DRAWing: *
>* *
>* Icode = 10x: call from Kaskad *
>* 100: elastic interaction secondaries *
>* 101: inelastic interaction secondaries *
>* 102: particle decay secondaries *
>* 103: delta ray generation secondaries *
>* 104: pair production secondaries *
>* 105: bremsstrahlung secondaries *
>* 110: decay products *
>* Icode = 20x: call from Emfsco *
>* 208: bremsstrahlung secondaries *
>* 210: Moller secondaries *
>* 212: Bhabha secondaries *
>* 214: in-flight annihilation secondaries *
>* 215: annihilation at rest secondaries *
>* 217: pair production secondaries *
>* 219: Compton scattering secondaries *
>* 221: photoelectric secondaries *
>* 225: Rayleigh scattering secondaries *
>* Icode = 30x: call from Kasneu *
>* 300: interaction secondaries *
>* Icode = 40x: call from Kashea *
>* 400: delta ray generation secondaries *
>* For all interactions secondaries are put on GENSTK common (kp=1,np) *
>* but for KASHEA delta ray generation where only the secondary elec- *
>* tron is present and stacked on FLKSTK common for kp=npflka *
>* *
>*======================================================================*
>*
> ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
> IF ( .NOT. LFCOPE ) THEN
> LFCOPE = .TRUE.
> IF ( KOMPUT .EQ. 2 ) THEN
> FILNAM = '/'//CFDRAW(1:8)//' DUMP A'
> ELSE
> FILNAM = CFDRAW
> END IF
> OPEN ( UNIT = IODRAW, FILE = FILNAM, STATUS = 'NEW', FORM =
> & 'FORMATTED' )
> END IF
>* No output by default:
>*----------------------------------------------------------------------*
>
>
> IF (MREG .EQ. 3) THEN
> IF (ICODE .EQ. 101) THEN
> WRITE(70, *) NCASE, JTRACK, ICODE, ETRACK,
> & CXTRCK, CYTRCK, CZTRCK, WTRACK
> do 10 ip = 1, NP
> WRITE(70, *) ip, KPART(ip), TKI(ip), CXR(ip), CYR(ip),
> & CZR(ip)
> 10 continue
> END IF
>
> END IF
>
>*----------------------------------------------------------------------*
> RETURN
>*==== end of subrutine Mgdraw =========================================*
> END
>
--
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 Feb 08 2019 - 16:45:18 CET