Skip to content
Snippets Groups Projects
Commit 02520bd4 authored by Mark Jones's avatar Mark Jones Committed by Mark K Jones
Browse files

Update SHMS forward and recon matrix elements 

1) Create a new fox file for COSY based on worked
done by Aruni Nadeeshani which fine tuned the Q1,Q2 and Q3
settings to match focal plane xfp/yfp distributions.

2) Run cosy with new fox file to create new forward matrix file:
shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_forward.dat
and link file to shms_forward.dat

3) Create a new reconstruction matrix file:
shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_recon_fit.dat

which was created by fitting SHMS single arm MC data
using the new forward matrix element. It was found that
the reconstruction matrix produced by COSY had very bad
reconstruction of xptar in regions of xptar/delta phase space.
parent f90e33d2
No related branches found
No related tags found
No related merge requests found
Show whitespace changes
Inline Side-by-side
cosy/shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018.fox 0 → 100755
+ 328
− 0
View file @ 02520bd4
{------------------------------------------------------------------------------
Super-HMS spectrometer, with forward and reconstruction maps.
VERSION for Monte Carlo maps
2017 version updated the effective lenghts
------------------------------------------------------------------------------
}
INCLUDE 'JLAB' ;
PROCEDURE RUN ;
{=========================== Variable Declarations ===========================}
VARIABLE Z1C 1 ; VARIABLE L1 1 ; VARIABLE L1E 1 ; VARIABLE Q1MID 1 ;
VARIABLE Z2C 1 ; VARIABLE L2 1 ; VARIABLE L2E 1 ; VARIABLE Q2MID 1 ;
VARIABLE Z3C 1 ; VARIABLE L3 1 ; VARIABLE L3E 1 ; VARIABLE Q3MID 1;
VARIABLE ZDC 1 ; VARIABLE LD 1 ; VARIABLE LDE 1 ;
VARIABLE ZBC 1 ; VARIABLE LB 1; VARIABLE LBE 1;
VARIABLE DZD1 1 ; VARIABLE DZD2 1;
VARIABLE DZF1 1 ; VARIABLE DZF2 1;
VARIABLE XDIP 1;
VARIABLE DTILT 1;
VARIABLE FTILT1 1; VARIABLE FTILT2 1;
VARIABLE ZFM 1 ;
VARIABLE DR0 1;
VARIABLE DR1 1 ;
VARIABLE DR2 1 ;
VARIABLE DR3 1 ;
VARIABLE DR4 1 ;
VARIABLE DRF 1 ;
VARIABLE DRFOC 1;
VARIABLE DR0EN 1; VARIABLE DR0EX 1;
VARIABLE DR1EN 1; VARIABLE DR1EX 1;
VARIABLE DR2EN 1; VARIABLE DR2EX 1;
VARIABLE DR3EN 1; VARIABLE DR3EX 1;
VARIABLE DLEX 1;
VARIABLE Q1 1 ; { quadrupole strength }
VARIABLE Q2 1 ; { quadrupole strength }
VARIABLE Q3 1 ; { quadrupole strength }
VARIABLE Q3SCALE 1 ; { flag for scaling Q3 }
VARIABLE FILE 80 ; { filename}
VARIABLE QFILE 80 ; { filename}
VARIABLE TILT 1 ; { Focal plane tilt angle }
VARIABLE I 1 ;
VARIABLE TMP 1 ;
VARIABLE SNAME 4 ; {Spectrometer name}
VARIABLE FRMODE 1 ; {Fringe field mode }
VARIABLE MORDER 1 ; { Maximum order for transformations}
VARIABLE FORDER 1 ; {Order of forward maps }
VARIABLE RORDER 1 ; {Order of reconstruction map }
VARIABLE RES 1 4 ; {Resolution array }
VARIABLE TMPSPOS 1000 ; { Saved value of SPOS }
VARIABLE GMAP 2000 8 ; { Global transfer map }
VARIABLE RMAP 2000 5 ; { Reconstruction map }
{========================== Procedure declarations ===========================}
{ Output forward transformation data at an aperture and update internal maps }
PROCEDURE APERTURE LABEL REGION ;
WRITE 20 '!NAME: '&LABEL ;
WRITE 20 '!REGION: '&REGION ;
WRITE 20 '!OFFSET: '&SF(TMPSPOS,'(F25.16)')&' (in meters)' ;
WRITE 20 '!LENGTH: '&SF(SPOS,'(F25.16)')&' (canonical length in meters)' ;
PT 20 ; {Output data}
TMPSPOS := TMPSPOS + SPOS ;
ANM MAP GMAP GMAP ; UM ; {Update global map + reset}
WRITE 6 ' Aperture at S = '&S(TMPSPOS)&' '&LABEL ;
ENDPROCEDURE ;
{ Compute and output maps to inner aperture and exit of quad magnet. }
PROCEDURE MAGQUAD NAME LEFF Q A FRAC ;
VARIABLE DRIFT 1 ; VARIABLE LMID 1 ;
{ Compute the location of the inner aperture }
LMID := (LEFF*FRAC) ;
IF FRMODE#0 ; FR -1 ; MQ LEFF Q A ; ENDIF ; {Entrance fringe field}
FR 0 ; MQ LMID Q A ; {First part of quad}
APERTURE NAME&'_MID' 'FRONT' ;
MQ LEFF-LMID Q A ; {Second part of quad}
IF FRMODE#0 ; FR -2 ; MQ LEFF Q A ; ENDIF ; {Exit fringe field}
APERTURE NAME&'_MAGNETIC_EXIT' 'BACK' ;
FR FRMODE ; ENDPROCEDURE ; {Restore fringe field mode}
{ Q3 fieldmap }
PROCEDURE QUAD ;
VARIABLE IS 1 ;
LOOP IS 1 4 ; FIELDMAP IS 'SHMS_Q3' ; ENDLOOP ;
ENDPROCEDURE ;
{Dipole fieldmap}
PROCEDURE DIPOLE ;
VARIABLE IS 1 ;
LOOP IS 1 8 ; FIELDMAP IS 'SHMS_DIP' ; ENDLOOP ;
ENDPROCEDURE ;
{================================= Main Code =================================}
{ Beam specifications }
MORDER := 6 ; {Max order for calculations}
OV morder 3 0 ; {Order Value <order><phase dim><#par>}
RPM 11000.0 5.4858E-4 -1 ; {<momentum> <mass in amu> <charge>}
SB 0 .040 0 .050 .015 0 0 .15 0 0 0 ;
{ --> PX PA R12 PY PB R34 PT PD R56 PG PZ }
{ Spectrometer name }
SNAME := 'SHMS' ;
{ Fringe-field specifications.}
frmode := 3 ;
FR FRMODE ;
{ Magnet geometry }
ZBC := 1.76 ; {Center of bender}
Z1C := 4.140 ; {Center of Q1}
Z2C := 6.80 ; {Center of Q2}
Z3C := 9.45 ; {Center of Q3}
ZDC := 12.937; {Center of dipole - might not use for anything}
XDIP:= 0.26 ; {Dipole vertical offset}
DTILT:= 9.2*pi/180.0; {Dipole tilt}
DZF1:= sqrt(4.933^2+(1.632-XDIP)^2);
FTILT1:= atan((1.632-XDIP)/4.933);
FTILT2:= FTILT1-DTILT ;
DZF2:= DZF1*cos(FTILT2);
LB := 1.1522 ;
L1 := 2.44 ; {Physical length of Q1}
L2 := 2.369 ; {Physical length of Q2}
L3 := 2.369 ; {Physical length of Q3}
LD := 4.098 ; {Physical length of dipole}
LBE := 0.80 ; {Effective length of bender}
L1E := 1.873 ; {Effective length of Q1}
L2E := 1.587 ; {Effective length of Q2}
L3E := 1.587 ; {Effective length of Q3}
LDE := 2.90 ;
DRF:= (DZF2-LDE/2)/cos(DTILT);
{ Compute some drift distances }
DR0 := ZBC-(LB/2); {Drift from target to bender}
DR1 := (Z1C - ZBC) - ((LB + L1)/2) ; {Drift from Bender to Q1}
DR2 := (Z2C - Z1C) - ((L1 + L2)/2) ; {Drift from Q1 to Q2}
DR3 := (Z3C - Z2C) - ((L2 + L3)/2) ; {Drift from Q2 to Q3}
{ DR4 := (ZDC - Z3C) - ((L3 + LD)/2) ; {Drift from Q3 to dipole}}
DZD1:= XDIP/tan(DTILT);
DZD2:= (LDE/2-XDIP/sin(DTILT))/cos(DTILT);
DR4 := (ZDC-Z3C) - L3E/2 - DZD1 - DZD2 ; {Drift from Q3 to dipole}
DR0EN := (LB-LBE)/2; {Drift from bender mech to mag ent}
DR0EX := (LB-LBE)/2; {Drift from bender mag to mech ex}
DR1EN := (L1-L1E)/2; {Drift from Q1 mech to mag ent}
DR1EX := (L1-L1E)/2; {Drift from Q1 mag to mech ex}
DR2EN := (L2-L2E)/2; {Drift from Q2 mech to mag ent}
DR2EX := (L2-L2E)/2; {Drift from Q2 mag to mech ex}
DR3EN := (L3-L3E)/2; {Drift from Q3 mech to mag ent}
DR3EX := (L3-L3E)/2; {Drift from Q3 mag to mech ex}
DLEX := 0.599/cos(DTILT) ; {Drift from dipole mag to mech exit}
DRFOC := DRF-DLEX ;
{ QUAD STRENGTHS }
{COSY solution}
{Let Q1/Q2/Q3 all vary: ME(5,2) = -1.2}
Q1 := -2.136856 ;
Q2 := 4.5670 ;
Q3 := -2.954160 ;
{Let Q1/Q2/Q3 all vary: ME(5,2) = -0.8}
{ Q1 := -2.370969 ;}
{ Q2 := 4.330444 ;}
{ Q3 := -2.530578 ;}
{ Computation order specifications }
ASK forder 'Order for maps (1-'&SF(morder,'(i1)')&'): ' 1 morder ;
CO FORDER ; RORDER := FORDER ;
{ Begin composing forward maps }
qfile := SNAME&'_Q3S.dat' ; {Quad scaling symplectic map}
file := SNAME&'_fr'&sf(frmode,'(i1)')&'_for_maps.dat' ;
openf 20 file 'UNKNOWN' ;
FR FRMODE ;
UM ; SM GMAP ; TMPSPOS := SPOS ; {Initialize global map}
{ Insert splitter magnet }
DL DR0 ;
APERTURE 'BENDER_MECHANICAL_ENTRANCE' 'FULL';
DL DR0EN ;
APERTURE 'BENDER_MAGNETIC_ENTRANCE' 'FULL';
RA -90 ; DP 14.34959928 3.0 0.225 ; RA 90 ; {Bender}
APERTURE 'BENDER_MAGNETIC_EXIT' 'FULL' ;
DL DR0EX ;
APERTURE 'BENDER_MECHANICAL_EXIT' 'FULL' ;
{ QUAD 1 }
DL DR1 ; {Drift TO Q1}
APERTURE 'Q1_MECHANICAL_ENTRANCE' 'FULL' ;
DL DR1EN ;
APERTURE 'Q1_MAGNETIC_ENTRANCE' 'FULL' ;
MAGQUAD 'Q1' L1E Q1 0.25 1/2 ;
DL DR1EX ;
APERTURE 'Q1_MECHANICAL_EXIT' 'FULL' ;
{ QUAD 2 }
DL DR2 ; {Drift TO Q2}
APERTURE 'Q2_MECHANICAL_ENTRANCE' 'FULL' ;
DL DR2EN ;
APERTURE 'Q2_MAGNETIC_ENTRANCE' 'FULL' ;
MAGQUAD 'Q2' L2E Q2 0.35 1/2 ;
DL DR2EX ;
APERTURE 'Q2_MECHANICAL_EXIT' 'FULL' ;
{ QUAD 3 }
DL DR3 ; {Drift TO Q3}
APERTURE 'Q3_MECHANICAL_ENTRANCE' 'FULL' ;
DL DR3EN ;
APERTURE 'Q3_MAGNETIC_ENTRANCE' 'FULL' ;
MAGQUAD 'Q3' L3E Q3 0.35 1/2 ;
DL DR3EX ;
APERTURE 'Q3_MECHANICAL_EXIT' 'FULL' ;
DL 0.18 ; {dipole mechanical entrance}
APERTURE 'DIPOLE_MECHANICAL_ENTRANCE' 'D-FRONT' ;
DL 0.301 ; {dipole mechanical entrance}
APERTURE 'DIPOLE_FLARE' 'D-FRONT' ;
DL 0.394715 ; {drift to magnetic entrance}
APERTURE 'DIPOLE_MAGNETIC_ENTRANCE' 'D-FRONT' ;
{ Dipole }
{Break into 8 segments of 2.3 degrees = 18.4 deg total}
FR -1; DP 9.06924 18.4 0.35 ; {calculate entrance fringe field}
FR 0 ; DP 9.06924 2.3 0.35 ; {now use for first trans}
APERTURE 'DIPOLE_MID' 'dipole';
FR 0; {turn all fringe fields off}
LOOP I 1 6;
DP 9.06924 2.3 0.35 ;
APERTURE 'DIPOLE_MID' 'dipole' ;
ENDLOOP ;
FR -2; DP 9.06924 18.4 0.35 ; {calculate exit fringe field}
FR 0; DP 9.06924 2.3 0.35 ;
APERTURE 'DIPOLE_MAGNETIC_EXIT' 'dipole' ;
FR FRMODE ;
DL DLEX ;
APERTURE 'DIPOLE_MECHANICAL_EXIT' 'drift' ; {dipole mechanical exit}
{ Focal Plane }
DL DRFOC ; {Drift to detector focal plane}
APERTURE 'FOCAL_PLANE' 'drift' ;
closef 20 ; {Done outputting forward maps}
{ Output global map in transport notation }
file := SNAME&'_fr'&sf(frmode,'(i1)')&'_for_global.dat' ;
openf 20 file 'UNKNOWN' ; UM ; AM GMAP ; PT 20 ; closef 20 ; UM ;
{ Some diagnostics }
WRITE 6 ' Q1-3, D1-F, Z1-M, LENGTH = ' (Q1)&(Q2)&(Q3) ;
WRITE 6 DR1&DR2&DR3&DR4&DRF ;
WRITE 6 Z1C&Z2C&Z3C&ZFM&SPOS ;
{ Convert GMAP from canonical to TRANSPORT form}
CATR GMAP MAP TMPSPOS ; spos := tmpspos ;
WRITE 6 ' Total length = '&S(CONS(SPOS))&' meters' ;
{ Output diagnostic matrix elements }
WRITE 6 ' ME(5,2) = '&S(ME(5,2)) ;
WRITE 6 ' ME(1,2) = '&S(ME(1,2)) ;
WRITE 6 ' ME(3,4) = '&S(ME(3,4)) ;
{ Compute focal plane tilt angle }
TILT := (180*atan(ME(1,1)*ME(1,26)/ME(1,6))/pi)+90 ;
write 6 ' Focal plane tilt ='&SF(tilt,'(f8.3)')&' degrees' ;
{ Compute RECON map }
RR1 MAP 0 0 .04 0 .04 .15 0 0 10000 RORDER RMAP RES ;
write 6 ' RECON resolutions = ' RES(1)&RES(2)&RES(3)&RES(4) ;
file := SNAME&'_fr'&sf(frmode,'(i1)')&'_rec.dat' ;
openf 20 file 'UNKNOWN' ; pmr rmap 20 ; closef 20 ;
ENDPROCEDURE ;
RUN ;
END ;
shms/shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_forward.dat 0 → 100644
+ 2798
− 0
View file @ 02520bd4
Source diff could not be displayed: it is too large. Options to address this: view the blob.
shms/shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_recon_fit.dat 0 → 100644
+ 253
− 0
View file @ 02520bd4
0.000157023759 -2.11258058e-05 2.04962709e-05 0.00041841748 00000
0.289198426 -0.0588920715 0.0329057746 0.624183647 10000
-1.49489069 0.011666898 -0.00512945882 -0.110042128 01000
-0.000595996697 -0.793476279 0.302743193 -0.000941021081 00100
-0.00145396264 1.61701836 -1.90481581 -0.00243013071 00010
0.3755404 0.122977509 0.0101196802 0.957059556 00001
-0.614481569 0.189556215 -0.102259765 -1.40055036 20000
3.40887041 -0.989043689 0.571432623 9.84852943 11000
-1.36980683 0.191156861 -0.261473671 -3.98501635 02000
0.015260013 4.14917758 -0.807983529 0.0215125593 10100
-0.0886873031 -2.93450285 0.613842529 -0.134246215 01100
-0.0723254746 0.034373833 -0.0148245201 -0.115391982 00200
-0.0699263244 -41.205406 7.83233525 -0.11298568 10010
0.249382753 20.6545304 -3.71633084 0.379380533 01010
0.668249482 -0.16315596 0.109913592 1.21803519 00110
-2.71378223 6.27894223 -1.47767539 -4.46884607 10001
6.65207426 -27.5931841 6.12912236 16.4771463 01001
0.246651983 4.95645263 -0.948531566 0.435528043 00101
-2.17901123 0.133147919 -0.248363525 -4.32403867 00020
-1.84820477 -57.6992061 10.9461338 -1.02785953 00011
-26.9127344 32.9722604 -8.87451743 -45.6307228 00002
5.27642668 -0.377828084 0.554621911 9.8586213 30000
-49.4146844 3.49373613 -5.29526282 -100.697611 21000
137.824707 -9.18516534 14.8994468 312.118977 12000
-114.04636 9.95706854 -12.3848982 -244.310636 03000
0.213582279 -10.489066 2.06031194 0.298108433 20100
-0.836832155 46.7643012 -8.65906229 -0.942032639 11100
1.64257592 6.20967313 -2.47672234 1.49060224 02100
3.15206851 -0.26032761 0.319003031 5.4683025 10200
-11.0588218 2.91606685 -1.47125442 -20.5599679 01200
-0.231471953 -1.35749093 0.223211763 -0.350555723 00300
0.250347902 160.796395 -31.4037597 0.561434603 20010
-2.65943719 -914.921674 173.994775 -6.82953474 11010
-2.14888634 913.851705 -170.896063 2.1990872 02010
-37.2601415 6.15542771 -4.61708235 -74.0305521 10110
137.384762 -40.2930456 19.8454183 295.889831 01110
2.05677332 63.1019531 -11.6180822 2.81343548 00210
13.1186877 -86.6458723 18.1309676 20.5254858 20001
-88.6139664 588.037602 -122.710403 -158.900748 11001
122.922626 -888.003421 182.656084 281.750849 02001
-22.4235377 5.84400617 -3.09158607 -40.6703478 10101
114.567201 -66.4444283 23.6079085 196.424555 01101
11.1250277 -10.1463012 3.10074817 14.9287864 00201
132.078977 -19.5584392 16.5516681 298.477322 10020
-572.753002 148.223048 -81.4931495 -1370.40404 01020
-8.07457474 -761.131906 143.745275 -10.6321416 00120
-22.1272267 278.572419 -56.8901202 -13.6869497 10011
-170.799566 -807.564424 140.866159 -240.162628 01011
-69.3852595 76.2469848 -22.4686587 -95.6114342 00111
87.7710533 25.3995902 1.80038297 104.984648 10002
135.946832 74.7980762 -4.43412528 202.247868 01002
-94.2717801 521.995577 -111.848853 -156.68086 00102
20.2835063 2859.38301 -545.648091 29.7849224 00030
92.1970079 -1057.31757 215.840605 155.775516 00021
59.8675201 254.12426 -44.6444676 114.043952 00012
5.50861319 -42.7894502 8.85190492 8.9512421 00003
-30.9616255 -3.84655634 -1.91976343 -51.4055197 40000
371.047672 42.6774928 24.123319 636.563743 31000
-1461.54837 -154.625136 -99.0510636 -2580.44147 22000
1914.12438 192.126521 129.882072 3291.87506 13000
-73.9734998 -40.1868672 13.2458276 673.737179 04000
-4.03094022 -41.8485896 7.82448031 -6.28511032 30100
29.5423343 550.667339 -105.821595 46.7758515 21100
-64.8454851 -2363.92979 458.259165 -104.287958 12100
40.0409674 2952.58751 -562.633002 68.573733 03100
-27.3631913 1.97189975 -2.59678171 -42.3007327 20200
150.24486 -23.8847051 17.1519095 249.03585 11200
-131.366426 -19.1362323 -8.57893436 -273.114466 02200
2.39139049 -13.6984598 2.99152019 3.65633225 10300
-5.83168341 24.8745947 -5.79774987 -8.48522606 01300
-2.93237517 -0.0837280145 -0.212468335 -3.79972426 00400
2.44757157 -292.290306 58.5455436 5.04577987 30010
52.9369293 2146.99692 -414.014652 71.2142442 21010
-357.105949 -2067.61461 346.984958 -519.90538 12010
611.471227 -3947.86834 781.739969 849.404813 03010
236.435768 -19.4588502 22.7385868 353.054003 20110
-1061.9135 323.748754 -148.137931 -1672.4002 11110
285.276263 -78.5060826 43.3269672 767.324311 02110
-28.7236718 106.981075 -24.6123993 -43.0649608 10210
53.2629597 -109.358149 31.3182634 77.1651724 01210
50.7646622 5.06590698 3.11372035 72.7942289 00310
-15.8348585 108.347464 -22.3914958 -18.3757188 30001
183.547886 -1099.8034 230.1153 278.446052 21001
-354.833943 2292.41128 -476.610172 -646.784391 12001
-251.408266 1627.91933 -339.386258 -457.862399 03001
22.2505342 -75.2116033 16.7765575 40.0749352 20101
-349.165204 382.919861 -105.024196 -549.239078 11101
361.193371 95.2309395 10.283055 649.340971 02101
-84.265558 134.998123 -33.8943749 -129.750947 10201
243.115564 -420.575844 103.514303 425.321245 01201
26.6531667 -7.45882716 3.78100514 39.6073507 00301
-604.254678 -30.83159 -44.1796765 -860.767029 20020
2337.52536 -146.800697 214.732775 3303.27118 11020
1407.89605 -204.135877 153.793545 2022.5552 02020
84.6390156 -443.965079 94.3367711 107.249756 10120
-128.611502 2160.35449 -439.002347 -213.521415 01120
-285.4419 -126.448527 0.0587564477 -494.97619 00220
250.813526 -126.21976 46.1512958 321.288908 20011
146.068431 -50.8256543 21.7426274 231.876936 11011
150.39686 -23.0155735 16.6654918 236.013638 02011
90.3779362 -643.419325 132.036018 50.4406637 10111
272.98514 26.5942137 16.6015942 412.365853 01111
-339.232904 -131.437872 -3.66482886 -545.957505 00211
54.3166997 -81.5966119 25.2476867 245.364548 20002
57.0395384 8.44425268 4.36290146 127.658232 11002
23.6304845 10.5016877 0.230343479 46.4415885 02002
-187.646049 -239.273977 31.3657588 -379.729313 10102
-47.9340015 -33.09271 2.65180286 -96.703276 01102
-38.318044 45.2804743 -10.5616202 7.29192368 00202
21.0355574 -1202.41739 246.274058 113.069634 10030
-127.074791 531.711394 -111.9034 -168.071246 01030
487.016213 535.454642 -58.5889368 1078.73824 00130
86.210337 -190.113647 44.0731232 119.670706 10021
60.8280963 -32.4235011 11.4366001 94.778876 01021
-80.8990211 -23.9120045 -2.44913469 -133.731202 00121
-12.6700537 -7.94442544 0.602130242 -27.1994929 10012
-1.73198269 -0.949051214 0.0552687794 -4.56273394 01012
-6.64231164 5.03503777 -1.39701642 -3.38533848 00112
-1.66187478 9.0832072 -1.89489734 -2.64760277 10003
-0.399730115 1.98193231 -0.41677759 -0.667011242 01003
1.6887521 0.723361092 -0.0027023929 2.90429852 00103
170.064887 242.929805 -31.1599599 392.198155 00040
-23.1716949 -4.49503506 -1.11782808 -37.0490619 00031
-0.425420535 0.245401268 -0.0680297587 0.0667210074 00022
0.352229619 0.117375474 0.00605258235 0.566369764 00013
0.00847464136 -0.000763546654 0.000849084623 0.013736421 00004
52.7162774 10.0938835 2.53799861 83.8162156 50000
-656.80092 -121.240526 -33.4746575 -1076.0814 41000
2459.24459 459.582459 129.571484 4197.78364 32000
-1840.62693 -445.055341 -87.1881093 -3455.88591 23000
-2684.65046 -377.474582 -172.202456 -4764.16615 14000
-1792.83886 -164.209501 -129.785486 -3095.23136 05000
10.2199141 144.441766 -27.2830173 16.87881 40100
-79.3145968 -1631.9658 311.711282 -137.924867 31100
160.21563 5525.52502 -1062.42814 306.553582 22100
-58.1982227 -3885.92544 736.928877 -135.262348 13100
-49.594135 -4830.09834 923.604553 -136.403727 04100
61.5688836 -7.5630023 6.49004807 93.9001377 30200
-414.777581 113.498485 -57.457561 -670.293041 21200
608.167346 -386.424939 132.229561 1092.76907 12200
-325.411197 891.938551 -203.664561 -470.525664 03200
-4.78744392 70.4288473 -14.4118354 -6.82713289 20300
14.9990476 -285.143247 58.0413393 17.6954881 11300
6.49036545 296.60285 -58.2670736 16.9852755 02300
15.4181152 -3.57630364 1.94425019 21.8411906 10400
-20.7949732 21.2603895 -6.02706518 -36.2964973 01400
0.101297017 5.25506705 -1.04962535 0.187268694 00500
-38.1991914 -135.249601 20.4425477 -54.324444 40010
181.575653 1133.81361 -189.523145 261.556554 31010
-13.5311549 -1786.37419 345.950198 -2.84942603 22010
15.2126742 -2200.75349 424.081176 21.6202644 13010
36.2776051 -1414.77191 274.108199 46.2996312 04010
-347.367368 -32.1097691 -19.8481201 -448.757328 30110
1439.91268 -566.828065 220.657493 1875.1512 21110
648.468225 22.7367248 45.4752471 898.349909 12110
151.805402 246.084635 -36.6731717 212.166983 03110
38.0904707 -793.244299 158.399113 57.3571543 20210
-104.069618 1264.32628 -253.595682 -137.277028 11210
114.291581 1804.7126 -341.778286 179.841918 02210
-293.193853 -20.5311446 -19.7884294 -435.159858 10310
401.600365 -7.88699764 37.3382854 734.339838 01310
2.99360148 -112.543668 22.8882914 5.77981322 00410
-34.164077 222.287894 -47.4541035 -57.0716554 40001
195.574151 -1064.59179 229.753039 318.297302 31001
-4.35760809 -34.566028 8.92239598 -22.3662993 22001
-30.3889434 129.515038 -27.1289965 -60.8684666 13001
-16.932109 81.1418768 -17.1515204 -33.2588595 04001
283.140813 134.603033 -3.1505734 471.579747 30101
-979.191653 -588.71229 34.1929749 -1666.3796 21101
-366.382695 -332.633311 34.4254703 -623.891185 12101
-90.155689 -132.347089 18.2111008 -152.325768 03101
158.726292 164.77323 -20.5656842 243.340054 20201
-644.240009 -891.695328 124.115636 -993.501223 11201
-308.065419 -357.580505 45.9038732 -479.145493 02201
-20.6379986 -44.7150178 6.57453153 -32.3169874 10301
34.3436519 36.7996083 -5.32654989 2.48657831 01301
-44.1240782 66.9703256 -17.2634717 -73.2051022 00401
157.567058 719.062373 -128.876905 230.943608 30020
701.087502 81.8758263 38.794821 979.553249 21020
274.121913 59.7653386 9.76283489 376.506311 12020
79.3327747 49.1533899 -3.4194422 104.603002 03020
-96.5407548 1669.28153 -335.553359 -178.636344 20120
-37.173743 1231.89881 -243.899884 -59.8756082 11120
18.9193242 677.87014 -130.571116 30.5114953 02120
827.754584 74.1841322 50.1154698 1115.96359 10220
386.417627 -131.511472 57.1816898 574.348828 01220
-14.3997089 495.93297 -100.733318 -33.0482214 00320
-335.067423 -61.1885164 -15.5575849 -571.092986 30011
-141.012367 -79.6576921 3.93909126 -238.774398 21011
-40.5956085 -40.6085128 4.55378006 -70.1060089 12011
-9.61175917 -16.4628569 2.40198587 -17.1028707 03011
-135.582711 -224.769833 33.0518928 -215.758705 20111
-83.0751452 -99.4016381 12.8439318 -131.479573 11111
-27.9365238 -36.7235203 4.95506284 -44.2788661 02111
-34.4446694 -18.1458275 0.583708426 -58.7816269 10211
4.52962792 -3.71395976 1.00373335 1.10239295 01211
11.3062521 44.9090338 -7.35413565 26.6700946 00311
-60.9127786 -27.8910625 0.0476036121 -112.252783 30002
-17.1621734 -6.04045209 -0.360067521 -31.4753166 21002
-4.33881232 -0.529213718 -0.288076558 -7.83394703 12002
-0.847869219 0.327736636 -0.142067962 -1.50982205 03002
-13.6199311 -57.6243144 9.70912606 -19.2451109 20102
-4.85724972 -22.1486755 3.80885028 -7.65122318 11102
-1.67244889 -6.67319714 1.13490937 -2.8145588 02102
22.7185586 -8.25065047 3.31457354 26.6742311 10202
4.57599014 -0.331108444 0.40355797 5.71352379 01202
4.82892984 56.9790652 -10.8581986 10.1110942 00302
-22.0502577 481.799669 -95.0884671 -28.8811187 20030
-12.7572862 284.295673 -56.1375502 -16.815359 11030
-0.537550357 139.859471 -27.1099338 0.807755421 02030
308.493215 31.6576716 18.5018416 449.175104 10130
106.861275 -25.7292896 13.7325215 161.454254 01130
-27.3250741 60.403444 -15.0160464 -50.3326469 00230
-21.3083409 -17.5611962 1.75167592 -30.7284198 20021
-12.1931393 -11.0580629 1.2053295 -18.5432804 11021
-3.7827337 -5.2285266 0.724244316 -5.78946788 02021
0.399231817 -0.427289634 0.103318644 -0.318453409 10121
2.14605053 0.393427211 0.0876658872 2.39219477 01121
-2.67831817 9.85904433 -2.02397153 -1.31246518 00221
2.17348691 -2.71949722 0.664349335 4.86230714 20012
0.445088866 -1.93031695 0.39825623 1.00924102 11012
0.107826822 -0.6883375 0.138493389 0.231813003 02012
2.31365026 -1.15044515 0.374161344 1.87374956 10112
0.32585602 -0.227109098 0.061631153 0.166885059 01112
2.42057973 10.7297904 -1.90360649 4.64642829 00212
-0.0714229003 0.517632773 -0.104852572 -0.180726801 20003
-0.014753887 0.172750623 -0.034577206 -0.0448782761 11003
-0.000801963554 0.034690687 -0.00675636864 -0.00723120057 02003
-0.0268653279 -0.281351771 0.0514839842 -0.0769479387 10103
0.0079909334 -0.074911399 0.0149977661 0.00968523611 01103
0.0172303567 -0.0260853991 0.00697449781 -0.00593119924 00203
65.8689638 7.1828781 3.80673702 93.9724548 10040
20.8463901 -2.89561802 2.24650456 31.0841762 01040
-7.33306744 6.73666362 -2.18412322 -14.062807 00140
0.565828761 0.350983293 -0.019873045 0.828627523 10031
0.469000208 0.161283914 0.00563730692 0.604843125 01031
-1.66924761 1.65833653 -0.4362618 -2.03067165 00131
0.233619729 -0.240828688 0.0600196565 0.107053831 10022
0.0260439455 -0.059850412 0.012350628 -0.0128937727 01022
0.39954702 1.62751519 -0.285315519 0.760969995 00122
-0.0133441634 -0.0518328827 0.00885884163 -0.0284698313 10013
-0.00392536519 -0.0136589522 0.00230768841 -0.00725693148 01013
0.013525213 -0.0283700482 0.00660441717 0.0177359374 00113
0.000478508388 0.00527210373 -0.00100903784 0.000226191496 10004
0.000747528015 0.00142585037 -0.00022575333 0.00105517062 01004
0.000154281238 0.00219326893 -0.000426864024 0.000310920444 00104
-1.05368225 1.92538786 -0.523450156 -2.3184317 00050
-0.416617471 -0.0132738096 -0.0266684046 -0.554346008 00041
0.0597422499 0.258764722 -0.0455480693 0.113793274 00032
0.00237047206 -0.00881769056 0.00191019693 0.00336063576 00023
0.000275758827 0.000723512982 -0.00012021094 0.000500819536 00014
4.9019448e-05 -0.000290312675 6.08284584e-05 7.66299186e-05 00005
---------------------------------------------
shms/shms_forward.dat
+ 1
− 1
View file @ 02520bd4
shms-2017-26cm-monte_quads_p18_forward.dat shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_forward.dat
\ No newline at end of file \ No newline at end of file
shms/shms_recon.dat
+ 1
− 1
View file @ 02520bd4
shms-2017-26cm-monte_quads_p18_recon.dat shms-2017-26cm-monte_q1_1018_q2_1027_q3_1018_recon_fit.dat
\ No newline at end of file \ No newline at end of file
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment