Re: [fluka-discuss]: magnetic field

From: Alfredo Ferrari <alfredo.ferrari_at_cern.ch>
Date: Fri, 14 Feb 2014 08:16:55 +0100

  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 - 09:12:14 CET