PycWB 0.29.0__tar.gz → 0.30.0__tar.gz

{pycwb-0.29.0 → pycwb-0.30.0}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PycWB
-Version: 0.29.0
+Version: 0.30.0
 Summary: A Python package for coherent gravitational wave burst analysis
 Home-page: https://git.ligo.org/yumeng.xu/pycwb
 Author-email: The PycWB team <yumeng.xu@ligo.org>
@@ -32,6 +32,11 @@ Requires-Dist: dask
 Requires-Dist: dask_jobqueue
 Requires-Dist: htcondor
 Requires-Dist: psutil
+Requires-Dist: memspectrum
+Requires-Dist: python-ligo-lw<2.0.0
+Requires-Dist: scikit-learn
+Requires-Dist: jinja2
+Requires-Dist: plotly
 Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page

{pycwb-0.29.0 → pycwb-0.30.0}/PycWB.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PycWB
-Version: 0.29.0
+Version: 0.30.0
 Summary: A Python package for coherent gravitational wave burst analysis
 Home-page: https://git.ligo.org/yumeng.xu/pycwb
 Author-email: The PycWB team <yumeng.xu@ligo.org>
@@ -32,6 +32,11 @@ Requires-Dist: dask
 Requires-Dist: dask_jobqueue
 Requires-Dist: htcondor
 Requires-Dist: psutil
+Requires-Dist: memspectrum
+Requires-Dist: python-ligo-lw<2.0.0
+Requires-Dist: scikit-learn
+Requires-Dist: jinja2
+Requires-Dist: plotly
 Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page

{pycwb-0.29.0 → pycwb-0.30.0}/PycWB.egg-info/SOURCES.txt

@@ -19,7 +19,7 @@ benchmark/io_performance/parquet_perf.py
 benchmark/io_performance/parquet_read_perf.py
 benchmark/io_performance/read_pref.py
 benchmark/likelihood/README.md
-benchmark/likelihood/generate_data_for_likelihood.py
+benchmark/likelihood/data_generator.py
 benchmark/likelihood/perf_cuda.py
 benchmark/likelihood/performance_test_dpf.py
 benchmark/likelihood/performance_test_opt_sky.py
@@ -159,6 +159,8 @@ docs/source/_templates/.keep
 envs/Dockerfile
 envs/environment.yml
 envs/mamba/Dockerfile
+examples/GW190521_search/README.md
+examples/GW190521_search/user_parameters.yaml
 examples/autoencoder/pycwb_injection.ipynb
 examples/autoencoder/user_parameters_injection.yaml
 examples/batch_injection/generate_parameters.py
@@ -193,11 +195,35 @@ examples/lvk_sep_2023/demo.ipynb
 examples/lvk_sep_2023/user_parameters.yaml
 examples/multiple_injection/pycwb_injection.py
 examples/multiple_injection/user_parameters_injection.yaml
+examples/new_injection_infra_with_LHV/README.md
+examples/new_injection_infra_with_LHV/analyze_events.ipynb
+examples/new_injection_infra_with_LHV/injection_parameters.py
+examples/new_injection_infra_with_LHV/user_parameters.yaml
+examples/new_injection_infra_with_LHV/input/LIGO_zero_det_HP.txt
+examples/new_injection_infra_with_LHV/input/aligo_O4high.txt
+examples/new_injection_infra_with_LHV/input/avirgo_O4high_NEW.txt
+examples/new_injection_infra_with_gaussian_noise/README.md
+examples/new_injection_infra_with_gaussian_noise/injection_parameters.py
+examples/new_injection_infra_with_gaussian_noise/user_parameters.yaml
+examples/new_injection_infra_with_real_data/README.md
+examples/new_injection_infra_with_real_data/injection_parameters.py
+examples/new_injection_infra_with_real_data/user_parameters.yaml
+examples/new_injection_infra_with_sky_patch/README.md
+examples/new_injection_infra_with_sky_patch/analyze_events.ipynb
+examples/new_injection_infra_with_sky_patch/injection_parameters.py
+examples/new_injection_infra_with_sky_patch/user_parameters.yaml
+examples/new_injection_infra_with_sky_patch/input/aligo_O4high.txt
 examples/pyseobnr_injection/generate_parameters.py
 examples/pyseobnr_injection/user_parameters_injection.yaml
 examples/pyseobnr_injection/waveform_model/__init__.py
 examples/pyseobnr_injection/waveform_model/waveform.py
