| Last change
 on this file since 3762 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 | 
      
      
| Rev | Line |  | 
|---|
| [3683] | 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.