Re: [fluka-discuss]: magnetic field

From: Mina Nozar <nozarm_at_triumf.ca>
Date: Thu, 13 Feb 2014 19:02:23 -0800

Hello everyone,

How do I set the magnetic field direction if it is in the -y direction?

Uniform dipole field along +y:
          BTX = ZERZER
          BTY = ONEONE <====
          BTZ = ZERZER
          B = DIPFLD

What should BTY be set to?

Thanks,
Mina


On 14-02-12 06:04 PM, Mina Nozar wrote:
> Ok, I was bitten by the 8 character region name feature!
>
>
> What doesn't make sense is that when printing region numbers out for
> "rRMVVac" and "rRMHVac", I saw the correct region number for "rRMHVac"
> region but -1 for "rRMVVac". What doesn't make sense is that both
> region names are 7 characters and I was not putting a space after the
> "c" in the region name for either but was seeing the correct region
> number for "rRMHVac"!
>
> Mina
> =============================
> * Set dipole filed in tesla (B field is 250 Gauss)
> DIPFLD = 0.0250D+00
>
> * Some useful printout
> WRITE (LUNOUT,*) ''
> WRITE (LUNOUT,*) ' magfld.f called for EHDT Raster Magnets'
> WRITE (LUNOUT,*) ' Dipole field: ',DIPFLD,' T'
> WRITE (LUNOUT,*) ''
>
> * Retrieve the indices of the magnetic regions starting from their
> * names, and store them
>
> CALL GEON2R ("rRMVVac ",NRMDV, IERR)
> CALL GEON2R ("rRMHVac ",NRMDH, IERR)
>
> WRITE (LUNOUT,*) ''
> WRITE (LUNOUT,*) ' "rRMVVac" region number is: ', NRMDV
> WRITE (LUNOUT,*) ' "rRMHVac" region number is: ', NRMDH
> WRITE (LUNOUT,*) ''
>
> LFIRST=.FALSE.
> END IF
>
>
>
> On 14-02-12 02:25 PM, Mina Nozar wrote:
>> Hello everyone,
>>
>>
>> How does one know whether a magnetic field has actually been applied
>> to a region? I can't see anything in the output files that conveys
>> this message.
>>
>> In my first attempt to set the magnetic field for a dipole magnet, I
>> modified magfld.f, assigned the appropriate B field to a dipole
>> sensitive region (vacuum) and in the ASSIGNMA card, and set magnetic
>> field on. I compiled magfld and built the executable linking to
>> magfld.o. The magnetic field is 250 G, so I put in 0.025 T for the
>> intensity in magfld.f. The results didn't show any deflections in
>> the backscattered e+/e- in the beam pipe from a 75 MeV beam dump.
>> Most of these e+/e- have energies less than 10 MeV (with the peak
>> around 1-2 MeV).
>>
>> Reading through the manual, I realized I forgot to include the
>> MGNFIELD card in the input file. So I included this card, recompiled
>> and ran again. And I can't see any difference.
>>
>> MGNFIELD 0.0 0.0 0.0
>>
>>
>> Thanks,
>> Mina
>>
>
Received on Fri Feb 14 2014 - 05:22:56 CET