TOSCA simulation for Pol He3 Field
From HallCWiki
Revision as of 12:35, 9 September 2020 by Jixie (Talk | contribs) (→How to run these scripts to combine maps)
By Jixie Zhang, Sep 9th,2020
Contents
Introduction
The TOSCA simulation geometry files were prepared by Steven Lassiter.
The simulation were performed by both Steven and Jixie.
In D2N, the target base was not rotated. The axis of horizontal large coils (HL) is along Lab X axis.
The axis of horizontal small coils (HS) is along Lab Z axis.
In A1N, the target base was rotated clockwise (top view) about Y axis by 45 degrees.
Please keep in mind that the coordinate system in the TOSCA simulation is defined in this way:
X axis goes to beam right, y axis goes up, Z axis goes to upstream.
We have rotated the output about Y axis by 180 degrees to make it in the Lab frame.
Where is the maps located
The rotated TOSCA field maps for each single set of coils are stored in group dir:
ifarm1901:/u/group/c-polhe3/Users/jixie/TOSCA/A1NMap
ifarm1901:/u/group/c-polhe3/Users/jixie/TOSCA/D2NMap
Scripts to combine maps
Jixie provides python code to combine these 7 rotated individual maps (HB, HL, HS, VL, VS, HLCC and HSCC)
into one map with each map scaled to the given current. These scrips are stored at
ifarm1901:/u/group/c-polhe3/Users/jixie/TOSCA/Scripts
How to run these scripts to combine maps
- ) configure the 7 component maps in 'Field map block', which is just a list.
The structure of this list is: Path_of_table, Current_of_table, Target_current. Jixie will provide the Current_of_table for each map. What you need to change is the "Path_of_table" and "Target_current" for each of them. ("Target_current" is the current go through the coils, not the current you read back in the GUI monitor.)
- ) Keep in mind that in D2N experiment, HS coil current in read back is in opposite
sign, VL and VS are switched. In A1N experiment, HS coil current in read back is also in opposite sign. When you provide "Target_current", you really need to be careful.
- ) This script can make some plot using functions. For example "PlotGradientVsZ(xx,yy,figname)"
will plot the 9 gradients as a function of z at line (x=xx,y=yy). User need to specify xx and yy locations and output figure name.