source: BAORadio/AmasNancay/transfert/transfert/Abell85/etude_ds.pic @ 673

#exemple exec etude_ds 0 1411.8251 0.5

exec ds


set numcycle $1
# options de scan: num cycle, directory, first file, last file+1, central freq., bandwidth
set f0 $2
set bw $3

set spectype "medfiltmtx"
set namcycle "onoffscanmedcycle"
#set spectype "mspecmtx"
#set namcycle "nocutcycle"

set firstFile 0
set lastp1File 46


scan $numcycle $namcycle  $firstFile $lastp1File $f0 $bw $spectype


#vecfrascii dataCh0${numcycle} dataCh0${numcycle}.dat
#vecfrascii dataCh1${numcycle} dataCh1${numcycle}.dat

del data${numcycle}

c++exec TVector<r_8> data${numcycle} = ((r_8) 0.5)*(dataCh0${numcycle} + dataCh1${numcycle}); KeepObj(data${numcycle});




#Nb: 7/2/11 pour transformer n et sec, le coef = taille fenetre en nbre de paquets/taux de trigger en Hz
newwin 1 1
plot2d dataCh0${numcycle} n*5120/8500 val 1 "cpts blue nsta notit"
settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle}, Ch 0)"
setaxelabels "T (sec)" "Intensity (Jy)"

newwin 1 1
plot2d dataCh1${numcycle} n*5120/8500 val 1 "cpts blue nsta notit"
settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle}, Ch 1)"      
setaxelabels "T (sec)" "Intensity (Jy)"

newwin 1 1
plot2d data${numcycle} n*5120/8500 val 1 "cpts blue nsta notit"
settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle})"    
setaxelabels "T (sec)" "Mean Intensity (Jy)"

newwin 1 1
plot2d data${numcycle} n val 1 "cpts blue nsta notit"
settitle "Drift scan  at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle})"
setaxelabels "Window num." "Mean Intensity (Jy)"