+examples/sine_gaussian_injection/README.md
+examples/sine_gaussian_injection/injection_parameters.py
+examples/sine_gaussian_injection/user_parameters.yaml
 examples/waveform_reconstruction/reconstruct_waveform_from_cluster.ipynb
+examples/white_noise_burst_injection/README.md
+examples/white_noise_burst_injection/injection_parameters.py
+examples/white_noise_burst_injection/user_parameters.yaml
 prototypes/search.yaml
 prototypes/cwb_core/coherence.cc
 prototypes/cwb_core/coherence.hh
@@ -253,6 +279,7 @@ pycwb/cli/batch_setup.py
 pycwb/cli/flow.py
 pycwb/cli/get_external_modules.py
 pycwb/cli/gwosc.py
+pycwb/cli/gwosc_data.py
 pycwb/cli/merge_catalog.py
 pycwb/cli/post_process.py
 pycwb/cli/run.py
@@ -282,6 +309,7 @@ pycwb/modules/condor/condor.py
 pycwb/modules/condor/module.yaml
 pycwb/modules/cwb_conversions/__init__.py
 pycwb/modules/cwb_conversions/cluster.py
+pycwb/modules/cwb_conversions/detector.py
 pycwb/modules/cwb_conversions/module.yaml
 pycwb/modules/cwb_conversions/pixel.py
 pycwb/modules/cwb_conversions/series.py
@@ -313,8 +341,10 @@ pycwb/modules/gracedb/module.yaml
 pycwb/modules/gracedb/test_gracedb.py
 pycwb/modules/gwosc/__init__.py
 pycwb/modules/gwosc/gwosc.py
+pycwb/modules/gwosc/utils.py
 pycwb/modules/injection/__init__.py
 pycwb/modules/injection/distribution_utils.py
+pycwb/modules/injection/inj_generators.py
 pycwb/modules/injection/injection.py
 pycwb/modules/injection/module.yaml
 pycwb/modules/injection/par_generator.py
@@ -340,6 +370,7 @@ pycwb/modules/likelihoodWP/dpf.py
 pycwb/modules/likelihoodWP/dpf_cython.pyx
 pycwb/modules/likelihoodWP/likelihood.py
 pycwb/modules/likelihoodWP/sky_stat.py
+pycwb/modules/likelihoodWP/typing.py
 pycwb/modules/likelihoodWP/utils.py
 pycwb/modules/logger/__init__.py
 pycwb/modules/logger/logger.py
@@ -409,6 +440,7 @@ pycwb/modules/workflow_utils/job_setup.py
 pycwb/modules/xtalk/__init__.py
 pycwb/modules/xtalk/monster.py
 pycwb/modules/xtalk/monster_old.py
+pycwb/modules/xtalk/type.py
 pycwb/modules/xtalk/xtalk_data.py
 pycwb/post_production/__init__.py
 pycwb/post_production/workflow.py
@@ -426,6 +458,7 @@ pycwb/types/network_event.py
 pycwb/types/network_pixel.py
 pycwb/types/sparse_series.py
 pycwb/types/time_frequency_series.py
+pycwb/types/time_series.py
 pycwb/types/wdm.py
 pycwb/types/wdm_xtalk.py
 pycwb/utils/__init__.py

{pycwb-0.29.0 → pycwb-0.30.0}/PycWB.egg-info/requires.txt

@@ -22,3 +22,8 @@ dask
 dask_jobqueue
 htcondor
 psutil
+memspectrum
+python-ligo-lw<2.0.0
+scikit-learn
+jinja2
+plotly

{pycwb-0.29.0 → pycwb-0.30.0}/benchmark/likelihood/README.md

@@ -4,10 +4,10 @@ This benchmark is designed to test the performance of the likelihood in Python
 
 ## Prepare the test data for the script
 
-Run the script `generate_data_for_likelihood.py` to generate the cluster data for the likelihood benchmark.
+Run the user parameter file with the customized processor `data_generator.py`
 
 ```bash
-python generate_data_for_likelihood.py
+pycwb run user_parameters_injection.yaml
 ```
 
 A `test_data.pkl` file will be generated for all the Python variables required for likelihood benchmark.

pycwb-0.30.0/benchmark/likelihood/data_generator.py

@@ -0,0 +1,154 @@
+import logging
+import os
+import psutil
+from copy import copy
+from pycwb.config import Config
+from pycwb.modules.super_cluster.super_cluster import supercluster_wrapper
+from pycwb.modules.super_cluster.supercluster import supercluster
+from pycwb.modules.xtalk.monster import load_catalog
+from pycwb.modules.coherence.coherence import coherence
+from pycwb.modules.read_data import generate_injections, generate_noise_for_job_seg, read_from_job_segment, check_and_resample
+from pycwb.modules.data_conditioning import data_conditioning
+from pycwb.types.job import WaveSegment
+from pycwb.types.network import Network
+from pycwb.modules.workflow_utils.job_setup import print_job_info
+
+logger = logging.getLogger(__name__)
+
+
+def process_job_segment(working_dir: str, config: Config, job_seg: WaveSegment, compress_json: bool = True,
+                        catalog_file: str = None, queue=None, production_mode: bool = False, skip_lags: list = None):
+    """
+    The core workflow to process single job segment with trails or lags. 
+
+    Parameters
+    ----------
+    working_dir : str
+        The working directory for the run
+    config : Config
+        The configuration object
+    job_seg : WaveSegment
+        The job segment to process
+    compress_json : bool
+        Whether to compress the json files
+    catalog_file : str
+        The catalog file to save the triggers
+    queue : Queue
+        The queue to send the triggers to the collector for saving in production mode
+    production_mode : bool
+        Whether to run in production mode, if True, the triggers will be sent to the queue instead of saving them in this function
+    skip_lags : list
+        The options to skip certain lags. It is used for resuming the processing after a crash
+    
+    """
+    print_job_info(job_seg)
+
+    if not job_seg.frames and not job_seg.noise and not job_seg.injections:
+        raise ValueError("No data to process")
+
+    base_data = None
+
+    if job_seg.frames:
+        base_data = read_from_job_segment(config, job_seg)
+    if job_seg.noise:
+        base_data = generate_noise_for_job_seg(job_seg, config.inRate, data=base_data)
+    
+    data = base_data
+    
+    # get all the trail_idx from the injections, if there is no injections, use 0
+    trail_idxs = {0}
+    if job_seg.injections:
+        trail_idxs = set([injection.get('trail_idx', 0) for injection in job_seg.injections])
+        
+    # loop over all the trail_idx, if there is no injections or only one trail, only loop once for trail_idx=0
+    for trail_idx in trail_idxs:
+        if job_seg.injections:
+            # use sub_job_seg for each trail_idx to avoid passing the trail_idx to the following functions,
+            sub_job_seg = copy(job_seg)
+            sub_job_seg.injections = [injection for injection in job_seg.injections if injection.get('trail_idx', 0) == trail_idx]
+            logger.info(f"Processing trail_idx: {trail_idx} with {len(sub_job_seg.injections)} injections: {sub_job_seg.injections}")
+            data = generate_injections(config, sub_job_seg, base_data)
+        else:
+            logger.info(f"Processing trail_idx: {trail_idx} without injections")
+            sub_job_seg = job_seg
+        # add trail_idx to the sub_job_seg
+        sub_job_seg.trail_idx = trail_idx
+
+        # check and resample the data
+        data = [check_and_resample(data[i], config, i) for i in range(len(job_seg.ifos))]
+        logger.info("Memory usage: %f.2 MB", psutil.Process().memory_info().rss / 1024 / 1024)
+
+        # data conditioning
+        tf_maps, nRMS_list = data_conditioning(config, data)
+        logger.info("Memory usage: %f.2 MB", psutil.Process().memory_info().rss / 1024 / 1024)
+
+        # initialize network object 
+        network = Network(config, tf_maps, nRMS_list)
+
+        # coherence
+        fragment_clusters = coherence(config, tf_maps, nRMS_list, net=network)
+        logger.info("Memory usage: %f.2 MB", psutil.Process().memory_info().rss / 1024 / 1024)
+
+        # supercluster
+        if config.use_root_supercluster:
+            super_fragment_clusters = supercluster(config, network, fragment_clusters, tf_maps)
+        else:
+            xtalk_coeff, xtalk_lookup_table, layers, _ = load_catalog(config.MRAcatalog)
+            super_fragment_clusters = supercluster_wrapper(config, network, fragment_clusters, tf_maps,
+                                                        xtalk_coeff, xtalk_lookup_table, layers)
+        logger.info("Memory usage: %f.2 MB", psutil.Process().memory_info().rss / 1024 / 1024)
+
+        # only dump the first supercluster
+        d = super_fragment_clusters[0]
+
+        from pycwb.modules.cwb_conversions import convert_fragment_clusters_to_netcluster, convert_netcluster_to_fragment_clusters
+
+        pwc = network.get_cluster(0)
+        wdm_list = network.get_wdm_list()
+        for wdm in wdm_list:
+            wdm.setTDFilter(config.TDSize, config.upTDF)
+
+        # load delay index
+        network.set_delay_index(config.TDRate)
+
+        # load time delay data
+        pwc.cpf(convert_fragment_clusters_to_netcluster(d.dump_cluster(0)), False)
+        pwc.setcore(False, 1)
+        pwc.loadTDampSSE(network.net, 'a', config.BATCH, config.BATCH)
+
+
+        from pycwb.modules.likelihoodWP.likelihood import load_data_from_ifo
+        import numpy as np
+
+        acor = network.net.acor
+        network_energy_threshold = 2 * acor * acor * config.nIFO
+        gamma_regulator = network.net.gamma * network.net.gamma * 2 / 3
+        delta_regulator = abs(network.net.delta) if abs(network.net.delta) < 1 else 1
+        REG = [delta_regulator * np.sqrt(2), 0, 0]
+        netEC_threshold = network.net.netRHO * network.net.netRHO * 2
+
+        n_sky = network.net.index.size()
+
+        ml, FP, FX = load_data_from_ifo(network, config.nIFO)
+
+        cluster_test = convert_netcluster_to_fragment_clusters(pwc)
+        pixels = cluster_test.clusters[0].pixels
+
+        # save FP, FX, rms, n_sky, gamma_regulator, network_energy_threshold to pickle
+        test_data = {
+            'FP': FP,
+            'FX': FX,
+            'cluster': cluster_test.clusters[0],
+            'pixels': pixels,
+            'n_ifo': config.nIFO,
+            'ml': ml,
+            'n_sky': n_sky,
+            'delta_regulator': delta_regulator,
+            'gamma_regulator': gamma_regulator,
+            'netEC_threshold': netEC_threshold,
+            'network_energy_threshold': network_energy_threshold,
+            'netCC': network.net.netCC,
+        }
+        import pickle
+        with open('test_data.pkl', 'wb') as f:
+            pickle.dump(test_data, f)

