Created magnetic field used mgfield.f routine

Physics, Transport and Magnetic Fields
1

Dear FLUKA expert,
I am trying to add a solenoid magnet to my simulation․To approximate the magnetic field of the solenoid, the simulation uses the Biot-Savart law. I did this using a simplified expression.

Picture1
Picture11490×296 65.9 KB

As a result of FLUKA, we get the following magnetic field.

Picture2
Picture21192×450 299 KB

If we compare the FLUKA results for the magnetic field with the expression, we see that they agree well.
But I can’t understand the results that get from USRBIN [BEAMPART], which you can see below. In the region of the magnetic field, the intensity is incorrect.

Picture3
Picture3742×450 180 KB

Please can you help me understand what I am missing or doing wrong?
Also, you could find input files.

magfld_Solenoid.f (3.1 KB)
solenoid.flair (3.4 KB)
solenoid.inp (1.5 KB)

Thanks in advance for any help!
Best regards
Zohrab

3

Dear @zamirkhanyan ,

the direction cosines of the magnetic field must be normalised as indicated in the manual.

A modified magnetic field routine is here
magfld_Solenoid.f (3.4 KB)
where I implemented the normalisation.

I didn’t check the correctness of your code, but it is recommend to use everywhere DOUBLE PRECISION variables, i.e. 0.5 should be replaced by 0.5D0

Best, Luigi

1 Like
4

Dear Luigi,
Thank you VERY much for your help!
Best regards
Zohrab