RE: lattice and magnetic field

From: Chrysostomos Valderanis <Chrysostomos.Valderanis_at_cern.ch>
Date: Tue, 3 Aug 2010 15:57:58 +0200

So. just to report my investigation on the problem. It seems that calls
to GEON2L do not work as I expected. If I remove all the GEON2L calls
and hardcode their outcome into the magfld.f file it works, at least
partially. If I include even one call to GEON2L the program crashes.
I really would like to avoid hardcoding all the numbers in the routine.
I some on could explain to me why GEON2L doesn't work the way I expect
it to work I would be grateful.

Thank you,
Makis

________________________________

From: owner-fluka-discuss_at_mi.infn.it on behalf of Chrysostomos
Valderanis
Sent: Sat 24/07/2010 16:28
To: msantana_at_slac.stanford.edu; fluka-discuss_at_fluka.org
Subject: RE: lattice and magnetic field

Hi again,

I implemented your suggestion. To be fair, I implemented the way I
understood it. But now fluka crashes. It crashes somewhere in the
initialization phase. I tried my input file without magnetic field and
it works. I also tried to comment out different sections of the magfld
routine but it makes no difference. Fluka crashes to the same point. I
scan the out file for signs of problems but I couldn't find any. It
crashes to the point where it tries to give a summary of the lattice
regions.
Any hint abount what is going on is very welcomed.

Thank you,
Makis

________________________________

From: owner-fluka-discuss_at_mi.infn.it on behalf of Santana, Mario
Sent: Fri 23/07/2010 19:20
To: Chrysostomos Valderanis; fluka-discuss_at_fluka.org
Subject: RE: lattice and magnetic field

Hi Makis,

What you are asking to do is possible and it is actually quite common.

You need the following instruction in your magfld.f routine

      CALL GEOL2N(MLATTC,LATNAM,IRTLAT,Ierr)

And then a code that looks like in the example below:

[...]
        IF (LATNAM.eq."QU01L ".or.LATNAM.eq."QU03L ".or.
     $ LATNAM.eq."QU05L ".or.LATNAM.eq."QU07L ".or.
     $ LATNAM.eq."QU09L ".or.LATNAM.eq."QU11L ".or.
     $ LATNAM.eq."QU13L ".or.LATNAM.eq."QU15L ".or.
     $ LATNAM.eq."QU17L ".or.LATNAM.eq."QU19L ".or.
     $ LATNAM.eq."QU21L ".or.LATNAM.eq."QU23L ".or.
     $ LATNAM.eq."QU25L ".or.LATNAM.eq."QU27L ".or.
     $ LATNAM.eq."QU29L ".or.LATNAM.eq."QU31L ".or.
     $ LATNAM.eq."QU33L ") THEN
           ck =3D 3.85D-01 * 7.8/7.32 ! effective length=3D7.8 cm
        ELSE IF (LATNAM.eq."QU02L ".or.LATNAM.eq."QU04L ".or.
     $ LATNAM.eq."QU06L ".or.LATNAM.eq."QU08L ".or.
     $ LATNAM.eq."QU10L ".or.LATNAM.eq."QU12L ".or.
     $ LATNAM.eq."QU14L ".or.LATNAM.eq."QU16L ".or.
     $ LATNAM.eq."QU18L ".or.LATNAM.eq."QU20L ".or.
     $ LATNAM.eq."QU22L ".or.LATNAM.eq."QU24L ".or.
     $ LATNAM.eq."QU26L ".or.LATNAM.eq."QU28L ".or.
     $ LATNAM.eq."QU30L ".or.LATNAM.eq."QU32L ") THEN
           ck =3D -3.85D-01 * 7.8/7.32 ! effective length=3D7.8 cm
        END IF

[...]

Hope this helps,

Mario

I am trying to build a geometry with lattice and magnetic fields and I
am f=
acing some problems.

First question. My reading of the online courses says that the magnetic
field should only be set for the prototype region. Since I have more
than one magnetic field, I am using a routine (magfld.f). But what
should be written in this file?
-Should I find the prototype regions with GEON2R and set the magnetic
fields for them?
-Should I find the lattice regions with GEON2R and set the magnetic
field for the prototype regions they correspond to?
-Should I find the lattice regions and set the magnetic field as it
should be in the replica area?
For the moment I am using option2 but it doesn't seem to work. Anyway I
wou, ld prefer to calculate
the magnetic field just for the prototypes as it is a much simpler
expression.

Second question. I really would like to be able to set two replicas of
the same prototype to different magnetic fields. Is this possible? How?

Thank you in advance,
Makis
Received on Tue Aug 03 2010 - 20:27:38 CEST

This archive was generated by hypermail 2.2.0 : Tue Aug 03 2010 - 20:27:44 CEST