lalsuite 7.26.2.dev20251206__cp312-cp312-manylinux_2_28_x86_64.whl

lalapps/data/lib_pipe_example.ini

@@ -0,0 +1,303 @@
+###############################################################################
+####### lib (lalinference burst) example initialisation file #########################
+####### For more information, see docs page at        #########################
+# https://wiki.ligo.org/viewauth/Bursts/LIB//
+###############################################################################
+
+[analysis]
+# IFOs to use for this analysis
+ifos=['H1','L1','V1']
+
+# select the engine to use, only tested with lalinferencenest
+engine=lalinferencenest
+
+# Number of parallel runs for each event. Using parallel runs increases the number of sample, giving a better exploration of the parameter space. Should be set to 2 (or 3) unless you know what you are doing.
+nparallel=3
+
+# Turn the coherence test (will result in Nifos+1 runs)?
+coherence-test=True
+
+# Upload results to the gravitational-wave events database?
+upload-to-gracedb=False
+
+#chain seed (incremented by one for each parallel chain).
+seed=4321
+
+# noise generator seed (incremented by eventnum for each event). If commented the cpu time will be used to initialize the chain (that may make it problematic to reproduce results)
+dataseed=1234
+
+# Option to specify text file for dumping time slide information
+#time-slide-dump=path/to/file
+
+#####################################################################################
+[paths]
+# the following can be given with command line (see below)
+#webdir is the base output dir for results pages
+webdir=/home/albert.einstein/public_html/lalinferenceresults/
+
+#baseurl is the www address of the above path, e.g:
+baseurl=https://atlas1.atlas.aei.uni-hannover.de/~albert.einstein/lalinferenceresults/
+
+#####################################################################################
+[input]
+# User-specified length of the psd. if not specified, will be automatically calculated from segment availability
+# psd-length=1024
+
+# User-specified psd start time (used for manual over-ride, if not given pipeline determines this automatically)
+# psd-start-time= 968621415
+
+# Maximum length to use for automatically-determined psdlength options
+max-psd-length=100
+
+# spacing between trigger times and start of psd estimation
+padding=10
+
+# all will analyze all the events in the injetion table. Ranges and/or single events can also be specified. E.g. events=[0,4,5,10-12] will analyze the event 0,4,5,10,11,12.
+events=all
+
+# Overlap between segments when analysing a continuous span of time
+# segment-overlap = 32
+
+# Can manually over-ride time limits here - must use these when analysing
+# a stretch of time, i.e. when no other input files are given
+#gps-start-time=
+#gps-end-time=
+analyse-all-time=False
+
+# Can manually specify input files here or over-ride on the command line
+#gps-time-file=
+#injection-file=
+#sngl-inspiral-file=
+#coinc-inspiral-file=
+#pipedown-db=
+
+# If this option is set "true" and a pipedown database is used as input, the pipeline will analyse the time slide events
+timeslides=false
+
+# Uncomment the following line to ignore science segments. Useful when using fake noise
+#ignore-science-segments=True
+
+#####################################################################################
+[condor]
+# Replace with your own executable locations
+
+datafind=/home/albert.einstein/bin/ligo_data_find
+mergescript=/home/albert.einstein/bin/lalapps_nest2pos
+resultspage=/home/albert.einstein/bin/cbcBayesBurstPostProc
+segfind=/home/albert.einstein/bin/ligolw_segment_query
+ligolw_print=/home/albert.einstein/bin/ligolw_print
+coherencetest=/home/albert.einstein/bin/lalapps_coherence_test
+lalinferencenest=/home/albert.einstein/bin/lalinference_nest
+lalinferencemcmc=/home/albert.einstein/bin/lalinference_mcmc
+lalinferencebambi=/home/albert.einstein/bin/lalinference_bambi
+
+# MPIrun only needed for MCMC, if not available use /bin/true
+mpirun=/home/albert.einstein/bin/mpirun
+
+# GraceDB only needed for online analysis, if not available use /bin/true
+gracedb=/home/albert.einstein/bin/gracedb
+
+# Queue information if needed to run on specific nodes
+#queue=
+# Samples to injection file (only for prior PP test)
+pos_to_sim_inspiral=/home/albert.einstein/bin/cbcBayesPosToSimBurst
+ppanalysis=/home/albert.einstein/bin/cbcBayesBurstPPAnalysis
+
+#####################################################################################
+# Section used by the datafind jobs (not used with simulated noise)
+[datafind]
+# Specifies we look for data in the local filesystem (as opposed to remote data)
+url-type=file
+
+# "type" of data to look for with ligo_data_find
+types={'H1':'H1_LDAS_C02_L2','L1':'L1_LDAS_C02_L2','V1':'HrecOnline'}
+
+#####################################################################################
+[data]
+# Channels to use to read h(t) from frame files
+# S5 has LSC-STRAIN, S6 has LDAS-STRAIN.
+channels={'H1':'H1:LDAS-STRAIN','L1':'L1:LDAS-STRAIN','V1':'V1:h_16384Hz'}
+
+#####################################################################################
+[lalinference]
+
+# fake-cache is used to specify the fake names for simulated data
+# For simulated LIGO noise use 'LALLIGO', Virgo noise: 'LALVirgo'
+# For simulated advanced detectors use 'LALSimAdLIGO' and 'LALSimAdVirgo'
+# If not enabled, the pipeline will look for real data
+#fake-cache={'H1':'LALLIGO','L1':'LALLIGO','V1':'LALVirgo'}
+#fake-channels={'H1':'Science','L1':'Science','V1':'Science'}
+
+# fLow can be used to specify the low frequency cutoff for each detector
+# flow = {'H1':40,'L1':40,'V1':40}
+
+#####################################################################################
+[engine]
+# All options in this section are passed to the sampler.
+#Some useful ones are noted below.
+# Options passed to a sampler which does not know them will simply be ignored
+
+# REQUIRED SETTINGS:
+# Nlive specifies the number of live points for each job
+nlive=1024
+
+# Sampling rate for data
+srate=2048
+
+# Segment length to use for analysis (should be long enough for whole template
+seglen=4
+
+# OPTIONAL SETTINGS:
+
+# Use lalinference_nest's resume ability if the run is interrupted
+resume=
+
+# approx can be used to manually specify an approximant
+# If this is not given, the code will use SineGaussianF. Available approx are SineGaussian(F) and Gaussian(F) as well as all the usual CBC approx.
+#approx=SineGaussianF
+
+# maxmcmc set the maximum chain length for the nested sampling sub-chains. Default 5000
+#maxmcmc=5000 # Auto determination is on, but the length cannot be longer than that.
+
+# Priors
+# For all parameters known to lalinference, the min and max default prior can be overwritten with
+#parname-min=MIN
+#parname-max=MAX
+
+# The starting point for the MCMC chain(s) can be specified with
+#parname=START
+
+# Parameters can be fixed to some value with
+#fix-parname=FIXVALUE
+
+#currently known parameters, together with default [min-max] are (radiants for angle, Mpc for distance, Msun for masses)
+
+#time                         Waveform time [trigtime-0.1-trigtime+0.1]
+#rightascension               Rightascension [0,2Pi]
+#declination                  Declination [-Pi/2,Pi/2]
+#polarisation                 Polarisation angle [0,Pi]
+
+# -- CBC template parameters -- #
+
+#chirpmss                     Chirpmass [1,15.3]
+#eta                          Symmetric massratio (needs --use-eta) [0,0.0312]
+#q                            Asymmetric massratio (a.k.a. q=m2/m1 with m1>m2) [0.033,1]
+#theta_jn                     Angle between J and line of sight [0,Pi]
+#distance                     Distance [1,2000]
+#logdistance                  Log Distance (requires --use-logdistance) [log(1),log(2000)]
+
+#Spin Parameters:
+#a1                           Spin1 magnitude (for precessing spins) [0,1]
+#a_spin1                      Spin1 magnitude (for aligned spins) [-1,1]
+#a2                           Spin2 magnitude (for precessing spins) [0,1]
+#a_spin2                      Spin2 magnitude (for aligned spins) [-1,1]
+#tilt_spin1                   Angle between spin1 and orbital angular momentum [0,Pi]
+#tilt_spin2                   Angle between spin2 and orbital angular momentum  [0, Pi]
+#phi_12                       Difference between spins' azimuthal angles [0,2Pi]
+#phi_jl                       Difference between total and orbital angular momentum azimuthal angles [0,2Pi]
+
+#Equation of State parameters (requires --use-tidal or --use-tidalT):
+#lambda1                      lambda1 [0,3000]
+#lambda2                      lambda2 [0,3000]
+#lambdaT                      lambdaT [0,3000]
+#dLambdaT                     dLambdaT [-500,500]
+
+# Settings for allowed component masses in Solar Masses, with default values
+#component-max=30.0
+#component-min=1.0
+
+# Allowed total masses in Msun (note, used together with component masses, mc,q,eta priors may lead to inconsistencies. Be careful!)
+#mtotal-max=35
+#mtotal-min=2
+
+# -- Burst template parameters --#
+
+#frequency                     Central Frequency [40,1024]Hz
+#quality                       Quality [3, 100]
+#duration                      Duration (for Gaussian template) [1e-4,0.5]sec
+#hrss                          Hrss [1e-23,1e-15] Must specify use-hrss, otherwise loghrss is used
+#loghrss                       Natural log of hrss [-52.959,-34.539]
+#alpha                         Ellipticity parameter [0,2Pi]rads
+
+# Setting time prior [secs]
+#dt=0.1
+
+# The following three options control various marginalised likelihoods. Use at most one.
+# Analytically marginalise over phase (only for Newtonian amplitude CBC orders - not usable for bursts)
+#margphi=
+# Semi-analytically marginalise over time (both for burst and CBC. Need significant sampling rate)
+#margtime=
+# Semi-analytically marginalise over time and phase (only for Newtonian amplitude CBC orders -- not usable for bursts)
+#margtimephi=
+
+# By default the CBC code use spins if the choosen approximant supports spin. NOTE that this include TaylorF2, which will be run with aligned spins.
+# Several options, here below,  allows the user to choose a particular setting:
+
+#Disable spin for waveforms which would normally have spins
+#disable-spin=
+
+# Only allow aligned spins
+#aligned-spin=
+
+# Only allow the heavier mass to spin (can be used together with aligned-spin)
+#singleSpin=
+
+# Print progress information throughout the run
+progress=
+
+
+# Sample in log(distance) for improved sampling
+#use-logdistance=
+# Sample in eta instead than q=m2/m1
+#use-eta=
+
+# --- Miscellanea Burst --- #
+
+# Sample in hrss (default is loghrss)
+#use-hrss
+
+# assume cross-only polarization
+#cross_only=
+# assume plus-only polarization
+#plus_only=
+
+
+#####################################################################################
+[mpi]
+# Settings when running with MPI for lalinference_mcmc or lalinference_bambi
+
+# number of CPUs desired and how much memory on each (MB)
+machine-count=8
+machine-memory=512
+
+#####################################################################################
+[resultspage]
+# Settings for the results pages (see cbcBayesPostProc --help for more)
+
+# Size of the side of the sky bins in degrees
+skyres=0.5
+
+#disable production of tar.zg archival data of postproc
+archive=None
+# Do no create 2D plots, which take a while to generate
+#no2D=
+
+# Send an email linking to the final page
+#email=albert.einstein@ligo.org
+
+#####################################################################################
+[segfind]
+# URL of the segment database
+segment-url=https://segdb.ligo.caltech.edu
+
+#####################################################################################
+[segments]
+# Names of science segments to use
+l1-analyze = L1:DMT-SCIENCE:4
+h1-analyze = H1:DMT-SCIENCE:4
+v1-analyze = V1:ITF_SCIENCEMODE:7
+
+[ppanalysis]
+# This section only used with lalinference_pp_pipe
+# Webdir is the output path for the PP check page
+webdir=/home/albert.einstein/public_html/lalinference/review/PPcheck/

