fwhmdeg = 3./60 fm_stage4, dum, dum, lcenter, data, error = CMBspectra.fisher_analysis(pars, fsky, mukarcminT, fwhmdeg, basic_spectra=basic_spec, der=der, deltal=25, varnames=varnames, noT=noT, noE=noE, noTE=noTE) clf() a0,leg0 = FisherNew.plot_fisher(fm_stage4, pars,'k', varnames=varnames, onesigma=True) s,c = FisherNew.give_sigmas(fm_stage4, pars) ###### calcul plus réaliste des paramètres NET_T = 40 #muK.sqrt(s) Optimistic ? nbols = 10000 fsky = 1. fwhmdeg = 3./60 duration = 1. CMBspectra.muKarcmin(NET_T, nbols, fsky, fwhmdeg) ### explore fsky nbols = 10000 NET_T = 40. allfsky = linspace(0.01, 1,100) pars['lensing_residual'] = 0.1 fwhmdeg = 30./60 allsvals = np.zeros((len(pars), len(allfsky))) for i in np.arange(len(allfsky)): print(i) mukarcminT = CMBspectra.muKarcmin(NET_T, nbols, allfsky[i], fwhmdeg) fm, dum, dum, lcenter, data, error = CMBspectra.fisher_analysis(pars, allfsky[i], mukarcminT, fwhmdeg, basic_spectra=basic_spec, der=der, deltal=25, varnames=varnames, noT=noT, noE=noE, noTE=noTE) s, c = FisherNew.give_sigmas(fm, pars) allsvals[:,i] = s