コード例 #1
0
ファイル: LogErrorTest.py プロジェクト: mscrim/LogErrorTest
def zeroVelocityMocks():

    sigmastar2 = 250.0 ** 2

    # -------------------- Read in 6dFGSv data ----------------------

    Dz, xerr, x_arr, y_arr, z_arr, xobs, r_arr, ngal = readin6dFGSvData()

    # Calculate velocity uncertainties sigma_v from Eq. (10) of Scrimgeour et al. (2015)
    sigma_v = 0.324 * 100.0 * Dz

    # Calculate Dec (np array)
    Dec = getDec(x_arr, y_arr, z_arr)
    # Index of galaxies in Great Circle.
    inC = np.where(Dec > -20.0)[0]

    # ---------------------- Convert Dz to z -------------------------

    zfile = "/Users/Morag/Cosmology_codes/redshift_OM0.3175.dat"
    zo = cc.z_x_arr(Dz, zfile)

    # -------- Assign each galaxy zero x, perturbed by xerr ----------
    # Give all galaxies v=0 but perturb x by xerr
    # Generate array of num random deltax values with std dev = xerr[gal]

    BFArr = []

    for i in xrange(100):
        print i

        x = xerr * np.random.standard_normal(ngal)

        ## Calibrate zeropoint to 0
        # Mean x in Great Circle
        # meanx = np.mean(x[inC])
        # Subtract this from all galaxies
        # x -= meanx

        # Calculate derived Dr
        Dr_der = Dz / 10 ** x

        # Calculate new derived redshifts
        zr_der = cc.z_x_arr(Dr_der, zfile)

        vp_rand = speedlight * ((zo - zr_der) / (1.0 + zr_der))

        # ------------------------- Calculate BF ---------------------------

        # u = MLE_BF(x_arr, y_arr, z_arr, r_arr, vp_rand, sigma_v, sigmastar2)
        u = MV_BF_6dFGSv(vp_rand)
        BFArr.append((u[0], u[1], u[2]))

    print ""

    meanBF, stdBF = meanVector(BFArr)
    print "Mean spurious BF: (%d pm %d, %d pm %d, %d pm %d)" % (
        meanBF[0],
        stdBF[0],
        meanBF[1],
        stdBF[1],
        meanBF[2],
        stdBF[2],
    )

    return
コード例 #2
0
def x_to_v(x, z_obs, Dz, zfile):
    Dr_der = Dz / 10.**x
    zr_der = cc.z_x_arr(Dr_der, zfile)
    v = speedlight * ((z_obs - zr_der)/(1.+zr_der))
    return v