lalapps/data/power_pipe.ini

@@ -0,0 +1,129 @@
+; excess power offline pipeline construction script configuration file
+;
+
+[condor]
+universe = standard
+datafind  = $ENV(HOME)/local/bin/LSCdataFind
+lalapps_binj = $ENV(HOME)/local/bin/lalapps_binj
+lalapps_power = $ENV(HOME)/local/bin/lalapps_power
+ligolw_add = $ENV(HOME)/local/bin/ligolw_add
+lalburst_binjfind = $ENV(HOME)/local/bin/lalburst_injfind
+lalburst_cluster = $ENV(HOME)/local/bin/lalburst_cluster
+lalburst_coinc = $ENV(HOME)/local/bin/lalburst_coinc
+lalburst_cut = $ENV(HOME)/local/bin/lalburst_cut
+ligolw_sqlite = $ENV(HOME)/local/bin/ligolw_sqlite
+lalburst_power_meas_likelihood = $ENV(HOME)/local/bin/lalburst_power_meas_likelihood
+
+[pipeline]
+user_tag = PLAYGROUND
+out_dir = logs
+cache_dir = cache
+ifos = H1 H2 L1
+seglist_H1 = H1.segwizard
+seglist_H2 = H2.segwizard
+seglist_L1 = L1.segwizard
+;tiling_phase_H2 = 0.0055769230769230765
+psds_per_power = 16
+psds_per_injection = 16
+injection_bands = 4
+enable_clustering = True
+files_per_bucluster = 10
+files_per_bucut = 30
+files_per_burca = 1
+files_per_binjfind = 30
+likelihood_data_cache_base = power_burca2_likelihood
+
+[datafind]
+gaps =
+url-type = file
+type_H1 = H1_RDS_C04_LX
+type_H2 = H2_RDS_C04_LX
+type_L1 = L1_RDS_C04_LX
+
+[lalapps_binj]
+min-frequency = $(macroflow)
+max-frequency = $(macrofhigh)
+min-bandwidth = 4.0
+max-bandwidth = 256.0
+min-duration = 0.00390625
+max-duration = 0.25
+min-e-over-r2 = 1e-20
+max-e-over-r2 = 1e-13
+population = targeted
+; galactic core
+; ra = 27940.04 s = 7 h 45 m 40.04 s
+; dec = -29o 00' 28.1"
+ra-dec = 2.0318570464121519,-0.50628171572274738
+seed = $(macroseed)
+; one injection every time-step seconds
+time-step = 35.014087480216972
+
+[input]
+; Power jobs are subclassed from AnalysisJob, which gets the channel name
+; from the [input] section
+channel = LSC-STRAIN
+
+[lalapps_power]
+bandwidth = 2048
+filter-corruption = 4096
+; 0.0 Hz == disable high-pass filter
+high-pass = 0.0
+low-freq-cutoff = 70.0
+max-event-rate = 10000
+max-tile-bandwidth = 256.0
+max-tile-duration = 0.25
+psd-average-points = 2097152
+ram-limit = 1024
+resample-rate = 8192
+tile-stride-fraction = 0.5
+window-length = 65536
+;gaussian-noise-rms = 4e-20
+
+[lalapps_power_H1]
+;calibration-cache = /home/kipp/scratch/work/calibration/S4/H1-CAL-V04-793152133-795986713.cache
+confidence-threshold = 8.5
+
+[lalapps_power_H2]
+;calibration-cache = /home/kipp/scratch/work/calibration/S4/H2-CAL-V04-793132513-796025053.cache
+confidence-threshold = 7.3
+
+[lalapps_power_L1]
+;calibration-cache = /home/kipp/scratch/work/calibration/S4/L1-CAL-V04-793128493-795929533.cache
+confidence-threshold = 7.8
+
+[lalburst_cut]
+comment = $(macrocomment)
+program = lalapps_power
+inj-made-only =
+
+[lalburst_cluster]
+comment = $(macrocomment)
+cluster-algorithm = excesspower
+
+[ligolw_add]
+
+[lalburst_binjfind]
+comment = $(macrocomment)
+match-algorithm = excesspower
+
+[lalburst_coinc]
+comment = $(macrocomment)
+program = lalapps_power
+coincidence-algorithm = excesspower
+
+[ligolw_sqlite]
+replace =
+preserve-ids =
+database = $(macrodatabase)
+tmp-space = $ENV(TMPDIR)
+
+[lalburst_power_meas_likelihood]
+verbose =
+tmp-space = $ENV(TMPDIR)
+
+[lalburst_coinc2]
+verbose =
+comment = $(macrocomment)
+program = lalapps_power
+coincidence-algorithm = excesspower2
+tmp-space = $ENV(TMPDIR)

