;$Id: pspecm_scl5_hel_sig.pro,v 1.6 2020/11/01 15:10:52 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=12,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; ; sig grav1mm pol1 pol2 hel1 hel2 run ;-1.000 2.64e-06 -0.88109 -0.41315 -0.97212 -0.24926 F1152sigm1a ;-0.100 1.88e-06 -0.16662 -0.03742 -0.31665 -0.02251 F1152sigm01a ; 0.100 1.90e-06 0.22322 0.04499 0.31759 0.02252 F1152sig01a ; 1.000 2.67e-06 0.89864 0.39000 0.96944 0.26901 F1152d2_sig1_t11_M2c_double ; dir0='F1152sigm1a' dir1='F1152sig07c' dir2='F1152sig05a' ;dir2='M1152e_exp6k4_M4b_rerun' dir3='F1152sig03a' ;dir3='M1152e_exp6k4_M4b_rerun_sig05' dir4='F1152sig001a' ;dir4='M1152e_exp6k4_M4b_rerun_sig01' ; yr_pspecm_scl3=[2d-12,5d-9] yr_pspecm_scl3=[2d-20,5d-8] ; thick4=5 siz=1.1 !p.charsize=1.8 !x.margin=[8.0,1.3] !y.margin=[3.2,0.2] !p.multi=0 fac2=2. half='!s!u 1!n!r!s-!r!d 2!n' specm='/specm.sav' label='' ; ; determine normalization ; default,maglabel,'mag' default,plot_polarization,1 ;ytickf='logticks_exp' yr=yr_pspecm_scl3 ; Hstar=2.066d10 a0=1.2545d15 freq=Hstar/(2.*!pi*a0) H0=3.24d-18 rhocrit_fac=3./(8.*!pi) fac_h=sqrt(3./2.)*H0/!pi fac=(Hstar/H0)^2/rhocrit_fac*(1./a0)^4 facGW=1./(16.*!pi) xr=[2.4e-4,1.2e-1] circ_sym,0.7,1 energyfac=1e2 mixed=0. ; !p.title='!6';+run !x.title='!8f!6 [Hz]' !y.title='!6!8h!6!s!d0!n!r!u2!n!7X!6!dGW!n(!8f!6)' ; restore,'../'+dir0+specm grav_tot=grav1m+grah1m/tm^2+mixed ghel_tot=grav2m+grah2m/tm^2+mixed OmegaGW=facGW*fac*k*grav_tot XiGW=facGW*fac*k*ghel_tot plot_oo,freq*k,energyfac*OmegaGW,xr=xr,yr=yr,ytickf=ytickf ; restore,'../'+dir1+specm grav_tot=grav1m+grah1m/tm^2+mixed ghel_tot=grav2m+grah2m/tm^2+mixed OmegaGW=facGW*fac*k*grav_tot XiGW=facGW*fac*k*ghel_tot oplot,freq*k,energyfac*OmegaGW,col=55 ; restore,'../'+dir2+specm grav_tot=grav1m+grah1m/tm^2+mixed OmegaGW=facGW*fac*k*grav_tot loadct,6 oplot,freq*k,energyfac*OmegaGW,col=122 loadct,5 ; restore,'../'+dir3+specm grav_tot=grav1m+grah1m/tm^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,energyfac*OmegaGW,col=155 ; restore,'../'+dir4+specm grav_tot=grav1m+grah1m/tm^2+mixed OmegaGW=facGW*fac*k*grav_tot oplot,freq*k,energyfac*OmegaGW,col=122 ; ; Caprini+16, Cornish, Maggiore00 points ; loadct,6 dir='~/tex/roper/PencilGW/data/' file=dir+'Cornish.rtf' c1=rtable(file,2,head=1) ; file=dir+'Cornish20Yrs.dat' c2=rtable(file,2,head=1) ; file=dir+'Maggiore.rtf' m1=rtable(file,2,head=1) ; file=dir+'Caprini-C3.rtf' a3=rtable(file,2,head=1) ; file=dir+'Caprini-C1.rtf' a1=rtable(file,2,head=1) ; file=dir+'Robson_Cornish.rtf' n1=rtable(file,2,head=1) oplot,10^n1[0,*],((10^n1[1,*]/fac_h)*10^n1[0,*])^2,li=3,col=122,thick=thick4 loadct,5 ; ; overplot ; ;xx=[5e-3,3.7e-2] & oplot,li=3,xx,1e-15/xx^(8./3.),thick=6 ;xx=[5e-3,3.7e-2] & oplot,li=3,xx,1e-15/xx^(7./3.),thick=6 ;xx=[5e-3,3.7e-2] & oplot,li=3,xx,1e-14/xx^(7./3.),thick=6 ; xx=[3e-4,8e-4] & oplot,li=3,xx,1e-6*xx ;xx=[7e-3,1.7e-2] & oplot,li=3,xx,1e-19/xx^(11./3.) ;xx=[7e-3,1.7e-2] & oplot,li=3,xx,5e-24/xx^(14./3.) xx=[7e-3,2.0e-2] & oplot,li=3,xx,2.2e-20/xx^4. xx=[7e-3,1.7e-2] & oplot,li=3,xx,5.3e-25/xx^5. xyouts,4.0e-4,2.0e-9,'!9A!8f!6',siz=siz*1.3 ;xyouts,1.7e-2,3.4e-13,'!9A!8f!6!u-11/3!n',siz=siz*1.5 ;xyouts,0.8e-2,6.4e-17,'!9A!8f!6!u-14/3!n',siz=siz*1.5 xyouts,1.7e-2,6.4e-13,'!9A!8f!6!u-4!n',siz=siz*1.5 xyouts,0.8e-2,1.0e-16,'!9A!8f!6!u-5!n',siz=siz*1.5 ; ; legends ; xx=4e-4 & dx=xx*2 legend,xx,dx,1e-14,0,'!7r!6=0.01',col=122 legend,xx,dx,1e-15,0,'!7r!6=0.3',col=155 & loadct,6 legend,xx,dx,1e-16,0,'!7r!6=0.5',col=122 & loadct,5 legend,xx,dx,1e-17,0,'!7r!6=0.7',col=55 legend,xx,dx,1e-18,0,'!7r!6=1' ; print,'$mv idl.ps ~/tex/roper/helical/fig/pspecm_scl5_hel_sig.ps' END