problems for neutrino energies above 50 TeV in charm in charged current interactions

From: Salvatore Mangano <manganos_at_ific.uv.es>
Date: Tue, 03 Jul 2012 18:32:40 +0200

Dear Fluka experts,

I am using Fluka (Fluka 2011.2.11) for studies
of charm in charged current interactions.
Everything seem to work fine for neutrino energies below 10 TeV.
Unfortunately for neutrino energies above around 50 TeV the run crashes.
I tried several different neutrino energies, but cannot make a clearer
statement.
I am using the input card given in the attachment and the crash occurs
after 329 events.
Can you reproduce the problem?
Is 10 TeV the maximal neutrino energy for dimuon production or should it
work?
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 5.2E+04 NEUTRIM
BEAMPOS 0.0 0.0 4.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+07
* 2 vacuum
   SPH 2 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00
            1.000000000000000D+07
* sea: radius 9 km
   SPH 104 0.000000000000000D+00 0.000000000000000D+00 0.000000000000000D+00
            9.000000000000000D+06
* plane to sample the passing-through particles
   XYP plane1 0.000000000000000D+05
   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
*ASSIGNMAT LEAD 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 1500.0
STOP
Received on Wed Jul 04 2012 - 11:05:09 CEST

This archive was generated by hypermail 2.2.0 : Wed Jul 04 2012 - 11:05:10 CEST