diff --git a/.viewana.prm b/.viewana.prm index 3cd5ff857919716330fb685e331c3c057246f90c..31646ca55ea0bcea003b4cd23e49e64e09a53ee0 100644 --- a/.viewana.prm +++ b/.viewana.prm @@ -21,116 +21,6 @@ - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - @2D @@ -138,13 +28,3 @@ - - - - - - - - - - diff --git a/anapaw_last.kumacold b/anapaw_last.kumacold index cae4026f9ace690fa02e2b986f77a1ceb8625eaa..5f08b5de5ca8cba182c10fb1454e06c07e9128e6 100644 --- a/anapaw_last.kumacold +++ b/anapaw_last.kumacold @@ -1,8 +1,4 @@ -*** Begin of history file: Wed Jan 26 21:20:24 2022 -a/l -exe kumac/check.exe -exe kumac/check.kumac -exe kumac/check.kumac -a/l -exit -*** End of history file: Wed Jan 26 22:29:32 2022 +*** Begin of history file: Thu Jan 27 00:10:29 2022 +a/l ridf/high0220.ridf +q +*** End of history file: Thu Jan 27 00:12:42 2022 diff --git a/kumac/mpos_dump.kumac b/kumac/mpos_dump.kumac new file mode 100644 index 0000000000000000000000000000000000000000..64480a1f1039cd095f4ea93a2232ad45d22d2ac0 --- /dev/null +++ b/kumac/mpos_dump.kumac @@ -0,0 +1,22 @@ +macro mpos_dump + inu = 1 + hist/plot [1] + vlocate vx vy + + nname = blow $htitle([1]) + + vx1 = vx(1) + vx2 = vx(2) + + if [vx1] > [vx2] then + vxtmp = [vx2] + vx2 = [vx1] + vx1 = [vxtmp] + endif + + xmin = [vx1] + xmax = [vx2] + + hist/fit [1]([xmin]:[xmax]) g + vec/del v* +exitm diff --git a/kumac/pos_cal.kumac b/kumac/pos_cal.kumac new file mode 100644 index 0000000000000000000000000000000000000000..f1ffcdac384355acdee9bea8db2c344a915fa944 --- /dev/null +++ b/kumac/pos_cal.kumac @@ -0,0 +1,3 @@ +lgy +ht 102 +exec kumac/mpos_dump.kumac diff --git a/kumac/read_cid.kumac b/kumac/read_cid.kumac new file mode 100644 index 0000000000000000000000000000000000000000..6cfbd5444ac09a0874963443aaf7ff5dcce2c695 --- /dev/null +++ b/kumac/read_cid.kumac @@ -0,0 +1,6 @@ + macro read_cid + cid + v/del Current_HID + v/cre Current_HID(1) 'I' + v/read Current_HID(1) ./CurrentHID.txt + return Current_HID(1) diff --git a/kumac/xfit_le.kumac b/kumac/xfit_le.kumac new file mode 100644 index 0000000000000000000000000000000000000000..1812480d09f3c00450adab52260f36818c2d6634 --- /dev/null +++ b/kumac/xfit_le.kumac @@ -0,0 +1,77 @@ +macro xfit_le hid=0 textfile=src/mpos_prm_dump.txt chopt= + + exec def_func + + IF ([hid] = 0) THEN + hid = $exec(kumac/read_cid) + ENDIF + 1Dflag = $HINFO([hid],'1DHISTO') + IF ([1Dflag] = 0) THEN + mess ' ' + mess 'This Histogam is not 1D Histogram. Input 1D Historam ID! ' + mess ' ' + stopm + ENDIF + + hi/pl [hid] + vec/del vx + vec/del vy + vlocate vx vy + + np = $vlen(vx,1) + do i=1,[np] + vx[i]= vx([i]) + vy[i]= vy([i]) + enddo + + lg_yi = $sigma(log([vy1])) + lg_yf = $sigma(log([vy[np]])) + p4 = ([vx1]*[lg_yf]-[vx[np]]*[lg_yi])/([vx1]-[vx[np]]) + p5 = ([lg_yi]-[lg_yf])/([vx1]-[vx[np]]) + + cons = $sigma([vy2]-exp([p4]+([p5])*([vx2]))) + mean = [vx2] + width = $sigma(abs(([vx3]-[vx1])/5)) + + npar = 5 + + vec/cre par([npar]) r + vec/cre epar([npar]) r + vec/input par(1:5) [cons] [mean] [width] [p4] [p5] + vec/pri par + + hi/fit [hid]([vx1]:[vx[np]]) l_exp e 5 par ! ! ! epar + + IF ( [chopt] .eq.'d') THEN + hid2 = 101 + while $hexist([hid2]) .eq. 1 do + hid2 = [hid2] + 1 + endwhile + + p4 = par(4) + p5 = par(5) + fun1 [hid2] exp([p4]+([p5])*x) 100 [vx1] [vx[np]] s + hi/del [hid2] + ENDIF + + mes ----------------------------------- + mes P1 : $sigma(par(1)) + mes P2 : $sigma(abs(par(2))) + mes P3 : $sigma(abs(par(3))) + mes P4 : $sigma(par(4)) + mes P5 : $sigma(par(5)) + mes ----------------------------------- + mes peak : $sigma(abs(par(2))) + sh echo [HID] $sigma(abs(par(2))) >> [textfile] + return + + + macro def_func + app comis quit + REAL FUNCTION l_exp(x) + COMMON/PAWPAR/PAR(5) + l_exp = par(1)*denlan((x-par(2))/abs(par(3))) + & + exp(par(4) + par(5)*x) + END + quit + return diff --git a/paw.metafile b/paw.metafile new file mode 100644 index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 diff --git a/src/mpos_prm_dump.txt b/src/mpos_prm_dump.txt new file mode 100644 index 0000000000000000000000000000000000000000..139597f9cb07c5d48bed18984ec4747f4b4f3438 --- /dev/null +++ b/src/mpos_prm_dump.txt @@ -0,0 +1,2 @@ + + diff --git a/src/mssd.prm b/src/mssd.prm index 538969b7c3fc072c35352ebe20361cc335de9ef8..907f1bd460d948070fc4a499b5fa09bc5b09c3bf 100644 --- a/src/mssd.prm +++ b/src/mssd.prm @@ -103,7 +103,7 @@ c365.0, 365.0, 711.0, 711.0, 885.0, 885.0, 1301.0, 1301.0, 1456.0, 1456.0, 1654. c high206 c362.3, 362.3, 709.0, 709.0, 882.0, 882.0, 1303.0, 1303.0, 1440.0, 1440.0, 1640.0, 1640.0 c high219 -386.3, 386.3, 709.0, 709.0, 882.0, 882.0, 1303.0, 1303.0, 1440.0, 1440.0, 1640.0, 1640.0 +c386.3, 386.3, 709.0, 709.0, 882.0, 882.0, 1303.0, 1303.0, 1440.0, 1440.0, 1640.0, 1640.0 c high221 386.7, 386.7, 710.3, 710.3, 886.7, 886.7, 1309.8, 1309.8, 1460.5, 1460.5, 1661.7, 1661.7 c