drltools.py 25.6 KB
Newer Older
Yan's avatar
Yan committed
1
2
3
4
5
#!/usr/bin/env python3
from matplotlib.backends.backend_qt5agg import FigureCanvasQTAgg as FigureCanvas
from matplotlib.figure import Figure
from PyQt5 import QtCore
from PyQt5 import QtWidgets
6
from PyQt5 import QtGui
Yan's avatar
Yan committed
7
from PyQt5 import QtPrintSupport
Yan's avatar
Yan committed
8
9
import matplotlib
import numpy as np
Yan's avatar
Yan committed
10
11
import prasopes.datatools as dt
import prasopes.graphtools as gt
Yan's avatar
Yan committed
12
import prasopes.filetools as ft
Yan's avatar
Yan committed
13
import os.path
Yan's avatar
Yan committed
14
import logging
Yan's avatar
Yan committed
15
16
matplotlib.use("Qt5Agg")

Yan's avatar
Yan committed
17

Yan's avatar
Yan committed
18
19
20
logger = logging.getLogger('drlLogger')


Yan's avatar
Yan committed
21
class HBar(QtWidgets.QFrame):
22
    """horizontal bar class"""
Yan's avatar
Yan committed
23
24
25
26
27
28
    def __init__(self):
        super(HBar, self).__init__()
        self._main = QtWidgets.QWidget()
        self.setFrameShape(QtWidgets.QFrame.HLine)


Yan's avatar
Yan committed
29
def floatize(table, row, column, nonneg=True):
30
31
32
33
    """grabs the tableWidgetItem and transforms its text safely to
    float, if the text is not acceptable as float, returns zero"""
    imptext = table.item(row, column).text()
    validator = QtGui.QDoubleValidator()
Yan's avatar
Yan committed
34
35
    if nonneg == True:
        validator.setBottom(0)
36
37
38
39
40
41
42
43
44
    status = validator.validate(imptext, 0)[0]
    if status == QtGui.QValidator.Acceptable:
        outfloat = float(imptext)
    else:
        outfloat = 0
    return outfloat


def update_profile(pt, row, dataset):
Yan's avatar
Yan committed
45
    """parent table profile spectrum updating procedure"""
Yan's avatar
Yan committed
46
    # Dont do anything to graph when the spectrum is not populated
Yan's avatar
Yan committed
47
    logger.debug("updating parent table row {} profile".format(row))
Yan's avatar
Yan committed
48
    if isinstance(dataset['masses'], type(None)):
49
50
        return

51
    spectrum = pt.cellWidget(row, 3).figure.get_axes()[0]
Yan's avatar
Yan committed
52
    masses = dataset['masses']
53
54
    massargs = dt.argsubselect(
        masses, floatize(pt, row, 1), floatize(pt, row, 2))
55
    yshape = np.mean(dataset['matrix'], axis=0)
Yan's avatar
Yan committed
56
    spectrum.clear()
57
58
59
    spectrum.plot(masses, yshape, ':', color='gray')
    spectrum.plot(masses[massargs], yshape[massargs], 'r')
    xex = max((masses[massargs[-1]]-masses[massargs[0]])*0.25, 0.20)
Yan's avatar
Yan committed
60
61
62
    spectrum.set_xlim(masses[massargs[0]]-xex,
                      masses[massargs[-1]]+xex)
    ymax = max(yshape[massargs])
63
64
65
    spectrum.set_ylim(ymax*-0.1, ymax*1.2)
    spectrum.figure.canvas.draw()

Yan's avatar
Yan committed
66

3Yan's avatar
3Yan committed
67
68
69
def get_intensity(row, ds, drls):
    startm = floatize(drls['pt'], row, 1)
    endm = floatize(drls['pt'], row, 2)
Yan's avatar
Yan committed
70
    massargs = dt.argsubselect(ds['masses'], startm, endm)
Yan's avatar
Yan committed
71
72
73
74
    #prevent division by 0
    intensity = np.divide(np.sum( ds['matrix'].T[massargs].T, axis=1),
            np.clip(ds['chrom_dat'][1], np.nextafter(
                0, 1, dtype=np.float32), None))
Yan's avatar
Yan committed
75
76
77
    return intensity


3Yan's avatar
3Yan committed
78
def get_daughterset(ds, drls):
Yan's avatar
Yan committed
79
80
    """Fuction to acquire the curves of the daugher ions"""
    logger.info("getting set of the daughter ions")
Yan's avatar
Yan committed
81
    # TODO: write a less resources demanding function - probably "per-line"
Yan's avatar
Yan committed
82
    names = []
3Yan's avatar
3Yan committed
83
    times = ds['chrom_dat'][0, :] - drls['tshift'].value()
Yan's avatar
Yan committed
84
    intensities = []
3Yan's avatar
3Yan committed
85
86
87
    for row in range(drls['dt'].rowCount()):
        if drls['dt'].cellWidget(row, 0).checkState() == 2:
            intensity = get_intensity(row, ds, drls)
