|
Last change
on this file since 3805 was 3683, checked in by ansari, 16 years ago |
Mise a jour et ajout de fichier pour taritement multifibres apres
prise de donnees de Nov2009 a Pittsburgh
- Introduction des classes BRMultiFitsReader et BRBaseProcessor
Reza, 27/11/2009
|
|
File size:
1.9 KB
|
| Line | |
|---|
| 1 | # Script affichage spectre moyennes fait par Monitoring
|
|---|
| 2 | if ( $# < 2 ) then
|
|---|
| 3 | echo ' dspms.pic/Erreur - Usage dspms Directory NumFile [BaseDir] '
|
|---|
| 4 | return
|
|---|
| 5 | endif
|
|---|
| 6 | set dir $1
|
|---|
| 7 | set numf $2
|
|---|
| 8 |
|
|---|
| 9 | set bdir /Raid/PittsNov09
|
|---|
| 10 | if ( $# > 2 ) then
|
|---|
| 11 | set bdir $3
|
|---|
| 12 | endif
|
|---|
| 13 |
|
|---|
| 14 | delobjs *
|
|---|
| 15 | # openppf /Raid/PittsNov09/$dir/meanspec$numf.ppf
|
|---|
| 16 | set flnm $bdir/$dir/meanspec$numf.ppf
|
|---|
| 17 | echo '------ Opening file ' $flnm
|
|---|
| 18 | openppf $flnm
|
|---|
| 19 | del myspec
|
|---|
| 20 | mv spectre myspec
|
|---|
| 21 | # mv /home/meanspec$numf myspec
|
|---|
| 22 | set nchan ${myspec.nrows}
|
|---|
| 23 | for i 0:$nchan
|
|---|
| 24 | del smV$i
|
|---|
| 25 | objaoper myspec row $i smV$i
|
|---|
| 26 | end
|
|---|
| 27 |
|
|---|
| 28 | setaxesatt 'font=helvetica,bold,16 fixedfontsize minorticks'
|
|---|
| 29 | zone 1
|
|---|
| 30 |
|
|---|
| 31 | set exfreq (n+2212)*250/4096+1250
|
|---|
| 32 | # set exfreq (n*250/4096+1250)
|
|---|
| 33 | plot2d smV0 $exfreq 10*log10(val) val>1.e-6 'black nsta cpts notit'
|
|---|
| 34 | setaxelabels ' Freq (MHz)' 'Spectra - dB' 'font=helvetica,bolditalic,16'
|
|---|
| 35 |
|
|---|
| 36 | set cols ( black darkviolet blue navyblue red orange yellow green )
|
|---|
| 37 |
|
|---|
| 38 | for i 1:${myspec.nrows}
|
|---|
| 39 | plot2d smV$i $exfreq 10*log10(val) val>1.e-6 "$cols[i] cpts nsta same"
|
|---|
| 40 | end
|
|---|
| 41 |
|
|---|
| 42 | if ( $nchan < 3 ) then
|
|---|
| 43 | zny = $nchan
|
|---|
| 44 | znx = 1
|
|---|
| 45 | else
|
|---|
| 46 | if ( $nchan < 7 ) then
|
|---|
| 47 | zny = $nchan/2
|
|---|
| 48 | znx = 2
|
|---|
| 49 | else
|
|---|
| 50 | zny = 3
|
|---|
| 51 | znx = 3
|
|---|
| 52 | endif
|
|---|
| 53 | endif
|
|---|
| 54 |
|
|---|
| 55 | defscript dspec
|
|---|
| 56 | zone $znx $zny
|
|---|
| 57 | for i 0:$nchan
|
|---|
| 58 | disp smV$i "$cols[i]"
|
|---|
| 59 | end
|
|---|
| 60 | endscript
|
|---|
| 61 |
|
|---|
| 62 | defscript dspeclg
|
|---|
| 63 | zone $znx $zny
|
|---|
| 64 | for i 0:$nchan
|
|---|
| 65 | disp smV$i "logy $cols[i]"
|
|---|
| 66 | end
|
|---|
| 67 | endscript
|
|---|
| 68 |
|
|---|
| 69 | defscript dhval
|
|---|
| 70 | zone $znx $zny
|
|---|
| 71 | for i 0:$nchan
|
|---|
| 72 | disp hvalV$i "$cols[i]"
|
|---|
| 73 | end
|
|---|
| 74 | endscript
|
|---|
| 75 |
|
|---|
| 76 | defscript dh2
|
|---|
| 77 | zone $znx $zny
|
|---|
| 78 | for i 0:$nchan
|
|---|
| 79 | plot2d h2dV$i x y val>0 "notit nsta $cols[i]"
|
|---|
| 80 | end
|
|---|
| 81 | endscript
|
|---|
| 82 |
|
|---|
| 83 |
|
|---|
| 84 | set lines ( B1.V1 B1.V2 B1.V3 B1.V4 B2.V1 B2.V2 B2.V3 B2.V4 )
|
|---|
| 85 | textdrawer lines cols 'font=helvetica,bold,16'
|
|---|
| 86 |
|
|---|
| 87 | # settitle 'Spectre 8 voies FFT - ScalingFactor=9' ' ' 'font=helvetica,bold,18'
|
|---|
| 88 |
|
|---|
| 89 | newwin
|
|---|
| 90 | dspeclg
|
|---|
| 91 | newwin
|
|---|
| 92 | dspec
|
|---|
| 93 | newwin
|
|---|
| 94 | dhval
|
|---|
| 95 | newwin
|
|---|
| 96 | dh2
|
|---|
| 97 |
|
|---|
| 98 | echo '----> use dspec dhval and dh2 scripts to check value histograms'
|
|---|
| 99 |
|
|---|
| 100 |
|
|---|
Note:
See
TracBrowser
for help on using the repository browser.