drltools.py 29.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
3Yan's avatar
3Yan committed
8
from copy import copy
Yan's avatar
Yan committed
9
import gc
Yan's avatar
Yan committed
10
11
import matplotlib
import numpy as np
Yan's avatar
Yan committed
12
13
import prasopes.datatools as dt
import prasopes.graphtools as gt
Yan's avatar
Yan committed
14
import prasopes.filetools as ft
Yan's avatar
Yan committed
15
import os.path
Yan's avatar
Yan committed
16
import logging
Yan's avatar
Yan committed
17
18
matplotlib.use("Qt5Agg")

Yan's avatar
Yan committed
19

Yan's avatar
Yan committed
20
21
22
logger = logging.getLogger('drlLogger')


Yan's avatar
Yan committed
23
24
class StretchWidget(QtWidgets.QWidget):
    """horizontal stretch class"""
Yan's avatar
Yan committed
25
    def __init__(self):
Yan's avatar
Yan committed
26
        super(StretchWidget, self).__init__()
Yan's avatar
Yan committed
27
        self._main = QtWidgets.QWidget()
Yan's avatar
Yan committed
28
29
        self.setSizePolicy(QtWidgets.QSizePolicy.Expanding
                          ,QtWidgets.QSizePolicy.Preferred)
Yan's avatar
Yan committed
30
31


