from DTCurrentPlot import * import numpy as np # import data fill1 = DTCurrentPlot('fills/4364/') fill2 = DTCurrentPlot('fills/4381/') # combine luminosity and currents fill2.data.luminosity = np.vstack((fill1.data.luminosity.reshape((-1,1)), fill2.data.luminosity.reshape((-1,1)))).flatten() fill2.data.currents = np.ma.concatenate([fill1.data.currents, fill2.data.currents], axis=6) fill2.data.background = np.ma.concatenate([fill1.data.background, fill2.data.background], axis=6) # rename fill number (fill1+fill2) for plots fill2.data.fill = '{} + {}'.format(fill1.data.fill, fill2.data.fill) # save preconfigured plots fill2.plot_data()
from DTCurrentPlot import * # import data plot = DTCurrentPlot('fills/3992/') # draw plots plot.draw(y='current', x='luminosity', series='superlayer', wheel=1, station=4, sector=4, format=None) plot.draw(y='slope', x='wheel', series='superlayer', station=4, format=None) plot.draw(y='maxcurrent', x='sector', station=4, format=None) plot.draw_slope_2d(station=4, format=None) # save in pdf plot.draw(y='slope', x='wheel', format='svg') # save preconfigured plots plot.plot_data()
from DTCurrentPlot import * import numpy as np # import data fill1 = DTCurrentPlot('fills/4364/') fill2 = DTCurrentPlot('fills/4381/') # combine luminosity and currents fill2.data.luminosity = np.vstack((fill1.data.luminosity.reshape( (-1, 1)), fill2.data.luminosity.reshape((-1, 1)))).flatten() fill2.data.currents = np.ma.concatenate( [fill1.data.currents, fill2.data.currents], axis=6) fill2.data.background = np.ma.concatenate( [fill1.data.background, fill2.data.background], axis=6) # rename fill number (fill1+fill2) for plots fill2.data.fill = '{} + {}'.format(fill1.data.fill, fill2.data.fill) # save preconfigured plots fill2.plot_data()