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)