| 
            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
           | 
        
      
      
| 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.