3Yan's avatar
3Yan committed
88
89
90
91
            corlist = []
            for i in range(drls['cors']):
                cor = drls['dt'].cellWidget(row, 1+i*2).currentIndex() - 1
                if cor not in (-2, -1):
Yan's avatar
Yan committed
92
                    factor = floatize(drls['dt'], row, 2+i*2, False)
3Yan's avatar
3Yan committed
93
94
95
96
97
98
                    correction = get_intensity(cor, ds, drls) * factor
                    intensity = intensity - correction
                    corlist.append("{} * {}".format(
                        drls['dt'].item(row,2+i*2).text(),
                        drls['dt'].cellWidget(row,1+i*2).currentText()))
            cortext = " + ".join(corlist)
Yan's avatar
Yan committed
99
            intensities.append(intensity)
3Yan's avatar
3Yan committed
100
            names.append("{} - ({})".format(
3Yan's avatar
3Yan committed
101
                drls['dt'].item(row, 0).text(),
3Yan's avatar
3Yan committed
102
                cortext))
Yan's avatar
Yan committed
103
104
105
    return names, times, intensities


3Yan's avatar
3Yan committed
106
def get_parentset(ds, drls):
107
108
109
110
111
    names = []
    times = ds['chrom_dat'][0, :]
    intensities = []
    # TODO: resolve intensities trouble
    rowlist = []
3Yan's avatar
3Yan committed
112
113
    for row in range(drls['dt'].rowCount()):
        if drls['dt'].cellWidget(row, 0).checkState() == 2:
114
            rowlist.append(row)
3Yan's avatar
3Yan committed
115
116
            if drls['dt'].cellWidget(row, 1).currentIndex() > 0\
                    and floatize(drls['dt'], row, 2) != 0:
117
                rowlist.append(
3Yan's avatar
3Yan committed
118
                    drls['dt'].cellWidget(row, 1).currentIndex()-1)
119
    for row in set(rowlist):
3Yan's avatar
3Yan committed
120
            intensity = get_intensity(row, ds, drls)
121
            intensities.append(intensity)
3Yan's avatar
3Yan committed
122
            names.append(drls['dt'].item(row, 0).text())
123
124
125
    return names, times, intensities


3Yan's avatar
3Yan committed
126
def update_drlspectrum(ds, drls, drlspectrum):
Yan's avatar
Yan committed
127
128
    """Generic DRL spectrum updating procedure"""
    logger.info("updating DRL spectrum")
129
    # Dont do anything when the dataset is not populated
Yan's avatar
Yan committed
130
    if isinstance(ds['masses'], type(None)):
Yan's avatar
Yan committed
131
        return
132
    colors = np.array([[0, 0, 0], [255, 0, 0], [0, 255, 0], [0, 0, 255],
133
134
135
                       [0, 200, 255], [255, 200, 0], [255, 100, 0],
                       [200, 50, 0], [255, 0, 200], [0, 100, 0],
                       [0, 100, 255], [100, 100, 100]])
Yan's avatar
Yan committed
136

137
    # TODO: write a less resources demanding function
3Yan's avatar
3Yan committed
138
    names, times, intensities = get_daughterset(ds, drls)
Yan's avatar
Yan committed
139
140
    for i in range(len(drlspectrum.lines)):
        drlspectrum.lines[0].remove()
141

3Yan's avatar
3Yan committed
142
    drlspectrum.axvline(0, color="#FF000088", linestyle=":")
3Yan's avatar
3Yan committed
143

144
    i = 0
3Yan's avatar
3Yan committed
145
146
147
148
    for row in range(drls['dt'].rowCount()):
        if drls['dt'].cellWidget(row, 0).checkState() == 2:
            drls['dt'].blockSignals(True)
            drls['dt'].item(row, 0).setBackground(QtGui.QBrush(
149
                QtGui.QColor(*colors[row % len(colors)], alpha=50)))
3Yan's avatar
3Yan committed
150
151
            drls['dt'].blockSignals(False)
            label = " {}".format(drls['pt'].item(row, 0).text())
3Yan's avatar
3Yan committed
152
153
            intensity = intensities[i]
            if drls['rel'].checkState() == 2:
Yan's avatar
Yan committed
154
155
156
157
                #Do not divide by 0
                intensity = np.divide(intensity, np.clip(np.sum(
                    intensities, 0), np.nextafter(
                        0, 1, dtype=np.float32), None))
3Yan's avatar
3Yan committed
158
            drlspectrum.plot(times, intensity, label=label,
159
160
                             color=(colors[row % len(colors)] / 255))
            i += 1
161
        else:
3Yan's avatar
3Yan committed
162
            drls['dt'].item(row, 0).setBackground(QtGui.QBrush())
163

164
    if len(names) != 0:
3Yan's avatar
3Yan committed
165
166
167
168
        gmax = np.amax(intensities)
        if drls['rel'].checkState() == 2:
            gmax = 1
        drlspectrum.set_ylim(top=gmax*1.1, bottom=gmax*-0.01)
169
        drlspectrum.legend(loc=2)
