From: Lindley Winslow (lawinslow@lbl.gov)
Date: Wed Nov 15 2006 - 20:48:51 CET
Hello all,
I am running 300 GeV muons through liquid scintillator. I wanted to
output a text file of secondaries and some information about the
parent particle (i.d. and energy) and the interaction that created
them. I am particularly interested residual nuclei. I modified
USDRAW as follows. I had two questions:
1) How do you determine if a reaction resulted in a RESNUC?
2) Is there a more correct way to access this information?
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 =
& 'UNFORMATTED' )
END IF
* No output by default:
* +-------------------------------------------------------------------*
* Output to text file.
WRITE( 71.0, * ) 'Routine USDRAW called'
WRITE( 71.0, * ) NCASE, NPFLKA, JTRACK, ICODE, NP
* Loop over Secondaries.
do 10 ip = 1, NP
WRITE( 71.0, *) ip, KPART(ip), TKI(ip)
10 continue
* Now the heavy secondaries and possibly resnuc
WRITE( 71.0, * ) NPHEAV, ICRES, IBRES
do 20 ip = 1, NPHEAV
WRITE( 71.0, *) ip, KHEAVY(ip), TKHEAV(ip)
20 continue
* ENDIF
RETURN
*=== End of subrutine Mgdraw ==========================================*
END
Thank you!
Lindley
This archive was generated by hypermail 2.1.6 : Thu Nov 16 2006 - 09:26:22 CET