Differences between revisions 23 and 68 (spanning 45 versions)
Revision 23 as of 2019-07-13 10:38:16
Size: 1447
Comment:
Revision 68 as of 2019-08-02 14:46:56
Size: 4381
Comment:
Deletions are marked like this. Additions are marked like this.
Line 1: Line 1:
= Engineering Run 2 (2019/7/13) = = Engineering Run 2 (JST 2019/7/13) =
<<TableOfContents(2)>>
Line 3: Line 4:
 * JST: Jul. 13, 2019, 12:00 - 18:00 JST
 * GPS: 1247022018 - 1247043618
 * JST: Jul. 13, 2019, --(12:00)-- 12:37 - 18:00 JST
Line 6: Line 6:
=== task list ===  * GPS: --(1247022018)-- 1247024238 - 1247043618

== Data related ==
Line 10: Line 12:
Line 11: Line 14:
Line 12: Line 16:
Line 13: Line 18:
Line 16: Line 22:
Line 17: Line 24:
  * K1:CAL-CS_PROC_MICH_DISPLACEMENT_DQ   * K1:CAL-CS_PROC_MICH_DISPLACEMENT_DQ (We used this channel.)
Line 26: Line 33:
  * K1:GRD-LSC_LOCK_OK = 1
* K1:GRD-LSC_LOCK_STATE_N = 157
Line 32: Line 40:
 * code development
  * generate cache file
  * read the data
  * PSD estimation
 * CBC lock segments
  * We checked lock flag of each data and collected the locked data. We considered lock segments by file units for convenience. If compared with [[http://gwwiki.icrr.u-tokyo.ac.jp/JGWwiki/KAGRA/Subgroups/DET/RUN#ER_with_MICH_on_Jul._13 | lock segment(DetChar)]], this lock segment is small slightly.
Line 37: Line 43:
 * search
 * how to output the result?
  * histogram of SNR
  * glitch list(GPS, JST, mass, ISCO frequency, duration)
 *
  || '''#''' || '''start GPS''' || '''stop GPS''' || '''duration''' ||
  || 1 || 1247024288 || 1247043616 || 19328 ||

== Off line analysis ==

=== Check point ===

 * Check the available locked data '''(done)'''
  * Summarize the data length '''(done)'''

 * Estimate power spectrum density '''(done)'''
  * Summarize information in pdf '''(done)'''

 * Generate lock list and segment list '''(done)'''

 * Generate template bank '''(skipped)'''
  * How many waveforms we need? '''(skipped)'''

 * Search signals by matched filter and chi^2 '''(done)'''

 * Summary the output of search '''(done)'''
  * histogram of SNR '''(done)'''
  * glitch list(GPS, JST, mass, ISCO frequency, duration) '''(done)'''

=== Output of search ===

 * Summary of locked data [[http://gwwiki.icrr.u-tokyo.ac.jp/JGWwiki/KAGRA/DAWG/CBC/ER2#Data_related | CBC locked segment]]

 * Averaged PSD [[https://gwdoc.icrr.u-tokyo.ac.jp/cgi-bin/private/DocDB/ShowDocument?docid=10448 | JGWdoc]]

  * Power spectrum density (PSD) is averaged over 1 sec chunks by taking median.

 * Calculation time to search for signals

  * Read 500 sec data : < 1 min

  * Estimate PSD : 2~3 min.

  * matched filter+chi^2 with single waveform (not bank) : < 1 min
Line 44: Line 84:
=== Related page ===  * Histogram of detection statistics [[https://gwdoc.icrr.u-tokyo.ac.jp/cgi-bin/private/DocDB/ShowDocument?docid=10448 | JGWdoc]]
  * If data is Gaussian, <snr>=1.
Line 46: Line 87:
[[http://gwwiki.icrr.u-tokyo.ac.jp/JGWwiki/KAGRA/Subgroups/DET/RUN#ER_with_MICH_on_Jul._13 | ER2 (DetChar)]]  * Time series of detection statistics [[https://gwdoc.icrr.u-tokyo.ac.jp/cgi-bin/private/DocDB/ShowDocument?docid=10448 | JGWdoc]]
  * waveform : 1.4 - 1.4 Msun
  * Beginning and ending 50 sec data were removed from the result.
  * We uploaded the 3 kinds of time series results which were maximized over rho(matched filter), chi square, rhohat(detection statistics) during 1 second window.
   We added time series of rho and chi square, because time series of rhohat don't contain the information of the glitch.
  * File names are
    * "cbc_max_rho_ER2.txt" (this data indicates the time when data contains high power)
    * "cbc_max_chi2_ER2.txt" (this data indicates the time when data contains high power in specific frequency band)
    * "cbc_max_rhohat_ER2.txt" (this data indicates the time when data "might" contains gravitational wave signal)

== Online analysis ==

 * Online matched filtering
  * July 13, 17:08 - 18:00 JST
  * Using low latency data files (1 sec) dumped from shared memory
  * About 1 minutes deley
  * CH: K1:DCS-CALIB_STRAIN_C10
  * Template: (30,30)Msun, no spin, Waveform: IMRPhenomD
  * Used pre-computed PSD

 * Plots

== Link ==

 * [[http://gwwiki.icrr.u-tokyo.ac.jp/JGWwiki/KAGRA/DAWG/CBC/ER1 | CBC search in ER1]]

 * [[http://gwwiki.icrr.u-tokyo.ac.jp/JGWwiki/KAGRA/Subgroups/DET/RUN#ER_with_MICH_on_Jul._13 | ER2 (DetChar)]]

Engineering Run 2 (JST 2019/7/13)

  • JST: Jul. 13, 2019, 12:00 12:37 - 18:00 JST

  • GPS: 1247022018 1247024238 - 1247043618

  • KAGRA full data
    • m31:/data/full : Full KAGRA data
  • Low latency data
    • m31:/data/ll_merged/ : Low latency data packed in 4096 seconds
    • hydra:/data/ll_merged/ : Low latency data packed in 4096 seconds
    • hydra: /data/ll-sharing/shm/llhoft/ : Low latency data (1 second) dumped from shared memory (retained only for 8192 seconds)
  • channel name
    • From ER2 (DetChar)

    • K1:CAL-CS_PROC_MICH_DISPLACEMENT_DQ (We used this channel.)
    • K1:CAL-CS_PROC_IMC_FREQUENCY_DQ
  • Low latency data channels
    • K1:DCS-CALIB_STRAIN_C10 16384
    • K1:DCS-CALIB_STATE_VECTOR_C10 16
    • K1:CAL-CS_PROC_IMC_FREQUENCY_DQ 16384
  • lock flag
    • K1:GRD-LSC_LOCK_STATE_N = 157
  • data quality vector
    • definition

    • We cbc group should check the second bit is 1.
  • CBC lock segments
    • We checked lock flag of each data and collected the locked data. We considered lock segments by file units for convenience. If compared with lock segment(DetChar), this lock segment is small slightly.

      #

      start GPS

      stop GPS

      duration

      1

      1247024288

      1247043616

      19328

Off line analysis

Check point

  • Check the available locked data (done)

    • Summarize the data length (done)

  • Estimate power spectrum density (done)

    • Summarize information in pdf (done)

  • Generate lock list and segment list (done)

  • Generate template bank (skipped)

    • How many waveforms we need? (skipped)

  • Search signals by matched filter and chi^2 (done)

  • Summary the output of search (done)

    • histogram of SNR (done)

    • glitch list(GPS, JST, mass, ISCO frequency, duration) (done)

  • Summary of locked data CBC locked segment

  • Averaged PSD JGWdoc

    • Power spectrum density (PSD) is averaged over 1 sec chunks by taking median.
  • Calculation time to search for signals
    • Read 500 sec data : < 1 min

    • Estimate PSD : 2~3 min.
    • matched filter+chi^2 with single waveform (not bank) : < 1 min

  • Histogram of detection statistics JGWdoc

    • If data is Gaussian, <snr>=1.

  • Time series of detection statistics JGWdoc

    • waveform : 1.4 - 1.4 Msun
    • Beginning and ending 50 sec data were removed from the result.
    • We uploaded the 3 kinds of time series results which were maximized over rho(matched filter), chi square, rhohat(detection statistics) during 1 second window.
      • We added time series of rho and chi square, because time series of rhohat don't contain the information of the glitch.
    • File names are
      • "cbc_max_rho_ER2.txt" (this data indicates the time when data contains high power)
      • "cbc_max_chi2_ER2.txt" (this data indicates the time when data contains high power in specific frequency band)
      • "cbc_max_rhohat_ER2.txt" (this data indicates the time when data "might" contains gravitational wave signal)

Online analysis

  • Online matched filtering
    • July 13, 17:08 - 18:00 JST
    • Using low latency data files (1 sec) dumped from shared memory
    • About 1 minutes deley
    • CH: K1:DCS-CALIB_STRAIN_C10
    • Template: (30,30)Msun, no spin, Waveform: IMRPhenomD
    • Used pre-computed PSD
  • Plots

KAGRA/DAWG/CBC/ER2 (last edited 2019-08-02 14:46:56 by HirotakaYuzurihara)