Dear FLUKA experts,
I am running a simulation where I irradiate a PbWO4 crystal with protons to understand the hadron damage produced. The features of the experimentally observed damage hint at lattice disorder that might be caused by fragments of heavy elements that have a very short track length and a high energy deposition in the crystal. Therefore, I am trying to retrieve the length track and energy deposition of the heavy fission fragments with the mgdraw.f user routine, entry USRDRAW. Find attached my mgdraw.f file and the output file (pwo_fragments001_fort.44) I get.
I obtain the information of the heavy fragments produced after each interaction the following way:
do 20 ip = 1, NPHEAV
WRITE( 44.0, '(A2,I5.5,1X,I3,1X,A2,1X,I2,1X,I2,
&1X,F10.6,1X,F10.6,1X,F10.6,1X,F10.6,1X,I2,1X,I2,1X,I2,1X,I2,1X,
&F10.6,1X,F10.6,1X,F10.6)')'H ',NCASE,ICODE,'-',ip,KHEAVY(ip),
&TKHEAV(ip)+AMNHEA(KHEAVY(ip)),CXHEAV(ip),
&CYHEAV(ip),CZHEAV(ip),NTRACK,MTRACK,ICHEAV(KHEAVY(ip)),
&IBHEAV(KHEAVY(ip)),AMNHEA(KHEAVY(ip))
20 continue
and I get an output of the type:
H 00001 101 - 1 6 3.738171 -0.910575 -0.127528 0.393178 0 0 2 4 3.727380
i.e. NTRACK=0 and MTRACK=0. From the manual, I understand this means that the energy deposition is point-like, so it's not recorded in (TRACKR). I have seen other methods to get the energy deposition, but I see some problems and I wonder which one would suit best my goals:
(1) Scoring the energy deposition with USRBIN, using AUXSCORE to filter the heavy fragments: this way I cannot identify what energy deposit belongs to what fragment, right?
(2) Using the comscw.f routine: I could get the generalized particle id (IJ), the deposited energy (RULL) and the particle generation (LLO), but how do I know the ID in the case of the heavy fragments? Should I just add the (FHEAVY) include file and get the information from there?
I would really appreciate if you could give me a hint as to what path to follow for my study.
Best regards,
Cristina
__________________________________________________________________________
You can manage unsubscription from this mailing list at
https://www.fluka.org/fluka.php?id=acc_info
Received on Tue May 09 2017 - 14:06:57 CEST