Changeset 3756 in Sophya


Ignore:
Timestamp:
Apr 28, 2010, 3:41:15 PM (15 years ago)
Author:
ansari
Message:

Ajout des programmes de calcul de la sensibilite de l'interfero (plan(u,v), PNoise(k)) , Reza 28/04/2010

Location:
trunk/Cosmo/RadioBeam
Files:
7 added
3 edited

Legend:

Unmodified
Added
Removed

  • trunk/Cosmo/RadioBeam/README

    r3489 r3756  
    33  Simulation programs for radio lobe reconstruction, noise estimate etc ...
    44 
    5   2006-2008  R. Ansari (UPS/LAL-IN2P3)    C. Magneville (DAPNIA/SPP)
     5  2006-2010  R. Ansari (UPS/LAL-IN2P3)    C. Magneville (DAPNIA/SPP)
    66              (C) IN2P3/CNRS              (C)  DAPNIA/CEA
    77 ===============================================================================
    88
    9 
    10 A/ 3D reconstruction (alpha,delta,freq) simulation program
    11 ===========================================================
    12 
    139Note :  Uses SOPHYA class library  (http://www.sophya.org )
    1410
    15 List of files and classes :
    16 ---------------------------
    17 
    18 - treccyl.cc : main program (for tests)
    19 - multicyl.h .cc : Class MultiCylindres , 3D (freq,AngX,AngY) signal simulation
    20   and reconstruction. Uses class MultiBeamCyl
    21 - mbeamcyl.h .cc : Class MultiBeamCyl 2D (freq, AngX) signal simulation and reconstruction
    22 - brsrc.h .cc : Class BRSourceGen for managing list of sources
    23   (frequencies and angular position)
    24 - makefile
    25 
    26 - recbeam.pic : script with example commands for display with piapp
    27 
    28 
     11A/ (u-v) plane coverage and Noise-P(k) computation for different interferometer configuration
     12==============================================================================================
     13List of files :
     14  - specpk.h specpk.cc : 3D power spectrum computation classes 
     15  - mdish.h mdish.cc : Classes for representing Multi-Dish interferometer configurations
     16  - pknoise.cc : main program , makefile
     17  - plpkn.pic : spiapp script for plotting pknoise.cc program output
     18   
    2919B/ P(k) in temperature plot + "cosmic variance" + noise ...
    3020=============================================================
    … …  
    4030C.2/ pa.pic + pa.cc + pa.h : Lobe synthesis (dipoles sum)
    4131
     32D/ 3D reconstruction (alpha,delta,freq) simulation program
     33===========================================================
     34
     35List of files and classes :
     36---------------------------
     37
     38- treccyl.cc : main program (for tests)
     39- multicyl.h .cc : Class MultiCylindres , 3D (freq,AngX,AngY) signal simulation
     40  and reconstruction. Uses class MultiBeamCyl
     41- mbeamcyl.h .cc : Class MultiBeamCyl 2D (freq, AngX) signal simulation and reconstruction
     42- brsrc.h .cc : Class BRSourceGen for managing list of sources
     43  (frequencies and angular position)
     44- makerecyl
     45
     46- recbeam.pic : script with example commands for display with piapp

  • trunk/Cosmo/RadioBeam/anapkn.pic

    r3489 r3756  
    5454
    5555##  ---- VOL1,VOL2,VOL3
    56 ##  VOL1: 10deg x 10deg x Delta_z=0.2 : 335 x 335 x 680 Mpc^3 = 76 10^6 Mpc^3
    57 ##  VOL2: 20deg x 20deg x Delta_z=0.2 : VOL2 = 300 10^6 Mpc^3
    58 ##  VOL3: 60deg x 60deg x Delta_z=0.2 : VOL4 = 2700 10^6 Mpc^3
     56##  VOL0: 10deg x 10deg x Delta_z=0.2 : 335 x 335 x 680 Mpc^3 = 76 10^6 Mpc^3
     57##  VOL1: 20deg x 20deg x Delta_z=0.2 : VOL1 = 300 10^6 Mpc^3
     58##  VOL2: 360deg x 10deg x Delta_z=0.2 : VOL2 = 2700 10^6 Mpc^3
     59##  VOL3: 360deg x 30deg x Delta_z=0.2 : VOL3 = 10800 10^6 Mpc^3
    5960##  VOL4: 40deg x 40deg x Delta_z=0.2 : VOL3 = 1200 10^6 Mpc^3
    60 set angcov ( '10x10 deg^2' '20x20 deg^2' '60x60 deg^2' '40x40 deg^2' )
    61 VOL1 = 76e6
    62 VOL2 = 300e6
    63 VOL3 = 2700e6
     61set angcov ( '20x20 deg^2' '360x10 deg^2' '360x30 deg^2' '40x40 deg^2' )
     62VOL0 = 76e6
     63VOL1 = 300e6
     64VOL2 = 2700e6
     65VOL3 = 10800e6
    6466VOL4 = 1200e6
    6567
    … …  
    104106endscript
    105107
     108#  Setup 64 x D=5m sur 100mx100m , 10000 m^2
     109defscript setup5
     110  vpixz05 = 100
     111  vpix = 100
     112  kappa = 100*100/1200
     113  kappa = $kappa*$kappa
     114  echo "setup-5: VPix=$vpix  Kappa=$kappa"
     115endscript
    106116
    107117
    … …  
    112122
    113123#  --- Temps d'integration  a z=0.5 avec FOV=4x4 degres^2  (D=5m)
    114 #  120 jours pour VOl1 en ~ 7 passes (7 couvertures elementaires)
    115 #  120 jours pour VOl2 en ~ 4x7 passes (28 couvertures elementaires)
    116 #  120 jours pour VOl3 en ~ 16x7 passes (225 couvertures elementaires)
    117 set ndays ( '120 days' '120 days'  '120 days' '120 days' )
    118 TINT1 = 1.5e6
    119 TINT2 = 3.7e5
    120 TINT3 = 4.5e4
     124#  120 jours pour VOl0 en ~ 7 passes (7 couvertures elementaires)
     125#  360 jours pour VOl2 en ~ 16x7 passes (225 couvertures elementaires)
     126#  360 jours pour VOl1 en ~ 4x7 passes (28 couvertures elementaires)
     127#  360 jours pour VOl2 en ~ 90x3 passes (270 couvertures elementaires)
     128#  360 jours pour VOl3 en ~ 90x9 passes (810 couvertures elementaires)
     129
     130set ndays ( '360 days' '360 days'  '360 days' '360 days' )
     131TINT1 = 1.1e6
     132TINT2 = 1.1e5
     133TINT3 = 3.8e4
    121134TINT4 = 1.e5
    122135
    123136#  Bruit Tsys = 100 K = 10^5 mK
    124137TSYS = 1e5
     138#  Bruit Tsys = 50 K = 5.10^4 mK
     139TSYS = 5e4
    125140
    126141DELNU = 1.e6
    … …  
    151166set vols ( $VOL1 $VOL2 $VOL3 )
    152167set tints ( $TINT1 $TINT2 $TINT3 )
    153 set cols ( green blue red )
     168set cols ( red orange blue )
    154169
    155170#############################################################
    … …  
    239254  xyl = "xylimits=0.01,0.5,1.e-3,0.4 logx logy minorticks"
    240255  set vols ( $VOL1 $VOL2 $VOL3 )
    241   set cols ( green blue red )
     256  set cols ( red orange blue )
    242257  csigpk $vols[0]
    243258  n/plot hpkz.$CSPK/$kk%$kk ! ! "notit nsta connectpoints $cols[0] $xyl line=solid,2"

  • trunk/Cosmo/RadioBeam/makefile

    r3192 r3756  
     1############################################################
     2##  makefile for interferometer response
     3##  R.Ansari - Avril 2010 
     4############################################################
     5
    16include $(SOPHYABASE)/include/sophyamake.inc
    27
    3 all : treccyl
     8all : pknoise 
    49
    5 clean :
    6         rm Objs/*
    7 depend :
    8         mkdir Objs
     10clean :
     11        rm Objs/*
    912
    10 treccyl : Objs/treccyl
    11         echo 'makefile : treccyl made'
     13pknoise : Objs/pknoise
     14        echo 'makefile : pknoise made'
    1215
    13 testFFT : Objs/testFFT.o
    14 ######
     16Objs/pknoise : Objs/pknoise.o Objs/specpk.o Objs/mdish.o
     17        $(CXXLINK) -o Objs/pknoise Objs/pknoise.o Objs/specpk.o Objs/mdish.o $(SOPHYAEXTSLBLIST)
    1518
    16 Objs/treccyl : Objs/treccyl.o Objs/multicyl.o Objs/mbeamcyl.o  Objs/brsrc.o
    17         $(CXXLINK) -o Objs/treccyl \
    18         Objs/treccyl.o Objs/multicyl.o Objs/mbeamcyl.o Objs/brsrc.o \
    19         $(SOPHYAEXTSLBLIST)
     19Objs/pknoise.o : pknoise.cc specpk.h  mdish.h
     20        $(CXXCOMPILE) -o Objs/pknoise.o pknoise.cc
    2021
    21 Objs/treccyl.o : treccyl.cc  multicyl.h mbeamcyl.h brsrc.h
    22         $(CXXCOMPILE) -o Objs/treccyl.o  treccyl.cc
     22Objs/specpk.o : specpk.cc specpk.h  mdish.h
     23        $(CXXCOMPILE) -o Objs/specpk.o specpk.cc
    2324
    24 Objs/brsrc.o : brsrc.cc  brsrc.h
    25         $(CXXCOMPILE) -o Objs/brsrc.o  brsrc.cc
    26 
    27 Objs/mbeamcyl.o : mbeamcyl.cc  mbeamcyl.h brsrc.h
    28         $(CXXCOMPILE) -o Objs/mbeamcyl.o  mbeamcyl.cc
    29 
    30 Objs/multicyl.o : multicyl.cc  multicyl.h mbeamcyl.h brsrc.h
    31         $(CXXCOMPILE) -o Objs/multicyl.o  multicyl.cc
    32 
    33 testFFT :  Objs/testFFT.o
    34         $(CXXLINK) -o testFFT \
    35         Objs/testFFT.o \
    36         $(SOPHYAEXTSLBLIST)
    37 
    38 Objs/testFFT.o : testFFT.cc
    39         $(CXXCOMPILE) -o Objs/testFFT.o testFFT.cc
     25Objs/mdish.o : mdish.cc mdish.h
     26        $(CXXCOMPILE) -o Objs/mdish.o mdish.cc
Note: See TracChangeset for help on using the changeset viewer.