170
171
    drlspectrum.figure.canvas.draw()

Yan's avatar
Yan committed
172

Yan's avatar
Yan committed
173
def gettableitemlist(ptable):
174
    ion_list = []
Yan's avatar
Yan committed
175
176
177
    for row in range(ptable.rowCount()):
        text = []
        for i in range(3):
178
179
            if not isinstance(ptable.item(row, i), type(None)):
                frg = ptable.item(row, i).text()
Yan's avatar
Yan committed
180
181
182
183
184
            else:
                frg = ""
            text.append(frg)
        line = "{} ({}-{})".format(*text)
        ion_list.append(line)
185
186
187
    return ion_list


3Yan's avatar
3Yan committed
188
def update_corrfors(drls):
189
    """update corrections selection layout of the daughter table"""
3Yan's avatar
3Yan committed
190
191
    ionlist = gettableitemlist(drls['pt'])
    for row in range(drls['dt'].rowCount()):
3Yan's avatar
3Yan committed
192
193
194
195
196
197
198
199
200
        for i in range(drls['cors']):
            corfor = drls['dt'].cellWidget(row, 1+i*2)
            index = corfor.currentIndex()
            corfor.blockSignals(True)
            corfor.clear()
            corfor.addItem("")
            corfor.addItems(ionlist)
            corfor.setCurrentIndex(index)
            corfor.blockSignals(False)
Yan's avatar
Yan committed
201

202

3Yan's avatar
3Yan committed
203
def ptable_changed(row, column, ds, drls, drlspectrum):
204
    """routine called by change of the ptable spectra"""
Yan's avatar
Yan committed
205
    logger.debug("ptable changed routine called")
3Yan's avatar
3Yan committed
206
207
    update_corrfors(drls)
    drls['dt'].item(row, 0).setText(gettableitemlist(drls['pt'])[row])
208
    if column in (1, 2):
3Yan's avatar
3Yan committed
209
        update_profile(drls['pt'], row, ds)
210

211

3Yan's avatar
3Yan committed
212
def dtable_changed(row, column, ds, drls, drlspectrum):
213
    """routine called by change of the dtable spectra"""
Yan's avatar
Yan committed
214
    logger.info("Change in the daughter ion table detected")
3Yan's avatar
3Yan committed
215
    #TODO: multiple corrections not incorporated!!!
3Yan's avatar
3Yan committed
216
    if drls['dt'].cellWidget(row, 0).checkState() == 2:
Yan's avatar
Yan committed
217
218
219
        if (column == 0) or (column % 2 == 0 and drls['dt'].cellWidget(
            row, column-1).currentIndex() != 0):
            update_drlspectrum(ds, drls, drlspectrum)
220

Yan's avatar
Yan committed
221

3Yan's avatar
3Yan committed
222
def corr_changed(row, ds, drls, drlspectrum):
223
    """routine called by change of correction for ion"""
Yan's avatar
Yan committed
224
    logger.info('''Change of the "correct to" detected''')
3Yan's avatar
3Yan committed
225
    #TODO: multiple corrections not incorporated!!!
3Yan's avatar
3Yan committed
226
    if (drls['dt'].cellWidget(row, 0).checkState() == 2
Yan's avatar
Yan committed
227
       and floatize(drls['dt'], row, 2, False) != 0):
3Yan's avatar
3Yan committed
228
        update_drlspectrum(ds, drls, drlspectrum)
229
230


3Yan's avatar
3Yan committed
231
def remove_rows(ds, drls, drlspectrum, rows=None):
Yan's avatar
Yan committed
232
    # TODO: maybe nicer selection in future, but this works for now
233
234
    if rows == None:
        rows = reversed(list(set(
3Yan's avatar
3Yan committed
235
            map(lambda x: x.row(), drls['pt'].selectedIndexes()))))
236
    for row in rows:
3Yan's avatar
3Yan committed
237
238
239
240
        drls['dt'].cellWidget(row,0).setCheckState(0)
        drls['dt'].removeRow(row)
        drls['pt'].removeRow(row)
        for i in range(drls['dt'].rowCount()):
3Yan's avatar
3Yan committed
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
            for cornum in range(drls['cors']):
                corfor = drls['dt'].cellWidget(i, 1+cornum*2)
                corfor.disconnect()
                index = corfor.currentIndex()
                corfor.clear()
                corfor.addItem("")
                corfor.addItems(gettableitemlist(drls['pt']))
                if index == row+1:
                    corfor.setCurrentIndex(0)
                    corr_changed(i, ds, drls, drlspectrum)
                elif index > row+1:
                    corfor.setCurrentIndex(index-1)
                else:
                    corfor.setCurrentIndex(index)
                    corfor.currentIndexChanged.connect(lambda:
                        corr_changed(i, ds, drls, drlspectrum))
257
258


3Yan's avatar
3Yan committed
259
def add_line(ds, drls, drlspectrum):
260
    """add parent ion to the table"""
Yan's avatar
Yan committed
261
    logger.debug("adding line")
3Yan's avatar
3Yan committed
262
    newrow = drls['pt'].rowCount()
