Skip to content

GitLab

Commit 33aa5123 authored by Ignacio Magana's avatar Ignacio Magana 💬
Browse files
Options

Merge branch 'O3a' into 'O3a' 

Analysis breakdown

See merge request lscsoft/gwcosmo!45
parents a9185a0f 2816c188
Pipeline #127586 passed with stage
    in 19 seconds
    Hide whitespace changes
    Inline Side-by-side
    Showing with 3 additions and 1 deletion
    gwcosmo/gwcosmo.py
    View file @ 33aa5123
    ...@@ -720,12 +720,14 @@ class gwcosmoLikelihood(object): ...@@ -720,12 +720,14 @@ class gwcosmoLikelihood(object):
    likelihood = self.pGD*(pxG/self.pDG) + self.pnGD*(pxnG/self.pDnG) # Eq 3 likelihood = self.pGD*(pxG/self.pDG) + self.pnGD*(pxnG/self.pDnG) # Eq 3
    np.savez('likelihood_breakdown_'+str(self.Lambda)+'.npz',[H0,likelihood, pxG,self.pDG,self.pGD, pxnG,self.pDnG,self.pnGD])
    if self.whole_cat == False: if self.whole_cat == False:
    pDnG_rest_of_sky = self.pD_H0nG(H0,allsky=False) pDnG_rest_of_sky = self.pD_H0nG(H0,allsky=False)
    pxnG_rest_of_sky = self.px_H0nG(H0,allsky=False) pxnG_rest_of_sky = self.px_H0nG(H0,allsky=False)
    likelihood = likelihood*self.catalog_fraction + (pxnG_rest_of_sky/pDnG_rest_of_sky)*self.rest_fraction # Eq 4 likelihood = likelihood*self.catalog_fraction + (pxnG_rest_of_sky/pDnG_rest_of_sky)*self.rest_fraction # Eq 4
    np.savez('likelihood_breakdown_'+str(self.Lambda)+'.npz',[H0,likelihood, pxG,self.pDG,self.pGD*self.catalog_fraction, pxnG,self.pDnG,self.pnGD*self.catalog_fraction, pxnG_rest_of_sky,pDnG_rest_of_sky,self.rest_fraction])
    return likelihood return likelihood
    ......
    Markdown is supported
    0% or .
    You are about to add 0 people to the discussion. Proceed with caution.
    Finish editing this message first!
    Please register or to comment