{pycwb-0.29.0 → pycwb-0.30.0}/benchmark/likelihood/performance_test_opt_sky.py

@@ -1,13 +1,15 @@
 # add .. to path[0]
 import sys, os
-sys.path.insert(0, '..')
+sys.path.insert(0, '../..')
 
 import pickle
 from math import sqrt
 import numpy as np
 from pycwb.modules.likelihoodWP.likelihood import calculate_dpf, find_optimal_sky_localization, calculate_sky_statistics
-from pycwb.modules.likelihoodWP.likelihood import load_data_from_pixels
+from pycwb.modules.likelihoodWP.likelihood import load_data_from_pixels, threshold_cut, fill_detection_statistic
 from pycwb.modules.xtalk.monster import load_catalog, getXTalk_pixels
+from pycwb.modules.xtalk.type import XTalk
+from pycwb.modules.likelihoodWP.typing import SkyMapStatistics
 import time
 
 #################################################################################
@@ -20,12 +22,14 @@ FP = test_data['FP']
 FX = test_data['FX']
 ml = test_data['ml']
 n_sky = test_data['n_sky']
+cluster = test_data['cluster']
 pixels = test_data['pixels']
 gamma_regulator = test_data['gamma_regulator']
 delta_regulator = test_data['delta_regulator']
 network_energy_threshold = test_data['network_energy_threshold']
 netCC = test_data['netCC']
 n_ifo = test_data['n_ifo']
