Commit e78492f2 authored by Thomas Planche's avatar Thomas Planche

starting to write report

parent 5765b6bd
**HRS FIT TOSCA !!Simplest HRS (i.e. magnifying section quads OFF) using a TOSCA map
**HRS FIT TOSCA
'OBJET' 1
'OBJET'
546.198 !60 keV U1+ = 546.198
2 !2: All the initial coordinates must be entered explicitly
1 1 !total number of particles; number of distinct momenta
0. 0. 0. 0. 0. 1. ’o’ !Y; T; Z; P; S; D; 'marker'
1 !1 or -9 (-9 disables the tracking of this particle)
'TOSCA' 2
0 2 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
'TOSCA'
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling (adjusted so that R at the center=120cm+60cm*sqrt(2), 120cm is reference radius of the HRS, 60cm is the amound by which the map has been enlarge to include fringe fields); X coordinate scaling; Y; Z
HEADER_2 ! Title
181 47 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
......@@ -17,13 +17,13 @@ HEADER_2 ! Title
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
180.1396693 0. 180.1396693 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
200.0 0. 200.0 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'FIT' 4
'FIT'
1 !number of variables
2 70 2.072 0.5 !element#2 (TOSCA); parameter #70 (RE); correlation set so that RE=RS (which is parameters#72 of element #2)
1 !number of constraints
3 1 2 2 0.0 1. 0 !3:particle coordinate; 1: particle #1; 2: Y(i.e. R); where?:at the end of element #2; wanted value=0.0; weight; 0:no additional parameters
'END' 6
'END'
This diff is collapsed.
**HRS FIT TOSCA !!Simplest HRS (i.e. magnifying section quads OFF) using a TOSCA map
**HRS FIT TOSCA
'OBJET' 1
546.198 !60 keV U1+ = 546.198
......@@ -8,7 +8,7 @@
1 !1 or -9 (-9 disables the tracking of this particle)
'TOSCA' 2
0 2 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
0 0 !print the map (no print = 0); output along partile(s) trajectory(ies) = 2
1.0365 1. 1. 1. !Magnetic field scaling (adjusted so that R at the center=120cm+60cm*sqrt(2), 120cm is reference radius of the HRS, 60cm is the amound by which the map has been enlarge to include fringe fields); X coordinate scaling; Y; Z
HEADER_2 ! Title
181 47 1 22. ! Nb of nodes in Theta; R; Z (=1 for 2D map); MOD(see manual, =0 for Cartesian, with mid-plane symmetry)
......@@ -17,7 +17,7 @@ HEADER_2 ! Title
4 !IORDRE(=2, 4 or 25 if 2D, unused if 3D)
0.1 !Integration step size[cm]
2 !KPOS, normally=2
180.1396693 0. 180.1396693 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
200.0 0. 200.0 0. !RE has been fitted to be =RS; TE=0 because the map opening = the bending angle = 90 deg; RS=RE; TS=-TE
'FIT' 3
......@@ -26,7 +26,7 @@ HEADER_2 ! Title
1 !number of constraints
3 1 2 2 0.0 1. 0 !3:particle coordinate; 1: particle #1; 2: Y(i.e. R); where?:at the end of element #2; wanted value=0.0; weight; 0:no additional parameters
'END' 6 4
'END' 4
************************************************************************************************************************************
1 Keyword, label(s) : OBJET
......@@ -43,10 +43,6 @@ HEADER_2 ! Title
2 Keyword, label(s) : TOSCA
OPEN FILE zgoubi.plt
FOR PRINTING TRAJECTORIES
NDIM = 2 ; Number of data file sets used is 1 ; Stored in field array # IMAP = 1 ;
Value of MOD.MOD2 is 22.1
......@@ -69,14 +65,14 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
Field map limits, angle, min, max, max-min (rad) : -0.785398E+00 0.785398E+00 0.157080E+01
Field map limits, radius, min, max, max-min (cm) : 0.165853E+03 0.234853E+03 0.690000E+02
Min / max fields drawn from map data : -1.90 / 4.24
@ x-node, y-node, z-node : 0.314 232. 0.00 / 0.777 185. 0.00
Normalisation coeff. BNORM : 1.036
Field min/max normalised : -1.97 / 4.39
Nbre of nodes in x/y/z : 181/ 47/ 1
Node distance in x/y/z : 8.7266E-03/ 1.500 / 0.000
Field map limits, angle : min, max, max-min (rad) : -0.78539816E+00 0.78539816E+00 0.15707963E+01
Field map limits, radius : min, max, max-min (cm) : 0.16585281E+03 0.23485281E+03 0.69000000E+02
Min / max fields drawn from map data : -1.8960 / 4.2372
@ X-node, Y-node, z-node : 0.314 232. 0.00 / 0.777 185. 0.00
Normalisation coeff. BNORM : 1.03650
Field min/max normalised : -1.9652 / 4.3919
Nbre of nodes in X/Y/Z : 181/ 47/ 1
Node distance in X/Y/Z : 8.726646E-03/ 1.50000 / 0.00000
Option for interpolation : 4
Smoothing using 25 point
......@@ -84,20 +80,10 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
Integration step : 0.1000 cm (i.e., 4.9912E-04 rad at mean radius RM = 200.4 )
KPOS = 2. Position of reference orbit on mechanical faces
at entrance RE = 180.140 cm TE = 0.00000 rad
at exit RS = 180.140 cm TS = 0.00000 rad
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 180.140 -0.000 0.000 0.000 1
CONDITIONS DE MAXWELL ( 3086. PAS ) :
DIV(B) LAPLACIEN(B) ROTATIONNEL(B)
0.000 0.000 0.000
0.000 0.000
0.000 0.000
LAPLACIEN SCALAIRE = 0.000
at entrance RE = 200.000 cm TE = 0.00000 rad
at exit RS = 200.000 cm TS = 0.00000 rad
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 182.857 -0.118 0.000 0.000 1
Cumulative length of optical axis = 0.00000000 m ; Time (for ref. rigidity & particle) = 0.00000 s
......@@ -114,18 +100,18 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
constraint # 1 IR = 2 , ok.
constraint # 1 I = 1 , ok.
FIT variables in good order, FIT will proceed.
FIT variables and constraints in good order, FIT will proceed.
Final FIT status will NOT be saved. For so, use the 'save [FileName]' command
STATUS OF VARIABLES (Iteration # 0 / 90 max.)
LMNT VAR PARAM MINIMUM INITIAL FINAL MAXIMUM STEP NAME LBL1 LBL2
2 1 70 90.1 180. 180.1396693 270. 0.600 TOSCA - -
2 1 72 90.1 180. 180.1396693 270. 0.600
STATUS OF VARIABLES (Iteration # 0 / 999 max.)
LMNT VAR PARAM MINIMUM INITIAL FINAL MAXIMUM STEP NAME LBL1 LBL2
2 1 70 100. 180. 180.13968 300. 1.129E-05 TOSCA - -
2 1 72 100. 180. 180.13968 300. 1.129E-05
STATUS OF CONSTRAINTS (Target penalty = 1.0000E-10)
TYPE I J LMNT# DESIRED WEIGHT REACHED KI2 NAME LBL1 LBL2 Parameter(s)
3 1 2 2 0.0000000E+00 1.0000E+00 8.8493140E-06 1.0000E+00 TOSCA - - 0 :
Fit reached penalty value 7.8310E-11
TYPE I J LMNT# DESIRED WEIGHT REACHED KI2 NAME LBL1 LBL2 Nb param. [value]
3 1 2 2 0.000000E+00 1.000E+00 -9.469218E-07 1.00E+00 TOSCA - - 0
Fit reached penalty value 8.9666E-13
************************************************************************************************************************************
......@@ -146,9 +132,6 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
2 Keyword, label(s) : TOSCA
zgoubi.plt
already open...
NDIM = 2 ; Number of data file sets used is 1 ; Stored in field array # IMAP = 1 ;
Value of MOD.MOD2 is 22.1
......@@ -159,14 +142,14 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
Pgm toscap, restored mesh coordinates for field map # 1, name : ../TOSCAMOD22.dat
Field map limits, angle, min, max, max-min (rad) : -0.785398E+00 0.785398E+00 0.157080E+01
Field map limits, radius, min, max, max-min (cm) : 0.165853E+03 0.234853E+03 0.690000E+02
Min / max fields drawn from map data : -1.90 / 4.24
@ x-node, y-node, z-node : 0.314 232. 0.00 / 0.777 185. 0.00
Normalisation coeff. BNORM : 1.036
Field min/max normalised : -1.97 / 4.39
Nbre of nodes in x/y/z : 181/ 47/ 1
Node distance in x/y/z : 8.7266E-03/ 1.500 / 0.000
Field map limits, angle : min, max, max-min (rad) : -0.78539816E+00 0.78539816E+00 0.15707963E+01
Field map limits, radius : min, max, max-min (cm) : 0.16585281E+03 0.23485281E+03 0.69000000E+02
Min / max fields drawn from map data : -1.8960 / 4.2372
@ X-node, Y-node, z-node : 0.314 232. 0.00 / 0.777 185. 0.00
Normalisation coeff. BNORM : 1.03650
Field min/max normalised : -1.9652 / 4.3919
Nbre of nodes in X/Y/Z : 181/ 47/ 1
Node distance in X/Y/Z : 8.726646E-03/ 1.50000 / 0.00000
Option for interpolation : 4
Smoothing using 25 point
......@@ -179,16 +162,6 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
A 1 1.0000 0.000 0.000 0.000 0.000 0.785 180.140 -0.000 0.000 0.000 1
CONDITIONS DE MAXWELL ( 3086. PAS ) :
DIV(B) LAPLACIEN(B) ROTATIONNEL(B)
0.000 0.000 0.000
0.000 0.000
0.000 0.000
LAPLACIEN SCALAIRE = 0.000
Cumulative length of optical axis = 0.00000000 m ; Time (for ref. rigidity & particle) = 0.00000 s
************************************************************************************************************************************
......@@ -206,7 +179,15 @@ R_min (cm), DR (cm), DTTA (deg), DZ (cm) : 1.658528E+02 1.500000E+00 5.00000
Made it to the end : 1
************************************************************************************************************************************
ZGOUBI RUN COMPLETED.
Zgoubi, author's dvlpmnt version.
Job started on 08-08-2017, at 16:06:17
JOB ENDED ON 08-08-2017, AT 16:06:18
CPU time, total : 0.23859899999999998
An updated version of the input data file, with variables in FIT'ed state, has been saved in zgoubi.FIT.out.dat.
Pgm zgoubi : Execution ended upon key END
************************************************************************************************************************************
%!PS-Adobe-2.0 EPSF-2.0
%%Title: layout.eps
%%Creator: gnuplot 4.6 patchlevel 4
%%CreationDate: Wed Aug 9 10:01:24 2017
%%CreationDate: Wed Aug 9 10:11:14 2017
%%DocumentFonts: (atend)
%%BoundingBox: 50 50 410 302
%%EndComments
......@@ -461,7 +461,7 @@ SDict begin [
/Author (tplanche)
% /Producer (gnuplot)
% /Keywords ()
/CreationDate (Wed Aug 9 10:01:24 2017)
/CreationDate (Wed Aug 9 10:11:14 2017)
/DOCINFO pdfmark
end
} ifelse
......
\documentclass[usletter, 11pt]{article}
\usepackage{graphicx}
\usepackage{amsmath,amssymb}
\DeclareMathAlphabet\mathbfcal{OMS}{cmsy}{b}{n}% for \mathbfcal to work
\usepackage{amsfonts}
\usepackage{color}
\definecolor{purple}{rgb}{1,0,1}
\definecolor{turquoise}{rgb}{0,1,1}
\usepackage[center, small, sl]{caption}
\usepackage{hyperref}%to include hyper links
\hypersetup{colorlinks,linkcolor=blue,citecolor=blue}%links
\definecolor{listinggray}{gray}{0.9}
\usepackage{listings}%to include lines of code
\usepackage[section]{placeins}% to be able to use \FloatBarrier
\usepackage{cleveref}%to use \cref
\usepackage{longtable}
\addtolength{\textwidth}{1.5 cm}
\addtolength{\hoffset}{-0.8 cm}
%\addtolength{\topmargin}{-2 cm}
%\addtolength{\textheight}{3. cm}
%opening
\title{Zgoubi, would you do for me\\TRIUMF's High Resolution Spectrometer?}
\author{Thomas Planche}
%\footnote{tplanche@triumf.ca}\\~\\ TRIUMF\thanks{This work has been supported by the Natural Sciences and Engineering Research Council of Canada. TRIUMF also receives federal funding via a contribution agreement through the National Research Council of Canada.}}
\begin{document}
\maketitle
\tableofcontents
% \section{Introduction}
% \begin{figure}[htb]
% \centering
% \includegraphics*[width=80mm]{figures/TUO1A03_Fig3}
% \caption{}\label{fig:fold}
% \end{figure}
\FloatBarrier
% \begin{equation}
% f(x,y,z)= \left\{
% \begin{array}{lc}
% \dfrac{Q}{\epsilon_0 4\pi r} & r \ge R~,\\
% ~\\
% \dfrac{Q}{\epsilon_0 8 \pi R} \left(3-\dfrac{r^2}{R^2}\right) & r < R~,\\
% \end{array}
% \right.
% \end{equation}
% \begin{equation}
% \begin{aligned}
% \delta E_k =& q V_{\rm rf} \sin{\phi}\\
% \delta z' =& K \, \frac{1 - a}{2} z\\
% \delta r' =& K \, \frac{1 + a}{2} (r-r_0)\\
% \end{aligned}
% \end{equation}
% \begin{figure}[!ht]
% \begin{center}
% \includegraphics[width=14cm]{figures/turn300}
% \caption{Turn \#300 shown together with neighbouring turns. The Blue lines is an arc of circle centered on the machine center.}\label{fig:turn30_50k_500uA}
% \end{center}
% \end{figure}
% \begin{figure}[!ht]
% \begin{minipage}{.5\textwidth}
% \includegraphics[width=7cm]{figures/0b9.png}
% \end{minipage}
% \begin{minipage}{.5\textwidth}
% \includegraphics[width=7cm]{figures/410b9.png}
% \end{minipage}
% \caption{}\label{fig:spunch}
% \end{figure}
\bibliographystyle{/home/tplanche/utils/latex/formating/elsarticle-num}
\bibliography{/home/tplanche/utils/latex/Mybib,/home/tplanche/utils/latex/AllDN}
\end{document}
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment