put back the neutron star options

bin/manifold_cbc_bank.py

@@ -68,11 +68,11 @@ def parse_args(args=None):
 	#p.add_option("", "--max-M", type="float", default=numpy.inf, help="maximum total mass. default inf")
 	p.add_option("", "--min-q", type="float",  default=1., help="minimum q. default 1")
 	p.add_option("", "--max-q", type="float",  default=100., help="maximum q. default 100")
-	#p.add_option("", "--min-ns-mass", type=float, default=0.1, help="minimum mass to consider an NS. default 0.1")
-	#p.add_option("", "--max-ns-mass", type=float, default=0.1, help="maximum mass to consider an NS. default 0.1")
+	p.add_option("", "--min-ns-mass", type=float, default=0.0, help="minimum mass to consider an NS. default 0")
+	p.add_option("", "--max-ns-mass", type=float, default=3.0, help="maximum mass to consider an NS. default 3.0")
 	#p.add_option("", "--min-ns-s1z", type=float, default=-0.05, help="minimum spin to allow for an NS. default -0.05")
 	#p.add_option("", "--max-ns-s1z", type=float, default=0.05, help="maximum spin to allow for an NS. default 0.05")
-	#p.add_option("", "--max-ns-abs-s1z", type=float, default=0.05, help="maximum abs z spin to allow for an NS. default 0.05")
+	p.add_option("", "--max-ns-abs-s1z", type=float, default=0.05, help="maximum abs z spin to allow for an NS. default 0.05")
 	p.add_option("", "--psd-xml", default=(manifold.utilities.data.DATA_ROOT / 'psd_for_treebank.xml.gz').as_posix(), help="Specify the xml file where the psd is stored.")
 	p.add_option("", "--instrument", default="H1", help="Specify which instrument the bank is for. Default H1")
 	p.add_option("", "--freq-low", type="float", default=15., help="flow of the serach. Default = 15.")
@@ -88,11 +88,6 @@ def parse_args(args=None):
 
 	(opt, args) = p.parse_args(args)
 
-	# These options are not currently supported but could be
-	opt.max_ns_abs_s1z = 1.0
-	opt.min_ns_mass = 1e-6
-	opt.max_ns_mass = 1e-6	
-
 	opt.min_M, opt.max_M = min_max_M_from_mcs_qs([opt.min_mc, opt.max_mc], [opt.min_q, opt.max_q])
 	opt.min_m1, opt.max_m1 = min_max_mi_from_Ms_qs([opt.min_M, opt.max_M],[opt.min_q, opt.max_q])
 	if not opt.min_m2: