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

From: Francesco Cerutti <Francesco.Cerutti_at_cern.ch>
Date: Thu, 8 Aug 2019 19:35:30 +0200

Dear Riya Dey,

- TIMED OUT is not necessarily an error tag. In principle the run might be
still ongoing without writing since a while. More often, it actually hides
an abort, but to find the reasons one has to look into the log, err, out
files in the fluka_# temporary directory (TIMED OUT as such does not give
any clue).

- the parameters of the special USRBIN you input look fine at first glance

- I do not fully understand the meaning of the region 1003, which may
imply a geometry mistake. If your 1000 spheres are defined inside the 'BOX
0' box and are replicated (by proper box transformations) in the other 124
boxes, then 'BOX 0' (as the full container body) cannot be a region
itself, since it would overlap with the contained spheres (one should
rather define the region obtained by the subtraction of all spheres from
the prototype container).

Kind regards

Francesco

**************************************************
Francesco Cerutti
CERN-EN/STI
CH-1211 Geneva 23
Switzerland
tel. +41 22 7678962

On Thu, 8 Aug 2019, riya dey wrote:

>
> Dear Luigi,
>
>  
>
> Thank you for explanation. But I am getting “TIMED OUT” error by using the following conditions in an
> input. Can you please tell me whether I1, I2 and F3 are properly assigned?
>
>  
>
> In the input,
>
>  
>
> Region 1 : BLKBDY
>
> Region2 : Void
>
> Region 3 to Region 1002 : spheres (with index 1 to 1000)
>
> Region 1003: BOX 0 (proto type box)
>
>  
>
> There are 124 lattices i.e. replications of proto type box.
>
> Region 1004 to Region 1127 : BOX 1 to BOX 124.
>
>  
>
>  
>
> I want to score energy deposition in spheres only. So I put
>
>  
>
> I1 min (WHAT 7) = 3
>
> I1 max (WHAT 4) = 1002
>
>  
>
>  
>
> I2 min (WHAT 8) = 0
>
> I2 max (WHAT 5) = 130
>
>  
>
> F3 min (WHAT 9) = -10000
>
> F3 max (WHAT 6) = 10000
>
>  
>
>  
>
>  
>
>  
>
>  
>
> Thanks and Regards
>
>  
>
> ~~~~~~~~~~~~~~~~~~~~~~~~~
>
> Riya Dey
>
>  
>
> From: 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 17:05
> To: riya dey <riyadey_at_barc.gov.in>
> Cc: fluka discussion <fluka-discuss_at_fluka.org>; riyadey0293_at_gmail.com
> Subject: Re: [fluka-discuss]: REPLICATION USING LATTICE CARD
>
>  
>
> 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> 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] On Behalf Of Luigi
> Salvatore Esposito
> Sent: 02 August 2019 03:39
> To: riya dey <riyadey_at_barc.gov.in>
> Cc: fluka discussion <fluka-discuss_at_fluka.org>; 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> 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] On Behalf
> Of Luigi Salvatore Esposito
> Sent: 31 July 2019 18:58
> To: riya dey <riyadey_at_barc.gov.in>
> Cc: fluka-discuss_at_fluka.org; 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> 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 Thu Aug 08 2019 - 20:43:43 CEST

This archive was generated by hypermail 2.3.0 : Thu Aug 08 2019 - 21:26:15 CEST