Dear Cristina,
in USDRAW your intercept your heavy fragments when they are generated, not
during their transport, which a track length is associated to.
To your purposes (damage, lattice disorder study), there is a much more
direct and proper way to proceed, that is calculation of DPAs (it has been
implemented in the code exactly to those purposes). You shall score
DPA-SCO with USRBIN, after activating ion transport and lowering hadron
thresholds as much as possible (as I explained last week). This way, also
the DPA contribution of low energy neutron products - not generated
explicitly - is taken into account. You will eventually get a spatial DPA
distribution, giving a picture of damage distribution and level.
Cheers
Francesco
**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. ++41 22 7678962
fax ++41 22 7668854
On Tue, 9 May 2017, Cristina Martin Perez wrote:
> 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 - 17:09:54 CEST