+netEC_threshold = test_data['netEC_threshold']
 REG = np.array([delta_regulator * sqrt(2), 0., 0.])
 
 rms, td00, td90, td_energy = load_data_from_pixels(pixels, n_ifo)
@@ -34,20 +38,17 @@ n_pix = rms.shape[1]
 #####################
 # Load xtalk catalog
 #####################
-if not os.environ.get('HOME_WAT_FILTERS'):
-    print('Please set HOME_WAT_FILTERS to the directory of WAT filters')
-    exit(1)
-
-fn = f"{os.environ.get('HOME_WAT_FILTERS')}/wdmXTalk/OverlapCatalog16-1024.bin"
+fn = f"./wdmXTalk/OverlapCatalog16-1024.bin"
 
 start_time = time.time()
-xtalk_coeff, xtalk_lookup_table, layers, nRes = load_catalog(fn)
+# xtalk_coeff, xtalk_lookup_table, layers, nRes = load_catalog(fn)
+xtalk = XTalk.load(fn)
 print(f"Time for load_catalog: {time.time() - start_time} s")
 
 #######################
 # Get xtalk for pixels
 #######################
-cluster_xtalk_lookup, cluster_xtalk = getXTalk_pixels(pixels, True, layers, xtalk_coeff, xtalk_lookup_table)
+cluster_xtalk_lookup, cluster_xtalk = xtalk.get_xtalk_pixels(pixels, True)
 
 #############################################
 # Convert FP, FX, rms, td00, td90 to float32
@@ -63,20 +64,35 @@ rms = rms.T.astype(np.float32)
 #############################################
 REG[1] = calculate_dpf(FP, FX, rms, n_sky, n_ifo, gamma_regulator, network_energy_threshold)
 
