Diff for "KAGRA/Subgroups/DET/gwpy/coupling"

Differences between revisions 5 and 11 (spanning 6 versions)

Coupling plot tool

Tool direction

In k1det1

[detchar@k1det1 ~]$ cd /home/detchar/gwpy-tools/kagra_coupling_plot/
[detchar@k1det1 kagra_coupling_plot]$ ls
kagra_coupling_plot.py  plot  run_plot.sh

In seikai

[user@seikai ~]$ cd /home/detchar/gwpy-tools/kagra_coupling_plot
[user@seikai kagra_coupling_plot]$ ls
kagra_coupling_plot.py  plot  run_plot.sh

How to use

The run_plot.sh file contains usage examples and argument explanations.

[detchar@k1det1 kagra_coupling_plot]$ cat run_plot.sh 

python kagra_coupling_plot.py -ach K1:VIS-EY_TM_OPLEV_SUM_OUT_DQ -bch K1:LSC-MICH_CTRL_CAL_OUT_DQ -s 1145621548.0 -d 600 -a 100 -f0 10 -f1 500 -p 1 

## [EX] python kagra_coupling_plot.py -ach K1:VIS-EY_TM_OPLEV_SUM_OUT_DQ -bch K1:LSC-MICH_CTRL_CAL_OUT_DQ -s 1199756064.0 -d 600 -a 100 -f 103.5 -p 1
## [1] -ach     :: Enter first channel name to compare
## [2] -bch     :: Enter second channel name to compare
## [3] -s       :: Enter start GPS time
## [4] -d       :: Enter time from star GPS time
## [5] -a       :: Enter Averaging Number, CAN NOT be used with -df at the same time
## [6] -df      :: Enter Frequency interval, df, CAN NOT be used with -a at the same time.
## [7] -f0      :: Enter start frequency to obtain amplitude value list
## [8] -f1      :: Enter end frequency to obtain amplitude value list
## [9] -p       :: Enter 0 or 1 to determine whether to save the plot. 0 or None is not drawn. 1 is saved.

and When you run,

[detchar@k1det1 kagra_coupling_plot]$ python kagra_coupling_plot.py -ach K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ -bch K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ -s 1145621548.0 -d 736 -df 0.1 -f0 10 -f1 100 -p 1

Loading Time series cache datas and time segment...

---- Entered A Chanel : K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ
---- Entered B Chanel : K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ
---- Entered Start GPS time : 1145621548.0
---- Entered End GPS time : 1145622284.0
---- Entered Time duration : 736
---- Entered Average : 73.6
---- Entered df : 0.1 Hz

Checking each channel sample rate...

---- A Channel Sample Rate : 2048 Hz
PASS
---- B Channel Sample Rate : 2048 Hz
PASS


### Calcurate coupling function during 1145621548.0 to 1145622284.0 ###

    Saved Coupling plot in ./plot/1145621548.0-736/K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ-to-K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ

    Amplitude list saveed in ./plot/1145621548.0-736/K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ-to-K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ direction

    Do you want to print the saved list? Enter the Numer (Yes = 1 / No = 0) : 

１

10 Hz : 2.65813895339e-07
11 Hz : 4.86833812374e-07
12 Hz : 3.2527766507e-06 
.
.
.
97 Hz : 3.06364533742e-08
98 Hz : 3.87633960486e-08
99 Hz : 3.20635358264e-08
100 Hz : 4.65378668935e-08


##### Completed #####

Slide

Ratio plot of two channels Report.001.png

-  ⇤ ← Revision 5 as of 2018-01-11 13:25:41 → 
  Size: 2220
  Editor: PJJung
  Comment:
+   ← Revision 11 as of 2018-02-20 21:18:36 → ⇥
  Size: 3545
  Editor: PJJung
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 Line 5:
+ * In k1det1
-Line 6:
+Line 8:
-[detchar@k1det1 ~]$ cd gwpy-tools/kagra_coupling_plot/
+[detchar@k1det1 ~]$ cd /home/detchar/gwpy-tools/kagra_coupling_plot/
-Line 8:
+Line 10:
+kagra_coupling_plot.py  plot  run_plot.sh
}}}


 * In seikai

{{{
[user@seikai ~]$ cd /home/detchar/gwpy-tools/kagra_coupling_plot
[user@seikai kagra_coupling_plot]$ ls
-Line 18:
+Line 29:
-python kagra_coupling_plot.py -ach K1:VIS-EY_TM_OPLEV_SUM_OUT_DQ -bch K1:LSC-MICH_CTRL_CAL_OUT_DQ -s 1145621548.0 -d 600 -a 100 -f 103.5 -p 1
+python kagra_coupling_plot.py -ach K1:VIS-EY_TM_OPLEV_SUM_OUT_DQ -bch K1:LSC-MICH_CTRL_CAL_OUT_DQ -s 1145621548.0 -d 600 -a 100 -f0 10 -f1 500 -p 1
-Line 20:
+Line 31:
+## [EX] python kagra_coupling_plot.py -ach K1:VIS-EY_TM_OPLEV_SUM_OUT_DQ -bch K1:LSC-MICH_CTRL_CAL_OUT_DQ -s 1199756064.0 -d 600 -a 100 -f 103.5 -p 1
-Line 24:
+Line 36:
-## [5] -a       :: Averaging Number. default avg. is 100.
## [6] -f       :: Enter frequency to obtain amplitude value
## [7] -p       :: Enter 0 or 1 to determine whether to save the plot. 0 or None is not drawn. 1 is saved.
+## [5] -a       :: Enter Averaging Number, CAN NOT be used with -df at the same time
## [6] -df      :: Enter Frequency interval, df, CAN NOT be used with -a at the same time.
## [7] -f0      :: Enter start frequency to obtain amplitude value list
## [8] -f1      :: Enter end frequency to obtain amplitude value list
## [9] -p       :: Enter 0 or 1 to determine whether to save the plot. 0 or None is not drawn. 1 is saved.
-Line 33:
+Line 47:
-[detchar@k1det1 kagra_coupling_plot]$ python kagra_coupling_plot.py -ach K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ -bch K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ -s 1145621548.0 -d 736 -a 1 -f 103.5 -p 1
+[detchar@k1det1 kagra_coupling_plot]$ python kagra_coupling_plot.py -ach K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ -bch K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ -s 1145621548.0 -d 736 -df 0.1 -f0 10 -f1 100 -p 1
-Line 37:
+Line 51:
-### Plotting PSD Ratio of K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ to K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ ###
+---- Entered A Chanel : K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ
---- Entered B Chanel : K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ
---- Entered Start GPS time : 1145621548.0
---- Entered End GPS time : 1145622284.0
---- Entered Time duration : 736
---- Entered Average : 73.6
---- Entered df : 0.1 Hz
-Line 39:
+Line 59:
- Saved Interploation Ratio plot in ./plot/1145621548.0-736/K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ-to-K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ direction
+Checking each channel sample rate...
-Line 41:
+Line 61:
-    Ratio for selected frequency 103.5 Hz : 5.95881874822
+---- A Channel Sample Rate : 2048 Hz
PASS
---- B Channel Sample Rate : 2048 Hz
PASS


### Calcurate coupling function during 1145621548.0 to 1145622284.0 ###

    Saved Coupling plot in ./plot/1145621548.0-736/K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ-to-K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ

    Amplitude list saveed in ./plot/1145621548.0-736/K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ-to-K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ direction

    Do you want to print the saved list? Enter the Numer (Yes = 1 / No = 0) : 

１

10 Hz : 2.65813895339e-07
11 Hz : 4.86833812374e-07
12 Hz : 3.2527766507e-06 
.
.
.
97 Hz : 3.06364533742e-08
98 Hz : 3.87633960486e-08
99 Hz : 3.20635358264e-08
100 Hz : 4.65378668935e-08