def compute(self):
self.rmsAnalysis = rmsIterator(self.fileInput,self.framesPerChunk)
self.rmsAnalysis.iterate(self.howManyChunks)
print "I'm trying to compute."
fig = plt.figure()
fig.set_size_inches(8,8)
fig.subplots_adjust(hspace=0.25)
ax1 = fig.add_subplot(221)
ax1.set_xlabel('Time ('+str("%.6f" % (self.framesPerChunk*self.howManyChunks*4096/19.6e6))+' Seconds total)')
ax1.set_ylabel('ADC Magnitude')
ax2 = fig.add_subplot(222)
ax2.set_xlabel('Time ('+str("%.6f" % (self.framesPerChunk*self.howManyChunks*4096/19.6e6))+' Seconds total)')
ax2.set_ylabel('ADC Magnitude')
ax3 = fig.add_subplot(223)
ax3.set_xlabel('Time ('+str("%.6f" % (self.framesPerChunk*self.howManyChunks*4096/19.6e6))+' Seconds total)')
ax3.set_ylabel('ADC Magnitude')
ax4 = fig.add_subplot(224)
ax4.set_xlabel('Time ('+str("%.6f" % (self.framesPerChunk*self.howManyChunks*4096/19.6e6))+' Seconds total)')
ax4.set_ylabel('ADC Magnitude')
ax1.plot(self.rmsAnalysis.returnDataT1P0[:],label='B'+str(self.rmsAnalysis.beam)+':T'+str("%.1f" % (self.rmsAnalysis.tune1/1e6))+':P0')
ax2.plot(self.rmsAnalysis.returnDataT1P1[:],label='B'+str(self.rmsAnalysis.beam)+':T'+str("%.1f" % (self.rmsAnalysis.tune1/1e6))+':P1')
ax3.plot(self.rmsAnalysis.returnDataT2P0[:],label='B'+str(self.rmsAnalysis.beam)+':T'+str("%.1f" % (self.rmsAnalysis.tune2/1e6))+':P0')
ax4.plot(self.rmsAnalysis.returnDataT2P1[:],label='B'+str(self.rmsAnalysis.beam)+':T'+str("%.1f" % (self.rmsAnalysis.tune2/1e6))+':P1')
ax1.legend(loc='best')
ax2.legend(loc='best')
ax3.legend(loc='best')
ax4.legend(loc='best')
fig.savefig(self.fileOutput)
def compute(self):
self.rmsAnalysis = rmsIterator(self.fileInput, self.framesPerChunk)
self.rmsAnalysis.iterate(self.howManyChunks)
print "I'm trying to compute."
fig = plt.figure()
fig.set_size_inches(8, 8)
fig.subplots_adjust(hspace=0.25)
ax1 = fig.add_subplot(221)
ax1.set_xlabel(
'Time (' +
str("%.6f" %
(self.framesPerChunk * self.howManyChunks * 4096 / 19.6e6)) +
' Seconds total)')
ax1.set_ylabel('ADC Magnitude')
ax2 = fig.add_subplot(222)
ax2.set_xlabel(
'Time (' +
str("%.6f" %
(self.framesPerChunk * self.howManyChunks * 4096 / 19.6e6)) +
' Seconds total)')
ax2.set_ylabel('ADC Magnitude')
ax3 = fig.add_subplot(223)
ax3.set_xlabel(
'Time (' +
str("%.6f" %
(self.framesPerChunk * self.howManyChunks * 4096 / 19.6e6)) +
' Seconds total)')
ax3.set_ylabel('ADC Magnitude')
ax4 = fig.add_subplot(224)
ax4.set_xlabel(
'Time (' +
str("%.6f" %
(self.framesPerChunk * self.howManyChunks * 4096 / 19.6e6)) +
' Seconds total)')
ax4.set_ylabel('ADC Magnitude')
ax1.plot(self.rmsAnalysis.returnDataT1P0[:],
label='B' + str(self.rmsAnalysis.beam) + ':T' +
str("%.1f" % (self.rmsAnalysis.tune1 / 1e6)) + ':P0')
ax2.plot(self.rmsAnalysis.returnDataT1P1[:],
label='B' + str(self.rmsAnalysis.beam) + ':T' +
str("%.1f" % (self.rmsAnalysis.tune1 / 1e6)) + ':P1')
ax3.plot(self.rmsAnalysis.returnDataT2P0[:],
label='B' + str(self.rmsAnalysis.beam) + ':T' +
str("%.1f" % (self.rmsAnalysis.tune2 / 1e6)) + ':P0')
ax4.plot(self.rmsAnalysis.returnDataT2P1[:],
label='B' + str(self.rmsAnalysis.beam) + ':T' +
str("%.1f" % (self.rmsAnalysis.tune2 / 1e6)) + ':P1')
ax1.legend(loc='best')
ax2.legend(loc='best')
ax3.legend(loc='best')
ax4.legend(loc='best')
fig.savefig(self.fileOutput)
def computeHist(self):
self.rmsAnalysis = rmsIterator(self.fileInput,self.framesPerChunk)
self.rmsAnalysis.iterate(self.howManyChunks)
print "<br>I'm trying to compute."
fig = plt.figure()
fig.set_size_inches(8,8)
fig.subplots_adjust(wspace=0.5)
fig.subplots_adjust(hspace=0.35)
#fig.tight_layout()
ax1 = fig.add_subplot(221)
fileAx1,bins1,trash1 = plt.hist(self.rmsAnalysis.returnDataT1P0[:],bins=8,histtype='step',normed=True,label = 'B'+str(self.rmsAnalysis.beam)+':T'+str("%.1f" % (self.rmsAnalysis.tune1/1e6))+":P0")
ax1.set_xticks(bins1)
ax1.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax1.get_xticklabels(), rotation=90)
ax1.set_xlabel('Step')
ax1.set_ylabel('Count')
ax1.legend(loc='best')
ax2 = fig.add_subplot(222)
fileAx2,bins2,trash2 = plt.hist(self.rmsAnalysis.returnDataT1P1[:],bins=8,histtype='step',normed=True,label = "B"+str(self.rmsAnalysis.beam)+":T"+str("%.1f" % (self.rmsAnalysis.tune1/1e6))+":P1")
ax2.set_xticks(bins2)
ax2.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax2.get_xticklabels(), rotation=90)
ax2.set_xlabel('Step')
ax2.set_ylabel('Count')
ax2.legend(loc='best')
ax3 = fig.add_subplot(223)
fileAx3,bins3,trash3 = plt.hist(self.rmsAnalysis.returnDataT2P0[:],bins=8,histtype='step',normed=True,label = "B"+str(self.rmsAnalysis.beam)+":T"+str("%.1f" % (self.rmsAnalysis.tune2/1e6))+":P0")
ax3.set_xticks(bins3)
ax3.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax3.get_xticklabels(), rotation=90)
ax3.set_xlabel('Step')
ax3.set_ylabel('Count')
ax3.legend(loc='best')
ax4 = fig.add_subplot(224)
fileAx4,bins4,trash4 = plt.hist(self.rmsAnalysis.returnDataT2P1[:],bins=8,histtype='step',normed=True,label = "B"+str(self.rmsAnalysis.beam)+":T"+str("%.1f" % (self.rmsAnalysis.tune2/1e6))+":P1")
ax4.set_xticks(bins4)
ax4.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax4.get_xticklabels(), rotation=90)
ax4.set_xlabel('Step')
ax4.set_ylabel('Count')
ax4.legend(loc='best')
fig.savefig(self.fileOutput+".png")
to_json(self.rmsAnalysis.jsonT1P0,self.fileOutput+"_t1p0.json")
to_json(self.rmsAnalysis.jsonT1P1,self.fileOutput+"_t1p1.json")
to_json(self.rmsAnalysis.jsonT2P0,self.fileOutput+"_t2p0.json")
to_json(self.rmsAnalysis.jsonT2P1,self.fileOutput+"_t2p1.json")
#a numpy array for 1 frame is 32,768 bytes. Fill CPU cache with frames for optimal performance. Saxon head node cache is 1024 KB.
#therefore 31.25 frames fit in the cache. Heavy lifting in thsi case is dot product of identical array, so assume ideal frame size of 15, with wiggle room built in.
chunks = 15
# how many chunks to plot
howManyChunks = (60*10)/0.003135
#fileArr is a list of four files, with the following order. HFON->HFOFF->LFON->LFOff
fileArr = [ '/u/data/leap/observations/056954_000126527', '/u/data/leap/observations/056954_000126528' ]
fig = plt.figure()
ax = fig.gca()
ax.set_xlabel('Chunk # ('+str(0.000209*chunks)+' Seconds)')
ax.set_ylabel('ADC Magnitude')
hf_on = rmsIterator(fileArr[0],chunks)
hf_on.iterate(howManyChunks)
hf_off = rmsIterator(fileArr[1],chunks)
hf_off.iterate(howManyChunks)
#lf_on = rmsIterator(fileArr[2],chunks)
#lf_on.iterate(howManyChunks)
#lf_off = rmsIterator(fileArr[3],chunks)
#lf_off.iterate(howManyChunks)
ax.plot(hf_on.returnData[:],label='HF On')
ax.plot(hf_off.returnData[:],label='HF Off')
#ax.plot(lf_on.returnData[:],label='LF On')
#ax.plot(lf_off.returnData[:],label='LF Off')
def computeHist(self):
self.rmsAnalysis = rmsIterator(self.fileInput, self.framesPerChunk)
self.rmsAnalysis.iterate(self.howManyChunks)
print "I'm trying to compute."
fig = plt.figure()
fig.set_size_inches(8, 8)
fig.subplots_adjust(hspace=0.25)
ax1 = fig.add_subplot(221)
fileAx1, bins1, trash1 = plt.hist(
self.rmsAnalysis.returnDataT1P0[:],
bins=8,
histtype='step',
normed=True,
label='B' + str(self.rmsAnalysis.beam) + ':T' +
str("%.1f" % (self.rmsAnalysis.tune1 / 1e6)) + ":P0")
ax1.set_xticks(bins1)
ax1.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax1.get_xticklabels(), rotation=90)
ax1.set_xlabel('Step')
ax1.set_ylabel('Count')
ax1.legend(loc='best')
ax2 = fig.add_subplot(222)
fileAx2, bins2, trash2 = plt.hist(
self.rmsAnalysis.returnDataT1P1[:],
bins=8,
histtype='step',
normed=True,
label="B" + str(self.rmsAnalysis.beam) + ":T" +
str("%.1f" % (self.rmsAnalysis.tune1 / 1e6)) + ":P1")
ax2.set_xticks(bins2)
ax2.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax2.get_xticklabels(), rotation=90)
ax2.set_xlabel('Step')
ax2.set_ylabel('Count')
ax2.legend(loc='best')
ax3 = fig.add_subplot(223)
fileAx3, bins3, trash3 = plt.hist(
self.rmsAnalysis.returnDataT2P0[:],
bins=8,
histtype='step',
normed=True,
label="B" + str(self.rmsAnalysis.beam) + ":T" +
str("%.1f" % (self.rmsAnalysis.tune2 / 1e6)) + ":P0")
ax3.set_xticks(bins3)
ax3.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax3.get_xticklabels(), rotation=90)
ax3.set_xlabel('Step')
ax3.set_ylabel('Count')
ax3.legend(loc='best')
ax4 = fig.add_subplot(224)
fileAx4, bins4, trash4 = plt.hist(
self.rmsAnalysis.returnDataT2P1[:],
bins=8,
histtype='step',
normed=True,
label="B" + str(self.rmsAnalysis.beam) + ":T" +
str("%.1f" % (self.rmsAnalysis.tune2 / 1e6)) + ":P1")
ax4.set_xticks(bins4)
ax4.xaxis.set_major_formatter(FormatStrFormatter('%0.2f'))
plt.setp(ax4.get_xticklabels(), rotation=90)
ax4.set_xlabel('Step')
ax4.set_ylabel('Count')
ax4.legend(loc='best')
fig.savefig(self.fileOutput)