lalapps/data/unimodalMeans.csv

@@ -0,0 +1,2 @@
+m1,  m2,  costheta_jn,  phi_orb,  psi,  ra,  dec,  logdistance,  time,  a1,  a2,  tilt1,  tilt2,  phi12,  phi_jl
+16.0, 7.0, 0.0, 3.141592653589793, 3.141592653589793, 3.141592653589793, 0.0, 4.0, 0.0, 0.5, 0.5, 1.5707963267948966, 1.5707963267948966, 3.141592653589793, 3.141592653589793

lalapps/git_version.py

@@ -0,0 +1,64 @@
+# -*- coding: utf-8 -*-
+# git_version.py - vcs information module
+#
+# Copyright (C) 2010 Nickolas Fotopoulos
+# Copyright (C) 2012-2013 Adam Mercer
+# Copyright (C) 2016 Leo Singer
+#
+# This program is free software; you can redistribute it and/or modify
+# it under the terms of the GNU General Public License as published by
+# the Free Software Foundation; either version 2 of the License, or (at
+# your option) any later version.
+#
+# This program is distributed in the hope that it will be useful, but
+# WITHOUT ANY WARRANTY; without even the implied warranty of
+# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
+# General Public License for more details.
+#
+# You should have received a copy of the GNU General Public License
+# along with with program; see the file COPYING. If not, write to the
+# Free Software Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston,
+# MA 02110-1301 USA
+
+id = "eeff03c6df54cc1b9aa861018eac9ca5072a0a77"
+date = "2025-12-4 04:45:18 +0000"
+branch = "None"
+tag = "None"
+if tag == "None":
+    tag = None
+author = "Karl Wette <karl.wette@ligo.org>"
+builder = "Unknown User <>"
+committer = "Karl Wette <karl.wette@ligo.org>"
+status = "CLEAN: All modifications committed"
+version = id
+verbose_msg = """Branch: None
+Tag: None
+Id: eeff03c6df54cc1b9aa861018eac9ca5072a0a77
+
+Builder: Unknown User <>
+Repository status: CLEAN: All modifications committed"""
+
+import warnings
+
+class VersionMismatchError(ValueError):
+    pass
+
+def check_match(foreign_id, onmismatch="raise"):
+    """
+    If foreign_id != id, perform an action specified by the onmismatch
+    kwarg. This can be useful for validating input files.
+
+    onmismatch actions:
+      "raise": raise a VersionMismatchError, stating both versions involved
+      "warn": emit a warning, stating both versions involved
+    """
+    if onmismatch not in ("raise", "warn"):
+        raise ValueError(onmismatch + " is an unrecognized value of onmismatch")
+    if foreign_id == "eeff03c6df54cc1b9aa861018eac9ca5072a0a77":
+        return
+    msg = "Program id (eeff03c6df54cc1b9aa861018eac9ca5072a0a77) does not match given id (%s)." % foreign_id
+    if onmismatch == "raise":
+        raise VersionMismatchError(msg)
+
+    # in the backtrace, show calling code
+    warnings.warn(msg, UserWarning)