| 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"] |
| 164 |
|
os.chdir(dirname1) |
| 165 |
|
|
| 166 |
|
|
| 167 |
< |
for isp in xrange(7,13): |
| 167 |
> |
for isp in xrange(1,7): |
| 168 |
|
# dirname2="SP"+str(SPNum) |
| 169 |
|
# os.mkdir(dirname2) |
| 170 |
|
# os.chdir(dirname2) |
| 220 |
|
ewc_content[2].append(ewc[nn][2]) |
| 221 |
|
|
| 222 |
|
|
| 223 |
< |
aa=subplot(3,2,isp-6) |
| 224 |
< |
text(180,32.0,"Sector %d" %(isp-6)) |
| 225 |
< |
if isp<5+6: |
| 223 |
> |
aa=subplot(3,2,isp) |
| 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,20.0,"Sector %d" %(isp)) |
| 231 |
> |
if isp<5: |
| 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) |
| 238 |
|
|
| 239 |
< |
if isp>4+6: |
| 239 |
> |
if isp>4: |
| 240 |
|
xticks(ind_tick,time_stamp[::30],rotation=90) |
| 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: |
| 261 |
|
else: |
| 262 |
|
axis(ymax=SUP_YMAX) |
| 263 |
|
|
| 264 |
< |
suptitle('history of halo rate(/10s) in ME- trigger sectors, Sep 11') |
| 265 |
< |
savefig('HaloRateHistory_MEMinus.png') |
| 266 |
< |
savefig('HaloRateHistory_MEMinus.svg') |
| 267 |
< |
savefig('HaloRateHistory_MEMinus.eps') |
| 264 |
> |
suptitle('history of halo rate(/10s) in ME+ trigger sectors, Sep 11') |
| 265 |
> |
savefig('HaloRateHistory_MEPlus.png') |
| 266 |
> |
savefig('HaloRateHistory_MEPlus.svg') |
| 267 |
> |
savefig('HaloRateHistory_MEPlus.eps') |
| 268 |
|
|
| 269 |
|
|
| 270 |
|
|
