Dear Alberto Fasso' :
I name the modified routine MGDRAW as mgdrawF18.f and put it in the
directory of /fluka/flutil.
( is it correct or should I put it in directory usermax ?)
I complied it and got error as below.
lee_at_lee-X402CA:~/physics/fluka/flutil$ ./fff mgdrawF18.f
awk: line 0: regular expression compile failed (missing '('))
./fff: 92: [: -le: unexpected operator
./fff: 109: ./fff: [[: not found
Please advise me
Thanks in advance
Lee
2014-01-22 Alberto Fasso' <fasso_at_slac.stanford.edu>
> Dear Lee,
>
> you get a big output file because you are using the default version of
> MGDRAW. It is better to follow the instructions below (that you find also
> in
> the manual description of USERDUMP):
>
> The default version of the routine writes a complete dump
> (unformatted) of
> one or more of the following:
> - all source particles
> - all particle trajectories
> - all local (pointlike) energy deposition events
> - all continuous energy deposition events
> - user-defined events
> Users can modify the routine by removing the existing lines of code
> and by
> writing their own code under one or more of the entries.
>
> That's it: remove all the lines of code (except the INCLUDE statements at
> the beginning) and open your own files and write your own logics under the
> ENTRYes you are interested in.
>
> In a similar way,
>
> WHAT(2) : if the default version of MGDRAW is used, ....
> ... forget this
>
> if a user version is used: ...
> ... follow these instructions for the user version
>
> Alberto
>
>
>
> On Sat, 18 Jan 2014, Chenyen Lee wrote:
>
> Dear Anton and Forum :
>>
>> Pardon me, I would like to ask a very elementary question.
>> I do not know how to use the USDRAW in the MGDRAW to modify to as Anton's
>> instuction.
>> I add the USERDUMP, but just get a big output file of DUMP.
>> I find the mgdraw.f file but do not know how to do it correctly.
>>
>> I read the manual and the UserRoutines.pdf
>> www.fluka.org/content/course/NEA/.../UserRoutines.pdf
>> Is there any more elementary documents or instructions?
>>
>> Thanks in advance
>>
>> Lee
>>
>>
>>
>>
>> 2013/11/30 Anton Lechner <Anton.Lechner_at_cern.ch>
>> Dear Lee,
>>
>> You can interpret the values you specify (from the summary file
>> created with USYSUW) as follows: for example, the differential
>> yield of 0.1 GeV^-1 you get in the energy interval 0-0.01 GeV
>> simply means that 0.1 GeV^-1*0.01GeV=1E-03 neutrons with E<10
>> MeV were produced in inelastic hadronic interactions within
>> your entire geometry (also in the surrounding air). Of course
>> this is consistent with the value you got for the larger energy
>> interval, since most neutrons are below 10 MeV.
>>
>> From the link you included in your file, I guess that you are
>> only interested in neutrons produced in following reaction:
>> p + O18 -> F18 + n
>>
>> If this is the case, you should not forget that a
>> non-negligible fraction of neutrons also emerge from other
>> processes, e.g.:
>> p + O18 -> O17 + n + p
>> p + O18 -> N14 + n + alpha
>>
>> USRYIELD will count all neutrons, also those from the latter
>> interactions.
>>
>> If you are only interested in the first reaction, you could for
>> example use USDRAW in the MGDRAW user routine. By including
>> '(RESNUC)', you have access to the atomic and mass number of
>> the residual nucleus and you can filter neutrons accordingly. I
>> included some sample code below.
>>
>> Cheers, Anton
>>
>>
>> * Check that incoming particle is a proton and that interaction
>> is inelastic:
>> IF ( JTRACK .EQ. 1 .AND. ICODE .EQ. 101 ) THEN
>>
>> * Check that residual nucleus is F18:
>> IF( ICRES .EQ.9 .AND. IBRES .EQ.18 ) THEN
>>
>> * Loop over all secondaries
>> DO 99 INMB = 1, NP
>>
>> * Check that secondary particle is a neutron
>> IF ( KPART( INMB ) .EQ. 8 ) THEN
>> * Here you might want to write to a file e.g. the neutron
>> energy TKI (INMB ) or simply count neutrons
>> ENDIF
>>
>> 99 CONTINUE
>>
>> ENDIF
>> ENDIF
>>
>>
>>
>> Thanks in advance
>>
>>
>> Lee
>>
>>
>> --
>> 03-3281200-3591
>> 0972164977
>> Chenyen Lee
>> Department of Nuclear Medicine
>> Chang Gang Memorial Hospital
>> Linko,Taiwan
>>
>>
>>
>>
>>
> --
> Alberto Fasso`
> SLAC-RP, MS 48, 2575 Sand Hill Road, Menlo Park CA 94025
> Phone: (1 650) 926 4762 Fax: (1 650) 926 3569
> fasso_at_slac.stanford.edu
Received on Sun Jan 26 2014 - 07:28:51 CET