263

3Yan's avatar
3Yan committed
264
265
    drls['pt'].blockSignals(True)
    drls['dt'].blockSignals(True)
266

3Yan's avatar
3Yan committed
267
    drls['pt'].setRowCount(newrow + 1)
268
    for i in range(3):
3Yan's avatar
3Yan committed
269
        drls['pt'].setItem(newrow, i, QtWidgets.QTableWidgetItem())
270
        if newrow is not 0:
3Yan's avatar
3Yan committed
271
272
            drls['pt'].item(newrow, i).setText(str(floatize(
                    drls['pt'], newrow-1, i)+1))
273

274
    ion_graph = Figure(figsize=(3, 1.5), dpi=100, facecolor="None")
275
276
    ion_graph.add_subplot(111, facecolor=(1, 1, 1, 0.8),
                          position=(-0.01, -0.01, 1.02, 1.02))
277
278
279
    graph_canvas = FigureCanvas(ion_graph)
    graph_canvas.setStyleSheet("background-color:transparent;")
    graph_canvas.setAutoFillBackground(False)
3Yan's avatar
3Yan committed
280
    drls['pt'].setCellWidget(newrow, 3, graph_canvas)
281

3Yan's avatar
3Yan committed
282
    drls['dt'].setRowCount(newrow + 1)
Yan's avatar
Yan committed
283
    checkbox = QtWidgets.QCheckBox()
3Yan's avatar
3Yan committed
284
    checkbox.setFocusProxy(drls['dt'])
285
286
    dname = QtWidgets.QTableWidgetItem()
    dname.setFlags(dname.flags() & ~QtCore.Qt.ItemIsEditable)
Yan's avatar
Yan committed
287
    dname.setTextAlignment(QtCore.Qt.AlignRight)
3Yan's avatar
3Yan committed
288
289
    drls['dt'].setItem(newrow, 0, dname)
    drls['dt'].setCellWidget(newrow, 0, checkbox)
3Yan's avatar
3Yan committed
290
291
292

    for i in range(drls['cors']):
        corfor = QtWidgets.QComboBox()
Yan's avatar
Yan committed
293
294
        #TODO: verify on windows that this fixed the weird selection issues
        #corfor.setFocusProxy(drls['dt'])
3Yan's avatar
3Yan committed
295
296
297
298
299
        corfor.setFrame(False)
        drls['dt'].setCellWidget(newrow, 1+i*2, corfor)
        drls['dt'].setItem(newrow, 2+i*2, QtWidgets.QTableWidgetItem())
        corfor.currentIndexChanged.connect(lambda: corr_changed(
            newrow, ds, drls, drlspectrum))
300

3Yan's avatar
3Yan committed
301
302
    drls['pt'].blockSignals(False)
    drls['dt'].blockSignals(False)
303

3Yan's avatar
3Yan committed
304
    ptable_changed(newrow, 1, ds, drls, drlspectrum)
305

Yan's avatar
Yan committed
306
    checkbox.stateChanged.connect(lambda: update_drlspectrum(
3Yan's avatar
3Yan committed
307
        ds, drls, drlspectrum))
308
309


Yan's avatar
Yan committed
310
def iontable(labels):
311
312
    """creates a table for ions"""
    table = QtWidgets.QTableWidget(columnCount=len(labels))
313
    table.setSizePolicy(QtWidgets.QSizePolicy.Expanding,
314
315
                        QtWidgets.QSizePolicy.Expanding)
    table.setHorizontalHeaderLabels(labels)
316
317
318
319
320
    table.setVerticalScrollBarPolicy(QtCore.Qt.ScrollBarAlwaysOn)

    for n in range(table.columnCount()):
        table.horizontalHeader().setSectionResizeMode(
            n, QtWidgets.QHeaderView.Stretch)
Yan's avatar
Yan committed
321
    table.setMinimumSize(600, 0)
322

323
    return table
324
325


3Yan's avatar
3Yan committed
326
def load_drltables(parent, dataset, drls, drlspectrum):
Yan's avatar
Yan committed
327
328
329
330
331
332
333
334
    filename = QtWidgets.QFileDialog.getOpenFileName(
            caption="Load DRL config tables",
            filter="comma-separated values (*.csv)")[0]
    if filename is not '':
        names = []
        start_masses = []
        end_masses = []
        states = []
3Yan's avatar
3Yan committed
335
        corrections = []
Yan's avatar
Yan committed
336

3Yan's avatar
3Yan committed
337

Yan's avatar
Yan committed
338
339
        with open(filename, 'r') as cfile:
            rawdata = cfile.read().splitlines()
3Yan's avatar
3Yan committed
340
341
        for i in range(len(rawdata[0].split(","))-4):
            corrections.append([])
Yan's avatar
Yan committed
342
343
        for i in range(1, len(rawdata)):
            rawline = rawdata[i].split(",")
344
            # TODO: rawline[4] can be -1 - decide if accept this behaviour
