RE: [fluka-discuss]: numerically counting electrons

From: Anton Lechner <Anton.Lechner_at_cern.ch>
Date: Tue, 7 Feb 2017 11:00:28 +0000

Dear Elif,

Another way to count electrons is with the MGDRAW routine. In this way you can identify the process in which electrons were generated (e.g. pair production, compton), i.e. you can distinguish electrons produced by photons from delta-electrons produced by electrons/positrons. MGDRAW has an entry (USDRAW) which is called whenever you have an interaction (you need to activate it by means of the USERDUMP card). The type of interaction can be distinguished with ICODE (which is e.g. 217 for pair production, 219 for Compton and 221 for Photo-electric effect). In addition, you can for example select only interactions of primary photons if you wish so (using the generation number LTRACK, which is 1 for primaries).

Here is a code snippet (to be inserted into USRDRAW) which dumps all electrons produced by primary photons to a file (fortran unit 50).
* Interaction by a primary photon (assuming you have
* defined a photon beam)
      IF ( LTRACK .EQ. 1) THEN
*
* Loop over all secondaries and dump kinetic
* energy, statistical weight and direction cosines if
* the particle is an electron
         DO 10 INP = 1, NP
            IF( KPART(INP) .EQ. 3 ) THEN
            WRITE( 50,* )
     & ' E=',TKI(INP),' Weight=',WEI(INP),
     & ' Tx=',CXR(INP), ' Ty=',CYR(INP),' Tz=',CZR(INP)
             END IF
  10 CONTINUE
*
      END IF

Just to reiterate on a very important remark already made by Andrea: the number of particles produced of course depends on the production thresholds which you select in your simulation (EMFCUT card).

Cheers, Anton



________________________________________
From: owner-fluka-discuss_at_mi.infn.it [owner-fluka-discuss_at_mi.infn.it] on behalf of Andrea Fontana [andrea.fontana_at_pv.infn.it]
Sent: 07 February 2017 10:56
To: 91120020413_at_ogrenci.ege.edu.tr; fluka-discuss_at_fluka.org
Subject: Re: [fluka-discuss]: numerically counting electrons

Dear Elif,
     many thanks for you inquiry to the FLUKA forum.

In FLUKA there are different methods to count secondary particles and
perhaps the simplest way in your case is to use a USRBDX score
with a linear one-way current from the region of your target to the
region surrounding it:

USRBDX 1. ELECTRON -21. TARGET VOID ecount1
USRBDX 1E-3 1E-6 1. &

You can integrate on the energy bins to get the total number,
as in this example, or you can specify more bins to see the
electrons energy distribution.

Please bear in mind that the result depends on the thresholds that
you have set in your problem (electron transport cutoff to be set
with EMFCUT) and also that the result is always normalized per primary
particle, in your case per incident photon.

Also looking in the output file, as suggested by Joachim, is always
very useful. However, being yours a purely electro-magnetic problem, the
table he indicated is not produced since it refers to the products of
nuclear interactions only.

Kind regards,
Andrea


Il 07/02/2017 07:41, 91120020413_at_ogrenci.ege.edu.tr ha scritto:
> Dear FLUKA users,
>
> How can I numerically detect and count electrons which are generated by
> photons?
> For example I sent 100 photons, how many electrons were generated?
> Thank you for your help.
>
> Elif Ermis
>
> __________________________________________________________________________
> You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
>

--
========================================================================
Dr. Andrea Fontana                    tel: +39 0382 987991
Istituto Nazionale                    fax: +39 0382 423241
di Fisica Nucleare
Sezione di Pavia                      e-mail: andrea.fontana_at_pv.infn.it
Via Bassi 6                           web   : www.pv.infn.it/~fontana
27100 PAVIA, Italy
========================================================================
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info
__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id¬c_info
Received on Tue Feb 07 2017 - 13:45:13 CET

This archive was generated by hypermail 2.3.0 : Tue Feb 07 2017 - 13:45:15 CET