[673] | 1 | #exemple exec etude_ds 0 1411.8251 0.5 |
---|
| 2 | |
---|
| 3 | exec ds |
---|
| 4 | |
---|
| 5 | |
---|
| 6 | set numcycle $1 |
---|
| 7 | # options de scan: num cycle, directory, first file, last file+1, central freq., bandwidth |
---|
| 8 | set f0 $2 |
---|
| 9 | set bw $3 |
---|
| 10 | |
---|
| 11 | set spectype "medfiltmtx" |
---|
| 12 | set namcycle "onoffscanmedcycle" |
---|
| 13 | #set spectype "mspecmtx" |
---|
| 14 | #set namcycle "nocutcycle" |
---|
| 15 | |
---|
| 16 | set firstFile 0 |
---|
| 17 | set lastp1File 46 |
---|
| 18 | |
---|
| 19 | |
---|
| 20 | scan $numcycle $namcycle $firstFile $lastp1File $f0 $bw $spectype |
---|
| 21 | |
---|
| 22 | |
---|
| 23 | #vecfrascii dataCh0${numcycle} dataCh0${numcycle}.dat |
---|
| 24 | #vecfrascii dataCh1${numcycle} dataCh1${numcycle}.dat |
---|
| 25 | |
---|
| 26 | del data${numcycle} |
---|
| 27 | |
---|
| 28 | c++exec TVector<r_8> data${numcycle} = ((r_8) 0.5)*(dataCh0${numcycle} + dataCh1${numcycle}); KeepObj(data${numcycle}); |
---|
| 29 | |
---|
| 30 | |
---|
| 31 | |
---|
| 32 | |
---|
| 33 | #Nb: 7/2/11 pour transformer n et sec, le coef = taille fenetre en nbre de paquets/taux de trigger en Hz |
---|
| 34 | newwin 1 1 |
---|
| 35 | plot2d dataCh0${numcycle} n*5120/8500 val 1 "cpts blue nsta notit" |
---|
| 36 | settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle}, Ch 0)" |
---|
| 37 | setaxelabels "T (sec)" "Intensity (Jy)" |
---|
| 38 | |
---|
| 39 | newwin 1 1 |
---|
| 40 | plot2d dataCh1${numcycle} n*5120/8500 val 1 "cpts blue nsta notit" |
---|
| 41 | settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle}, Ch 1)" |
---|
| 42 | setaxelabels "T (sec)" "Intensity (Jy)" |
---|
| 43 | |
---|
| 44 | newwin 1 1 |
---|
| 45 | plot2d data${numcycle} n*5120/8500 val 1 "cpts blue nsta notit" |
---|
| 46 | settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle})" |
---|
| 47 | setaxelabels "T (sec)" "Mean Intensity (Jy)" |
---|
| 48 | |
---|
| 49 | newwin 1 1 |
---|
| 50 | plot2d data${numcycle} n val 1 "cpts blue nsta notit" |
---|
| 51 | settitle "Drift scan at ${f0} +/- [${bw}/2]MHz (cycle ${numcycle})" |
---|
| 52 | setaxelabels "Window num." "Mean Intensity (Jy)" |
---|