Yan's avatar
Yan committed
345
346
347
348
            n = len(rawline)
            if n < 6 or not n % 2 == 0 or (int(rawline[3]) not in
                range(3)) or not (set(map(int, rawline[4:n:2])) &
                set(range(-1, n))):
Yan's avatar
Yan committed
349
350
351
352
353
354
                QtWidgets.QMessageBox.warning(
                    parent, "Load DRL config tables",
                    "Config file corrupted on line {},"
                    " cancelling request".format(i+1))
                return
            for j, k in enumerate((names, start_masses, end_masses,
3Yan's avatar
3Yan committed
355
                                   states, *corrections)):
Yan's avatar
Yan committed
356
                k.append(rawline[j])
3Yan's avatar
3Yan committed
357
358
359
        for row in reversed(range(drls['pt'].rowCount())):
            drls['dt'].removeRow(row)
            drls['pt'].removeRow(row)
360
        # first populate only the parent table
Yan's avatar
Yan committed
361
        for i in range(len(names)):
3Yan's avatar
3Yan committed
362
363
364
365
            add_line(dataset, drls, drlspectrum)
            drls['pt'].item(i, 0).setText(names[i])
            drls['pt'].item(i, 1).setText(start_masses[i])
            drls['pt'].item(i, 2).setText(end_masses[i])
366
        # and after that the daughter table
Yan's avatar
Yan committed
367
        for i in range(len(names)):
3Yan's avatar
3Yan committed
368
            for j in range(int((len(rawline)-4)/2)):
3Yan's avatar
3Yan committed
369
                drls['dt'].cellWidget(i, 1+j*2).setCurrentIndex(
3Yan's avatar
3Yan committed
370
371
                    int(corrections[0+j*2][i]))
                drls['dt'].item(i, 2+j*2).setText(corrections[1+j*2][i])
3Yan's avatar
3Yan committed
372
            drls['dt'].cellWidget(i, 0).setCheckState(int(states[i]))
Yan's avatar
Yan committed
373
374


3Yan's avatar
3Yan committed
375
def save_drlconfig(drls, parent, exp_f_name=None):
Yan's avatar
Yan committed
376
    """safe DRL table layout so it can be summoned when needed"""
377
378
379
380
    if exp_f_name == None:
        exp_f_name = ft.get_save_filename(
            "Save DRL table layout", "comma-separated values (*.csv)",
            "csv", parent)
Yan's avatar
Yan committed
381
    if exp_f_name is not '':
3Yan's avatar
3Yan committed
382
383
384
385
386
        corlist = []
        for i in range(drls['cors']):
            corlist.append("corrected_to_{}, factor_{} ".format(
                i+1, i+1))
        cortext = ", ".join(corlist)
Yan's avatar
Yan committed
387
388
        expf = open(exp_f_name, 'w')
        expf.write("#ion_name, start m/z, end m/z, visible,"
3Yan's avatar
3Yan committed
389
                   "{}\n".format(cortext))
3Yan's avatar
3Yan committed
390
        for row in range(drls['pt'].rowCount()):
Yan's avatar
Yan committed
391
392
            vals = []
            for i in range(3):
3Yan's avatar
3Yan committed
393
394
                vals.append(drls['pt'].item(row, i).text())
            vals.append(drls['dt'].cellWidget(row, 0).checkState())
3Yan's avatar
3Yan committed
395
396
397
            for i in range(drls['cors']):
                vals.append(drls['dt'].cellWidget(row, 1+i*2).currentIndex())
                vals.append(drls['dt'].item(row, 2+i*2).text())
Yan's avatar
Yan committed
398
399
400
            expf.write("{}\n".format((",".join(map(str, vals)))))
        expf.close()

Yan's avatar
Yan committed
401

3Yan's avatar
3Yan committed
402
def export_drlspectrum(parent, fn, ds, drls):
Yan's avatar
Yan committed
403
404
    if fn[0] is None:
        QtWidgets.QMessageBox.warning(
405
            None, "Export DRL dataset",
Yan's avatar
Yan committed
406
407
            "Nothing to export, cancelling request")
        return
3Yan's avatar
3Yan committed
408
409
410
    names, times, intensities = get_daughterset(ds, drls)
    subset = np.where(times > 0)[0]
    times = times[subset]
Yan's avatar
Yan committed
411
    intensities = list(map(lambda x: x[subset], intensities))
3Yan's avatar
3Yan committed
412
413

    pnames, ptimes, pintensities = get_parentset(ds, drls)
414
415
416
417
418
419
420
421
422
423
424
425
426
427
    if names == []:
        QtWidgets.QMessageBox.warning(
            None, "Export DRL dataset",
            "Nothing to export, cancelling request")
        return
    fname = QtWidgets.QFileDialog.getSaveFileName(
            None,"Export DRL data",
        options=(QtWidgets.QFileDialog.DontConfirmOverwrite |
        QtWidgets.QFileDialog.HideNameFilterDetails))[0]
    if fname is '':
        return
    exp_f_name = list(map(lambda x: "{}/{}-{}.csv".format(
        fname, os.path.basename(fname), x),
        ["raw", "corrected", "input"]))
