RE: lattice with new format

From: Mario Santana Leitner (Mario.Santana@cern.ch)
Date: Tue Jan 23 2007 - 22:57:15 CET

  • Next message: foucher@subatech.in2p3.fr: "dose calculation"

    Hi Fabio,

    Just replace 'number' by 'NAME' in the FLUKA manual if you are working with names. If you don't need to go to fancy, maybe the best way is:
    1) to define your prototype and enclose it by a simple body, e.g. an RCC.

      RCC ustbk 0.0 10000. 47.02 0.0 0.0 340.
      ...
      ...
    2) to define the placeholders that will host the repetitions of your prototype, e.g.:
      RCC ustbk01 0.000 0.000 51745.66 0.000 0.000 340.00 25.00
      RCC ustbk02 0.000 0.000 52132.66 0.000 0.000 340.00 25.00
      RCC ustbk03 0.000 0.000 52519.66 0.000 0.000 340.00 25.00
      RCC ustbk04 0.000 0.000 52949.45 0.000 0.000 340.00 25.00
      RCC ustbk05 0.000 0.000 53336.45 0.000 0.000 340.00 25.00
      ...
    3) with the corresponding dummy regions, e.g.:
      USTBK01 25 +ustbk01
      USTBK02 25 +ustbk02
      USTBK03 25 +ustbk03
      USTBK04 25 +ustbk04
      USTBK05 25 +ustbk05
      ...
    4) and then to use the LATTICE card to cite the association between the placeholder and the prototype (coordinade transformation), e.g.
    LATTICE USTBK01 USTBK01L ROT#1
    LATTICE USTBK02 USTBK02L ROT#2
    LATTICE USTBK03 USTBK03L ROT#3
    LATTICE USTBK04 USTBK04L ROT#4
    LATTICE USTBK05 USTBK05L ROT#5
    ...
      - here WHAT(1) is the placeholder, WHAT(2) is empty (thus WHAT(3)) and USTBK01L is the name of the lattice (use whichever you like).
      - the transformation is dictated by the SDUM. If it were empty then you would need to write a LATTIC routine. Here it refers to a roto-translation which is defined below.
    5) Specify in which consists that association for each replica.
    ROT-DEFI 101.0 0.0 0.0 0.0 10000.0 -51698.64
    ROT-DEFI 102.0 0.0 0.0 0.0 10000.0 -52085.64
    ROT-DEFI 103.0 0.0 0.0 0.0 10000.0 -52472.64
    ROT-DEFI 104.0 0.0 0.0 0.0 10000.0 -52902.44
    ROT-DEFI 105.0 0.0 0.0 0.0 10000.0 -53289.44
    ...
    In this case there is no rotation (WHAT(2)=WHAT(3)=0.0) but there is a translation in Y and Z, dY=10000.0, dZ=-51698.64,.... Note that if you move the objects defined in (2) by dX,dY,dZ you will fall into the prototype defined in (1), which is what you need to match prototype and repetitions.

    Mario
     

    ---
         WHAT(1) = "Container-region" NAME of the first lattice cell
                   ("From region WHAT(1)...")
                   No default
         WHAT(2) = "Container-region" NAME of the last lattice cell
                   ("...to region WHAT(2)...")
                   Default = WHAT(1)
         WHAT(3) = step length in assigning "Container-region" numbers
                   ("...in steps of  WHAT(3)").
                   Default = 1.
         WHAT(4) = lattice NAME of the first lattice cell, assigned to
                   region WHAT(1)
                   No default
         WHAT(5) = lattice NAME of the last lattice cell, assigned to
                   region WHAT(2)
         WHAT(6) = step length in assigning cell numbers
                   ("...in steps of  WHAT(6)").
                   Default = 1.
         SDUM    = possible index of transformation associated with this lattice.
                   Exceptionally, here SDUM can contain an integer number, in free
                   format, following any of the strings "ROT#", "Rot#", "rot#",
                   "RO#", "Ro#", "ro#. If any one of such strings is present, an
                   integer identifying the associated roto-translation is read in
                   the following characters. If no such string is found, the LATTIC
                   user routine will beg called whenever a transformation is
                   required.
    

  • Next message: foucher@subatech.in2p3.fr: "dose calculation"

    This archive was generated by hypermail 2.1.6 : Wed Jan 24 2007 - 09:18:41 CET