Sphere vs. patches

Comparaison reconstruction directe et reconstruction par patches de la sphère

L'objectif est de comparer les performances d'une reconstruction des lentilles par un estimateur agissant directement sur la sphere à une méthode en deux etapes, consistant a decouper la sphere en patches puis appliquer sur les patches un estimateur "flat sky".

Tests et harmonisation des outils de simulation/reconstruction par patches

Bilan des outils disponibles pour la simulation et reconstruction sur les patches :

les outils de simulation :

  • [LP] package IDL :
    • code : td_sqmap_creator.pro dans ${RMPDIR}/rmp/v1/patch/idl/simu
    • commande : dans idl faire : td_sqmap_creator,paramfile,niter
    • exemple de fichier de parametres : /sps/planck/Users/perotto/lensdata/alexisse/alexis_simuparams.txt
  • [AL] package C++ :
    • code : LensedMapGenerator.cc dans ${RMPDIR}/rmp/v1/proj/src/exec
    • A COMPLETER
  • les différences identifiées
    • [AL] interpolation closest ou bilinear, [LP] interpolation bicubic
    • generer la carte de T non-lentillé jusqu'à lmax=2*lmax

les outils de reconstruction

  • [LP] : reconstruction en 3 etapes :
    1. estimation du spectre de puissance des cartes lentillees (partie IDL). 3 methodes :
      • raccoursis : on donne le spectre theorique (e.g. issu de CAMB)
      • estimation carte à carte :
        • Un spectre estimé pour chacune des cartes.
        • Méthode basée sur la routine map2cbin.pro.
        • Dans idl, lancer analyse_obsmaps,paramfile,niter, /interpol
      • estimation globale d'un spectre pour toutes les cartes
        • en particulier, preparation des cartes de c(k) pour la méthode d'estimation de Steph (MU=ZUP)
        • dans idl, pour calculer la carte de c(k) des cartes lentillées, lancer prepare_ckmap, paramfile, niter, /dolensed
    2. estimation du champ de deflexion non-normalisé : 2 méthodes
      • Avec la MU :
        • code : TT_MUestimator.cc dans ${RMPDIR}/rmp/v1/patch/src
        • commande à lancer : ./TT_MUestimator params.txt ${numero_de_la_carte}
        • exemple de fichier de params : /sps/planck/Users/perotto/lensdata/alexisse/alexis_MUsimuparams.txt
        • exemple de script :
      • Sans la MU (clTT estimé carte à carte ou ClTT theorique) :
        • code : TT_estimator.cc dans ${RMPDIR}/rmp/v1/patch/src
        • commande à lancer : ./TT_estimator params.txt ${numero_de_la_carte}
        • exemple de fichier de params : /sps/planck/Users/perotto/lensdata/alexisse/alexis_simuparams.txt
        • exemple de script : /afs/in2p3.fr/group/planck/perotto/SL4/rmp/v1/patch/scripts/laurence_estidtt.bash
        • NB : ~10sec de CPU consommé : lancer par paquet de 50 iterations sur anastasie !
    3. calcul de la normalisation de l'estimateur : 2 méthodes
      • Avec la MU :
        • code : TT_MUnormalisor.cc dans ${RMPDIR}/rmp/v1/patch/src
        • commande à lancer : ./TT_MUnormalisor params.txt
        • exemple de fichier de params : /sps/planck/Users/perotto/lensdata/alexisse/alexis_MUsimuparams.txt
        • NB : ne lancer qu'une fois pour toutes les cartes
      • Sans la MU (clTT estimé carte à carte ou ClTT theorique) :
        • code : TT_normalisor.cc dans ${RMPDIR}/rmp/v1/patch/src
        • commande à lancer : ./TT_normalisor params.txt ${numero_de_la_carte}
        • exemple de fichier de params : /sps/planck/Users/perotto/lensdata/alexisse/alexis_simuparams.txt
        • NB : ~10sec de CPU consommé : lancer par paquet de 50 iterations sur anastasie !
        • NB : si ClTT theorique, ne lancer qu'une fois ./TT_normalisor params.txt
  • [AL] : A COMPLETER
    1. estimation du spectre de puissance des cartes lentillees
      • raccoursis : on donne le spectre theorique (e.g. issu de CAMB)
      • estimation globale d'un spectre pour toutes les cartes :
        • code : PatchPSEstimator.cc
        • basé sur la classe Power_Spectrum_Accumulator , méthode addPowerSpectrumOf
        • exemple de fichier de params :
        • exemple de script :
    2. estimation du champ de deflexion :
      • code : PatchLensingPotentialEstimator.cc dans ${RMPDIR}/rmp/v1/proj/exec
      • basé sur la classe d_estimator , méthode estimate
      • NB : la normalisation est recalculée à chaque fois par l'appel à getNormalisation
      • commande à lancer : ./PatchLensingPotentialEstimator paramFile [index_begin] [index_end]
      • exemple de fichier de params :
      • exemple de script :
    3. calcul de la normalisation de l'estimateur :
      • code : PatchLensingNL0Estimator.cc dans ${RMPDIR}/rmp/v1/proj/exec
      • basé sur la classe Lensing_Noise_Estimator, méthode calculate
      • commande à lancer : PatchLensingNL0Estimator paramFile
      • exemple de fichier de params :

