| Last change
 on this file since 3871 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.