those scripts are to plot Aggregate plot for coolpup. the script for myself is in ~/script/plot_pup1.sh and ~/script/plot_pup.py
python ~/script/plot_pup1.sh cool1 cool2 /broad/boxialab/shawn/projects/J/self_data/ChIP-seq/12-22/1.mapping/HCT116_J_chip_85_87_sort_dedup_q005_peaks.narrowPeak 10-9 10-10 0.6 1.4 -0.2 0.2
python work.py --input aa.txt,bb.txt --vmin 0.6 --vmax 1.4 --vmin_diff -0.2 --vmax_diff 0.2
python work.py --input HEK_WT1_merge_on_hek_H7_J_on_H7_merge_q0.05_CTCF_cobound_descend_00_5k.txt --vmin 0.6 --vmax 1.4