| 164 |
|
os.chdir(dirname1) |
| 165 |
|
|
| 166 |
|
|
| 167 |
< |
for isp in xrange(1,7): |
| 167 |
> |
for isp in xrange(7,13): |
| 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) |
| 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,20.0,"Sector %d" %(isp)) |
| 231 |
< |
if isp<5: |
| 230 |
> |
text(180,20.0,"Sector %d" %(isp-6)) |
| 231 |
> |
if isp<5+6: |
| 232 |
|
subplots_adjust(hspace=0.0) |
| 233 |
|
## this turn to be the best option for logscale |
| 234 |
|
semilogy(ewc_content[2],linestyle='steps') |
| 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: |
| 239 |
> |
if isp>4+6: |
| 240 |
|
xticks(ind_tick,time_stamp[::30],rotation=90) |
| 241 |
|
else: |
| 242 |
|
xticks([],[]) |
| 261 |
|
else: |
| 262 |
|
axis(ymax=SUP_YMAX) |
| 263 |
|
|
| 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') |
| 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') |
| 268 |
|
|
| 269 |
|
|
| 270 |
|
|