Yan's avatar
Yan committed
32
def floatize(table, row, column, nonneg=True):
33
34
35
36
    """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
37
38
    if nonneg == True:
        validator.setBottom(0)
39
40
41
42
43
44
45
46
47
    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
48
    """parent table profile spectrum updating procedure"""
Yan's avatar
Yan committed
49
    # Dont do anything to graph when the spectrum is not populated
Yan's avatar
Yan committed
50
    logger.debug("updating parent table row {} profile".format(row))
Yan's avatar
Yan committed
51
    if isinstance(dataset['masses'], type(None)):
52
53
        return

54
    spectrum = pt.cellWidget(row, 3).figure.get_axes()[0]
Yan's avatar
Yan committed
55
    masses = dataset['masses']
56
57
    massargs = dt.argsubselect(
        masses, floatize(pt, row, 1), floatize(pt, row, 2))
58
    yshape = np.mean(dataset['matrix'], axis=0)
Yan's avatar
Yan committed
59
    spectrum.clear()
60
61
62
    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
63
64
65
    spectrum.set_xlim(masses[massargs[0]]-xex,
                      masses[massargs[-1]]+xex)
    ymax = max(yshape[massargs])
66
67
68
    spectrum.set_ylim(ymax*-0.1, ymax*1.2)
    spectrum.figure.canvas.draw()

Yan's avatar
Yan committed
69

3Yan's avatar
3Yan committed
70
71
72
def get_intensity(row, ds, drls):
    startm = floatize(drls['pt'], row, 1)
    endm = floatize(drls['pt'], row, 2)
Yan's avatar
Yan committed
73
    massargs = dt.argsubselect(ds['masses'], startm, endm)
Yan's avatar
Yan committed
74
75
76
77
    #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
78
79
80
    return intensity


3Yan's avatar
3Yan committed
81
def get_daughterset(ds, drls):
Yan's avatar
Yan committed
82
83
    """Fuction to acquire the curves of the daugher ions"""
    logger.info("getting set of the daughter ions")
Yan's avatar
Yan committed
84
    names = []
3Yan's avatar
3Yan committed
85
    times = ds['chrom_dat'][0, :] - drls['tshift'].value()
Yan's avatar
Yan committed
86
    intensities = []
3Yan's avatar
3Yan committed
87
88
89
    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
90
            corlist = []
Yan's avatar
Yan committed
91
            for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
92
93
                cor = drls['dt'].cellWidget(row, 1+i*2).currentIndex() - 1
                if cor not in (-2, -1):
Yan's avatar
Yan committed
94
                    factor = floatize(drls['dt'], row, 2+i*2, False)
3Yan's avatar
3Yan committed
95
96
97
98
99
100
                    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
101
            intensities.append(intensity)
3Yan's avatar
3Yan committed
102
            names.append("{} - ({})".format(
3Yan's avatar
3Yan committed
103
                drls['dt'].item(row, 0).text(),
3Yan's avatar
3Yan committed
104
                cortext))
Yan's avatar
Yan committed
105
106
107
    return names, times, intensities


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


3Yan's avatar
3Yan committed
127
def update_drlspectrum(ds, drls, drlspectrum):
Yan's avatar
Yan committed
128
129
    """Generic DRL spectrum updating procedure"""
    logger.info("updating DRL spectrum")
Yan's avatar
Yan committed
130
    # Do not do anything when the data set is not populated.
Yan's avatar
Yan committed
131
    if isinstance(ds['masses'], type(None)):
Yan's avatar
Yan committed
132
        return
133
    colors = np.array([[0, 0, 0], [255, 0, 0], [0, 255, 0], [0, 0, 255],
134
135
136
                       [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
137

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
                #Do not divide by 0
                intensity = np.divide(intensity, np.clip(np.sum(
                    intensities, 0), np.nextafter(
3Yan's avatar
3Yan committed
157
                        0, 1, dtype=np.float32), None), dtype=np.float64)
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()):
Yan's avatar
Yan committed
192
        for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
193
194
195
196
197
198
199
200
            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
    if drls['dt'].cellWidget(row, 0).checkState() == 2:
Yan's avatar
Yan committed
216
217
218
        if (column == 0) or (column % 2 == 0 and drls['dt'].cellWidget(
            row, column-1).currentIndex() != 0):
            update_drlspectrum(ds, drls, drlspectrum)
219

Yan's avatar
Yan committed
220

3Yan's avatar
3Yan committed
221
def corr_changed(correction, ds, drls, drlspectrum):
222
    """routine called by change of correction for ion"""
Yan's avatar
Yan committed
223
    logger.info('''Change of the "correct to" detected''')
3Yan's avatar
3Yan committed
224
225
226
227
    for i in range(drls['dt'].rowCount()):
        for j in range(drls['dt'].columnCount()):
            if correction == drls['dt'].cellWidget(i, 1+j*2):
                row, column = i, 1 + j*2
Yan's avatar
Yan committed
228
229
230
                logger.debug('''Change of the "correct to" on '''
                             '''row {}, column {}'''.format(
                                 row + 1, column + 1))
3Yan's avatar
3Yan committed
231
    if (drls['dt'].cellWidget(row, 0).checkState() == 2
3Yan's avatar
3Yan committed
232
       and floatize(drls['dt'], row, column+1, False) != 0):
3Yan's avatar
3Yan committed
233
        update_drlspectrum(ds, drls, drlspectrum)
234
235


Yan's avatar
Yan committed
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
def corcount_changed(ds, drls, drlspectrum):
    """routine called by change of the correction factors count"""
    logger.info('''change in the correction count detected''')
    diff = int(drls['cors'].value() - ((drls['dt'].columnCount() - 1) / 2))
    if diff == 0:
        return
    elif diff > 0:
        drls['dt'].setColumnCount(1 + (drls['cors'].value() * 2))
        newcors = []
        for col in range(drls['dt'].columnCount() - (diff * 2),
                       drls['dt'].columnCount(), 2):
            for row in range(drls['dt'].rowCount()):
                    drls['dt'].setCellWidget(
                        row, col, QtWidgets.QComboBox())
                    drls['dt'].cellWidget(row, col).setFrame(False)
                    drls['dt'].setItem(
                        row, col+1, QtWidgets.QTableWidgetItem())
                    newcors.append(drls['dt'].cellWidget(row, col))
        list(map(lambda x: x.currentIndexChanged.connect(lambda:
            corr_changed(x, ds, drls, drlspectrum)), newcors))
        update_corrfors(drls)
        dcolums = ["Name"]
        for i in range(drls['cors'].value()):
            dcolums.append("corrected for ({})".format(i+1))
            dcolums.append("factor ({})".format(i+1))
        drls['dt'].setHorizontalHeaderLabels(dcolums)
    else:
        drls['dt'].setColumnCount(1 + (drls['cors'].value() * 2))
        update_drlspectrum(ds, drls, drlspectrum)


3Yan's avatar
3Yan committed
267
def remove_rows(ds, drls, drlspectrum, rows=None):
3Yan's avatar
3Yan committed
268
    logger.info("remowing rows")
269
270
    if rows == None:
        rows = reversed(list(set(
3Yan's avatar
3Yan committed
271
            map(lambda x: x.row(), drls['pt'].selectedIndexes()))))
272
    cors = []
273
    for row in rows:
3Yan's avatar
3Yan committed
274
275
276
277
        drls['dt'].cellWidget(row,0).setCheckState(0)
        drls['dt'].removeRow(row)
        drls['pt'].removeRow(row)
        for i in range(drls['dt'].rowCount()):
Yan's avatar
Yan committed
278
            for cornum in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
279
                corfor = drls['dt'].cellWidget(i, 1+cornum*2)
3Yan's avatar
3Yan committed
280
                cors.append(corfor)
3Yan's avatar
3Yan committed
281
282
283
284
285
286
287
                corfor.disconnect()
                index = corfor.currentIndex()
                corfor.clear()
                corfor.addItem("")
                corfor.addItems(gettableitemlist(drls['pt']))
                if index == row+1:
                    corfor.setCurrentIndex(0)
3Yan's avatar
3Yan committed
288
                    corr_changed(corfor, ds, drls, drlspectrum)
3Yan's avatar
3Yan committed
289
290
291
292
                elif index > row+1:
                    corfor.setCurrentIndex(index-1)
                else:
                    corfor.setCurrentIndex(index)
Yan's avatar
Yan committed
293
294
        list(map(lambda x: x.currentIndexChanged.connect(lambda:
            corr_changed(x, ds, drls, drlspectrum)), cors))
295
296


3Yan's avatar
3Yan committed
297
def add_line(ds, drls, drlspectrum):
298
    """add parent ion to the table"""
Yan's avatar
Yan committed
299
    logger.debug("adding line")
3Yan's avatar
3Yan committed
300
    newrow = drls['pt'].rowCount()
301

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

3Yan's avatar
3Yan committed
305
    drls['pt'].setRowCount(newrow + 1)
306
    for i in range(3):
3Yan's avatar
3Yan committed
307
        drls['pt'].setItem(newrow, i, QtWidgets.QTableWidgetItem())
308
        if newrow is not 0:
3Yan's avatar
3Yan committed
309
310
            drls['pt'].item(newrow, i).setText(str(floatize(
                    drls['pt'], newrow-1, i)+1))
311

312
    ion_graph = Figure(figsize=(3, 1.5), dpi=100, facecolor="None")
313
314
    ion_graph.add_subplot(111, facecolor=(1, 1, 1, 0.8),
                          position=(-0.01, -0.01, 1.02, 1.02))
315
316
317
    graph_canvas = FigureCanvas(ion_graph)
    graph_canvas.setStyleSheet("background-color:transparent;")
    graph_canvas.setAutoFillBackground(False)
3Yan's avatar
3Yan committed
318
    drls['pt'].setCellWidget(newrow, 3, graph_canvas)
319

3Yan's avatar
3Yan committed
320
    drls['dt'].setRowCount(newrow + 1)
Yan's avatar
Yan committed
321
    checkbox = QtWidgets.QCheckBox()
3Yan's avatar
3Yan committed
322
    checkbox.setFocusProxy(drls['dt'])
323
324
    dname = QtWidgets.QTableWidgetItem()
    dname.setFlags(dname.flags() & ~QtCore.Qt.ItemIsEditable)
Yan's avatar
Yan committed
325
    dname.setTextAlignment(QtCore.Qt.AlignRight)
3Yan's avatar
3Yan committed
326
327
    drls['dt'].setItem(newrow, 0, dname)
    drls['dt'].setCellWidget(newrow, 0, checkbox)
3Yan's avatar
3Yan committed
328

Yan's avatar
Yan committed
329
    for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
330
331
332
333
        col = copy(1+i*2)
        drls['dt'].setCellWidget(newrow, col, QtWidgets.QComboBox())
        drls['dt'].cellWidget(newrow, col).setFrame(False)
        drls['dt'].setItem(newrow, col+1, QtWidgets.QTableWidgetItem())
334

3Yan's avatar
3Yan committed
335
336
    drls['pt'].blockSignals(False)
    drls['dt'].blockSignals(False)
337

3Yan's avatar
3Yan committed
338
    cors = list(map(lambda x: drls['dt'].cellWidget(newrow, 1+x*2),
Yan's avatar
Yan committed
339
                    range(drls['cors'].value())))
3Yan's avatar
3Yan committed
340
341
342
    list(map(lambda x: x.currentIndexChanged.connect(lambda:
        corr_changed(x, ds, drls, drlspectrum)), cors))

3Yan's avatar
3Yan committed
343
    ptable_changed(newrow, 1, ds, drls, drlspectrum)
344
    select_all_btn_up(ds, drls, drlspectrum)
345

Yan's avatar
Yan committed
346
    checkbox.stateChanged.connect(lambda: update_drlspectrum(
3Yan's avatar
3Yan committed
347
        ds, drls, drlspectrum))
348
349
    checkbox.stateChanged.connect(lambda: select_all_btn_up(
        ds, drls, drlspectrum))
350
351


Yan's avatar
Yan committed
352
def iontable(labels):
353
354
    """creates a table for ions"""
    table = QtWidgets.QTableWidget(columnCount=len(labels))
355
    table.setSizePolicy(QtWidgets.QSizePolicy.Expanding,
356
357
                        QtWidgets.QSizePolicy.Expanding)
    table.setHorizontalHeaderLabels(labels)
358
359
360
361
362
    table.setVerticalScrollBarPolicy(QtCore.Qt.ScrollBarAlwaysOn)

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

365
    return table
366
367


3Yan's avatar
3Yan committed
368
def load_drltables(parent, dataset, drls, drlspectrum):
Yan's avatar
Yan committed
369
370
371
372
373
374
375
376
    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
377
        corrections = []
Yan's avatar
Yan committed
378

3Yan's avatar
3Yan committed
379

Yan's avatar
Yan committed
380
381
        with open(filename, 'r') as cfile:
            rawdata = cfile.read().splitlines()
3Yan's avatar
3Yan committed
382
383
        for i in range(len(rawdata[0].split(","))-4):
            corrections.append([])
Yan's avatar
Yan committed
384
385
        for i in range(1, len(rawdata)):
            rawline = rawdata[i].split(",")
Yan's avatar
Yan committed
386
387
388
389
            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
390
391
392
393
394
395
                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
396
                                   states, *corrections)):
Yan's avatar
Yan committed
397
                k.append(rawline[j])
3Yan's avatar
3Yan committed
398
399
400
        for row in reversed(range(drls['pt'].rowCount())):
            drls['dt'].removeRow(row)
            drls['pt'].removeRow(row)
401
        # first populate only the parent table
Yan's avatar
Yan committed
402
        for i in range(len(names)):
3Yan's avatar
3Yan committed
403
404
405
406
            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])
407
        # and after that the daughter table
Yan's avatar
Yan committed
408
        drls['cors'].setValue(int(1 + (len(names) / 2)))
Yan's avatar
Yan committed
409
        for i in range(len(names)):
3Yan's avatar
3Yan committed
410
            for j in range(int((len(rawline)-4)/2)):
3Yan's avatar
3Yan committed
411
                drls['dt'].cellWidget(i, 1+j*2).setCurrentIndex(
3Yan's avatar
3Yan committed
412
413
                    int(corrections[0+j*2][i]))
                drls['dt'].item(i, 2+j*2).setText(corrections[1+j*2][i])
3Yan's avatar
3Yan committed
414
            drls['dt'].cellWidget(i, 0).setCheckState(int(states[i]))
Yan's avatar
Yan committed
415
416


3Yan's avatar
3Yan committed
417
def save_drlconfig(drls, parent, exp_f_name=None):
Yan's avatar
Yan committed
418
    """safe DRL table layout so it can be summoned when needed"""
419
420
421
422
    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
423
    if exp_f_name is not '':
3Yan's avatar
3Yan committed
424
        corlist = []
Yan's avatar
Yan committed
425
        for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
426
427
428
            corlist.append("corrected_to_{}, factor_{} ".format(
                i+1, i+1))
        cortext = ", ".join(corlist)
Yan's avatar
Yan committed
429
430
        expf = open(exp_f_name, 'w')
        expf.write("#ion_name, start m/z, end m/z, visible,"
3Yan's avatar
3Yan committed
431
                   "{}\n".format(cortext))
3Yan's avatar
3Yan committed
432
        for row in range(drls['pt'].rowCount()):
Yan's avatar
Yan committed
433
434
            vals = []
            for i in range(3):
3Yan's avatar
3Yan committed
435
436
                vals.append(drls['pt'].item(row, i).text())
            vals.append(drls['dt'].cellWidget(row, 0).checkState())
Yan's avatar
Yan committed
437
            for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
438
439
                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
440
441
442
            expf.write("{}\n".format((",".join(map(str, vals)))))
        expf.close()

Yan's avatar
Yan committed
443

3Yan's avatar
3Yan committed
444
def export_drlspectrum(parent, fn, ds, drls):
Yan's avatar
Yan committed
445
446
    if fn[0] is None:
        QtWidgets.QMessageBox.warning(
Yan's avatar
Yan committed
447
448
            None, "Export DRL dat aset",
            "No file opened. Nothing to export, canceling request")
Yan's avatar
Yan committed
449
        return
3Yan's avatar
3Yan committed
450
451
452
    names, times, intensities = get_daughterset(ds, drls)
    subset = np.where(times > 0)[0]
    times = times[subset]
Yan's avatar
Yan committed
453
    intensities = list(map(lambda x: x[subset], intensities))
3Yan's avatar
3Yan committed
454
455

    pnames, ptimes, pintensities = get_parentset(ds, drls)
456
457
    if names == []:
        QtWidgets.QMessageBox.warning(
Yan's avatar
Yan committed
458
459
460
            None, "Export DRL data set",
            "No rows in the Corrected ions table selected. "
            "Nothing to export, canceling request")
461
462
463
464
        return
    fname = QtWidgets.QFileDialog.getSaveFileName(
            None,"Export DRL data",
        options=(QtWidgets.QFileDialog.DontConfirmOverwrite |
Yan's avatar
Yan committed
465
466
        QtWidgets.QFileDialog.HideNameFilterDetails),
        directory=os.path.basename(fn[0][:-4]))[0]
467
468
469
470
471
    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
472
473
474
475
476
477
478
479
480
481
    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
482
483
484
485
486
487
    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]))))
488
        expf.write("timeshift = {}".format(drls['tshift'].value()))
489
        for j in range(len(times)):
Yan's avatar
Yan committed
490
            dataset = list()
491
492
493
            dataset.append(times[j])
            for intensity in table[1]:
                dataset.append(intensity[j])
Yan's avatar
Yan committed
494
495
            expf.write("{}\n".format((",".join(map(str, dataset)))))
        expf.close()
3Yan's avatar
3Yan committed
496
    save_drlconfig(drls, parent, exp_f_name[2])
Yan's avatar
Yan committed
497

Yan's avatar
Yan committed
498

3Yan's avatar
3Yan committed
499
def print_graph(labels, ds, drls):
Yan's avatar
Yan committed
500
    printfig = Figure(figsize=(5, 2), dpi=100)
Yan's avatar
Yan committed
501
502
503
    printplot = printfig.add_subplot(111)
    printcanvas = FigureCanvas(printfig)
    gt.pop_plot(printplot, labels)
3Yan's avatar
3Yan committed
504
    update_drlspectrum(ds, drls, printplot)
Yan's avatar
Yan committed
505
    widget = QtWidgets.QDialog(None, windowTitle='Print preview')
Yan's avatar
Yan committed
506
507
    layout = QtWidgets.QVBoxLayout(widget)
    layout.addWidget(printcanvas)
Yan's avatar
Yan committed
508
    widget.resize(600, 400)
Yan's avatar
Yan committed
509
510
511
512
513
514
515
    widget.show()
    dialog = QtPrintSupport.QPrintDialog()
    if dialog.exec_() == QtWidgets.QDialog.Accepted:
        printcanvas.render(dialog.printer())
    widget.close()


3Yan's avatar
3Yan committed
516
def clip_range(drls):
517
    activeObject = QtWidgets.QApplication.focusWidget()
3Yan's avatar
3Yan committed
518
    if activeObject in (drls['pt'], drls['dt']):
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
        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
544
def paste_clip(ds, drls, drlspectrum):
545
    activeObject = QtWidgets.QApplication.focusWidget()
3Yan's avatar
3Yan committed
546
    if activeObject in (drls['pt'], drls['dt'])\
547
548
549
550
551
552
553
554
555
556
557
558
559
560
        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
561
                if table == drls['pt'] and j < 3:
562
                    if i >= (table.rowCount()):
3Yan's avatar
3Yan committed
563
                        add_line(ds, drls, drlspectrum)
564
                    table.item(i, j).setText(col)
565
                if table == drls['dt'] and j != 0 and j % 2 == 0\
566
567
568
569
                     and i < table.rowCount():
                     table.item(i,j).setText(col)


3Yan's avatar
3Yan committed
570
def key_pressed(event, ds, drls, drlspectrum):
571
    if event.key() == QtCore.Qt.Key_Delete:
572
        rows = reversed(list(map(lambda x: x.row(),
3Yan's avatar
3Yan committed
573
574
            drls['pt'].selectionModel().selectedRows())))
        remove_rows(ds, drls, drlspectrum, rows)
Yan's avatar
Yan committed
575
    if event.key() == QtCore.Qt.Key_F5:
3Yan's avatar
3Yan committed
576
577
578
        update_drlspectrum(ds, drls, drlspectrum)
        for row in range(drls['pt'].rowCount()):
            update_profile(drls['pt'], row, ds)
579
580
    if event.key() == QtCore.Qt.Key_C\
        and event.modifiers().__int__() == QtCore.Qt.ControlModifier:
3Yan's avatar
3Yan committed
581
        clip_range(drls)
582
583
    if event.key() == QtCore.Qt.Key_V\
        and event.modifiers().__int__() == QtCore.Qt.ControlModifier:
3Yan's avatar
3Yan committed
584
        paste_clip(ds, drls, drlspectrum)
585

Yan's avatar
Yan committed
586

587
588
def select_all_btn_up(ds, drls, drlspectrum, state=None):
    logger.info('''Change of the select-all button detected''')
589
590
591
592
593
594
595
596
597
598
    if state in (0,2):
        for row in range(drls['dt'].rowCount()):
            drls['dt'].cellWidget(row, 0).blockSignals(True)
            drls['dt'].cellWidget(row, 0).setCheckState(state)
            drls['dt'].cellWidget(row, 0).blockSignals(False)
        update_drlspectrum(ds, drls, drlspectrum)
    elif state == 1:
        drls['checkAll'].setCheckState(2)
    else:
        drls['checkAll'].blockSignals(True)
599
600
601
        btns = []
        for row in range(drls['dt'].rowCount()):
            btns.append(drls['dt'].cellWidget(row, 0).checkState())
602
        if all(btns):
603
            drls['checkAll'].setCheckState(2)
604
605
        elif any(btns):
            drls['checkAll'].setCheckState(1)
606
607
        else:
            drls['checkAll'].setCheckState(0)
608
        drls['checkAll'].blockSignals(False)
609

610

Yan's avatar
Yan committed
611
def main_window(parent, ds, filename, cache):
612
    """constructs a dialog window"""
3Yan's avatar
3Yan committed
613
614
    def saveonclose(widget, event, buffer, drls, canvas):
        buffer[0], buffer[1] = drls, canvas
Yan's avatar
Yan committed
615
        QtWidgets.QMainWindow.closeEvent(widget, event)
Yan's avatar
Yan committed
616

3Yan's avatar
3Yan committed
617
618
    # pt = parenttable
    # dt = daughtertable
Yan's avatar
Yan committed
619
    drls = dict(pt=None, dt=None, tshift=None, cors=None, rel=None,
620
                checkAll=None)
3Yan's avatar
3Yan committed
621

Yan's avatar
Yan committed
622
    window = QtWidgets.QMainWindow(
Yan's avatar
Yan committed
623
        parent, windowTitle='Delayed reactant labeling')
Yan's avatar
Yan committed
624
625
    main_widget = QtWidgets.QWidget(window)
    window.setCentralWidget(main_widget)
Yan's avatar
Yan committed
626

Yan's avatar
Yan committed
627
    window.closeEvent = lambda event: saveonclose(
3Yan's avatar
3Yan committed
628
        window, event, cache, drls, graph_canvas)
Yan's avatar
Yan committed
629

Yan's avatar
Yan committed
630
631
    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
632
633
634
    graph_canvas = FigureCanvas(dial_graph)
    graph_canvas.setStyleSheet("background-color:transparent;")
    graph_canvas.setAutoFillBackground(False)
Yan's avatar
Yan committed
635
    graphlabels = dict(x=[0], y=[0], line=None, name="",
Yan's avatar
Yan committed
636
637
                       xlabel="time(min)",
                       ylabel="relative intensity")
Yan's avatar
Yan committed
638
639
640
641
    gt.pan_factory(chromplot)
    gt.zoom_factory(chromplot, 1.15)
    gt.pop_plot(chromplot, graphlabels)

3Yan's avatar
3Yan committed
642
643
644
    time_title = QtWidgets.QLabel("Time shift (min):")
    drls['tshift'] = QtWidgets.QDoubleSpinBox(
            minimum=-100, maximum=1440, decimals=3)
3Yan's avatar
3Yan committed
645
    drls['rel'] = QtWidgets.QCheckBox("Steady state approximation")
646
    drls['checkAll'] = QtWidgets.QCheckBox("Select all")
Yan's avatar
Yan committed
647
    drls['cors'] = QtWidgets.QSpinBox(minimum=0)
648

649
650
651
    drl_load = QtWidgets.QPushButton("&Load")
    drl_save = QtWidgets.QPushButton("&Save")
    drl_export = QtWidgets.QPushButton("&Export")
Yan's avatar
Yan committed
652
    drl_print = QtWidgets.QPushButton("&Print")
653
    close = QtWidgets.QPushButton("&Close")
Yan's avatar
Yan committed
654
    close.clicked.connect(window.close)
655

Yan's avatar
Yan committed
656
657
658
    btn_add = QtWidgets.QPushButton("&Add")
    btn_rem = QtWidgets.QPushButton("Remove")

3Yan's avatar
3Yan committed
659
    if cache == [None, None]:
3Yan's avatar
3Yan committed
660
        dcolums = ["Name"]
Yan's avatar
Yan committed
661
662
        drls['cors'].setValue(3)
        for i in range(drls['cors'].value()):
3Yan's avatar
3Yan committed
663
664
665
            dcolums.append("corrected for ({})".format(i+1))
            dcolums.append("factor ({})".format(i+1))
        drls['dt'] = iontable(dcolums)
3Yan's avatar
3Yan committed
666
667
        drls['pt'] = iontable(["Name", "start (m/z)", "end (m/z)",
                               "profile"])
3Yan's avatar
3Yan committed
668
669
670
671
        #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
672
        add_line(ds, drls, chromplot)
Yan's avatar
Yan committed
673
    else:
3Yan's avatar
3Yan committed
674
675
        drls = cache[0]
        graph_canvas = cache[1]
676

677
    window.keyPressEvent = lambda event: key_pressed(
3Yan's avatar
3Yan committed
678
        event, ds, drls, chromplot)
679

680
    btn_add.clicked.connect(lambda: add_line(
3Yan's avatar
3Yan committed
681
        ds, drls, chromplot))
682
    btn_rem.clicked.connect(lambda: remove_rows(
3Yan's avatar
3Yan committed
683
        ds, drls, chromplot))
Yan's avatar
Yan committed
684
    drl_load.clicked.connect(lambda: load_drltables(
3Yan's avatar
3Yan committed
685
        main_widget, ds, drls, chromplot))
Yan's avatar
Yan committed
686
    drl_save.clicked.connect(lambda: save_drlconfig(
3Yan's avatar
3Yan committed
687
        drls, main_widget))
Yan's avatar
Yan committed
688
    drl_print.clicked.connect(lambda: print_graph(
3Yan's avatar
3Yan committed
689
        graphlabels, ds, drls))
Yan's avatar
Yan committed
690
    drl_export.clicked.connect(lambda: export_drlspectrum(
3Yan's avatar
3Yan committed
691
692
693
694
695
696
697
698
        main_widget, filename, ds, drls))

    drls['pt'].itemChanged.connect(lambda item: ptable_changed(
        item.row(), item.column(), ds, drls, chromplot))
    drls['dt'].itemChanged.connect(lambda item: dtable_changed(
        item.row(), item.column(), ds, drls, chromplot))
    drls['tshift'].valueChanged.connect(lambda: update_drlspectrum(
        ds, drls, chromplot))
Yan's avatar
Yan committed
699
700
    drls['cors'].valueChanged.connect(lambda: corcount_changed(
        ds, drls, chromplot))
3Yan's avatar
3Yan committed
701
702
    drls['rel'].stateChanged.connect(lambda: update_drlspectrum(
        ds, drls, chromplot))
703
704
705
    drls['checkAll'].stateChanged.connect(lambda state:
        select_all_btn_up(ds, drls, chromplot, state))

706
707
708
709
710
711
712
713
714
    actionBar = QtWidgets.QToolBar(window)
    window.addToolBar(QtCore.Qt.BottomToolBarArea, actionBar)
    actionBar.setAllowedAreas(QtCore.Qt.BottomToolBarArea)
    actionBar.setFloatable(False)
    actionBar.setMovable(False)
    actionBar.setToolButtonStyle(QtCore.Qt.ToolButtonTextBesideIcon)
    actionBar.addWidget(drl_load)
    actionBar.addWidget(drl_save)
    actionBar.addWidget(drl_print)
Yan's avatar
Yan committed
715
    actionBar.addWidget(StretchWidget())
716
717
718
719
720
721
722
    actionBar.addWidget(drl_export)
    actionBar.addWidget(close)

    dtdock = QtWidgets.QDockWidget()
    dtdock.setWidget(QtWidgets.QWidget())
    dtdock_layout = QtWidgets.QVBoxLayout(dtdock.widget())
    dtdock.setWindowTitle("Corrected ions table")
Yan's avatar
Yan committed
723
724
725
726
727
728
729
    dt_butlayout = QtWidgets.QHBoxLayout()
    dt_butlayout.addWidget(drls['checkAll'])
    dt_butlayout.addWidget(StretchWidget())
    dt_butlayout.addWidget(QtWidgets.QLabel("Number of corrections:"))
    dt_butlayout.addWidget(drls['cors'])
    dt_butlayout.addWidget(StretchWidget())
    dtdock_layout.addLayout(dt_butlayout)
730
    dtdock_layout.addWidget(drls['dt'])
Yan's avatar
Yan committed
731
    window.addDockWidget(QtCore.Qt.BottomDockWidgetArea, dtdock)
732
733
734
735
736
737

    ptdock = QtWidgets.QDockWidget()
    ptdock.setWidget(QtWidgets.QWidget())
    ptdock_layout = QtWidgets.QVBoxLayout(ptdock.widget())
    ptdock.setWindowTitle("Raw ions table")
    ptdock_layout.addWidget(drls['pt'])
Yan's avatar
Yan committed
738
739
740
741
    pt_butlayout = QtWidgets.QHBoxLayout()
    pt_butlayout.addWidget(btn_add)
    pt_butlayout.addWidget(btn_rem)
    pt_butlayout.addStretch(0)
742
743
744
745
746
747
748
    ptdock_layout.addLayout(pt_butlayout)
    window.addDockWidget(QtCore.Qt.RightDockWidgetArea, ptdock)

    main_layout = QtWidgets.QVBoxLayout(main_widget)
    graphparams_layout = QtWidgets.QHBoxLayout()
    main_layout.addWidget(graph_canvas, stretch=1)
    main_layout.addLayout(graphparams_layout)
3Yan's avatar
3Yan committed
749
750
751
    graphparams_layout.addWidget(time_title)
    graphparams_layout.addWidget(drls['tshift'])
    graphparams_layout.addStretch(1)
3Yan's avatar
3Yan committed
752
753
    graphparams_layout.addWidget(drls['rel'])
    graphparams_layout.addStretch(1)
754
755


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