paramprompt off read data 1 off_ACIS-S.dat ascii 1 2 read data 2 off_ACIS-I.dat ascii 1 2 read data 3 off_HRC-S.dat ascii 1 2 read data 4 off_HRC-I.dat ascii 1 2 read data 5 off_ACIS-S.dat ascii 1 3 read data 6 off_ACIS-I.dat ascii 1 3 read data 7 off_HRC-S.dat ascii 1 3 read data 8 off_HRC-I.dat ascii 1 3 errors 1:8 = 0.2 ignore 1:8 filter :50 shexp[dyexp] shexp[dzexp] poly[asy] poly[asz] poly[aiy] poly[aiz] poly[hsy] poly[hsz] poly[hiy] poly[hiz] asy integrate off asz integrate off aiy integrate off aiz integrate off hsy integrate off hsz integrate off hiy integrate off hiz integrate off #freeze asy #freeze asz #freeze aiy #freeze aiz #freeze hsy #freeze hsz #freeze hiy #freeze hiz freeze dyexp freeze dzexp dyexp.offset = 0 dyexp.ampl = 3:2.0:4.0 dyexp.coeff = -0.003:-0.006:-0.001 dzexp.offset = 0 dzexp.ampl = 2:1.0:3.0 dzexp.coeff = -0.003:-0.006:-0.001 asy.c0 = -2.9:-10:10 asz.c0 = 2.0:-10:10 aiy.c0 = -3.2:-10:10 aiz.c0 = 2.0:-10:10 hsy.c0 = -2.5:-10:10 hsz.c0 = 2.0:-10:10 hiy.c0 = -1.5:-10:10 hiz.c0 = 1.5:-10:10 source 1 = dyexp + asy source 2 = dyexp + aiy source 3 = dyexp + hsy source 4 = dyexp + hiy source 5 = -dzexp + asz source 6 = -dzexp + aiz source 7 = -dzexp + hsz source 8 = -dzexp + hiz # fit 1,5