EMFCUT: production thresholds overrides transport thresholds

From: George Kharashvili <georgek_at_jlab.org>
Date: Wed, 20 Apr 2011 16:46:18 -0400

Dear FLUKA experts,

It appears that the production thresholds set by EMFCUT with PROD-CUT
sdum override the transport thresholds set by EMFCUT with blank sdum.
For example, if production cuts are set to 3 MeV for e-, e+ and photons,
and transport cuts are set to 2 MeV, the output file shows that Ae, Ap,
Ecut, and Pcut were ALL taken to be 3 MeV (as is happens in case of the
attached input).
However, the manual (EMFCUT) and slide 7 of the lecture on Ionization
and Transport of the last FLUKA course suggest that when the production
threshold is higher than the transport threshold, particles with
energies between the two thresholds are still transported and loose
energy by dE/dx.
Could this be a bug in the code?
I noticed that FLUKA-2008 was doing the same.

Best regards,
  - George

-- 
George Kharashvili                     | Tel: +1 (757) 269-6435
Radiation Control Dept., Jefferson Lab | Fax: +1 (757) 269-6050
12050 Jefferson Avenue, MS 6B          |
Newport News, Virginia 23606, USA      | E-mail: georgek_at_jlab.org
--------------040804000501030200090709
Content-Type: text/plain;
 name="emfcuttest.txt"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
 filename="emfcuttest.txt"
TITLE
GLOBAL                                         1.0       1.0
DEFAULTS                                                             
NEW-DEFA
BEAM            -5.0                                                 
ELECTRON
BEAMPOS
GEOBEGIN                                                             
COMBNAME
    0    0
* Black body
SPH blkbody    0.0 0.0 0.0 100000.0
* Void sphere
SPH void       0.0 0.0 0.0 10000.0
* Cylindrical target
RCC target     0.0 0.0 0.1 0.0 0.0 10.0 5.0
END
* Black hole
BLKBODY      5 +blkbody -void
* Void around
VOID         5 +void -target
* Target
TARGET       5 +target
END
GEOEND
* ..+....1....+....2....+....3....+....4....+....5....+....6....+....7..
ASSIGNMA    BLCKHOLE   BLKBODY
ASSIGNMA      VACUUM      VOID
ASSIGNMA      COPPER    TARGET
EMFCUT        -0.003     0.003            BLCKHOLE  @LASTMAT         
PROD-CUT
EMFCUT        -0.002     0.002             BLKBODY  @LASTREG
RANDOMIZ         1.0
START          100.0
STOP
--------------040804000501030200090709
Content-Type: text/plain;
 name="emfcuttest_out.txt"
Content-Transfer-Encoding: 7bit
Content-Disposition: attachment;
 filename="emfcuttest_out.txt"
1   
****************************************************************************************************
   FLUKA2011 Version 2.2  Apr-11 by A. Ferrari DATE:  4/19/11    TIME:
16:52: 8
     ****************************************************************************************************
 ***** Next control card *****   TITLE      0.000       0.000       0.000 
     0.000       0.000       0.000
 ***** Next control card *****   GLOBAL     0.000       0.000       0.000 
     1.000       1.000       0.000
 ***** Next control card *****   DEFAULTS   0.000       0.000       0.000 
     0.000       0.000       0.000      NEW-DEFA
 *----------------   
..+....1....+....2....+....3....+....4....+....5....+....6....+....7..   
       ----------------*
 ***** Next control card *****   BEAM      -5.000       0.000       0.000 
     0.000       0.000       0.000      ELECTRON
 **** Density scaling factors start at location        1 and end at    
4000 (I*4 addr.) ****
 ***** Next control card *****   BEAMPOS    0.000       0.000       0.000 
     0.000       0.000       0.000
 ***** Next control card *****   GEOBEGIN   0.000       0.000       0.000 
     0.000       0.000       0.000      COMBNAME
1
                    IVOPT = 0          IDBG = 0
                                                  Body data
 *
 *---------------    Black body                                           
                           ---------------*
 *
 *
 *---------------    Void sphere                                          
                           ---------------*
 *
  SPH blkbody          1  0.00000000E+00  0.00000000E+00  0.00000000E+00 
0.10000000E+06  0.00000000E+00  0.00000000E+00         5
 *
 *---------------    Cylindrical target                                   
                           ---------------*
 *
  SPH void             2  0.00000000E+00  0.00000000E+00  0.00000000E+00 
0.10000000E+05  0.00000000E+00  0.00000000E+00        15
  RCC target           3  0.00000000E+00  0.00000000E+00  0.10000000E+00 
0.00000000E+00  0.00000000E+00  0.10000000E+02        25
                          0.50000000E+01
  END                  4  0.00000000E+00  0.00000000E+00  0.00000000E+00 
0.00000000E+00  0.00000000E+00  0.00000000E+00        36
 Number of bodies         3
 Length of FPD-Array     41
  Accuracy parameter:  3.89E-04 suggested:  3.89E-04
  Detected max. geometry dimension:  1.00E+04 cm
                                                  Input region data
 *
 *---------------    Black hole                                           
                           ---------------*
 *
BLKBODY      5 +        blkbody  -        void
 *
 *---------------    Void around                                          
                           ---------------*
 *
VOID         5 +        void     -        target
 *
 *---------------    Target                                               
                           ---------------*
 *
TARGET       5 +        target
  END
 Number of input regions       3
 Number of  code zones         3
 Length of integer array      78
 CODE ZONE    INPUT REGION   ZONE DATA LOC.    NO. OF BODIES    REGION NO.
      1            1              25                 2               1
      2            2              34                 2               2
      3            3              43                 1               3
    I    IR1(I)    IR2(I)
    1       1         1
    2       2         2
    3       3         3
 Interpreted body echo
 Body n.        1  SPH  blkbody    Rot.     0
          0.000000              0.000000              0.000000            
 100000.0
 Body n.        2  SPH  void       Rot.     0
          0.000000              0.000000              0.000000            
 10000.00
 Body n.        3  RCC  target     Rot.     0
          0.000000              0.000000             0.1000000            
 0.000000              0.000000              10.00000
          5.000000
 Interpreted region echo
 Region n.        1 BLKBODY
                  1       -2
 Region n.        2 VOID
                  2       -3
 Region n.        3 TARGET
                  3
1 OPTION  0 WAS USED IN CALCULATING VOLUMES, FOR     3 REGIONS
 3: INPUT VOLUMES, ANYTHING ELSE: VOLUMES = 1.0
           VOLUMES (CM**3)
1   REG        1           2           3
 VOLUME   1.000E+00   1.000E+00   1.000E+00
0NGEOM=      4002,   NGLAST=      4218
 **** Stars/Energy accumulation arrays start at location     4221 and end
at     4244 (I*4 addr.) ****
 **** Minimum step size array start at location     4247 and end at    
4256 (I*4 addr.) ****
 **** Maximum step size array start at location     4259 and end at    
4268 (I*4 addr.) ****
GEOEND
 Total time used for geometry initialization:  9.999E-03 s
 ***** Next control card *****   RANDOMIZ   1.000       0.000       0.000 
     0.000       0.000       0.000
 RM64 INITIALIZED:       98765     12345           0         0
 ***** Next control card *****   ASSIGNMA   1.000       1.000       0.000 
     0.000       0.000       0.000
 ***** Next control card *****   ASSIGNMA   2.000       2.000       0.000 
     0.000       0.000       0.000
 ***** Next control card *****   ASSIGNMA   12.00       3.000       0.000 
     0.000       0.000       0.000
 ***** Next control card *****   EMFCUT   -3.0000E-03  3.0000E-03   0.000 
     1.000       710.0       0.000      PROD-CUT
 ***** Next control card *****   EMFCUT   -2.0000E-03  2.0000E-03   0.000 
     1.000       3.000       0.000
 ***** Next control card *****   START      100.0       0.000       0.000 
     0.000       0.000       0.000
 ***** EMF   tabulations in blank common start at  location10542885 *****
 ***** end at location10566262 (I*4 addresses) *****
1 Quantities/Biasing associated with each media:
 COPPER
      Rho =    8.96000     g/cm**3     Rlc=    1.43559     cm
      Ae =    3.51100     MeV    Ue =    13963.8     MeV
      Ap =    3.00000     MeV    Up =    13963.3     MeV
      dE/dx fluctuations activated for this medium, level  1
      below the threshold for explicit secondary electron production
     (up to  2I discrete levels, up to   2 K-edges)
1 Correspondence of regions and EMF-FLUKA material numbers and names:
       Region          EMF                     FLUKA
          1         0   VACUUM               1  BLCKHOLE
     Ecut =  0.0000E+00 MeV,    Pcut =  0.0000E+00 MeV,    BIAS = F,  
Ray. = F,   S(q,Z) = F,  Pz(q,Z) = F
          2         0   VACUUM               2  VACUUM
     Ecut =  0.0000E+00 MeV,    Pcut =  0.0000E+00 MeV,    BIAS = F,  
Ray. = F,   S(q,Z) = F,  Pz(q,Z) = F
          3         1   COPPER              12  COPPER
     Ecut =  3.5110E+00 MeV,    Pcut =  3.0000E+00 MeV,    BIAS = F,  
Ray. = F,   S(q,Z) = T,  Pz(q,Z) = F
--------------040804000501030200090709--
Received on Thu Apr 21 2011 - 00:05:39 CEST

This archive was generated by hypermail 2.2.0 : Thu Apr 21 2011 - 00:05:42 CEST