memo361 : WP-1019
Created Tue Oct 19 04:22:30 2010
Last Modified Wed Oct 20 14:01:31 2010
*
#1 Tue Oct 19 04:22:30 2010 / Tue Oct 19 04:22:30 2010
＜今日の目標＞
・ゲイン調整が終わったらいいなあ
・SOLDIRAsでVETOのピーク位置をシミュレートする。
*
#2 Tue Oct 19 14:50:33 2010 / Wed Oct 20 14:01:31 2010
<gaintuneのストラテジー>
1.移動平均+ピークディテクションでピークを出す
2.出したピークch-535chの値を計算する
3.符号が+ => 行き過ぎ、ゲインを下げる
HVはとりあえず25きざみにしておく。
なおVETOは後回し(シミュレーション待ち)
Initial 1st 2nd 3rd nextHV
U101 -297.81 -205.30 -67.99 -0.47 1725
U102 -280.08 -167.12 -34.35 -3.72 1610
U103 -317.79 -245.12 -132.05 -12.49 1705
U104 -218.38 -157.71 -58.33 -1.27 1725
U105 -195.32 -61.98 70.06 -4.18 1675
U106 -329.68 -281.09 -203.04 -92.96 1850
U107 -240.76 -124.62 6.11 1.78 1698
U108 -88.39 -34.21 18.02 3.22 1745
U109 -336.09 -284.74 -204.32 -85.51 1800
U110 49.94 -37.42 -10.15 0.46 1690
U111 -238.84 -113.56 25.44 19.27 1690
U112 -336.98 -277.10 -187.17 -61.90 1800
U113 -6.00 -16.87 -8.33 -4.95 1705
U114 7.33 2.19 -7.11 -10.32 1698
U115 -123.94 -1.00 -13.36 -10.08 1750
U201 -290.27 -208.12 -95.32 24.51 1790
U202 -321.94 -265.65 -183.39 -64.85 1725
U203 -320.43 -263.83 -181.18 -69.07 1750
U204 -285.42 -199.46 -104.75 9.04 1825
U205 -143.61 -51.82 2.66 9.71 1725
U206 -276.97 -181.17 -77.21 -8.34 1775
U207 -177.64 -72.58 -10.65 -6.24 1727
U208 -176.57 -88.62 -41.47 32.04 1735
U209 -150.90 -53.47 5.30 3.73 1725
U210 -178.78 -77.55 -8.34 -17.98 1740
U211 -111.59 8.72 3.00 4.49 1748
U212 -43.44 9.35 0.91 3.29 1723
U213 -27.85 43.74 -3.19 -0.87 1710
U214 -38.35 33.87 -22.76 -1.77 1690
U215 -177.18 -65.71 -12.29 12.99 1750
D101 -125.79 -9.20 -8.78 7.81 1555
D102 -229.88 -118.38 30.12 -8.24 1590
D103 -195.96 -104.42 -122.82 14.40 1645
D104 -353.90 -305.02 -212.57 -97.81 1700
D105 -4.31 -11.02 1.42 12.83 1652
D106 88.37 -39.82 -47.10 19.29 1605
D107 -125.68 -70.62 -6.93 -4.08 1650
D108 -179.60 -118.46 11.07 -0.73 1670
D109 -123.49 -112.64 9.13 -3.14 1675
D110 -268.51 -165.53 -45.78 20.79 1720
D111 -39.06 12.31 -9.72 -5.20 1670
D112 -29.11 -6.51 -17.58 4.65 1715
D113 -153.02 -38.66 20.54 -3.90 1740
D114 -204.99 -76.29 -21.53 12.49 1735
D115 -60.27 -3.41 -10.46 0.46 1700
D201 -251.50 -159.27 -37.84 -40.89 1785
D202 -305.83 -254.72 -164.61 -162.08 1775*
D203 -303.31 -240.67 -145.08 -157.65 1800*
D204 -210.72 -119.67 -14.80 -13.41 1710
D205 -313.27 -231.64 -185.66 -178.61 1800*
D206 -175.96 -72.31 45.31 46.11 1670
D207 -249.72 -145.00 -14.63 -22.79 1715
D208 -314.95 -241.23 -139.43 -139.72 1800*
D209 -266.27 -163.30 30.47 38.47 1675
D210 -156.93 -48.91 -164.91 -170.15 1750*
D211 -33.99 34.04 -41.82 -36.16 1700
D212 -157.00 -47.84 13.72 7.42 1720
D213 -69.40 -4.85 -4.94 -2.09 1675
D214 -136.70 -17.66 -13.43 0.08 1710
D215 -224.89 -125.89 -5.41 -3.07 1750
*印は飽和している可能性がある。
*
#3 Tue Oct 19 19:37:21 2010 / Tue Oct 19 19:50:35 2010
今日のまとめ
ゲインチューンを３回繰り返した。
↑
クリックで詳細
U側はほとんど調整できている感があるのに対し、D側はなんとも不揃い。
D側のいくつかは既にHVを上げてもゲインが変化しなくなっているので、妥協が必要?
次のiterationでかける予定の電圧は
memo361
のリストの左端のカラムに記述。
＜懸案事項＞
*ゲイン飽和
*VETOのゲイン調整法
*
#4 Tue Oct 19 19:53:19 2010 / Tue Oct 19 19:53:19 2010
ゲイン調整に使ったピーク検知ルーチン
---pkd48.kumac---
macro detect thresh=10 bin=4000 min=10 mm=1 max=4000 id smooth=40
v/del hist
v/cre xx([bin]) R
sigma xx=array([bin],1#[bin])
v/cre hist([bin]) R [bin]*0
v/cre par(4) R
v/cre imean([bin]) R
hi/get_vect/content [id] hist
i=[min]
j=0
zone 1 2
sum = 0
do s=0, [smooth]
sum = [sum] + hist([i]+[s])
enddo
* mess [sum]
lbin = [bin]-[smooth]
while [i]<[lbin] do
sum = [sum] + hist([i]+[smooth])
off = hist([i])
inpp = [sum]/([smooth]+2)
vec/inp hist([i]) [inpp]
sum = [sum] - [off]
i = [i] + 1
endwhile
graph [bin] xx hist
hi/pl [id]
i = [min]
maxv = $sigma(vmax(hist))
mess [maxv]
maxv = [maxv]*0.75
loop:
if [i]<[lbin] then
if hist([i])<[maxv] then
i = [i]+1
goto loop
endif
imin=[i]
while hist([i])>=[maxv] do
i = [i]+1
if [i]>=[bin] then
breakl
endif
endwhile
imax=[i]
endif
* mess [imin]
* mess [imax]
v/fit xx([imin]:[imax]) hist([imin]:[imax]) ! g QN 0 par
revised = par(2)
revised = [revised] + [smooth]/2
v/inp par(2) [revised]
v/write par(2) ! '(1x,f8.4)'
*
* v/fit xx([imin]:[imax]) hist([imin]:[imax]) ! g QN 0 par
* v/write par(2) ! '(f7.1)'
* par(3) = [imin]
* par(4) = [imax]
** v/input par(1) [imax]
return
macro many thresh=10 bin=4000 min=10 mm=1 max=4000 from num=1 smooth=40 center=323
vec/cre storepeak([num]) R
vec/cre storeoffset([num]) R
vec/cre ids([num]) R
do i=1,[num]
id=[from]+[i]-1
exec kumac/pkd48#detect [thresh] [bin] [min] [mm] [max] [id] [smooth]
vec/inp ids([i]) [id]
vec/inp storepeak([i]) par(2)
cent = par(2)-[center]
vec/inp storeoffset([i]) [cent]
wait
enddo
vec/write ids,storepeak,storeoffset ! '(f4, 1x, f8.2, 1x, f8.2)'
vec/write ids,storepeak ./gaintune.txt '(f4, 1x, f8.2)'
return
*
Link to:
*
Linked from:
memo361 : WP-1019
*#1 Tue Oct 19 04:22:30 2010 / Tue Oct 19 04:22:30 2010
＜今日の目標＞
・ゲイン調整が終わったらいいなあ
・SOLDIRAsでVETOのピーク位置をシミュレートする。
*#2 Tue Oct 19 14:50:33 2010 / Wed Oct 20 14:01:31 2010
1.移動平均+ピークディテクションでピークを出す
2.出したピークch-535chの値を計算する
3.符号が+ => 行き過ぎ、ゲインを下げる
HVはとりあえず25きざみにしておく。
なおVETOは後回し(シミュレーション待ち)
Initial 1st 2nd 3rd nextHV
U101 -297.81 -205.30 -67.99 -0.47 1725
U102 -280.08 -167.12 -34.35 -3.72 1610
U103 -317.79 -245.12 -132.05 -12.49 1705
U104 -218.38 -157.71 -58.33 -1.27 1725
U105 -195.32 -61.98 70.06 -4.18 1675
U106 -329.68 -281.09 -203.04 -92.96 1850
U107 -240.76 -124.62 6.11 1.78 1698
U108 -88.39 -34.21 18.02 3.22 1745
U109 -336.09 -284.74 -204.32 -85.51 1800
U110 49.94 -37.42 -10.15 0.46 1690
U111 -238.84 -113.56 25.44 19.27 1690
U112 -336.98 -277.10 -187.17 -61.90 1800
U113 -6.00 -16.87 -8.33 -4.95 1705
U114 7.33 2.19 -7.11 -10.32 1698
U115 -123.94 -1.00 -13.36 -10.08 1750
U201 -290.27 -208.12 -95.32 24.51 1790
U202 -321.94 -265.65 -183.39 -64.85 1725
U203 -320.43 -263.83 -181.18 -69.07 1750
U204 -285.42 -199.46 -104.75 9.04 1825
U205 -143.61 -51.82 2.66 9.71 1725
U206 -276.97 -181.17 -77.21 -8.34 1775
U207 -177.64 -72.58 -10.65 -6.24 1727
U208 -176.57 -88.62 -41.47 32.04 1735
U209 -150.90 -53.47 5.30 3.73 1725
U210 -178.78 -77.55 -8.34 -17.98 1740
U211 -111.59 8.72 3.00 4.49 1748
U212 -43.44 9.35 0.91 3.29 1723
U213 -27.85 43.74 -3.19 -0.87 1710
U214 -38.35 33.87 -22.76 -1.77 1690
U215 -177.18 -65.71 -12.29 12.99 1750
D101 -125.79 -9.20 -8.78 7.81 1555
D102 -229.88 -118.38 30.12 -8.24 1590
D103 -195.96 -104.42 -122.82 14.40 1645
D104 -353.90 -305.02 -212.57 -97.81 1700
D105 -4.31 -11.02 1.42 12.83 1652
D106 88.37 -39.82 -47.10 19.29 1605
D107 -125.68 -70.62 -6.93 -4.08 1650
D108 -179.60 -118.46 11.07 -0.73 1670
D109 -123.49 -112.64 9.13 -3.14 1675
D110 -268.51 -165.53 -45.78 20.79 1720
D111 -39.06 12.31 -9.72 -5.20 1670
D112 -29.11 -6.51 -17.58 4.65 1715
D113 -153.02 -38.66 20.54 -3.90 1740
D114 -204.99 -76.29 -21.53 12.49 1735
D115 -60.27 -3.41 -10.46 0.46 1700
D201 -251.50 -159.27 -37.84 -40.89 1785
D202 -305.83 -254.72 -164.61 -162.08 1775*
D203 -303.31 -240.67 -145.08 -157.65 1800*
D204 -210.72 -119.67 -14.80 -13.41 1710
D205 -313.27 -231.64 -185.66 -178.61 1800*
D206 -175.96 -72.31 45.31 46.11 1670
D207 -249.72 -145.00 -14.63 -22.79 1715
D208 -314.95 -241.23 -139.43 -139.72 1800*
D209 -266.27 -163.30 30.47 38.47 1675
D210 -156.93 -48.91 -164.91 -170.15 1750*
D211 -33.99 34.04 -41.82 -36.16 1700
D212 -157.00 -47.84 13.72 7.42 1720
D213 -69.40 -4.85 -4.94 -2.09 1675
D214 -136.70 -17.66 -13.43 0.08 1710
D215 -224.89 -125.89 -5.41 -3.07 1750
*印は飽和している可能性がある。
*#3 Tue Oct 19 19:37:21 2010 / Tue Oct 19 19:50:35 2010
ゲインチューンを３回繰り返した。
↑
クリックで詳細
U側はほとんど調整できている感があるのに対し、D側はなんとも不揃い。
D側のいくつかは既にHVを上げてもゲインが変化しなくなっているので、妥協が必要?
次のiterationでかける予定の電圧はmemo361のリストの左端のカラムに記述。
＜懸案事項＞
*ゲイン飽和
*VETOのゲイン調整法
*#4 Tue Oct 19 19:53:19 2010 / Tue Oct 19 19:53:19 2010
---pkd48.kumac---
macro detect thresh=10 bin=4000 min=10 mm=1 max=4000 id smooth=40
v/del hist
v/cre xx([bin]) R
sigma xx=array([bin],1#[bin])
v/cre hist([bin]) R [bin]*0
v/cre par(4) R
v/cre imean([bin]) R
hi/get_vect/content [id] hist
i=[min]
j=0
zone 1 2
sum = 0
do s=0, [smooth]
sum = [sum] + hist([i]+[s])
enddo
* mess [sum]
lbin = [bin]-[smooth]
while [i]<[lbin] do
sum = [sum] + hist([i]+[smooth])
off = hist([i])
inpp = [sum]/([smooth]+2)
vec/inp hist([i]) [inpp]
sum = [sum] - [off]
i = [i] + 1
endwhile
graph [bin] xx hist
hi/pl [id]
i = [min]
maxv = $sigma(vmax(hist))
mess [maxv]
maxv = [maxv]*0.75
loop:
if [i]<[lbin] then
if hist([i])<[maxv] then
i = [i]+1
goto loop
endif
imin=[i]
while hist([i])>=[maxv] do
i = [i]+1
if [i]>=[bin] then
breakl
endif
endwhile
imax=[i]
endif
* mess [imin]
* mess [imax]
v/fit xx([imin]:[imax]) hist([imin]:[imax]) ! g QN 0 par
revised = par(2)
revised = [revised] + [smooth]/2
v/inp par(2) [revised]
v/write par(2) ! '(1x,f8.4)'
*
* v/fit xx([imin]:[imax]) hist([imin]:[imax]) ! g QN 0 par
* v/write par(2) ! '(f7.1)'
* par(3) = [imin]
* par(4) = [imax]
** v/input par(1) [imax]
return
macro many thresh=10 bin=4000 min=10 mm=1 max=4000 from num=1 smooth=40 center=323
vec/cre storepeak([num]) R
vec/cre storeoffset([num]) R
vec/cre ids([num]) R
do i=1,[num]
id=[from]+[i]-1
exec kumac/pkd48#detect [thresh] [bin] [min] [mm] [max] [id] [smooth]
vec/inp ids([i]) [id]
vec/inp storepeak([i]) par(2)
cent = par(2)-[center]
vec/inp storeoffset([i]) [cent]
wait
enddo
vec/write ids,storepeak,storeoffset ! '(f4, 1x, f8.2, 1x, f8.2)'
vec/write ids,storepeak ./gaintune.txt '(f4, 1x, f8.2)'
return
*Link to: