In [1]:

    

import numpy as np
import os, sys
path = os.path.abspath('../library/')
if path not in sys.path:
    sys.path.append(path)
from cumhpx import do_likelihood, local_normalisation
name = 'raster_0.7_1.2'


    

In [ ]:

    

print(local_normalisation(name, verbose = True))
outputDir = '../output/%s_0.001' %name
fSample = 0.001
print(do_likelihood(outputDir,fSample= fSample))


    

    




(25.910147711235833, 4.467184215036124, 313826)
GDR2mock
nstar before cleaning: 1506672
nstar after g-band cleaning: 1161038
nstar before cleaning: 1161038
nstar after rp-band cleaning: 844364
844364
after smc cut: 844040
after lmc cut: 842305
GDR2
nstar before cleaning: 1253908
nstar after g-band cleaning: 1106539
nstar before cleaning: 1106539
nstar after rp-band cleaning: 930786
930786
after smc cut: 929405
after lmc cut: 919862
hpx 2055 is empty
hpx 2064 is empty
hpx 2065 is empty
hpx 2066 is empty
hpx 2067 is empty
hpx 2068 is empty
hpx 2069 is empty
hpx 2070 is empty
hpx 2071 is empty
hpx 2072 is empty
hpx 2073 is empty
hpx 2080 is empty
hpx 2081 is empty
hpx 2082 is empty
hpx 2088 is empty
hpx 2112 is empty
hpx 2336 is empty
hpx 2338 is empty
hpx 2344 is empty
hpx 2346 is empty
likelihood:  -328108.99258733296
total diff:  131017
too many: 26730, too little: -104287
hpx 2055 is empty
hpx 2064 is empty
hpx 2065 is empty
hpx 2066 is empty
hpx 2067 is empty
hpx 2068 is empty
hpx 2069 is empty
hpx 2070 is empty
hpx 2071 is empty
hpx 2072 is empty
hpx 2073 is empty
hpx 2080 is empty
hpx 2081 is empty
hpx 2082 is empty
hpx 2088 is empty
hpx 2112 is empty
hpx 2336 is empty
hpx 2338 is empty
hpx 2344 is empty
hpx 2346 is empty
hpx 2055 is empty
hpx 2064 is empty
hpx 2065 is empty
hpx 2066 is empty
hpx 2067 is empty
hpx 2068 is empty
hpx 2069 is empty
hpx 2070 is empty
hpx 2071 is empty
hpx 2072 is empty
hpx 2073 is empty
hpx 2080 is empty
hpx 2081 is empty
hpx 2082 is empty
hpx 2088 is empty
hpx 2112 is empty
hpx 2336 is empty
hpx 2338 is empty
hpx 2344 is empty
hpx 2346 is empty
842305 919862






    




/home/rybizki/anaconda3/lib/python3.6/site-packages/healpy/pixelfunc.py:304: RuntimeWarning: invalid value encountered in less_equal
  return np.absolute(m - badval) <= atol + rtol * np.absolute(badval)
/home/rybizki/anaconda3/lib/python3.6/site-packages/healpy/projaxes.py:1053: RuntimeWarning: invalid value encountered in less
  result.data[result.data<0]=0.0
/home/rybizki/anaconda3/lib/python3.6/site-packages/healpy/projaxes.py:1054: RuntimeWarning: invalid value encountered in greater
  result.data[result.data>1]=1.0






    




(842305, -328108.99258733296, 26730, 104287)






    




/home/rybizki/anaconda3/lib/python3.6/site-packages/matplotlib/colors.py:1031: RuntimeWarning: invalid value encountered in less_equal
  mask |= resdat <= 0

In [ ]:

    

name_comp = 'raster_0.7_1.0'
compDir = '../output/%s_0.001/GDR2mock_20.7Gmag.fits' %name_comp
fSample = 0.0025
do_likelihood(outputDir,fSample= fSample, comp_file = compDir, savefile=False)


    

    




GDR2mock

In [ ]: