Skip to content
Snippets Groups Projects
Commit 9583ce3a authored by Kipp Cannon's avatar Kipp Cannon
Browse files

gstlal-inspiral: add O2 hyperbank mass model

- switch ranking statistic to this unconditionally for the October 2018
  open public alert test run.  after the test run and before O3 we will
  introduce more user-friendly plumbing allowing this to work with other
  template banks.
- NOTE:  git-lfs is used for any file ending in .hdf5 in the
  population_models/ tree
parent dbbf803b
No related branches found
No related tags found
No related merge requests found
......@@ -72,6 +72,10 @@ __all__ = [
class UniformInTemplatePopulationModel(object):
# FIXME: this needs to learn the template chirp masses so it can
# take out the factor of mchirp**(5./3.) factor on which the
# horizon distance depends (which appears in a volume factor in the
# numerator of the likelihood ratio)
def __init__(self, template_ids):
"""
Assumes uniform in template population model, no
......@@ -87,7 +91,8 @@ class UniformInTemplatePopulationModel(object):
class SourcePopulationModel(object):
DEFAULT_FILENAME = os.path.join(gstlal_config_paths["pkgdatadir"], "lnP_template_signal_BNS_gaussian_lowspin_Ozel.hdf5")
POPULATION_MODELS_PATH = os.path.join(gstlal_config_paths["pkgdatadir"], "population_models")
DEFAULT_FILENAME = os.path.join(POPULATION_MODELS_PATH, "O2/lnP_template_signal_BBH_logm_reweighted_mchirp.hdf5")
def __init__(self, template_ids, filename = None):
......
......@@ -261,14 +261,8 @@ class LnSignalDensity(LnLRDensity):
self.horizon_history = horizonhistory.HorizonHistories((instrument, horizonhistory.NearestLeafTree()) for instrument in self.instruments)
# source population model
if self.template_ids:
self.population_model = inspiral_intrinsics.UniformInTemplatePopulationModel(self.template_ids)
# FIXME: switch to this when a model file becomes
# available
#self.population_model = inspiral_intrinsics.SourcePopulationModel(self.template_ids)
else:
# default lnP = 1/len(templates) = 0
self.population_model = inspiral_intrinsics.UniformInTemplatePopulationModel([0])
# FIXME: introduce a mechanism for selecting the file
self.population_model = inspiral_intrinsics.SourcePopulationModel(self.template_ids)
self.InspiralExtrinsics = inspiral_extrinsics.InspiralExtrinsics(self.min_instruments)
......
SUBDIRS = O2
population_modelsdir = $(pkgdatadir)/$(subdir:share/%=%)
population_modelsdir = $(pkgdatadir)/$(subdir:share/%=%)
population_models_DATA = \
lnP_template_signal_BBH_logm_reweighted_mchirp.hdf5
File added
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment