Dear Fluka experts,
In continuation of the problem I came up with couple of days ago I want to add few more things. After removing “SPECSOUR” card from my input file while keeping DPMJET physics process and executing FLUKA after linking DPMJET to it , “core dump” problem got removed. So basically the input file is with just unidirectional Au beam with momentum of 100GeV/nucleon. For p+p collision by keeping “SPECSOUR” card it takes forever to finish even 1 FLUKA cycle. Is it some limitation regarding SPECSOUR card for invoking heavy ion collision or something related to it is missing in my input file ? I am hereby attaching my input file. The collision vertex is defined in vacuum within +/- 0.001 cm along X-Y-Z coordinate and it has been divided into upstream and downstream along Z axis by XYP cut plane.
Thanks
-Sourav
On Nov 9, 2014, at 12:44 PM, Sourav Tarafdar <Sourav.Tarafdar_at_weizmann.ac.il<mailto:Sourav.Tarafdar_at_weizmann.ac.il>> wrote:
Dear Francesco,
Thanks for the suggestions. I implemented your suggestions as far as defining the collision vertex is concerned , i.e., inclusion of BEAM card and HI-PROPE card on top of SPECSOUR card. Also I removed Au material region relevant for my collision vertex. However I got the error on my terminal which looks like
======================= Running FLUKA for cycle # 1 =======================
$FLUPRO/flutil/rfluka: line 359: 13280 Floating point exception(core dumped) "${EXE}" < "$INPN" 2> "$LOGF" > "$LOGF"
Just wondering whether I am missing something in my input file ? My error file and output file are completely empty.
Thanks
-Sourav
On Nov 9, 2014, at 12:28 AM, Francesco Cerutti <Francesco.Cerutti_at_cern.ch<mailto:Francesco.Cerutti_at_cern.ch>> wrote:
Dear Sourav,
1) Shall I have to define Target region of Au material at collision vertex
?
not at all (unless your collisions are supposed to take place in gold instead of vacuum...). You defined the second beam species in the SPECSOUR card, where you have also said that your first beam is made by HEAVYIONs. So you need a HI-PROPE card to specify the HEAVYION nature. And you still need a BEAM card (where you shall put HEAVYION - making redudant your SPECSOUR WHAT(11) setting - and an energy per nucleon exceeding 100 GeV/n for transport initialization purposes).
2) Which scoring card should be used in order to retrieve kinematic
information of produced particles and nuclear remnants after H.I. collision
? So far I have used �USRYIELD� card.
USRYIELD does not work yet for the SPECSOUR collision event. You can inspect the properties of the products by customizing the usrein.f routine - which is automatically called before source particles start to be transported - in order to look at the FLKSTK content (to this purpose remember to add the
INCLUDE '(FLKSTK)'
statement).
Cheers
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
On Sat, 8 Nov 2014, Sourav Tarafdar wrote:
Dear Fluka users,
I have been trying to simulate Au+Au collision at C.M. energy of 200
GeV/nucleon using DPMJET3 within FLUKA framework. For defining the colliding
beam I used SPECSOUR card with PPSOURCE type. I am mostly interested in
getting the kinematic variables (px, py, pz, energy , pid etc) of the
particles after the collision and also the remnants of the nuclei after
collision. My questions are
1) Shall I have to define Target region of Au material at collision vertex
?
2) Which scoring card should be used in order to retrieve kinematic
information of produced particles and nuclear remnants after H.I. collision
? So far I have used �USRYIELD� card.
Any suggestions will be of great help.
Please find the attached input file for FLUKA and the output file after
executing FLUKA using the command $FLUPRO/flutil/rfluka -e flukadpm3 test
Any suggestions will be of great help.
Thanks
-Sourav
<fluka_dpm_wotarg.inp>
Received on Mon Nov 10 2014 - 02:40:10 CET