Differences between revisions 44 and 68 (spanning 24 versions)
Revision 44 as of 2019-07-16 12:04:38
Size: 2873
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 46: Line 46:
== Check point == == Off line analysis ==
Line 48: Line 48:
 * Check the available locked data(OK)
  * Summarize the data length
=== Check point ===
Line 51: Line 50:
 * Estimate power spectrum density(OK)
  * Sumarize information in pdf(OK?)
 * Check the available locked data '''(done)'''
  * Summarize the data length '''(done)'''
Line 54: Line 53:
 * Generate lock list and segment list(Skipped)  * Estimate power spectrum density '''(done)'''
  * Summarize information in pdf '''(done)'''
Line 56: Line 56:
 * Generate template bank(not yet)
  * How many waveforms we need?(not yet)
 * Generate lock list and segment list '''(done)'''
Line 59: Line 58:
 * Search signals by matched filter and chi^2(not yet)  * Generate template bank '''(skipped)'''
  * How many waveforms we need? '''(skipped)'''
Line 61: Line 61:
 * Summary the output of search (not yet)
  * histogram of SNR(not yet)
  * glitch list(GPS, JST, mass, ISCO frequency, duration)(not yet)
  * other?
 * Search signals by matched filter and chi^2 '''(done)'''
Line 66: Line 63:
 * Summary the output of search '''(done)'''
  * histogram of SNR '''(done)'''
  * glitch list(GPS, JST, mass, ISCO frequency, duration) '''(done)'''
Line 67: Line 67:

== Output of search ==
=== Output of search ===
Line 72: Line 71:
 * Averaged PSD [[https://gwdoc.icrr.u-tokyo.ac.jp/cgi-bin/private/DocDB/ShowDocument?docid=10448 | PSDs]]  * Averaged PSD [[https://gwdoc.icrr.u-tokyo.ac.jp/cgi-bin/private/DocDB/ShowDocument?docid=10448 | JGWdoc]]
Line 76: Line 75:
 * Calculation time to search for signals(PSD estimate, template bank generation, matched filter, chi^2 calculation)  * Calculation time to search for signals
Line 78: Line 77:
 * Histogram of detection statistics   * Read 500 sec data : < 1 min
Line 80: Line 79:
 * Time series of detection statistics   * Estimate PSD : 2~3 min.

  * matched filter+chi^2 with single waveform (not bank) : < 1 min
Line 83: Line 84:
=== Link ===  * 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 85: 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)