== 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 ===

{{attachment:Ratio plot of two channels Report.001.png||width=1200}}
{{attachment:Ratio plot of two channels Report.002.png||width=1200}}
{{attachment:Ratio plot of two channels Report.003.png||width=1200}}
{{attachment:Ratio plot of two channels Report.004.png||width=1200}}
{{attachment:Ratio plot of two channels Report.005.png||width=1200}}
{{attachment:Ratio plot of two channels Report.006.png||width=1200}}
{{attachment:Ratio plot of two channels Report.007.png||width=1200}}
{{attachment:Ratio plot of two channels Report.008.png||width=1200}}
{{attachment:Ratio plot of two channels Report.009.png||width=1200}}