Skip to content
Snippets Groups Projects
Commit 4cac23fe authored by Gregory Ashton's avatar Gregory Ashton
Browse files

Update to showing one and 2 glitch distribution

parent cb7c547e
No related branches found
No related tags found
No related merge requests found
...@@ -4,7 +4,7 @@ ...@@ -4,7 +4,7 @@
export PATH="/home/gregory.ashton/anaconda2/bin:$PATH" export PATH="/home/gregory.ashton/anaconda2/bin:$PATH"
export MPLCONFIGDIR=/home/gregory.ashton/.config/matplotlib export MPLCONFIGDIR=/home/gregory.ashton/.config/matplotlib
for ((n=0;n<90;n++)) for ((n=0;n<5;n++))
do do
/home/gregory.ashton/anaconda2/bin/python generate_data.py "$1" /local/user/gregory.ashton --no-interactive /home/gregory.ashton/anaconda2/bin/python generate_data.py "$1" /local/user/gregory.ashton --no-interactive
done done
......
...@@ -7,6 +7,7 @@ from tqdm import tqdm ...@@ -7,6 +7,7 @@ from tqdm import tqdm
from oct2py import octave from oct2py import octave
import glob import glob
from scipy.stats import rv_continuous, chi2 from scipy.stats import rv_continuous, chi2
from latex_macro_generator import write_to_macro
filenames = glob.glob("CollectedOutput/*.txt") filenames = glob.glob("CollectedOutput/*.txt")
...@@ -21,21 +22,27 @@ for fn in filenames: ...@@ -21,21 +22,27 @@ for fn in filenames:
df['CLUSTER_ID'] = fn.split('_')[1] df['CLUSTER_ID'] = fn.split('_')[1]
df_list.append(df) df_list.append(df)
df = pd.concat(df_list) df = pd.concat(df_list)
print 'Number of samples = ', len(df)
print 'Max twoF', df.twoF.max()
fig, ax = plt.subplots() fig, ax = plt.subplots()
for ng in np.unique(df.nglitches.values): for ng in np.unique(df.nglitches.values):
ax.hist(df[df.nglitches==ng].twoF, bins=20, histtype='step', normed=True, print 'ng={}'.format(ng)
Nsamples = len(df[df.nglitches == ng])
MaxtwoF = df[df.nglitches == ng].twoF.max()
print 'Number of samples = ', Nsamples
print 'Max twoF', MaxtwoF
ax.hist(df[df.nglitches == ng].twoF, bins=40, histtype='step', normed=True,
linewidth=1, label='$N_\mathrm{{glitches}}={}$'.format(ng)) linewidth=1, label='$N_\mathrm{{glitches}}={}$'.format(ng))
ax.set_xlabel('$\widehat{2\mathcal{F}}$ for 1 glitch') write_to_macro('DirectedMC{}GlitchNoiseOnlyMaximum'.format(ng),
'{:1.1f}'.format(MaxtwoF), '../macros.tex')
write_to_macro('DirectedMC{}GlitchNoiseN'.format(ng),
'{:1.0f}'.format(Nsamples), '../macros.tex')
ax.set_xlabel('$\widehat{2\mathcal{F}}$')
ax.set_xlim(0, 90) ax.set_xlim(0, 90)
ax.legend(frameon=False, fontsize=6) ax.legend(frameon=False, fontsize=6)
fig.tight_layout() fig.tight_layout()
fig.savefig('single_glitch_noise_twoF_histogram.png') fig.savefig('glitch_noise_twoF_histogram.png')
#from latex_macro_generator import write_to_macro
#write_to_macro('DirectedMCNoiseOnlyMaximum', '{:1.1f}'.format(np.max(df.twoF)),
# '../macros.tex')
#write_to_macro('DirectedMCNoiseN', len(df), '../macros.tex')
...@@ -2,11 +2,11 @@ Executable= SingleGlitchMCNoiseOnly.sh ...@@ -2,11 +2,11 @@ Executable= SingleGlitchMCNoiseOnly.sh
Arguments= $(Cluster)_$(Process) Arguments= $(Cluster)_$(Process)
Universe=vanilla Universe=vanilla
Input=/dev/null Input=/dev/null
accounting_group = ligo.dev.o2.cw.explore.test accounting_group=ligo.dev.o1.cw.directedisolatedother.semicoherent
Output=CollectedOutput/out.$(Cluster).$(Process) Output=CollectedOutput/out.$(Cluster).$(Process)
Error=CollectedOutput/err.$(Cluster).$(Process) Error=CollectedOutput/err.$(Cluster).$(Process)
Log=CollectedOutput/log.$(Cluster).$(Process) Log=CollectedOutput/log.$(Cluster).$(Process)
request_cpus = 1 request_cpus = 1
request_memory = 16 GB request_memory = 16 GB
Queue 100 Queue 170
Paper/glitch_noise_twoF_histogram.png

21.6 KiB

...@@ -12,8 +12,13 @@ ...@@ -12,8 +12,13 @@
\def\BasicExamplehzero{1.0{\times}10^{-24}} \def\BasicExamplehzero{1.0{\times}10^{-24}}
\def\BasicExamplenburn{50.0} \def\BasicExamplenburn{50.0}
\def\BasicExamplenprod{50.0} \def\BasicExamplenprod{50.0}
\def\DirectedMC{}GlitchNoiseN{970.0}
\def\DirectedMCNoiseN{1.0{\times}10^{4}} \def\DirectedMCNoiseN{1.0{\times}10^{4}}
\def\DirectedMCNoiseOnlyMaximum{52.4} \def\DirectedMCNoiseOnlyMaximum{52.4}
\def\DirectedMConeGlitchNoiseN{10000}
\def\DirectedMConeGlitchNoiseOnlyMaximum{82.6}
\def\DirectedMCtwoGlitchNoiseN{970}
\def\DirectedMCtwoGlitchNoiseOnlyMaximum{80.8}
\def\SingleGlitchDepth{10.0} \def\SingleGlitchDepth{10.0}
\def\SingleGlitchFCMismatch{0.7} \def\SingleGlitchFCMismatch{0.7}
\def\SingleGlitchFCtwoF{718.5} \def\SingleGlitchFCtwoF{718.5}
......
Paper/single_glitch_noise_twoF_histogram.png

21.1 KiB

0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment