charm in charged current: core dump for high beam energy

From: Salvatore Mangano <manganos_at_ific.uv.es>
Date: Mon, 07 Nov 2011 16:28:07 +0100

Dear Fluka experts,

I am using Fluka (Fluka 2011.2.6) for studies
of charm in charged current interactions.
The goal should be to select dimuon events.

My file.inp (see also attachment) does something
for a BEAM energy of 1.E+03 GeV, but gives me a core dump
if I increase the energy to 1.E+04 GeV.
Error message:
/sps/km3net/users/mangano/fluka/fluka_32/flutil/rfluka: line 358: 6653
Floating point exception(core dumped) "${EXE}" <"$INPN" 2>"$LOGF" >"$LOGF"

My file.inp file starts as follows:
BEAM 1.0E+04
NEUTRIM
BEAMPOS 0.0 0.0 1.E+05 NEGATIVE
DISCARD -27.0
PHYSICS 2.0 2.0 2.0 2.0 NEUTRINO

If you need some additional information please let me know.
Thank you for any suggestion.
Best regards,
  Salvatore Mangano

*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
TITLE
Tracciamento di particelle prodotte da interazione di neutrino in acqua
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* Units cm, GeV, sec, GeV/c
BEAM 1.0E+04 NEUTRIM
BEAMPOS 0.0 0.0 1.E+05 NEGATIVE
* YOU NEED THIS OTHERWISE NO NEUTRINO PROPAGATION
DISCARD -27.0
PHYSICS 2.0 2.0 2.0 2.0 NEUTRINO
*PHYSICS -2.0 -2.0 -2.0 2.0 NEUTRINO
USRICALL 0.0 SFERA
*SOURCE 0.0 0. 0. 0. 10. 200.
*
GEOBEGIN COMBNAME
     0 0 Only sea water
* 1 black hole
   SPH 1 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00
            5.000000000000000D+06
* 2 vacuum
   SPH 2 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00
            1.000000000000000D+06
* sea: radius 2 km
   SPH 104 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00
            2.000000000000000D+05
* plane to sample the passing-through particles
   XYP plane1 0.000000000000000D+03
   END
********** regions ****************
* reg. 1 external BH:
   EBH 5 1 -2
* reg. 2 external void:
   Vac 5 2 -104
* reg. 3 sea above the plane:
   Sea1 5 104 -plane1
* reg. 4 sea below the plane:
   Sea2 5 104 plane1
   END
GEOEND
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
ASSIGNMAT 1.0 1.0
ASSIGNMAT 2.0 2.0
* Water:
ASSIGNMAT WATER 3.0 4.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
* e+e- and gamma production threshold set at 1 GeV
*EMFCUT 900.0 900.0 1.0 3.0 29.0 PROD-CUT
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
USERDUMP 101.0 22.0 +2.0 Output
RANDOMIZE 1.0
*...+....1....+....2....+....3....+....4....+....5....+....6....+....7....+....8
*START 1000.0
START 100.0
STOP
Received on Tue Nov 08 2011 - 16:28:44 CET

This archive was generated by hypermail 2.2.0 : Tue Nov 08 2011 - 16:28:44 CET