ma non va bene fuori dal piano galattico. Ho covergenza solo sul piano e su pochi altri patches sparsi.
Non so se sto sbagliando qualcosa o se la cosa sia possibile, chiedero' bene ad Anna ma sembra che CCA non abbia abbastanza informazioni.
Se fisso l'indice del synchrotrone a 2.8 ho piu' covergenze ma non ho fatto un run completo per cui non so bene. Lo faro' sta sera. Comunque non vorrei fissare cosi' categoricamente questo numero.
Ricapitolando: facendo l'analisi in polarizzazione ho solo 7 canali (mancano i 2 ad alta frequenza) e in piu' per ora non sto introducendo nessuna prior perche' non credo che abbiamo una mappa credibile in polarizzazione per dust e synch. La cosa positiva e' che invece abbiamo solo 3 componenti (synch dust cmb) invece di 4 (synch ff dust cmb).
Quindi rispetto alla challenge in T riduciamo il numero di canali informativi (9ch+2prior) a 7(ch) e diminuiamo le componenti da stimare da 4 a 3.
Questo e' un motivo sufficiente a farmi diminuire i patches su cui ho convergenza o c'e qualcosa che non va?
Un'altra domanda veramente stupida e' che non mi ricordo se normalmente applichiamo la maschera solo ai canali per fare l'inversione o la maschera serve anche nella stima della matrice di mixing. In principio servirebbe da subito ma in polarizzazione per questa simulazione non ho sorgenti per cui la maschera forse non serve. vabbuoz
Wednesday, July 11, 2007
Monday, July 9, 2007
july 9
Aggiornamenti e questioni:
Julian ha runnato di nuovo le mappe per generare anche le relative mappe di rumore. Ore i canali hanno 9 fields che rappresentano le 3 mappe piu' i 6 elementi della matrice di mixing che determinano completamente le caratterisiche dell rumore (matrice diagonale). Io uso gli elementi diagonali (sigma^2) e costruisco le mappe del rumore per ogni canale.
Fatto questo le mappe del rumore vanno trattate come le mappe stesse (smoothing per il beam che manca a 33.+ cambio di coordinate). Questo e' stato fatto.
Poi ho runnato CCA con la configurazione dell'altra volta sembra funzionare bene.Parametri incogniti solo indice spettrale del synch e quello della dust mentre la T_dust e' fissata a 18K. Ho provato a fare un run variando anche la temperatura della polvere ma non ce l'ha fatta. Non converge. Adesso faro' un run completo su tutti i patches con la configurazione che ho visto essere abbastanza stabile e vediamo.
Nello stesso tempo sto cercando di capire come interfacciarmi con loro che devono runnare madmap.
Loro hanno bisogno di una matrice di mixing per ogni pixel. Per i pixel che non sono nell'intersezione delle pixel list dei vari patches non c'e' problema, ma per i pixel nelle intersezione ho vari dubbi. Non so se e' corretto fare la media delle matrici di mixing perche' ho paura di perdere la normalizzazione. Invece potrei calcolare per ogni pixel i valori corretti per i due parametri che stimo da CCA e poi calcolare la matrice di mixing per ogni pixel.
L'unico dubbio (probabilmente idiota) e' che questi numeri stanno all'esponente, tipo e^-(beta_synch) per cui non mi piace tanto l'idea di mediare i beta synch tra loro nelle varie regioni... ma forse e' solo una mia preoccupazione inutile.
Ancora non riesco a trattare bene le intersezioni delle pixel list. Volevo fare un plot come quello che ha fatto Anna nella figura 2 del paper ma non ci riesco :(
Per ora e' tutto
s
Julian ha runnato di nuovo le mappe per generare anche le relative mappe di rumore. Ore i canali hanno 9 fields che rappresentano le 3 mappe piu' i 6 elementi della matrice di mixing che determinano completamente le caratterisiche dell rumore (matrice diagonale). Io uso gli elementi diagonali (sigma^2) e costruisco le mappe del rumore per ogni canale.
Fatto questo le mappe del rumore vanno trattate come le mappe stesse (smoothing per il beam che manca a 33.+ cambio di coordinate). Questo e' stato fatto.
Poi ho runnato CCA con la configurazione dell'altra volta sembra funzionare bene.Parametri incogniti solo indice spettrale del synch e quello della dust mentre la T_dust e' fissata a 18K. Ho provato a fare un run variando anche la temperatura della polvere ma non ce l'ha fatta. Non converge. Adesso faro' un run completo su tutti i patches con la configurazione che ho visto essere abbastanza stabile e vediamo.
Nello stesso tempo sto cercando di capire come interfacciarmi con loro che devono runnare madmap.
Loro hanno bisogno di una matrice di mixing per ogni pixel. Per i pixel che non sono nell'intersezione delle pixel list dei vari patches non c'e' problema, ma per i pixel nelle intersezione ho vari dubbi. Non so se e' corretto fare la media delle matrici di mixing perche' ho paura di perdere la normalizzazione. Invece potrei calcolare per ogni pixel i valori corretti per i due parametri che stimo da CCA e poi calcolare la matrice di mixing per ogni pixel.
L'unico dubbio (probabilmente idiota) e' che questi numeri stanno all'esponente, tipo e^-(beta_synch) per cui non mi piace tanto l'idea di mediare i beta synch tra loro nelle varie regioni... ma forse e' solo una mia preoccupazione inutile.
Ancora non riesco a trattare bene le intersezioni delle pixel list. Volevo fare un plot come quello che ha fatto Anna nella figura 2 del paper ma non ci riesco :(
Per ora e' tutto
s
Thursday, June 28, 2007
bandwidht simulation (top hat)
100-1a 100.0e9 83.3e9 116.7e9
100-1b 100.0e9 83.3e9 116.7e9
100-2a 100.0e9 83.3e9 116.7e9
100-2b 100.0e9 83.3e9 116.7e9
100-3a 100.0e9 83.3e9 116.7e9
100-3b 100.0e9 83.3e9 116.7e9
100-4a 100.0e9 83.3e9 116.7e9
100-4b 100.0e9 83.3e9 116.7e9
143-1a 143.0e9 119.2e9 166.8e9
143-1b 143.0e9 119.2e9 166.8e9
143-2a 143.0e9 119.2e9 166.8e9
143-2b 143.0e9 119.2e9 166.8e9
143-3a 143.0e9 119.2e9 166.8e9
143-3b 143.0e9 119.2e9 166.8e9
143-4a 143.0e9 119.2e9 166.8e9
143-4b 143.0e9 119.2e9 166.8e9
143-5 143.0e9 119.2e9 166.8e9
143-6 143.0e9 119.2e9 166.8e9
143-7 143.0e9 119.2e9 166.8e9
143-8 143.0e9 119.2e9 166.8e9
217-1 217.0e9 180.8e9 253.2e9
217-2 217.0e9 180.8e9 253.2e9
217-3 217.0e9 180.8e9 253.2e9
217-4 217.0e9 180.8e9 253.2e9
217-5a 217.0e9 180.8e9 253.2e9
217-5b 217.0e9 180.8e9 253.2e9
217-6a 217.0e9 180.8e9 253.2e9
217-6b 217.0e9 180.8e9 253.2e9
217-7a 217.0e9 180.8e9 253.2e9
217-7b 217.0e9 180.8e9 253.2e9
217-8a 217.0e9 180.8e9 253.2e9
217-8b 217.0e9 180.8e9 253.2e9
353-1 353.0e9 294.2e9 411.8e9
353-2 353.0e9 294.2e9 411.8e9
353-3a 353.0e9 294.2e9 411.8e9
353-3b 353.0e9 294.2e9 411.8e9
353-4a 353.0e9 294.2e9 411.8e9
353-4b 353.0e9 294.2e9 411.8e9
353-5a 353.0e9 294.2e9 411.8e9
353-5b 353.0e9 294.2e9 411.8e9
353-6a 353.0e9 294.2e9 411.8e9
353-6b 353.0e9 294.2e9 411.8e9
353-7 353.0e9 294.2e9 411.8e9
353-8 353.0e9 294.2e9 411.8e9
545-1 545.0e9 454.3e9 635.8e9
545-2 545.0e9 454.3e9 635.8e9
545-3 545.0e9 454.3e9 635.8e9
545-4 545.0e9 454.3e9 635.8e9
857-1 857.0e9 714.2e9 999.9e9
857-2 857.0e9 714.2e9 999.9e9
857-3 857.0e9 714.2e9 999.9e9
857-4 857.0e9 714.2e9 999.9e9
LFI-18a 70.0e9 63.0e9 77.0e9
LFI-18b 70.0e9 63.0e9 77.0e9
LFI-19a 70.0e9 63.0e9 77.0e9
LFI-19b 70.0e9 63.0e9 77.0e9
LFI-20a 70.0e9 63.0e9 77.0e9
LFI-20b 70.0e9 63.0e9 77.0e9
LFI-21a 70.0e9 63.0e9 77.0e9
LFI-21b 70.0e9 63.0e9 77.0e9
LFI-22a 70.0e9 63.0e9 77.0e9
LFI-22b 70.0e9 63.0e9 77.0e9
LFI-23a 70.0e9 63.0e9 77.0e9
LFI-23b 70.0e9 63.0e9 77.0e9
LFI-24a 44.0e9 39.6e9 48.0e9
LFI-24b 44.0e9 39.6e9 48.0e9
LFI-25a 44.0e9 39.6e9 48.0e9
LFI-25b 44.0e9 39.6e9 48.0e9
LFI-26a 44.0e9 39.6e9 48.0e9
LFI-26b 44.0e9 39.6e9 48.0e9
LFI-27a 30.0e9 27.0e9 33.0e9
LFI-27b 30.0e9 27.0e9 33.0e9
LFI-28a 30.0e9 27.0e9 33.0e9
LFI-28b 30.0e9 27.0e9 33.0e9
June 28 11:35am
Ieri ho importato le mappe da seaborg sul mio laptop. Ho cominciato a scrivere il programma che le doveva portare a bassa risoluzione (30') come il canale a 33GHz.
Le mappe fatte da springtide per ora sono in coordinate Ecliptiche. Potremmo trovarci a dover guardare mappe in coordinate diverse dalle galattiche per cui un possibile improuvement per CCA potrebbe essere quello di lavorare con mappe in qualsiasi tipo di coordinate. Siccome CCA estae patches a partire da coordinate galattiche sipotrebbe mettere un if iniziale e nel caso coord ne 'G' fa la conversione.
Penso che sia piu' facile e corretto passare per gli alm. Dopo averli calcolati si riproietta la mappa nelle nuove coordinate. Non e' una procedura perfetta perche' bisognerebbe rifare la mappa da capo col mapmaking ma questo non e' sempre possibile.
Questi sarebbero gli steps:
calcolare gli alm con synfast o con map2alm
call coordsys2euler_zyz(2000.0_dp, 2000.0_dp, 'E', 'G', psi, theta, phi)
call rotate_alm(nside, alm_TGC, psi, theta, phi)
Oggi sto facendo la rotazione scrivendo un programmino IDL che mi cambia le coordinate in coda al programma che abbassa l;a risoluzione. Federico mi ha detto come lo fa di solito lui e sto tentando di riprodurre la sua strategia.
Le mappe fatte da springtide per ora sono in coordinate Ecliptiche. Potremmo trovarci a dover guardare mappe in coordinate diverse dalle galattiche per cui un possibile improuvement per CCA potrebbe essere quello di lavorare con mappe in qualsiasi tipo di coordinate. Siccome CCA estae patches a partire da coordinate galattiche sipotrebbe mettere un if iniziale e nel caso coord ne 'G' fa la conversione.
Penso che sia piu' facile e corretto passare per gli alm. Dopo averli calcolati si riproietta la mappa nelle nuove coordinate. Non e' una procedura perfetta perche' bisognerebbe rifare la mappa da capo col mapmaking ma questo non e' sempre possibile.
Questi sarebbero gli steps:
calcolare gli alm con synfast o con map2alm
call coordsys2euler_zyz(2000.0_dp, 2000.0_dp, 'E', 'G', psi, theta, phi)
call rotate_alm(nside, alm_TGC, psi, theta, phi)
Oggi sto facendo la rotazione scrivendo un programmino IDL che mi cambia le coordinate in coda al programma che abbassa l;a risoluzione. Federico mi ha detto come lo fa di solito lui e sto tentando di riprodurre la sua strategia.
weblog!
Ciao a tutti,
ho pensato che forse potrebbe essere comodo mettere tutte le info relative a questi miei tre mesi a Berkeley su un weblog cosi' io non mi perdo niente e tutte le persone interessate possono aggiornarsi sul mio lavoro.
E speriamo che va tutto bene.
e ora... al lavoro!
s
ho pensato che forse potrebbe essere comodo mettere tutte le info relative a questi miei tre mesi a Berkeley su un weblog cosi' io non mi perdo niente e tutte le persone interessate possono aggiornarsi sul mio lavoro.
E speriamo che va tutto bene.
e ora... al lavoro!
s
Subscribe to:
Posts (Atom)