tests

  • avec codes de simulation [LP], ont ete genere 768 patches dans /sps/planck/Users/perotto/lensdata/alexisse/ :
    • T non-lentille : TT_lp_pix350_ang0.175_lmax3000_*.fits
    • T lentille : LTT_lp_pix350_ang0.175_lmax3000_*.fits
    • Phi : PP_lp_pix700_ang0.175_lmax6000_*.fits
    • NB : la simulation de l'effet de lentille est correcte jusqu'a l=3000 (ensuite probable imprecision)
  • paramètres de l'estimation :
    • lmax_deflexion = 2600
    • ClTT et lensClTT théoriques
  • Estimation avec codes [LP]
    • code utilisé : TT_estimateur
    • deflexion estimée (non-normalisée) : bias_d_TT_lp_pix350_ang0.175_lmax2600_*.fits
    • normalisation : A_TT_lp_pix350_ang0.175_lmax2600.fits
  • Estimation avec codes [AL]
    • code utilisé :
    • deflexion estimée : RPP_al3_pix350_ang0.175_lmax3000_FDI_*.fits
  • plots
    • voir le script idl : ${RMPDIR}/rmp/v1/patch/idl/scripts/compare_sphere_patch.scr
    • NB : AL_cldd * k^2 * angsize^2

Tests de la méthode sphere by patch

description de la méthode

  • découpage de la spère en patches
    • code : SphereToPatch.cc dans ${RMPDIR}/rmp/v1/proj/exec
    • exemple de fichier de parametres :

tests

  • découpage de lensed_map-11187
    • output : 768 patches en Fourier : FDI_lensed_map-11187_Toep_pix200_ang0.175_*.fits
    • dir : /sps/planck/Users/lavabre/laurence/
  • reconstruction sur les patches
    • lmax_deflexion = 2600
    • ClTT et lensClTT théoriques
  • Estimation avec codes [LP]
    • code utilisé : TT_estimateur
    • deflexion estimée (non-normalisée) : bias_d_TT_lp2_pix350_ang0.175_lmax2600_*.fits
    • normalisation : A_TT_lp_pix350_ang0.175_lmax2600.fits
  • Estimation avec codes [AL]
    • code utilisé :
    • deflexion estimée : RPP_al_FDI_lensed_map-11187_Toep_pix200_ang0.175_*.fits
  • plots
    • voir le script idl : ${RMPDIR}/rmp/v1/patch/idl/tests/compare_sphere_patch.scr
    • voir fichier de parametres : /sps/planck/Users/lavabre/laurence/laurence_params.txt
    • NB : AL_cldd * k^2 * angsize^2
    • NB : les deux cldd reconstruits sont normalisés par angsize^5 :
      • les cartes sont a normaliser par angsize dans le sens : cl^TT(CAMB) = c^TT(estimé) * angsize^2
      • ce qui explique le facteur angsize^4 dans le cldd(estimé)
      • reste un facteur angsize supplementaire que je ne comprends pas.