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

Differences between revisions 8 and 9

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


### Plotting PSD Ratio of K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ to K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ ###

    coherence-plot Saved

    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? (Y / N)

Y

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 8 as of 2018-02-15 16:20:39 → 
  Size: 2437
  Editor: PJJung
  Comment:
+   ← Revision 9 as of 2018-02-20 19:26:20 → ⇥
  Size: 3445
  Editor: PJJung
  Comment:
-Deletions are marked like this.
+Additions are marked like this.
 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 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 34:
+Line 35:
-## [4] -d       :: Enter time from star GPS time duration
## [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.
+## [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.
-Line 44:
+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 48:
+Line 51:
+---- 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
-Line 50:
+Line 69:
- 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
+    coherence-plot Saved
-Line 52:
+Line 71:
-    Ratio for selected frequency 103.5 Hz : 5.95881874822
+    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? (Y / N)

Y

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