32 |
|
SUP_YMAX=55.0 # maximum y range for suppressed plots;mainly for vprates |
33 |
|
SUP_YMAX1=500.0 # maximum y range for tot:stub |
34 |
|
SUP_YMAX2=35.0 # maximum y range for tot:track |
35 |
< |
SUP_YMAX3=40.0 # maximum y range for tot:halo |
35 |
> |
SUP_YMAX3=50.0 # maximum y range for tot:halo |
36 |
|
isSummary=False # is this rate counter summary? |
37 |
|
isTotalrate=False |
38 |
|
ratelabel=["Stubs","Tracks","Halos"] |
221 |
|
|
222 |
|
|
223 |
|
aa=subplot(3,2,isp-6) |
224 |
+ |
# this equivalent to the best; but not so simple |
225 |
+ |
# yscale('log') |
226 |
+ |
# this doesn't scale well when using logscale |
227 |
+ |
# myax=gca() |
228 |
+ |
# myax.set_yscale('log') |
229 |
+ |
|
230 |
|
text(180,32.0,"Sector %d" %(isp-6)) |
231 |
|
if isp<5+6: |
232 |
|
subplots_adjust(hspace=0.0) |
233 |
< |
plot(ewc_content[2],linestyle='steps') |
233 |
> |
## this turn to be the best option for logscale |
234 |
> |
semilogy(ewc_content[2],linestyle='steps') |
235 |
|
ind_tick=xrange(0,len(time_stamp),30) |
236 |
|
# axvspan(50.0,100.0,facecolor='y',alpha=0.4) |
237 |
|
axvspan(105.0,175.0,facecolor='y',alpha=0.3) |
241 |
|
else: |
242 |
|
xticks([],[]) |
243 |
|
# xlabel('time stamp(/10s), from %s to %s' %(l_logtime[index_startstamp],l_logtime[-1])) |
244 |
< |
ylabel("halo rate(Hz)") |
244 |
> |
# due to logscale, remove the label from right subplots |
245 |
> |
if isp%2: |
246 |
> |
ylabel("halo rate(Hz)") |
247 |
|
linklabel='' |
248 |
|
linklabel="link_M%d" %(3) |
249 |
+ |
|
250 |
|
# savefig(pat_ewc+'_'+linklabel+'.png') |
251 |
|
|
252 |
|
# if isTotalrate: |