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