-l_max = find_optimal_sky_localization(n_ifo, n_pix, n_sky, FP, FX, rms, td00, td90, ml, REG, netCC, delta_regulator,
+skymap_statistics = find_optimal_sky_localization(n_ifo, n_pix, n_sky, FP, FX, rms, td00, td90, ml, REG, netCC, delta_regulator,
                                       network_energy_threshold)
-print(f"l_max: {l_max}")
+skymap_statistics = SkyMapStatistics.from_tuple(skymap_statistics)
 
 #############################################
 # Calculate sky statistics
 #############################################
-calculate_sky_statistics(l_max, n_ifo, n_pix, FP, FX, rms, td00, td90, ml, REG, network_energy_threshold, cluster_xtalk,
-                         cluster_xtalk_lookup)
-
+sky_statistics = calculate_sky_statistics(skymap_statistics.l_max, n_ifo, n_pix, FP, FX, rms, td00, td90, ml, REG, network_energy_threshold, cluster_xtalk,
+                         cluster_xtalk_lookup, DEBUG=True)
+
+rejected = threshold_cut(sky_statistics, network_energy_threshold, netEC_threshold)
+if rejected:
+    print(f"Cluster rejected due to threshold cuts: {rejected}")
+
+
+fill_detection_statistic(sky_statistics, skymap_statistics, cluster=cluster, 
+                            n_ifo=n_ifo, xtalk=xtalk,
+                            network_energy_threshold=network_energy_threshold)
+# Eo = 8255.51, Lo = 7805.59, Ep = 1552.52, Lp = 1467.77
+# Gn = 83.9425, Ec = 7571.68, Dc = 0.441406, Rc = 0.999999, Eh = 9.45466
+print(f"expected Eo: 8255.51, Lo: 7805.59, Ep: 1552.52, Lp: 1467.77")
+print(f"expected Gn: 83.9425, Ec: 7571.68, Dc: 0.441406, Rc: 0.999999, Eh: 9.45466, Es: 0, NC: 369, NS: 368")
+print(f"expected N = 83.0302, N2 = 83.9448")
+print(f"expected Np = 70.9542, Em = 1676.79, Lm = 1521.89, norm = 4.91776, Ec = 1539.66, Dc = 0.0897575, ch = 0.932774")
+print(f"expected cc = 1, rho = 27.7458, ecor = 1.4013e-45, penalty = 2.92971e-33, xrho = 4.2039e-45")
 #############################################
 # Performance test
 #############################################
-
+exit()
 total_time = 0
 # # convert FP, FX, rms to float32
 for i in range(10):

{pycwb-0.29.0/benchmark/supercluster → pycwb-0.30.0/benchmark/likelihood}/user_parameters_injection.yaml

@@ -1,7 +1,7 @@
 ###### project ######
 
 outputDir: "output"
-nproc: 6
+nproc: 1
 ###### search ######
 
 #calibVer: "C01"
@@ -69,6 +69,7 @@ iwindow: 100
 # simulation
 simulation: "all_inject_in_one_segment"
 nfactor: 1
+segment_processer: ./data_generator.process_job_segment
 
 injection:
   segment:
@@ -77,17 +78,6 @@ injection:
     noise:
       seeds: [1002, 1004]
   parameters:
-    - mass1: 20
-      mass2: 20
-      spin1z: 0
-      spin2z: 0
-      distance: 1000
-      inclination: 0
-      polarization: 0
-      gps_time: 1126259462.4
-      coa_phase: 0
-      ra: 0
-      dec: 0
     - mass1: 36
       mass2: 29
       spin1x: 0.3
@@ -100,5 +90,6 @@ injection:
       coa_phase: 0
       ra: 0
       dec: 0
+      pol: 0
+      f_lower: 16
       approximant: "IMRPhenomXPHM"
-  approximant: "IMRPhenomXHM"

{pycwb-0.29.0 → pycwb-0.30.0}/bin/pycwb

@@ -13,6 +13,7 @@ import pycwb.cli.xtalk as xtalk
 import pycwb.cli.merge_catalog as merge_catalog
 import pycwb.cli.post_process as post_process
 import pycwb.cli.gwosc as gwosc
+import pycwb.cli.gwosc_data as gwosc_data
 import pycwb.cli.get_external_modules as get_external_modules
 
 def main():
@@ -53,11 +54,16 @@ def main():
     post_process.init_parser(post_process_parser)
     post_process_parser.set_defaults(func=post_process.command)
 
-    #gwosc subcommand
+    # gwosc subcommand
     gwosc_parser = subparsers.add_parser('gwosc', help='Download files to run a pycwb analysis of a GW event of your choice')
     gwosc.init_parser(gwosc_parser)
     gwosc_parser.set_defaults(func=gwosc.command)
 
+    # gwosc_data subcommand
+    gwosc_data_parser = subparsers.add_parser('gwosc-data', help='Download data files to run a pycwb analysis with user parameters')
+    gwosc_data.init_parser(gwosc_data_parser)
+    gwosc_data_parser.set_defaults(func=gwosc_data.command)
+
     # get_external_modules subcommand
     get_external_modules_parser = subparsers.add_parser('get-external-modules', help='Get external modules')
     get_external_modules.init_parser(get_external_modules_parser)

{pycwb-0.29.0 → pycwb-0.30.0}/cwb-core/network.cc

@@ -69,6 +69,22 @@ network::network(const network& value)
 
 network::~network()
 {
+  // Delete detector pointers
+  for (detector* d : ifoList) {
+    delete d;
+  }
+
+  // Delete WDM pointers
+  for (WDM<double>* wdm : wdmList) {
+    if (wdm) {
+      delete wdm;    // delete WDM object
+    }
+  }
+
+  // Delete wc_List
+  for (size_t i = 0; i < wc_List.size(); ++i) {
+      wc_List[i].clear();
+  }
   return;
 }
 

{pycwb-0.29.0 → pycwb-0.30.0}/cwb-core/pycwb.hh

@@ -18,6 +18,8 @@ std::vector<double> inline pycwb_get_wavearray_data(wavearray<double> *wave) {
         data.push_back(wave->data[i]);
     }
 
+    // wave->resize(0); // Clear the wavearray to free memory
+
     return data;
 };
 
@@ -27,6 +29,8 @@ std::vector<short> inline pycwb_get_short_wavearray_data(wavearray<short> *wave)
         data.push_back(wave->data[i]);
     }
 
+    // wave->resize(0); // Clear the wavearray to free memory
+
     return data;
 };
 

{pycwb-0.29.0 → pycwb-0.30.0}/docs/source/mod_cwb.rst

@@ -62,3 +62,28 @@ Replace the following in `monster.cc`
 
 
 to let the output managed by the python code.
