Re: [fluka-discuss]: magnetic field

From: Mina Nozar <nozarm_at_triumf.ca>
Date: Fri, 14 Feb 2014 09:35:17 -0800

THANK YOU Alfredo and Mario!

Best wishes,
Mina

On 14-02-13 11:16 PM, Alfredo Ferrari wrote:
> BTY = -ONEONE
>
> will do the job. Of course the three cosines which define the magnetic
> field direction can have whichever value/sign, with the only
> constraint that
>
> BTX**2 + BTY**2 + BTZ**2 = ONEONE
>
>
>
>
> +----------------------------------------------------------------------+
> | Alfredo Ferrari || Tel.: +41.22.76.76119 |
> | CERN-EN/STI || Fax.: +41.22.76.69474 |
> | 1211 Geneva 23 || e-mail: Alfredo.Ferrari_at_cern.ch |
> | Switzerland || |
> +----------------------------------------------------------------------+
>
> On Thu, 13 Feb 2014, Mina Nozar wrote:
>
>> 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 - 19:38:00 CET