source: Sophya/trunk/Cosmo/RadioBeam/formnoise.pic

#######################################################
##### Script de trace de PNoise pour differentes configurations 
####### BAORadio --  R. Ansari - Avril 2008 
### ===> Fichiers formnoise.pic + vnpairs.cc + vnpairs.h

## --- Script d'execution 
defscript go 
  delobjs *
  c++execfrf vnpairs.cc vnpairs.h
endscript 
## ---- Script prjkt H2NPairsName ProfName Nbin 
defscript prjkt 'Usage: prjkt H2NPairsName ProfName Nbin' 
  if ( $# < 3 ) then
     echo 'Usage: prjkt H2NPairsName ProfName Nbin'
     return
  endif 
  set h2n $1
  set prn k_theta_$2
  set nbin $3
  # Longueur d'onde a z ~ 1 
  LAMBDA = 0.42
  # Facteur de distance transverse a z ~ 1 (Mpc)
  DCT = 3400.
  newprof $prn 0. 3000 $nbin 
  projprof $prn $h2n 2*M_PI*sqrt(x*x+y*y)/$LAMBDA val 
  set prn k_Mpc_$2
  newprof $prn 0. 1. $nbin 
  projprof $prn $h2n 2*M_PI*sqrt(x*x+y*y)/$LAMBDA/$DCT val 
  disp $prn 
endscript 

## --- Script de representation de config antennes
defscript drantennes
  zone 1 
  nt2d recp xr yr - - - - 'red marker=circle,13 nsta '
  settitle 'Configuration antennes - D=5m' ' ' 'font=times,bolditalic,18'
  setaxelabels 'x (metres)' 'y (metres)' 'font=times,bolditalic,18' 
endscript 


defscript toto 
setaxesatt 'font=times,bold,18 fixedfontsize minorticks'
# go
zone 1 3
prjkt h2_fa fa 100 
prjkt h2_fa_d2m fa_d2m 200 
prjkt h2hex hex 200 

# Config A : Packed 2 meter dishes 100m x 100m   TPose ~ 6.25 10^4 sec  (6.25 = (5/2)^2)
# Config B : Packed 5 meter dishes 100m x 100m   TPose ~ 10^4 sec
# Config C : 5 meter dishes - config hexa, cercle R ~ 100 m   TPose ~ 10^4 sec
# Config D : 5 meter dishes - config hexa, cercle R ~ 100 m  TPose ~ 50 10^4 sec

# Nb de paires A , B , C 
set npairs ( 2450 380 122 6000 )
set kMpc ( k_Mpc_fa_d2m k_Mpc_fa k_Mpc_hex k_Mpc_hex ) 
set maxprof ( ${k_Mpc_fa_d2m.vmax}  ${k_Mpc_fa.vmax} ${k_Mpc_hex.vmax} ${k_Mpc_hex.vmax} )  

echo $maxprof
set cols ( blue magenta red orange black )
zone 1
for i 0:4
  echo "----- I= $i " 
  ph = "$kMpc[i]"
  maxph = $maxprof[i]
  npair = $npairs[i]
  set gropt "nsta same connectpoints line=solid,2 $cols[i]" 
  if ( $i == 0 ) then 
    set gropt "xylimits=0.005,1.,1,5000 logx logy nsta connectpoints line=solid,2 $cols[i] notit" 
  endif
  plot2d $ph x 1e4/val/$npair*$maxph val>5 "$gropt" 
end
settitle 'PNoise (mK^2) Mpc^2  z~1' ' ' 'font=times,bolditalic,18'
setaxelabels 'k_Transverse Mpc^-1 ' ' PNoise (mK^2 Mpc^2)'  'font=times,bolditalic,18'

addline 0.04 1000 0.4 20 'line=dashed,2 black' 
set txt ( '(A) 100mx100m Packed 2m dishes, tint=6.2 10^4 s' '(B) 100mx100m Packed 5m dishes, tint=10^4 s' '(C) 400x5m dishes over 200m circle, , tint=10^4s' "(C') 400x5m dishes over 200m circle, , tint=50 10^4s" 'P(k) HI' ) 
textdrawer txt cols 'font=times,bold,16 frame'

endscript