Re: [fluka-discuss]: Neutrino primaries with source.f user routine

From: Paola Sala <paola.sala_at_mi.infn.it>
Date: Thu, 16 Mar 2017 17:03:34 +0100

Dear Andrea

the use of an user-written source for
neutrino interactions is possible, provided the following card
is added when a neutrino is loaded on the stack:
            LFRPHN (NPFLKA) = .TRUE.
 In this way, the neutrino will be forced to interact at the
     point specified as starting position in the source routine

 this information is present in the release notes..not propagated to the
user manual, sorry for the inconvenience).


As a consequence of this "forced interaction" behavior, if you need to
generate events with a given energy distribution, you have to sample from
a distribution in rate, not in flux.

Hope this helps
Paola
> Dear Fluka experts,
> I want to write a user routine to generate neutrino primaries within a
> volume, with given energy distribution, in order to study the neutrino
> interactions in that volume.
> I understood that fluka, whenever the primary particle ID is set to a
> neutrino one with the BEAM card, is NOT tracking the neutrino, rather it
> is forcing a neutrino interaction in the point where the primary is
> generated.
>
> If I just set the primary id to 27 (nu mu) in BEAM, and I use BEAMPOS to
> fix the neutrino position within my volume, I see that interactions are
> present.
>
> I see this by activating a data collision tape, with the USERDUMP card.
> My USERDRAW code is as follows:
>
> ENTRY USDRAW ( ICODE, MREG, XSCO, YSCO, ZSCO )
> if (ICODE.EQ.101) then
> if ((JTRACK .eq. 27).or.(JTRACK . eq. 28).or.(JTRACK.eq.5)
> & .or.(JTRACK.eq.6)) THEN
> print *,'NEUTRINO',JTRACK,ETRACK
> do IP=1,NP
> print *,'secondary',KPART(ip),Plr(ip
> end do
> end if
> end if
>
> If, however, I take the source.f file from $FLUPRO/usermvax (without
> modifying it), and I compile it and link to a fluka executable, then
> neutrino interactions are no longer present (i do not have any print in
> the output log) (I use the same input file in both cases, with BEAM and
> BEAMPOS cards)
>
> Is this normal? Is there any incompatibility between source.f user
> routine and primary neutrinos?
>
> Thanks
>
> Andrea
>
>
>


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 Mar 16 2017 - 18:16:10 CET

This archive was generated by hypermail 2.3.0 : Thu Mar 16 2017 - 18:16:13 CET