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