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