RE: lattice and magnetic field

Hi Makis

Please post input and user routines so that we can check.

Cheers
Stefan

On Tue, 3 Aug 2010, Chrysostomos Valderanis wrote:

> 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
>
>

-- 
___________________________________
Stefan Roesler
CERN, DGS-RP
CH-1211 Geneva 23
Switzerland
Phone:  +41-22-7679891
Fax:    +41-22-7669639
E-mail: Stefan.Roesler_at_cern.ch