Yan's avatar
Yan committed
428
429
430
431
432
433
434
435
436
437
    for name in exp_f_name:
        if os.path.isfile(name):
            quest = QtWidgets.QMessageBox.warning(
                parent, "Export DRL data",
                "File {} already exists in the filesystem.\n"
                "Do you want to overwrite it?"
                .format(os.path.basename(name)),
                QtWidgets.QMessageBox.Yes, QtWidgets.QMessageBox.No)
            if quest == QtWidgets.QMessageBox.No:
                return
438
439
440
441
442
443
444
    if not os.path.exists(fname):
        os.makedirs(fname)
    for i, table in enumerate([[pnames, pintensities],
                               [names, intensities]]):
        expf = open(exp_f_name[i], 'w')
        expf.write("times, {}\n".format((",".join(table[0]))))
        for j in range(len(times)):
Yan's avatar
Yan committed
445
            dataset = list()
446
447
448
            dataset.append(times[j])
            for intensity in table[1]:
                dataset.append(intensity[j])
Yan's avatar
Yan committed
449
450
            expf.write("{}\n".format((",".join(map(str, dataset)))))
        expf.close()
3Yan's avatar
3Yan committed
451
    save_drlconfig(drls, parent, exp_f_name[2])
Yan's avatar
Yan committed
452

Yan's avatar
Yan committed
453

3Yan's avatar
3Yan committed
454
def print_graph(labels, ds, drls):
Yan's avatar
Yan committed
455
    printfig = Figure(figsize=(5, 2), dpi=100)
Yan's avatar
Yan committed
456
457
458
    printplot = printfig.add_subplot(111)
    printcanvas = FigureCanvas(printfig)
    gt.pop_plot(printplot, labels)
3Yan's avatar
3Yan committed
459
    update_drlspectrum(ds, drls, printplot)
Yan's avatar
Yan committed
460
    widget = QtWidgets.QDialog(None, windowTitle='Print preview')
Yan's avatar
Yan committed
461
462
    layout = QtWidgets.QVBoxLayout(widget)
    layout.addWidget(printcanvas)
Yan's avatar
Yan committed
463
    widget.resize(600, 400)
Yan's avatar
Yan committed
464
465
466
467
468
469
470
    widget.show()
    dialog = QtPrintSupport.QPrintDialog()
    if dialog.exec_() == QtWidgets.QDialog.Accepted:
        printcanvas.render(dialog.printer())
    widget.close()


3Yan's avatar
3Yan committed
471
def clip_range(drls):
472
    activeObject = QtWidgets.QApplication.focusWidget()
3Yan's avatar
3Yan committed
473
    if activeObject in (drls['pt'], drls['dt']):
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
        table = activeObject
        if len(table.selectedRanges()) == 0:
            return
        if len(table.selectedRanges()) > 1:
            QtWidgets.QMessageBox.warning(
                None, "Operation not supported",
                "Operation not supported for multiple ranges,\n"
                "cancelling request")
            return
        sr = table.selectedRanges()[0]
        rows = []
        for row in range(sr.topRow(), sr.bottomRow() + 1):
            line = []
            for col in range(sr.leftColumn(), sr.rightColumn() + 1):
                if isinstance(table.cellWidget(row, col),
                            QtWidgets.QComboBox):
                    line.append(
                        table.cellWidget(row, col).currentText())
                elif isinstance(table.item(row, col),
                              QtWidgets.QTableWidgetItem):
                    line.append(table.item(row, col).text())
            rows.append(("\t").join(line))
        QtWidgets.QApplication.clipboard().setText(("\n").join(rows))


3Yan's avatar
3Yan committed
499
def paste_clip(ds, drls, drlspectrum):
500
    activeObject = QtWidgets.QApplication.focusWidget()
3Yan's avatar
3Yan committed
501
    if activeObject in (drls['pt'], drls['dt'])\
502
503
504
505
506
507
508
509
510
511
512
513
514
515
        and activeObject.selectedRanges() != []:
        table = activeObject
        if len(table.selectedRanges()) > 1:
            QtWidgets.QMessageBox.warning(
                None, "Operation not supported",
                "Operation not supported for multiple ranges,\n"
                "cancelling request")
        cliptext = QtWidgets.QApplication.clipboard().text()
        rows = cliptext.split("\n")
        startrow = table.selectedRanges()[0].topRow()
        startcol = table.selectedRanges()[0].leftColumn()
        for i, row in enumerate(rows, start=startrow):
            cols = row.split("\t")
            for j, col in enumerate(cols, start=startcol):
3Yan's avatar
3Yan committed
516
                if table == drls['pt'] and j < 3:
517
                    if i >= (table.rowCount()):
3Yan's avatar
3Yan committed
518
                        add_line(ds, drls, drlspectrum)
519
                    table.item(i, j).setText(col)
3Yan's avatar
3Yan committed
520
                if table == drls['dt'] and j == 2\
521
522
523
524
                     and i < table.rowCount():
                     table.item(i,j).setText(col)


