That works like charm! thanks.
saurabh
On Thu, Apr 2, 2020 at 10:18 PM Saurabh Mukherjee <sxm3816_at_gmail.com> wrote:
> 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 16 2020 - 14:02:23 CEST