Ignore:
Timestamp:
Apr 18, 2011, 5:30:44 PM (14 years ago)
Author:
ansari
Message:

Corrections diverses: choix lobe gaussien/triangle et specif DishDiameter au lieu de DoL ds applobe/calcpk2, possibilite application lobe freq.independante ds applobe, Reza 18/04/2011

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/RadioBeam/subtractradsrc.cmd

    r3830 r3973  
    1212csh> ~/Objs/exe/cmvginit3df -a -1 -2 -C -G 0. -F 0 -x 360,3 -y 360,3 -z 256,3 -Z 0.60 -8 1. -n 10000 -O 0,2 -o lssz060 -T 2
    1313
     14# 1.c/ To run SimLSS with GSM map parametersand 40x40 deg maps (DeltaFreq=500 MHz)
     15csh> ~/Objs/exe/cmvginit3df -a -1 -2 -C -G 0. -F 0 -x 600,1.9 -y 800,1.9 -z 256,2.8 -Z 0.60 -8 1. -n 10000 -O 0,2 -o lssz060 -T 2
     16
    1417# 1.c/ Change the X and Z axis of the cube to adapt it to RadioBeam package convention
    1518#  SimLSS output : the radial (redshift) direction along X axis of the cube (TArray)
    … …  
    1821
    1922csh> cat > racube.pic
    20 set f lssz056
     23set f lssz060
    2124readfits ${f}_r.fits
    2225rename ${f}_r map
    … …  
    3336print lsscube
    3437# expmeansig lsscube val
    35 saveppf lsscube lsscube.ppf
     38saveppf lsscube lsscubez060.ppf
    3639
    3740csh> spiapp -term -exec racube.pic
    3841
     42#### Cube LSS 40x30 deg (3')  @ z=0.6 ( lsscubez060.ppf )
     43#### -> Size= 122880000 Mean=-7.01664e-05 Sigma=2.53016 Min=-13.7439 Max=14.4648
    3944
    4045## Step 2/ Produce synchrotron and radio source sky cubes  (cube unit is Temparature- Kelvin)
    4146# 2.a/ Synchrotron map from HASLAM 400 MHz map
    42 csh> ./Objs/syncube syncmap_eq.fits syncube.ppf
     47csh> ./Objs/syncube syncmap_eq.fits syncube.ppf syncmap.ppf
    4348# 2.b/ radio source cube from NVSS catalog
    44 csh> ./Objs/srcat2cube nvss.fits nvsscube.ppf
     49csh> ./Objs/srcat2cube -nvss nvss.fits nvsscube.ppf nvssmap.ppf
     50# Or from the north20 catalog :
     51csh> ./Objs/srcat2cube -north20 north20cm.fits north20cube.ppf north20map.ppf
     52
    4553# 2.c/ Add the two cubes using the following spiapp script
    4654csh> cat > sumcubes.pic
    4755openppf syncube.ppf
    48 openppf srcnor3d.ppf
     56openppf radsrccube.ppf
    4957# expmeansig syncube val
    50 # expmeansig srcnor3d val
    51 c++exec TArray<r_4> fgndcube = syncube+srcnor3d; KeepObj(fgndcube);
     58# expmeansig nvsscube val
     59c++exec TArray<r_4> fgndcube = syncube+radsrccube; KeepObj(fgndcube);
    5260print fgndcube
    5361# expmeansig fgndcube val
    … …  
    5664csh> spiapp -term -exec sumcubes.pic
    5765
     66#### syncube:Mean= 1.8101  Sigma= 0.326538  Min= 0.857019 Max= 3.58987
     67#### nvsscube: Mean= 1.95073  Sigma= 1.68515  Min= 0.857019 Max= 428.398
     68#### fgndcube=syncube+nvsscube: Mean= 0.140623  Sigma= 1.65068  Min= 0 Max= 426.559
     69#### north20: fgndcube_north:
     70
     71## Step 2.b/ Produce foreground cube from GSM
     72csh> ./Objs/gsm2cube ../Catalogs/GSM/ 1 256 fgndcube_gsm.ppf
     73
    5874## Step 3/ Apply lobe (50 meter diameter array) effect on foreground cube and LSS cube
    59 csh> ./Objs/applobe 50. fgndcube.ppf fgndcube_lobe.ppf
    60 csh> ./Objs/applobe 50. lsscube.ppf lsscube_lobe.ppf
    61 ## Step 3.b/ Correct for the lobe effect by bringing all to the beam of Diam/Lambda = 130 
    62 csh> ./Objs/applobe 50. lsscube_lobe.ppf lsscube_corlobe.ppf 130
    63 csh> ./Objs/applobe 50. fgndcube_lobe.ppf fgndcube_corlobe.ppf 130
     75csh> set ddish=55.
     76csh> set ddishcor=50.
     77csh> ./Objs/applobe $ddish fgndcube.ppf fgndcube_lobe.ppf
     78csh> ./Objs/applobe -fib $ddish fgndcube.ppf fgndcube_flobe.ppf
     79csh> ./Objs/applobe $ddish lsscube.ppf lsscube_lobe.ppf
     80csh> ./Objs/applobe -fib $ddish lsscube.ppf lsscube_flobe.ppf
     81## Step 3.b/ Correct for the lobe effect by bringing all to the beam of Diam/Lambda = 150  (55 m @ z=0.7 - 820 MHz) 
     82csh> ./Objs/applobe $ddish lsscube_lobe.ppf lsscube_corlobe.ppf $ddishcor
     83csh> ./Objs/applobe $ddish fgndcube_lobe.ppf fgndcube_corlobe.ppf $ddishcor
     84
     85## Step 3.c/ Apply lobe (Filled 11x11 5m dishes array) effect on foreground cube and LSS cube
     86csh> ./Objs/applobe repf11x11.ppf fgndcube.ppf fgndcube_lobe.ppf
     87csh> ./Objs/applobe repf11x11.ppf lsscube.ppf lsscube_lobe.ppf
     88## Step 3.d/ Correct for the lobe effect by bringing all to the beam of Diam/Lambda = 150 (55 m @ z=0.7 - 820 MHz)
     89csh> ./Objs/applobe repf11x11.ppf lsscube_lobe.ppf lsscube_corlobe.ppf $ddishcor
     90csh> ./Objs/applobe repf11x11.ppf fgndcube_lobe.ppf fgndcube_corlobe.ppf $ddishcor
    6491
    6592### Step 4/ Compute power spectra
    66 ## mass to temperature converion factor   CT21 ~= 0.2 mK
     93## mass to temperature converion factor   CT21 ~= 0.21 mK for gHI=2% , 0.11 for gHI=1% , 0.13 for gHI=0.008x(1+0.6)
    6794## Foreground maps are in temperature
    6895## Noise fluctuations Sigma^2 ~ T_sys^2 / t_obs * DeltaFreq
    69 ## Tsys ~ 50 K , DeltaFreq ~ 0.275 MHz , t_obs ~ 1 day ~ 80 000 s.
    70 ## sigma_noise ~ 0.35 mK
     96## Tsys ~ 50 K , DeltaFreq ~ 0.5 MHz , t_obs ~ 1 day ~ 80 000 s.
     97## sigma_noise ~ 0.25 mK -> 3 mK
    7198# 4.a/ LSS power spectrum  without noise
    72 csh> ./Objs/calcpk lsscube.ppf lsspk.ppf 0.2
     99csh> ./Objs/calcpk lsscube.ppf lsspk.ppf 0.13
    73100# and with noise 
    74 csh> ./Objs/calcpk lsscube.ppf lsspkwn.ppf 0.2 0.35
     101csh> ./Objs/calcpk lsscube.ppf lsspkwn.ppf 0.13  3
    75102#  with the lobe effect
    76 csh> ./Objs/calcpk lsscube_lobe.ppf lsspklobe.ppf 0.2 
    77 csh> ./Objs/calcpk lsscube_lobe.ppf lsspklobewn.ppf 0.2 0.35 
    78 csh> ./Objs/calcpk lsscube_corlobe.ppf lsspkcorlobe.ppf 0.2 
     103csh> ./Objs/calcpk lsscube_lobe.ppf lsspklobe.ppf 0.13
     104csh> ./Objs/calcpk lsscube_flobe.ppf lsspkflobe.ppf 0.13
     105csh> ./Objs/calcpk lsscube_lobe.ppf lsspklobewn.ppf 0.13 3 
     106csh> ./Objs/calcpk lsscube_corlobe.ppf lsspkcorlobe.ppf 0.13 
    79107
    80108# 4.b/ Foreground power spectrum
    81109csh> ./Objs/calcpk fgndcube.ppf fgndpk.ppf 1000
    82110csh> ./Objs/calcpk fgndcube_lobe.ppf fgndpklobe.ppf 1000
     111csh> ./Objs/calcpk fgndcube_flobe.ppf fgndpkflobe.ppf 1000
    83112csh> ./Objs/calcpk fgndcube_corlobe.ppf fgndpkcorlobe.ppf 1000
    84113
    85114# 4.c/ Extract LSS P(k) from Foreground+LSS+noise , after cleaning/subtraction without beam
    86 csh> ./Objs/calcpk2 lsscube.ppf 0.2 fgndcube.ppf 1000 subpk.ppf 0.35 50. 0. 0. P2
     115csh> set beamdesc=repf11x11.ppf
     116csh> set ddishcor=50.
     117csh> set noiselev=1.
     118csh> ./Objs/calcpk2 lsscube.ppf 0.13 fgndcube.ppf 1000 subpk.ppf $noiselev $beamdesc 0. 0. P2
    87119# 4.d / Extract LSS P(k) from Foreground+LSS+noise and beam effect, without beam correction 
    88 csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.2 fgndcube_lobe.ppf 1000 subpklobe.ppf 0.35 50. 0. 3. P2
    89 # 4.e / Extract LSS P(k) from Foreground+LSS+noise and beam effect - correcting to a beam of Diam/Lambda = 130
    90 csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.2 fgndcube_lobe.ppf 1000 subpkcorlobe.ppf 0.35 50. 130. 3. P2
     120csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.13 fgndcube_lobe.ppf 1000 subpklobe.ppf $noiselev $beamdesc 0. 0. P2
     121# 4.e / Extract LSS P(k) from Foreground+LSS+noise and beam effect - correcting to a beam of Diam= $ddishcor
     122csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.13 fgndcube_lobe.ppf 1000 subpkcorlobe.ppf $noiselev $beamdesc $ddishcor 0. P2
    91123#  Or using a linear fit for foreground subtraction (old version)
    92 csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.2 fgndcube_lobe.ppf 1000 subpkcorlobep1.ppf 0.35 50. 130. 3. P1
     124csh> ./Objs/calcpk2 lsscube_lobe.ppf 0.13 fgndcube_lobe.ppf 1000 subpkcorlobep1.ppf $noiselev $beamdesc $ddishcor 0. P2
    93125# 4.f / Estimate residual noise from Foreground removal :
    94 csh> ./Objs/calcpk2 zerolss.ppf 0. fgndcube_lobe.ppf 1000 subpknolss.ppf 0.35 50. 130. 3. P2
    95 csh> ./Objs/calcpk2 zerolss.ppf 0. fgndcube_lobe.ppf 1000 subpknolssnocor.ppf 0.35 50. 0. 3.
     126csh> ./Objs/calcpk2 lsscube.ppf 0. fgndcube_lobe.ppf 1000 residcorlobe.ppf  $noiselev $beamdesc $ddishcor 0. P2
     127csh> ./Objs/calcpk2 lsscube.ppf 0. fgndcube_lobe.ppf 1000 residnocor.ppf $noiselev $beamdesc 0. 0. P2
    96128
    97129### Step 5 / Check the results using spiapp
    … …  
    100132openppf fgndpk.ppf
    101133openppf fgndpklobe.ppf
     134openppf fgndpkflobe.ppf
    102135openppf fgndpkcorlobe.ppf
    103136openppf lsspk.ppf
    104137openppf lsspklobe.ppf
     138openppf lsspkflobe.ppf
    105139openppf lsspklobewn.ppf
    106140openppf subpklobe.ppf
    … …  
    136170
    137171#  Calcul du volume total en Mpc^3
    138 set VOL 3*3*3*360*360*256
     172set VOL (1.9*1.9*2.8*800*600*256)
    139173plot2d lsspk x val*$VOL 1 'logx logy nsta xylimits=0.01,2.,10.,1e4 cpts marker=box,5 gold'
    140174plot2d lsspklobewn x val*$VOL 1 'same nsta cpts marker=box,5 siennared'
Note: See TracChangeset for help on using the changeset viewer.