HI
thank you very much for the reply. Let me try it and I will get back to you.
sincerely
saurabh
On Thu, Apr 2, 2020 at 9:15 PM Paola Sala <paola.sala_at_mi.infn.it> wrote:
> 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 - 20:08:26 CEST