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