3Yan's avatar
3Yan committed
525
def key_pressed(event, ds, drls, drlspectrum):
526
    if event.key() == QtCore.Qt.Key_Delete:
527
        rows = reversed(list(map(lambda x: x.row(),
3Yan's avatar
3Yan committed
528
529
            drls['pt'].selectionModel().selectedRows())))
        remove_rows(ds, drls, drlspectrum, rows)
Yan's avatar
Yan committed
530
    if event.key() == QtCore.Qt.Key_F5:
3Yan's avatar
3Yan committed
531
532
533
        update_drlspectrum(ds, drls, drlspectrum)
        for row in range(drls['pt'].rowCount()):
            update_profile(drls['pt'], row, ds)
534
535
    if event.key() == QtCore.Qt.Key_C\
        and event.modifiers().__int__() == QtCore.Qt.ControlModifier:
3Yan's avatar
3Yan committed
536
        clip_range(drls)
537
538
    if event.key() == QtCore.Qt.Key_V\
        and event.modifiers().__int__() == QtCore.Qt.ControlModifier:
3Yan's avatar
3Yan committed
539
        paste_clip(ds, drls, drlspectrum)
540
541


Yan's avatar
Yan committed
542
def main_window(parent, ds, filename, cache):
543
    """constructs a dialog window"""
3Yan's avatar
3Yan committed
544
545
    def saveonclose(widget, event, buffer, drls, canvas):
        buffer[0], buffer[1] = drls, canvas
Yan's avatar
Yan committed
546
        QtWidgets.QMainWindow.closeEvent(widget, event)
Yan's avatar
Yan committed
547

3Yan's avatar
3Yan committed
548
549
    # pt = parenttable
    # dt = daughtertable
3Yan's avatar
3Yan committed
550
    drls = dict(pt=None, dt=None, tshift=None, cors=3, rel=None)
3Yan's avatar
3Yan committed
551

Yan's avatar
Yan committed
552
    window = QtWidgets.QMainWindow(
Yan's avatar
Yan committed
553
        parent, windowTitle='Delayed reactant labelling')
Yan's avatar
Yan committed
554
555
    main_widget = QtWidgets.QWidget(window)
    window.setCentralWidget(main_widget)
Yan's avatar
Yan committed
556

Yan's avatar
Yan committed
557
    window.closeEvent = lambda event: saveonclose(
3Yan's avatar
3Yan committed
558
        window, event, cache, drls, graph_canvas)
Yan's avatar
Yan committed
559

Yan's avatar
Yan committed
560
561
    dial_graph = Figure(figsize=(5, 2), dpi=100, facecolor="None")
    chromplot = dial_graph.add_subplot(111, facecolor=(1, 1, 1, 0.8))
Yan's avatar
Yan committed
562
563
564
    graph_canvas = FigureCanvas(dial_graph)
    graph_canvas.setStyleSheet("background-color:transparent;")
    graph_canvas.setAutoFillBackground(False)
Yan's avatar
Yan committed
565
    graphlabels = dict(x=[0], y=[0], line=None, name="",
Yan's avatar
Yan committed
566
567
                       xlabel="time(min)",
                       ylabel="relative intensity")
Yan's avatar
Yan committed
568
569
570
571
    gt.pan_factory(chromplot)
    gt.zoom_factory(chromplot, 1.15)
    gt.pop_plot(chromplot, graphlabels)

3Yan's avatar
3Yan committed
572
573
574
575
    time_title = QtWidgets.QLabel("Time shift (min):")
    drls['tshift'] = QtWidgets.QDoubleSpinBox(
            minimum=-100, maximum=1440, decimals=3)

3Yan's avatar
3Yan committed
576
577
    drls['rel'] = QtWidgets.QCheckBox("Steady state approximation")

578
579
580
    drl_load = QtWidgets.QPushButton("&Load")
    drl_save = QtWidgets.QPushButton("&Save")
    drl_export = QtWidgets.QPushButton("&Export")
Yan's avatar
Yan committed
581
    drl_print = QtWidgets.QPushButton("&Print")
582
    close = QtWidgets.QPushButton("&Close")
Yan's avatar
Yan committed
583
    close.clicked.connect(window.close)
584

Yan's avatar
Yan committed
585
586
587
    btn_add = QtWidgets.QPushButton("&Add")
    btn_rem = QtWidgets.QPushButton("Remove")

3Yan's avatar
3Yan committed
588
    if cache == [None, None]:
3Yan's avatar
3Yan committed
589
590
591
592
593
        dcolums = ["Name"]
        for i in range(drls['cors']):
            dcolums.append("corrected for ({})".format(i+1))
            dcolums.append("factor ({})".format(i+1))
        drls['dt'] = iontable(dcolums)
3Yan's avatar
3Yan committed
594
595
        drls['pt'] = iontable(["Name", "start (m/z)", "end (m/z)",
                               "profile"])
