On 05.11.2014 15:31, VIRGINIA ALVAREZ SANCHEZ wrote:
>
> Dear fluka users,
>
> I am trying to plot the number of photons reaching a scintillator
> crystal versus time. Some time ago another person asked a similar
> question:
>
> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1964.html
>
> And I have thought that I could try the solution given in this answer:
>
> http://www.fluka.org/web_archive/earchive/new-fluka-discuss/1968.html
>
> But can you please tell me if there is another way in the lastest
> FLUKA version to do it?
>
> Thank you in advance
>
Dear, Virginia Alvarez Sanchez
I think the solution you have mentioned above is the easiest way to
obtain spectrum in time.
All you have to do is to add USERWEIG card, and add code to your
fluscw.f user routine.
Something like:
in your input file add:
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....
USERWEIG 1. 0.0 1.
in your fluscw_mod.f file:
INCLUDE '(SCOHLP)'
INCLUDE '(TRACKR)'
*
FLUSCW = ONEONE
LSCZER = .FALSE.
* spectrum in time from 10 to 25 nseconds
DTMIN = 10.D-09
DTMAX = 25.D-09
IF ( ATRACK .GE. DTMIN .AND. ATRACK .LT. DTMAX ) THEN
FLUSCW = ONEONE
LSCZER = .FALSE.
ELSE
FLUSCW = ZERZER
LSCZER = .TRUE.
END IF
RETURN
*=== End of function Fluscw ===========================================*
END
This will work for all USRTRACK and USRBDX cards.
Best regards,
Mikhail
Received on Thu Nov 06 2014 - 14:59:53 CET