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
Subscribe to:
Posts (Atom)