+
+Solved part of the memory leak
+--------------------------------
+
+In `network.cc`, the following line was added to the destructor of `Network` class:
+
+.. code-block:: c
+
+    // Delete detector pointers
+    for (detector* d : ifoList) {
+        delete d;
+    }
+
+    // Delete WDM pointers
+    for (WDM<double>* wdm : wdmList) {
+        if (wdm) {
+        delete wdm;    // delete WDM object
+        }
+    }
+
+    // Delete wc_List
+    for (size_t i = 0; i < wc_List.size(); ++i) {
+        wc_List[i].clear();
+    }
+    re

pycwb-0.30.0/examples/GW190521_search/README.md

@@ -0,0 +1,16 @@
+# An example to search for GW190521 
+
+
+## Fetch the required data
+
+```bash
+pycwb gwosc-data user_parameters.yaml
+```
+
+This command will read the `ifo`, `gps_start` and `gps_end` from the `user_parameters.yaml` file and download the frame files and data quality files from GWOSC.
+
+## Run the search
+
+```bash
+pycwb run user_parameters.yaml
+```

pycwb-0.30.0/examples/GW190521_search/user_parameters.yaml

@@ -0,0 +1,95 @@
+###### project ######
+
+outputDir: "output"
+nproc: 6
+###### search ######
+
+#calibVer: "C01"
+#channelName: "DCS-CALIB_STRAIN_CLEAN_C01" # //"DCS-CALIB_STRAIN_CLEAN_C01";  // C00 -> "GDS-CALIB_STRAIN", C01 -> "DCS-CALIB_STRAIN_C01", C02 -> "DCS-CALIB_STRAIN_C02"
+
+analysis: "2G"
+cfg_search: "r"
+
+plot_trigger: True
+save_waveform: True
+plot_waveform: True
+save_injection: True
+plot_injection: True
+compress_output_json: False
+
+optim: False
+
+###### network configuration ######
+ifo: ["L1","H1"]
+refIFO: "L1"
+
+# lags
+lagSize: 1
+lagStep: 1.
+lagOff: 0
+lagMax: 0
+
+# superlags
+slagSize: 1 # number of super lags (simulation=1) - if slagSize=0 -> Standard Segments
+slagMin: 0
+slagMax: 0
+slagOff: 0
+
+# job
+segLen: 1200
+segMLS: 600
+segTHR: 200
+segEdge: 10
+
+# frequency
+fLow: 16.
+fHigh: 1024.
+
+
+levelR: 3
+l_low: 4 # low frequency resolution level		// std (sthr = 2)
+l_high: 10 # high frequency resolution level	// std (sthr = 8)
+
+wdmXTalk: "wdmXTalk/OverlapCatalog16-1024.bin"
+
+healpix: 7
+
+###### cWB production thresholds & regulators ######
+
+bpp: 0.001
+subnet: 0.5
+subcut: -1.0
+netRHO: 5.0
+cedRHO: 5.0
+netCC: 0.5
+Acore: 1.7
+Tgap: 0.2
+Fgap: 128.0
+delta: 0.5
+cfg_gamma: -1.0
+LOUD: 300
+
+pattern: 10
+
+iwindow: 100
+
+# simulation
+simulation: "one_inject_in_one_segment"
+nfactor: 1
+
+
+gps_start: 1242439367.4
+gps_end: 1242446567.4
+
+channelNamesRaw: ['L1:GWOSC-16KHZ_R1_STRAIN', 'H1:GWOSC-16KHZ_R1_STRAIN']
+frFiles: ["input/L1_frames.in", "input/H1_frames.in"]
+
+DQF: [
+  [ "L1" ,"input/L1_cat0.txt", CWB_CAT0, 0., False, False ],
+  [ "H1" ,"input/H1_cat0.txt", CWB_CAT0, 0., False, False ],
+  [ "L1" ,"input/L1_cat1.txt", CWB_CAT1, 0., False, False ],
+  [ "H1" ,"input/H1_cat1.txt", CWB_CAT1, 0., False, False ],
+  [ "L1" ,"input/L1_cat2.txt", CWB_CAT2, 0., False, False ],
+  [ "H1" ,"input/H1_cat2.txt", CWB_CAT2, 0., False, False ],
+]
+