Re: [fluka-discuss]: Reading out heavy fragments in FLUKA Heavy Ion collision using DPMJET-III

From: Sourav Tarafdar <Sourav.Tarafdar_at_weizmann.ac.il>
Date: Mon, 24 Nov 2014 05:40:09 +0000

Hi Sunil,

Sorry it was typo on my part. Indeed I am creating executable with link to DPMJET-III by using $FLUPRO/flutil/ldpmqmd -e executable_name -m fluke name_userroutine.o and then rfluka script is being used to run the created executable over FLUKA input file. My apologies for typo in my previous message . However the issue related to reading the heavy fragments still remains.

Thanks
-Sourav


On Nov 24, 2014, at 5:38 AM, Sunil C <csunil11_at_gmail.com<mailto:csunil11_at_gmail.com>> wrote:


Hi

I assume you invoked ldpmqmd to create an executable. Your command shows rfluka, which is a run script.

It should thus be

$FLUPRO/flutil/ldpmqmd -e executable_name -m fluka name_userroutine.o

Sunil


On Nov 23, 2014 3:57 PM, "Sourav Tarafdar" <Sourav.Tarafdar_at_weizmann.ac.il<mailto:Sourav.Tarafdar_at_weizmann.ac.il>> wrote:
Dear Fluka experts,

I was trying to read out heavy fragments kinematic information in heavy ion collision . I defined a collision region in the form of sphere filled with vacuum and used SPECSOUR and HI-PROPE to define my colliding beam of Au+Au at C.M. energy of 200GeV/nucleon. I tried doing that in two ways.

1) Customizing usreou.f where I made a loop over NPHEAV and tried writing the kinematic variables in user defined output file. Somehow the output file shows NPHEAV = 0 . Attached is my usreou.f file along with the corresponding input file named “fluk_dpm_usreou.inp” .

2) Customizing mgdraw.f routine where again I was making loop over NPHEAV. The corresponding input file is modified to have “USERDUMP" as scoring card. However in this case fluke crashes. I am attaching the mgdraw.f file along with the corresponding input file named “fluk_dpm_mgdraw.inp”.

In both the above cases I first compiled customized user routine to generate the corresponding “.o” file by using
$FLUPRO/flutil/fff name_userroutine.
Then I created the executable with link to DPMJET-III by using $FLUPRO/flutil/rfluka -e executable_name -m fluke name_userroutine.o
The corresponding executable is then used to run FLUKA using the rfluka script and FLUKA input file.



Can anyone please provide information related to the correct way of reading the heavy fragments kinematic information ?

Thanks
-Sourav
Received on Mon Nov 24 2014 - 08:03:41 CET