Re: [fluka-discuss]: REPLICATION USING LATTICE CARD

From: Luigi Salvatore Esposito <luigi.salvatore.esposito_at_cern.ch>
Date: Fri, 2 Aug 2019 11:35:14 +0000

Dear Riya Dey,
please find enclosed a working example with an USRBIN scoring ENERGY.
Note that WHAT(5) of the USRBIN should be larger than the number of lattices you have in your geometry.

For reading the results, you need to converte the output into an ASCII file (with usbrea program present in flutil) and
then post-process this file in order to associate each value to the corresponding region and lattice.

After the header of the ASCII file, there is an array in a 10-column format, to be read row after row (a single
column has no specific meaning).
It contains the energy deposition matrix (in GeV per primary particle), element after element.

To properly interpret it, you have to consider that the matrix has 2 indices (the first for the regions and the
second for the lattices, as indicated in the header).
In the sequence of matrix elements, the first index runs first over its interval, then the second.
So the list is
E(1,0), E(2,0), … E(#reg,0)
E(1,1), E(2,1), … E(#reg,1)
….
E(1,#latt), E(2,#latt) …. E(#reg,#latt)
where #reg is the number of regions (value of WHAT(4)) and #latt is the number of lattices (value of WHAT(5)).

In other words, elements with the second index = 0 represent regions that to not belong to any lattice,
elements with second index =1 represent regions that belong to the first lattice, etc.

Hope this could help.

Best luigi




On 2 Aug 2019, at 10:36, riya dey <riyadey_at_barc.gov.in<mailto:riyadey_at_barc.gov.in>> wrote:

Dear Luigi,

I have gone through that link but the solution is not clear to me . If I use SPECIAL in TYPE in USRBIN card, what should I put in I1 min, I1 max, I2 min, I2max, F3 min, F3 max, N3 ? so that I can get energy scoring in all the spheres ?

Can you please give me one example by putting values in the previously sent input? It will be very helpful for me.

Thanks and Regards


Riya Dey
=========================

From: owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it> [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Luigi Salvatore Esposito
Sent: 02 August 2019 03:39
To: riya dey <riyadey_at_barc.gov.in<mailto:riyadey_at_barc.gov.in>>
Cc: fluka discussion <fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>>; riyadey0293_at_gmail.com<mailto:riyadey0293_at_gmail.com>
Subject: Re: [fluka-discuss]: REPLICATION USING LATTICE CARD

Dear Riya Dey,
for the scoring you can use the USRBIN with WHAT(1) = 8.
More details in the manual and from the slides about advanced geometry topics
https://indico.cern.ch/event/489973/contributions/2000434/
pages 51-52.

Hope this could help.
Best, luigi



On 1 Aug 2019, at 06:28, riya dey <riyadey_at_barc.gov.in<mailto:riyadey_at_barc.gov.in>> wrote:

Thank you for your help. The issue is solved. Now I can see spheres in the proto type as well as in the replicated regions. I have one more doubt regarding SCORING in the lattice region

Since, no region name has been assigned for spheres in the replicated region, now if I want to score energy or fluence in a sphere in the replicated region, how can I call it in the scoring card? Let’s say I want to score energy in all the spheres inside the replicated as well as proto type region. What should I put in WHAT (1) card (TYPE)?


Thanks and Regards


Riya Dey

From: owner-fluka-discuss_at_mi.infn.it<mailto:owner-fluka-discuss_at_mi.infn.it> [mailto:owner-fluka-discuss_at_mi.infn.it] On Behalf Of Luigi Salvatore Esposito
Sent: 31 July 2019 18:58
To: riya dey <riyadey_at_barc.gov.in<mailto:riyadey_at_barc.gov.in>>
Cc: fluka-discuss_at_fluka.org<mailto:fluka-discuss_at_fluka.org>; riyadey0293_at_gmail.com<mailto:riyadey0293_at_gmail.com>
Subject: Re: [fluka-discuss]: REPLICATION USING LATTICE CARD

Dear Riya Dey,

see my answers inline.

Best regards, luigi



On 31 Jul 2019, at 13:28, riya dey <riyadey_at_barc.gov.in<mailto:riyadey_at_barc.gov.in>> wrote:

Dear FLUKA Experts,

A sample input using LATTICE card is attached here. My intention is to replicate the region called “BOUND” (which is an RPP) in 6 directions with transformations delta x = + 0.5 , delta x = - 0.5; delta y = + 0.5 , delta y = - 0.5; delta z = + 0.5 , delta z = - 0.5 .

In the input I have made two replicas considering transformations delta x = + 0.5 and delta z = + 0.5. I have the following doubts:

1. For every replication, I am adding one ROT-DEFI card; start transform card and end transform card; within which I am adding the proto type regions and giving them different names for replication; next regions with the names of proto type and the replicated regions are defined in the region card; one lattice card is added for every replication with different transformation index. At the same time all these regions (proto type and the replicated regions) should be omitted from the void region.

Is this the correct way of lattice replication? Because if there are many prototype regions which need to be replicated 5 times in all 6 direction, then it will end up with a messy input with lot of different replicated region names. Is there any other way of replication by adding single transform card (at least along each axis) ?

For the replicas, you need only in external container bodies (REPBOUND and REPBOUN1 in your input.
No need to replicate the spheres inside. All geometry details need to be defined for the basic unit only.
This is actually the purpose of the lattice!

In addition, there was an error in the LATTICE card definition: you need to use the inverse transformation that is the one bringing from the replica to the basic unit.

See these modifications implemented in the attached input file.


2. While viewing the geometry, the material in the replicated region is showing vacuum although water is assigned in both proto type and replicated regions . Does this mean that the replication is not correctly done?

You need to switch to lattice layer to see inside the replica. See attached snapshot.





Thanks and Regards

Riya Dey

<geometry.PNG><lattice.inp>


<image001.png>




__________________________________________________________________________
You can manage unsubscription from this mailing list at https://www.fluka.org/fluka.php?id=acc_info

Received on Fri Aug 02 2019 - 15:06:46 CEST

This archive was generated by hypermail 2.3.0 : Fri Aug 02 2019 - 15:06:53 CEST