Hello
it works,
but there is a more direct way to do, without the need to set ad-hoc
thresolds: you can intercept pair production un the USDRAW entry of
mgdraw.f, selecting icode=217, and look at secondaries that are scored in
the GENSTCK common :
ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
IF ( .NOT. LFCOPE ) THEN
LFCOPE = .TRUE.
OPEN ( UNIT = 89, FILE = "collfile2", STATUS = 'UNKNOWN')
OPEN ( UNIT = 88, FILE = "collfile", STATUS = 'UNKNOWN')
END IF
CALL GEOR2N ( MREG, MRGNAM, IERR1 )
IF( ICODE .EQ. 217 .AND. MRGNAM .EQ. "WMesh1") THEN
DO KP = 1, NP
IF( KPART (KP) .EQ. 4 ) WRITE(89,'(1P,4G25.15)')
& XSCO, TKI (KP)
ENDDO
END IF
Paola
> Dear Users
> Let me explain the simulation I am trying to setup. A beam of thermal
> neutron is incident on cadmium. The resulting capture gammas are incident
> on a tungsten mesh. These gamma photons will cause pair production in the
> tungsten mesh.
> I am trying to record the kinetic energy of positrons right at the instant
> that they are born. To do this, I have set the transport threshold for
> positrons in the EMFCUT really high at 1 GeV. The production cut is set at
> 511 keV. Then in mgdraw file, I write the energy of the positron before it
> is killed (due to high transport threshold) using ENDDRAW. The code is-
> ENTRY ENDRAW ( ICODE, MREG, RULL, XSCO, YSCO, ZSCO )
> IF ( .NOT. LFCOPE ) THEN
> LFCOPE = .TRUE.
> OPEN ( UNIT = 88, FILE = "collfile", STATUS = 'UNKNOWN')
> END IF
> CALL GEOR2N ( MREG, MRGNAM, IERR1 )
> IF(MRGNAM .EQ. "WMesh1") THEN
> IF(JTRACK .EQ. 4) WRITE(88,'(1P,4G25.15)')
> & XSCO, RULL
> END IF
> RETURN
>
> The question is- Is this a correct way of recording the positron's KE
> right
> after its birth or is their another way?
>
> I am attaching the input file and mgdraw file.
> thanks for your attention
> sincerely
> saurabh mukherjee
>
Paola Sala
INFN Milano
tel. Milano +39-0250317374
tel. CERN +41-227679148
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Thu Apr 02 2020 - 19:20:04 CEST