from PLOT.PLOT import PLOT from RetrieveData.RetrieveDataForPlot import RetrieveDataForPlot from RetrieveData.Retriever import Retriever __author__ = 'mmateja' basic_path = "/home/mmateja/PycharmProjects/OperationPlots/" #nameFileEnd = "Mu" startfillnumber = 3819 endfillnumber = 4000 d = Retriever( basic_path, "lumi_LHCb", startfillnumber, endfillnumber) data = d.retrieve("LPC") #LPC, CONDITION r = RetrieveDataForPlot(data) datax,datay=r.retrieveData("fillnumber","max_lumi") for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): #-key2,0iteritems ze wzgledu na typy danych plot=PLOT(value1, value2) plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015") #key1 plot.setxlabel("LHC Fillnumber") plot.setylabel("Peak Inst Luminosity (10^24/cm^2*s^1)") plot.draw() datax,datay=r.retrieveData("fillnumber","average_lumi") for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): #-key2,0iteritems ze wzgledu na typy danych plot=PLOT(value1, value2) plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015") #key1 plot.setxlabel("LHC Fillnumber") plot.setylabel("Peak Inst Luminosity (10^24/cm^2*s^1)") plot.draw() #plot.savefig("OUTPUTPLOTS/2015PeakLumiFill.png") plot.show()
from PLOT.PLOT import PLOT, dateformat import numpy as np from datetime import datetime basic_path = "/home/mmateja/PycharmProjects/OperationPlots/" #nameFileEnd = "Mu" startfillnumber = 3819 endfillnumber = 4000 xaxis = "fillnumber" #TIME_DATE, FILLNUMBER yaxis = "max_lumi" d=PlotLPC(xaxis, yaxis, basic_path, "lumi_LHCb", startfillnumber, endfillnumber) #d=PlotCondition("time_year","value", basic_path, nameFileEnd) for (key1, value1),(key2, value2) in zip(d.datax.iteritems(), d.datay.iteritems()): #-key2,0iteritems ze wzgledu na typy danych #value1 = [datetime.fromtimestamp(v1) for v1 in value1] if key1 != key2 and ("time_sec" and "lumi" in xaxis and yaxis): print "ERROR DIFFERENT KEYS!!!!!!!" if "time_date" in xaxis and yaxis: value1 = dateformat(value1) plot=PLOT(value1, value2) plot.setTitle(key1) #key1 plot.setDate() plot.draw() else: plot=PLOT(value1, value2) plot.setTitle(key1) #key1 plot.draw() plot.show() raw_input("Press enter to continue")
startfillnumber = 3819 endfillnumber = 5000 # Start by retrieving data for LPC/lumi_LHCb files in the range of fillnumbers print 'Retrieving data from lumi_LHCb between fill# ' + str(startfillnumber) + ' and fill# ' + str(endfillnumber) d = Retriever( basic_path, "lumi_LHCb", startfillnumber, endfillnumber) lumi = d.retrieve("LPC") #LPC, CONDITION r = RetrieveDataForPlot(lumi) xaxis = "fillnumber" yaxis = "max_lumi_lumi" print 'Doing ' + yaxis + ' vs ' + xaxis datax,datay = r.retrieveData(xaxis,yaxis) for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): plot=PLOT() plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015") plot.setxlabel("LHC Fillnumber") plot.setylabel("Peak Inst Luminosity (Hz/ub)") plot.draw(value1, value2, 'b', 6) plot.savefig("OUTPUTPLOTS/2015PeakLumiFill.png") xaxis = "time_date" yaxis = "max_lumi_lumi" print 'Doing ' + yaxis + ' vs ' + xaxis datax,datay = r.retrieveData(xaxis,yaxis) for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): value1 = dateformat(value1) plot=PLOT() plot.setDate() plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015")
# plot.setTitle(key1) #key1 # plot.setDate() # plot.draw() # else: # plot=PLOT(value1, value2) # plot.setTitle(key1) #key1 # plot.draw() # ###################################################################################################################### d = Retriever( basic_path, "lumi_LHCb", startfillnumber, endfillnumber) data = d.retrieve("LPC") #LPC, CONDITION r = RetrieveDataForPlot(data) datax,datay=r.retrieveData("fillnumber","max_lumi") for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): #-key2,0iteritems ze wzgledu na typy danych plot=PLOT(value1, value2) plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015") #key1 plot.setxlabel("LHC Fillnumber") plot.setylabel("Peak Inst Luminosity (10^24/cm^2*s^1)") plot.draw() datax,datay=r.retrieveData("fillnumber","average_lumi") for (key1, value1),(key2, value2) in zip(datax.iteritems(), datay.iteritems()): #-key2,0iteritems ze wzgledu na typy danych plot=PLOT(value1, value2) plot.setTitle("LHCb Peak Instantaneous Lumi at p-p 6.5 TeV in 2015") #key1 plot.setxlabel("LHC Fillnumber") plot.setylabel("Peak Inst Luminosity (10^24/cm^2*s^1)") plot.draw() xaxis = "time_date" #TIME_DATE, FILLNUMBER yaxis = "average_lumi"