Procedure:

base = '/towerdata/meas/adsec_calib/if_functions/20130831_fittedKLs_04/'
m2c = readfits(base+'full_m2v.fits')
ff = adsec.ff_matrix
force = ff ## m2c
force = abs(force)
ring0 = adsec.mir_act2act[0:8]

print, intersection(wcl, ring0, ring0good)
out = adsec.mir_act2act[9:*]

print, intersection(wcl, out, outgood)
fr0 = dblarr(500)
for i = 0, 499 do fr0[i]=rms(force[i,ring0good])
fout = dblarr(500)
for i = 0, 499 do fout[i]=rms(force[i,outgood])

; it appears that the force on ring0 are much larger than on the other actuators

comparison with older m2v

base1 = '/towerdata/adsec_calib/M2C/KL_v7/'
m2c1 = readfits(base1+'m2c.fits')

ffbase ='/towerdata/meas/adsec_calib/if_functions/20100928_064906/adsec_save.sav' ; maybe correct

restore, ffbase

ff1 = adsec_save.ff_matrix

wcl1 = adsec_save.act_w_cl

force1 = ff1 ## m2c1
force1 = abs(force1)
ring0 = adsec.mir_act2act[0:8]

print, intersection(wcl1, ring0, ring0good1)
out = adsec.mir_act2act[9:*]
fr01 = dblarr(500)
for i = 0, 499 do fr01[i]=rms(force1[i,ring0])
fout1 = dblarr(500)
for i = 0, 499 do fout1[i]=rms(force1[i,out])

the comparison shows that with KLV7 the forces where basically similar

test of modal thresholding

base = '/towerdata/meas/adsec_calib/if_functions/20130831_fittedKLs_04/'
m2c = readfits(base+'full_m2v.fits')
ff = adsec.ff_matrix

wcl = adsec.act_w_cl

fhist = '/towerdata/meas/adsec_calib/dist_flao1_ts4_atm_s0.8_L040.0_v15.0_ovfreq1000.00_sd1983.fits

fhist = readfits(fhist) *2

flat = '/towerdata/meas/adsec_calib/flat/Applicable/default.sav'

restore, flat

fflat = transpose(flattened_status.current)

fflat = rebin(fflat, 4000, 672, /samp)

c2m = pseudo_invert(m2c, count_z =aa)

klc = c2m ## fhist

n2cut = [500,465,406,351,300]-2

maxforce = dblarr(n_elements(n2cut), 2)

klc1 = klc*0

for i = 0, n_elements(n2cut) -1 do begin

klc[*,n2cut[i]:*] = 0

;klc1[*,0:n2cut[i]-1] = klc[*,0:n2cut[i]-1]

mact = m2c ## klc

fact = ff ## (mact)+fflat

maxforce[i,*] = minmax(fact[*,wcl])

print, minmax(fact[*,wcl])

endfor
Topic revision: r4 - 10 Sep 2013, RunaBriguglio
This site is powered by FoswikiCopyright © by the contributing authors. All material on this collaboration platform is the property of the contributing authors.
Ideas, requests, problems regarding Foswiki? Send feedback