3Yan's avatar
3Yan committed
596
597
598
599
        #TODO: DIRTY, DIRTY, DIRTY !!! do it nicer when I'll know how i want it
        for n in range(drls['dt'].columnCount()):
            drls['dt'].horizontalHeader().setSectionResizeMode(
                n, QtWidgets.QHeaderView.Interactive)
3Yan's avatar
3Yan committed
600
        add_line(ds, drls, chromplot)
Yan's avatar
Yan committed
601
    else:
3Yan's avatar
3Yan committed
602
603
        drls = cache[0]
        graph_canvas = cache[1]
604

605
    window.keyPressEvent = lambda event: key_pressed(
3Yan's avatar
3Yan committed
606
        event, ds, drls, chromplot)
607

608
    btn_add.clicked.connect(lambda: add_line(
3Yan's avatar
3Yan committed
609
        ds, drls, chromplot))
610
    btn_rem.clicked.connect(lambda: remove_rows(
3Yan's avatar
3Yan committed
611
        ds, drls, chromplot))
Yan's avatar
Yan committed
612
    drl_load.clicked.connect(lambda: load_drltables(
3Yan's avatar
3Yan committed
613
        main_widget, ds, drls, chromplot))
Yan's avatar
Yan committed
614
    drl_save.clicked.connect(lambda: save_drlconfig(
3Yan's avatar
3Yan committed
615
        drls, main_widget))
Yan's avatar
Yan committed
616
    drl_print.clicked.connect(lambda: print_graph(
3Yan's avatar
3Yan committed
617
        graphlabels, ds, drls))
Yan's avatar
Yan committed
618
    drl_export.clicked.connect(lambda: export_drlspectrum(
3Yan's avatar
3Yan committed
619
620
621
622
        main_widget, filename, ds, drls))

    drls['pt'].itemChanged.connect(lambda item: ptable_changed(
        item.row(), item.column(), ds, drls, chromplot))
623

3Yan's avatar
3Yan committed
624
625
    drls['dt'].itemChanged.connect(lambda item: dtable_changed(
        item.row(), item.column(), ds, drls, chromplot))
626

3Yan's avatar
3Yan committed
627
628
    drls['tshift'].valueChanged.connect(lambda: update_drlspectrum(
        ds, drls, chromplot))
629

3Yan's avatar
3Yan committed
630
631
632
    drls['rel'].stateChanged.connect(lambda: update_drlspectrum(
        ds, drls, chromplot))

Yan's avatar
Yan committed
633
    main_layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
634
    sub_layout = QtWidgets.QHBoxLayout()
3Yan's avatar
3Yan committed
635
636
    graph_layout = QtWidgets.QVBoxLayout()
    graphparams_layout = QtWidgets.QHBoxLayout()
Yan's avatar
Yan committed
637
    tablelayout = QtWidgets.QVBoxLayout()
Yan's avatar
Yan committed
638
    pt_butlayout = QtWidgets.QHBoxLayout()
Yan's avatar
Yan committed
639
    main_butlayout = QtWidgets.QHBoxLayout()
Yan's avatar
Yan committed
640

Yan's avatar
Yan committed
641
642
643
644
    pt_butlayout.addWidget(btn_add)
    pt_butlayout.addWidget(btn_rem)
    pt_butlayout.addStretch(0)

645
646
    main_layout.addLayout(sub_layout)
    main_layout.addWidget(HBar())
Yan's avatar
Yan committed
647
    main_layout.addLayout(main_butlayout)
648

Yan's avatar
Yan committed
649
650
    main_butlayout.addWidget(drl_load)
    main_butlayout.addWidget(drl_save)
Yan's avatar
Yan committed
651
    main_butlayout.addWidget(drl_print)
Yan's avatar
Yan committed
652
653
    main_butlayout.addStretch(1)
    main_butlayout.addWidget(drl_export)
654
    main_butlayout.addWidget(close)
Yan's avatar
Yan committed
655

3Yan's avatar
3Yan committed
656
    sub_layout.addLayout(graph_layout, stretch=1)
Yan's avatar
Yan committed
657
    sub_layout.addLayout(tablelayout)
3Yan's avatar
3Yan committed
658
659
660
661
662
    graph_layout.addWidget(graph_canvas, stretch=1)
    graph_layout.addLayout(graphparams_layout)
    graphparams_layout.addWidget(time_title)
    graphparams_layout.addWidget(drls['tshift'])
    graphparams_layout.addStretch(1)
3Yan's avatar
3Yan committed
663
664
    graphparams_layout.addWidget(drls['rel'])
    graphparams_layout.addStretch(1)
665

Yan's avatar
Yan committed
666
    tablelayout.addWidget(QtWidgets.QLabel("Raw ions table:"))
3Yan's avatar
3Yan committed
667
    tablelayout.addWidget(drls['pt'])
Yan's avatar
Yan committed
668
669
    tablelayout.addLayout(pt_butlayout)
    tablelayout.addWidget(QtWidgets.QLabel("Corrected ions table:"))
3Yan's avatar
3Yan committed
670
    tablelayout.addWidget(drls['dt'])
671

Yan's avatar
Yan committed
672
    window.show()