Differences between revisions 9 and 10
Revision 9 as of 2018-02-20 19:26:20
Size: 3445
Editor: PJJung
Comment:
Revision 10 as of 2018-02-20 19:28:38
Size: 3513
Editor: PJJung
Comment:
Deletions are marked like this. Additions are marked like this.
Line 67: Line 67:
### Plotting PSD Ratio of K1:PEM-EX_SEIS_NS_BLRMS_MID_OUT_DQ to K1:PEM-EY_SEIS_WE_BLRMS_MID_OUT_DQ ### ### Calcurate coupling function during 1145621548.0 to 1145622284.0 ###
Line 69: Line 69:
    coherence-plot Saved     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

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? (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 Ratio plot of two channels Report.002.png Ratio plot of two channels Report.003.png Ratio plot of two channels Report.004.png Ratio plot of two channels Report.005.png Ratio plot of two channels Report.006.png Ratio plot of two channels Report.007.png Ratio plot of two channels Report.008.png Ratio plot of two channels Report.009.png

KAGRA/Subgroups/DET/gwpy/coupling (last edited 2018-02-20 21:18:36 by PJJung)