#!/usr/bin/python
# -*- coding: UTF-8 -*-
# BEGIN PYTHON 2/3 COMPATIBILITY BOILERPLATE
from __future__ import absolute_import
from __future__ import with_statement
from __future__ import division
from __future__ import nested_scopes
from __future__ import generators
from __future__ import unicode_literals
from __future__ import print_function
import sys
# more py2/3 compat
from neurotools.system import *
if sys.version_info<(3,):
from itertools import imap as map
# END PYTHON 2/3 COMPATIBILITY BOILERPLATE
import cgid.config
import numpy as np
import os
import pickle
############################################################################
# Gather usable unit info
# NOTE!! 1 INDEXED UNIT IDS !!
#
# defines a list of called units containing (session,area,uid) tuples
# execfile(expanduser('~/Workspace2/cgidanalysis_r2/grabunits.py'))
allunitsbysession = {
('SPK120918','PMv'): [165,138,155,226,48,108,118,53,101,51,37,
127,116,27,54,205,227,208,160,38,218,223,222,135,206,235,
50,168,19,184,174,198,39,4,188,185,22,187,171,114,87,150,
84,81,12,193,199,1,236,221,157,212,63],
('RUS120518','M1' ): [1,35,80,30,50,55,5,78,45,121,127,18,94,
9,122,106,118,64,65,99,102,75,107,110,68,98,88,91,73,
66,20,42,6],
('SPK120925','PMv'): [8,69,86,17,62,106,195,201,143,32,211,133,
58,212,196,76,83,210,26,103,190,139,171,163,160,122,18,
31,157,87,220,215,181,167,172,191,89,12,9,95,80,189,140,
25,22,54,94,178,108,1,85,41,179,159,52,35,82,113,166,
30,4,20],
('RUS120523','M1' ): [80,64,98,142,59,6,63,92,101,71,76,17,
146,74,106,105,81,119,118,43,87,128,124,129,38,133,123,
21,114,18,5,104,1,28,130],
('SPK120925','PMd'): [14,123,12,67,85,16,91,8,78,43,7,41,19,
106,109,64,11,30,104,86,34,98,57,112,70,62,113,102,103,
20,51,3,68,22,33,45,107,115,116,73,50,46,35],
('SPK120924','PMv'): [110,76,149,72,12,20,5,147,14,2,87,123,65,
31,102,193,17,114,39,42,135,108,150,10,104,51,36,55,120,
79,134,208,61,199,140,152,74,116,46,146,169,194,6,200,
214,155,188,50,1,90,22,162,130,26,18,175,159,221,166,184,
91,86,177,189,165,220,185,211,141],
('RUS120523','PMd'): [118,36,60,75,114,4,1,121,39,23,12,108,11,
117,59,22,119,81,71,37,99,56,34,104,125,30,54,92,112,65,
61,62,33,97,78,96,55],
('SPK120924','M1' ): [93,40,26,110,35,80,11,134,54,43,41,4,62,
100,90,6,107,112,18,92,16,83,49,30,82,61,87,56,25,53,
47,36,104],
('SPK120925','M1' ): [8,30,86,105,111,37,107,10,21,27,22,33,
51,95,85,59,46,52,19,98,63,13,15,90,56,81,109,80,77,
73,48,17,14,5,74,58,84],
('RUS120518','PMv'): [46,25,163,210,200,180,172,169,141,22,184,
100,36,90,34,18,161,145,21,136,6,2,111,29,49,165,62,50,
38,162,77,16,115,1,101,94,15],
('RUS120521','PMd'): [45,69,77,85,25,67,89,98,94,65,27,59,51,
26,42,30,19,99,31,13,62,93,74],
('SPK120924','PMd'): [25,115,23,92,20,99,100,90,7,81,8,111,
117,35,91,112,65,48,59,16,104,105,38,63,83,122,108,33,
107,75,72,53,3,11,69,114,19,46,45,97,102,120,80,41,42,
87,96,124],
('SPK120918','PMd'): [119,70,45,42,40,84,43,36,62,12,14,35,77,
49,4,126,74,98,109,69,50,22,122,54,115,82,88,55,32,97,
23,83,111,19,61,72,104,17,94,108,87,131,121],
('RUS120523','PMv'): [26,103,109,198,191,190,136,181,139,184,
202,199,65,36,32,174,149,92,59,23,172,41,80,165,5,164,
104,157,77,53,49,47,28,46,99,44,58,37,118,123,14,97,
148,17,86,22,94],
('RUS120521','PMv'): [38,43,196,191,155,166,185,4,186,147,161,
136,154,105,75,10,23,104,163,24,87,133,19,160,34,29,57,
2,158,47,56,15,110,81,45,41,35,156],
('RUS120521','M1' ): [105,57,35,55,51,42,100,64,93,62,102,86,
37,101,24,11,31,96,91,73,66,97,26,3,81,19,118,111],
('SPK120918','M1' ): [6,107,82,80,90,83,23,39,52,2,65,46,121,
79,96,102,118,21,49,86,93,73,19,55,59,70,33],
('RUS120518','PMd'): [95,98,8,68,110,94,44,105,86,30,69,75,21,
103,83,72,22,38,93,34,36,96,77,59,79,73,90,113,82,91,89]
}
units = []
for (session,area),uu in allunitsbysession.items():
units+=[(session,area,u) for u in uu]
# what script generated this?
[docs]def classification_results():
return pickle.load(open(cgid.config.CGID_PACKAGE + os.path.sep + 'classification_results.p','rb'))
colorcodes = np.array([[ 0. ,0.6,1.0],
[ 0.9,0. ,1.0]])
print('LOADED CLASSIFICATION RESULTS: CAUTION ZERO INDEXED (UNUSUAL,NOT LIKE MATLAB)')
acceptable = {('RUS120518','M1',65),
('RUS120518','M1',66),
('RUS120518','M1',68),
('RUS120518','M1',106),
('RUS120518','M1',107),
('RUS120518','M1',110),
('RUS120518','PMd',22),
('RUS120518','PMd',30),
('RUS120518','PMd',36),
('RUS120518','PMd',73),
('RUS120518','PMd',75),
('RUS120518','PMd',77),
('RUS120518','PMd',90),
('RUS120518','PMd',91),
('RUS120518','PMd',95),
('RUS120518','PMv',15),
('RUS120518','PMv',46),
('RUS120518','PMv',50),
('RUS120518','PMv',90),
('RUS120521','M1',24),
('RUS120521','M1',26),
('RUS120521','M1',31),
('RUS120521','M1',37),
('RUS120521','M1',62),
('RUS120521','M1',64),
('RUS120521','M1',66),
('RUS120521','M1',86),
('RUS120521','M1',91),
('RUS120521','M1',93),
('RUS120521','M1',96),
('RUS120521','M1',97),
('RUS120521','M1',101),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('RUS120521','PMd',27),
('RUS120521','PMd',31),
('RUS120521','PMd',42),
('RUS120521','PMd',45),
('RUS120521','PMd',59),
('RUS120521','PMd',62),
('RUS120521','PMd',67),
('RUS120521','PMd',85),
('RUS120521','PMd',98),
('RUS120521','PMd',99),
('RUS120521','PMv',23),
('RUS120521','PMv',34),
('RUS120521','PMv',57),
('RUS120521','PMv',87),
('RUS120523','M1',1),
('RUS120523','M1',28),
('RUS120523','M1',38),
('RUS120523','M1',43),
('RUS120523','M1',71),
('RUS120523','M1',74),
('RUS120523','M1',76),
('RUS120523','M1',81),
('RUS120523','M1',92),
('RUS120523','M1',106),
('RUS120523','M1',114),
('RUS120523','M1',118),
('RUS120523','M1',119),
('RUS120523','M1',123),
('RUS120523','M1',129),
('RUS120523','M1',130),
('RUS120523','M1',133),
('RUS120523','PMd',11),
('RUS120523','PMd',12),
('RUS120523','PMd',23),
('RUS120523','PMd',34),
('RUS120523','PMd',36),
('RUS120523','PMd',54),
('RUS120523','PMd',55),
('RUS120523','PMd',56),
('RUS120523','PMd',60),
('RUS120523','PMd',71),
('RUS120523','PMd',97),
('RUS120523','PMd',117),
('RUS120523','PMd',118),
('RUS120523','PMd',119),
('RUS120523','PMv',5),
('RUS120523','PMv',46),
('RUS120523','PMv',47),
('RUS120523','PMv',59),
('RUS120523','PMv',86),
('RUS120523','PMv',92),
('RUS120523','PMv',94),
('RUS120523','PMv',118),
('RUS120523','PMv',123),
('SPK120918','M1',2),
('SPK120918','M1',23),
('SPK120918','M1',80),
('SPK120918','PMd',12),
('SPK120918','PMd',23),
('SPK120918','PMd',35),
('SPK120918','PMd',36),
('SPK120918','PMd',43),
('SPK120918','PMd',50),
('SPK120918','PMd',55),
('SPK120918','PMd',62),
('SPK120918','PMd',70),
('SPK120918','PMd',84),
('SPK120918','PMd',88),
('SPK120918','PMd',109),
('SPK120918','PMd',122),
('SPK120918','PMv',51),
('SPK120918','PMv',53),
('SPK120918','PMv',84),
('SPK120918','PMv',116),
('SPK120918','PMv',185),
('SPK120918','PMv',188),
('SPK120918','PMv',199),
('SPK120918','PMv',205),
('SPK120918','PMv',208),
('SPK120918','PMv',221),
('SPK120918','PMv',222),
('SPK120918','PMv',223),
('SPK120918','PMv',236),
('SPK120924','M1',4),
('SPK120924','M1',25),
('SPK120924','M1',26),
('SPK120924','M1',43),
('SPK120924','M1',49),
('SPK120924','M1',53),
('SPK120924','M1',54),
('SPK120924','M1',56),
('SPK120924','M1',61),
('SPK120924','M1',62),
('SPK120924','M1',80),
('SPK120924','M1',83),
('SPK120924','M1',93),
('SPK120924','M1',100),
('SPK120924','M1',104),
('SPK120924','M1',110),
('SPK120924','PMd',7),
('SPK120924','PMd',16),
('SPK120924','PMd',20),
('SPK120924','PMd',23),
('SPK120924','PMd',25),
('SPK120924','PMd',33),
('SPK120924','PMd',35),
('SPK120924','PMd',41),
('SPK120924','PMd',45),
('SPK120924','PMd',53),
('SPK120924','PMd',59),
('SPK120924','PMd',69),
('SPK120924','PMd',72),
('SPK120924','PMd',75),
('SPK120924','PMd',90),
('SPK120924','PMd',92),
('SPK120924','PMd',104),
('SPK120924','PMd',105),
('SPK120924','PMd',112),
('SPK120924','PMd',114),
('SPK120924','PMd',115),
('SPK120924','PMv',31),
('SPK120924','PMv',50),
('SPK120924','PMv',51),
('SPK120924','PMv',87),
('SPK120924','PMv',91),
('SPK120924','PMv',102),
('SPK120924','PMv',110),
('SPK120924','PMv',140),
('SPK120924','PMv',159),
('SPK120924','PMv',166),
('SPK120924','PMv',185),
('SPK120924','PMv',188),
('SPK120924','PMv',220),
('SPK120924','PMv',221),
('SPK120925','M1',10),
('SPK120925','M1',13),
('SPK120925','M1',15),
('SPK120925','M1',19),
('SPK120925','M1',21),
('SPK120925','M1',30),
('SPK120925','M1',48),
('SPK120925','M1',56),
('SPK120925','M1',58),
('SPK120925','M1',59),
('SPK120925','M1',73),
('SPK120925','M1',74),
('SPK120925','M1',84),
('SPK120925','M1',85),
('SPK120925','M1',86),
('SPK120925','PMd',7),
('SPK120925','PMd',11),
('SPK120925','PMd',14),
('SPK120925','PMd',20),
('SPK120925','PMd',22),
('SPK120925','PMd',30),
('SPK120925','PMd',33),
('SPK120925','PMd',34),
('SPK120925','PMd',35),
('SPK120925','PMd',41),
('SPK120925','PMd',45),
('SPK120925','PMd',46),
('SPK120925','PMd',50),
('SPK120925','PMd',51),
('SPK120925','PMd',57),
('SPK120925','PMd',67),
('SPK120925','PMd',68),
('SPK120925','PMd',86),
('SPK120925','PMd',103),
('SPK120925','PMd',104),
('SPK120925','PMd',106),
('SPK120925','PMd',107),
('SPK120925','PMd',112),
('SPK120925','PMd',113),
('SPK120925','PMd',116),
('SPK120925','PMd',123),
('SPK120925','PMv',25),
('SPK120925','PMv',26),
('SPK120925','PMv',30),
('SPK120925','PMv',35),
('SPK120925','PMv',52),
('SPK120925','PMv',82),
('SPK120925','PMv',83),
('SPK120925','PMv',86),
('SPK120925','PMv',103),
('SPK120925','PMv',108),
('SPK120925','PMv',133),
('SPK120925','PMv',160),
('SPK120925','PMv',163),
('SPK120925','PMv',166),
('SPK120925','PMv',167),
('SPK120925','PMv',171),
('SPK120925','PMv',172),
('SPK120925','PMv',178),
('SPK120925','PMv',179),
('SPK120925','PMv',189),
('SPK120925','PMv',190),
('SPK120925','PMv',195),
('SPK120925','PMv',212),
('SPK120925','PMv',220)}
beta_rhythmic = {('RUS120518','M1',65),
('RUS120518','M1',66),
('RUS120518','M1',106),
('RUS120518','M1',107),
('RUS120518','M1',110),
('RUS120518','PMd',22),
('RUS120518','PMd',30),
('RUS120518','PMd',36),
('RUS120518','PMd',73),
('RUS120518','PMd',75),
('RUS120518','PMd',90),
('RUS120518','PMd',91),
('RUS120518','PMv',15),
('RUS120518','PMv',46),
('RUS120518','PMv',50),
('RUS120521','M1',24),
('RUS120521','M1',31),
('RUS120521','M1',37),
('RUS120521','M1',96),
('RUS120521','M1',101),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('RUS120521','PMd',31),
('RUS120521','PMd',42),
('RUS120521','PMd',45),
('RUS120521','PMd',59),
('RUS120521','PMd',85),
('RUS120521','PMd',99),
('RUS120521','PMv',34),
('RUS120521','PMv',57),
('RUS120523','M1',28),
('RUS120523','M1',106),
('RUS120523','M1',119),
('RUS120523','M1',123),
('RUS120523','M1',129),
('RUS120523','M1',130),
('RUS120523','PMd',12),
('RUS120523','PMd',23),
('RUS120523','PMd',54),
('RUS120523','PMd',55),
('RUS120523','PMd',56),
('RUS120523','PMd',71),
('RUS120523','PMd',117),
('RUS120523','PMd',118),
('RUS120523','PMv',5),
('RUS120523','PMv',46),
('RUS120523','PMv',47),
('RUS120523','PMv',59),
('RUS120523','PMv',118),
('SPK120918','M1',2),
('SPK120918','M1',80),
('SPK120918','PMd',23),
('SPK120918','PMd',35),
('SPK120918','PMd',43),
('SPK120918','PMd',50),
('SPK120918','PMd',55),
('SPK120918','PMd',62),
('SPK120918','PMd',84),
('SPK120918','PMd',88),
('SPK120918','PMd',109),
('SPK120918','PMd',122),
('SPK120918','PMv',51),
('SPK120918','PMv',188),
('SPK120918','PMv',199),
('SPK120918','PMv',236),
('SPK120924','M1',4),
('SPK120924','M1',56),
('SPK120924','M1',61),
('SPK120924','M1',62),
('SPK120924','M1',83),
('SPK120924','M1',110),
('SPK120924','PMd',7),
('SPK120924','PMd',20),
('SPK120924','PMd',33),
('SPK120924','PMd',35),
('SPK120924','PMd',53),
('SPK120924','PMd',59),
('SPK120924','PMd',72),
('SPK120924','PMd',92),
('SPK120924','PMd',104),
('SPK120924','PMd',105),
('SPK120924','PMd',112),
('SPK120924','PMv',87),
('SPK120924','PMv',185),
('SPK120924','PMv',220),
('SPK120924','PMv',221),
('SPK120925','M1',13),
('SPK120925','M1',15),
('SPK120925','M1',56),
('SPK120925','M1',58),
('SPK120925','M1',59),
('SPK120925','M1',74),
('SPK120925','M1',84),
('SPK120925','M1',85),
('SPK120925','M1',86),
('SPK120925','PMd',7),
('SPK120925','PMd',14),
('SPK120925','PMd',20),
('SPK120925','PMd',33),
('SPK120925','PMd',34),
('SPK120925','PMd',35),
('SPK120925','PMd',41),
('SPK120925','PMd',45),
('SPK120925','PMd',46),
('SPK120925','PMd',51),
('SPK120925','PMd',57),
('SPK120925','PMd',67),
('SPK120925','PMd',68),
('SPK120925','PMd',103),
('SPK120925','PMd',104),
('SPK120925','PMd',107),
('SPK120925','PMd',112),
('SPK120925','PMd',113),
('SPK120925','PMd',116),
('SPK120925','PMd',123),
('SPK120925','PMv',35),
('SPK120925','PMv',52),
('SPK120925','PMv',82),
('SPK120925','PMv',86),
('SPK120925','PMv',160),
('SPK120925','PMv',167),
('SPK120925','PMv',172),
('SPK120925','PMv',178),
('SPK120925','PMv',190),
('SPK120925','PMv',195)}
'''
Variable storing all units exhibiting β-LFP rhythmicity
'''
'''
# more units
beta_rhythmic = {('RUS120518','M1',65),
('RUS120518','M1',66),
('RUS120518','M1',75),
('RUS120518','M1',106),
('RUS120518','M1',107),
('RUS120518','M1',110),
('RUS120518','PMd',21),
('RUS120518','PMd',22),
('RUS120518','PMd',30),
('RUS120518','PMd',36),
('RUS120518','PMd',59),
('RUS120518','PMd',72),
('RUS120518','PMd',73),
('RUS120518','PMd',75),
('RUS120518','PMd',89),
('RUS120518','PMd',90),
('RUS120518','PMd',91),
('RUS120518','PMv',15),
('RUS120518','PMv',16),
('RUS120518','PMv',46),
('RUS120518','PMv',50),
('RUS120521','M1',24),
('RUS120521','M1',31),
('RUS120521','M1',37),
('RUS120521','M1',96),
('RUS120521','M1',101),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('RUS120521','PMd',31),
('RUS120521','PMd',42),
('RUS120521','PMd',59),
('RUS120521','PMd',85),
('RUS120521','PMd',99),
('RUS120521','PMv',34),
('RUS120521','PMv',57),
('RUS120523','M1',28),
('RUS120523','M1',106),
('RUS120523','M1',119),
('RUS120523','M1',123),
('RUS120523','M1',129),
('RUS120523','M1',130),
('RUS120523','PMd',12),
('RUS120523','PMd',23),
('RUS120523','PMd',54),
('RUS120523','PMd',55),
('RUS120523','PMd',56),
('RUS120523','PMd',62),
('RUS120523','PMd',71),
('RUS120523','PMd',117),
('RUS120523','PMd',118),
('RUS120523','PMv',5),
('RUS120523','PMv',36),
('RUS120523','PMv',46),
('RUS120523','PMv',47),
('RUS120523','PMv',58),
('RUS120523','PMv',59),
('RUS120523','PMv',118),
('RUS120523','PMv',149),
('SPK120918','M1',2),
('SPK120918','M1',80),
('SPK120918','PMd',23),
('SPK120918','PMd',35),
('SPK120918','PMd',43),
('SPK120918','PMd',50),
('SPK120918','PMd',55),
('SPK120918','PMd',62),
('SPK120918','PMd',83),
('SPK120918','PMd',84),
('SPK120918','PMd',88),
('SPK120918','PMd',109),
('SPK120918','PMd',122),
('SPK120918','PMv',38),
('SPK120918','PMv',51),
('SPK120918','PMv',188),
('SPK120918','PMv',199),
('SPK120918','PMv',236),
('SPK120924','M1',4),
('SPK120924','M1',16),
('SPK120924','M1',56),
('SPK120924','M1',61),
('SPK120924','M1',62),
('SPK120924','M1',83),
('SPK120924','M1',110),
('SPK120924','PMd',7),
('SPK120924','PMd',8),
('SPK120924','PMd',20),
('SPK120924','PMd',33),
('SPK120924','PMd',35),
('SPK120924','PMd',42),
('SPK120924','PMd',53),
('SPK120924','PMd',59),
('SPK120924','PMd',72),
('SPK120924','PMd',81),
('SPK120924','PMd',92),
('SPK120924','PMd',104),
('SPK120924','PMd',105),
('SPK120924','PMd',112),
('SPK120924','PMv',87),
('SPK120924','PMv',185),
('SPK120924','PMv',220),
('SPK120924','PMv',221),
('SPK120925','M1',13),
('SPK120925','M1',14),
('SPK120925','M1',15),
('SPK120925','M1',56),
('SPK120925','M1',58),
('SPK120925','M1',59),
('SPK120925','M1',74),
('SPK120925','M1',84),
('SPK120925','M1',85),
('SPK120925','M1',86),
('SPK120925','PMd',3),
('SPK120925','PMd',7),
('SPK120925','PMd',8),
('SPK120925','PMd',14),
('SPK120925','PMd',20),
('SPK120925','PMd',33),
('SPK120925','PMd',34),
('SPK120925','PMd',35),
('SPK120925','PMd',41),
('SPK120925','PMd',45),
('SPK120925','PMd',46),
('SPK120925','PMd',51),
('SPK120925','PMd',57),
('SPK120925','PMd',67),
('SPK120925','PMd',68),
('SPK120925','PMd',98),
('SPK120925','PMd',103),
('SPK120925','PMd',104),
('SPK120925','PMd',107),
('SPK120925','PMd',112),
('SPK120925','PMd',113),
('SPK120925','PMd',116),
('SPK120925','PMd',123),
('SPK120925','PMv',35),
('SPK120925','PMv',52),
('SPK120925','PMv',82),
('SPK120925','PMv',86),
('SPK120925','PMv',160),
('SPK120925','PMv',167),
('SPK120925','PMv',172),
('SPK120925','PMv',178),
('SPK120925','PMv',190),
('SPK120925','PMv',195)}
'''
unimodal = {('RUS120518','M1',65),
('RUS120518','M1',66),
('RUS120518','M1',68),
('RUS120518','M1',106),
('RUS120518','M1',107),
('RUS120518','M1',110),
('RUS120518','PMd',22),
('RUS120518','PMd',30),
('RUS120518','PMd',36),
('RUS120518','PMd',73),
('RUS120518','PMd',75),
('RUS120518','PMd',77),
('RUS120518','PMd',90),
('RUS120518','PMd',91),
('RUS120518','PMv',15),
('RUS120518','PMv',46),
('RUS120518','PMv',50),
('RUS120518','PMv',90),
('RUS120521','M1',24),
('RUS120521','M1',31),
('RUS120521','M1',37),
('RUS120521','M1',66),
('RUS120521','M1',86),
('RUS120521','M1',91),
('RUS120521','M1',93),
('RUS120521','M1',96),
('RUS120521','M1',101),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('RUS120521','PMd',31),
('RUS120521','PMd',42),
('RUS120521','PMd',45),
('RUS120521','PMd',59),
('RUS120521','PMd',62),
('RUS120521','PMd',67),
('RUS120521','PMd',85),
('RUS120521','PMd',98),
('RUS120521','PMd',99),
('RUS120521','PMv',23),
('RUS120521','PMv',34),
('RUS120521','PMv',57),
('RUS120521','PMv',87),
('RUS120523','M1',1),
('RUS120523','M1',28),
('RUS120523','M1',38),
('RUS120523','M1',43),
('RUS120523','M1',74),
('RUS120523','M1',81),
('RUS120523','M1',92),
('RUS120523','M1',106),
('RUS120523','M1',114),
('RUS120523','M1',118),
('RUS120523','M1',119),
('RUS120523','M1',123),
('RUS120523','M1',129),
('RUS120523','M1',130),
('RUS120523','M1',133),
('RUS120523','PMd',11),
('RUS120523','PMd',12),
('RUS120523','PMd',23),
('RUS120523','PMd',54),
('RUS120523','PMd',55),
('RUS120523','PMd',56),
('RUS120523','PMd',71),
('RUS120523','PMd',97),
('RUS120523','PMd',117),
('RUS120523','PMd',118),
('RUS120523','PMv',5),
('RUS120523','PMv',46),
('RUS120523','PMv',47),
('RUS120523','PMv',59),
('RUS120523','PMv',86),
('RUS120523','PMv',92),
('RUS120523','PMv',94),
('RUS120523','PMv',118),
('RUS120523','PMv',123),
('SPK120918','M1',2),
('SPK120918','M1',80),
('SPK120918','PMd',12),
('SPK120918','PMd',23),
('SPK120918','PMd',35),
('SPK120918','PMd',43),
('SPK120918','PMd',50),
('SPK120918','PMd',55),
('SPK120918','PMd',62),
('SPK120918','PMd',84),
('SPK120918','PMd',88),
('SPK120918','PMd',109),
('SPK120918','PMd',122),
('SPK120918','PMv',51),
('SPK120918','PMv',53),
('SPK120918','PMv',84),
('SPK120918','PMv',116),
('SPK120918','PMv',185),
('SPK120918','PMv',188),
('SPK120918','PMv',199),
('SPK120918','PMv',205),
('SPK120918','PMv',208),
('SPK120918','PMv',236),
('SPK120924','M1',4),
('SPK120924','M1',25),
('SPK120924','M1',26),
('SPK120924','M1',49),
('SPK120924','M1',53),
('SPK120924','M1',56),
('SPK120924','M1',61),
('SPK120924','M1',62),
('SPK120924','M1',83),
('SPK120924','M1',93),
('SPK120924','M1',100),
('SPK120924','M1',104),
('SPK120924','M1',110),
('SPK120924','PMd',7),
('SPK120924','PMd',16),
('SPK120924','PMd',20),
('SPK120924','PMd',33),
('SPK120924','PMd',35),
('SPK120924','PMd',41),
('SPK120924','PMd',45),
('SPK120924','PMd',53),
('SPK120924','PMd',59),
('SPK120924','PMd',72),
('SPK120924','PMd',75),
('SPK120924','PMd',92),
('SPK120924','PMd',104),
('SPK120924','PMd',105),
('SPK120924','PMd',112),
('SPK120924','PMd',114),
('SPK120924','PMv',31),
('SPK120924','PMv',50),
('SPK120924','PMv',51),
('SPK120924','PMv',87),
('SPK120924','PMv',102),
('SPK120924','PMv',110),
('SPK120924','PMv',140),
('SPK120924','PMv',166),
('SPK120924','PMv',185),
('SPK120924','PMv',188),
('SPK120924','PMv',220),
('SPK120924','PMv',221),
('SPK120925','M1',10),
('SPK120925','M1',13),
('SPK120925','M1',15),
('SPK120925','M1',21),
('SPK120925','M1',30),
('SPK120925','M1',48),
('SPK120925','M1',56),
('SPK120925','M1',58),
('SPK120925','M1',59),
('SPK120925','M1',73),
('SPK120925','M1',74),
('SPK120925','M1',84),
('SPK120925','M1',85),
('SPK120925','M1',86),
('SPK120925','PMd',7),
('SPK120925','PMd',11),
('SPK120925','PMd',14),
('SPK120925','PMd',20),
('SPK120925','PMd',22),
('SPK120925','PMd',30),
('SPK120925','PMd',33),
('SPK120925','PMd',34),
('SPK120925','PMd',35),
('SPK120925','PMd',41),
('SPK120925','PMd',45),
('SPK120925','PMd',46),
('SPK120925','PMd',50),
('SPK120925','PMd',51),
('SPK120925','PMd',57),
('SPK120925','PMd',67),
('SPK120925','PMd',68),
('SPK120925','PMd',86),
('SPK120925','PMd',103),
('SPK120925','PMd',104),
('SPK120925','PMd',107),
('SPK120925','PMd',112),
('SPK120925','PMd',113),
('SPK120925','PMd',116),
('SPK120925','PMd',123),
('SPK120925','PMv',25),
('SPK120925','PMv',30),
('SPK120925','PMv',35),
('SPK120925','PMv',52),
('SPK120925','PMv',82),
('SPK120925','PMv',86),
('SPK120925','PMv',103),
('SPK120925','PMv',108),
('SPK120925','PMv',133),
('SPK120925','PMv',160),
('SPK120925','PMv',163),
('SPK120925','PMv',166),
('SPK120925','PMv',167),
('SPK120925','PMv',171),
('SPK120925','PMv',172),
('SPK120925','PMv',178),
('SPK120925','PMv',179),
('SPK120925','PMv',189),
('SPK120925','PMv',190),
('SPK120925','PMv',195),
('SPK120925','PMv',220)}
allunits = {(s,a,u) for (s,a),us in allunitsbysession.items() for u in us
}
thinthicklabeled = {('SPK120924','PMd',20),
('SPK120924','PMd',81),
('RUS120523','M1',71),
('SPK120924','PMv',76),
('SPK120918','M1',6),
('SPK120918','M1',107),
('SPK120925','PMd',12),
('SPK120925','M1',111),
('RUS120523','PMd',56),
('RUS120523','M1',124),
('RUS120523','PMd',33),
('SPK120924','PMv',135),
('SPK120925','PMd',85),
('SPK120925','PMv',95),
('SPK120924','M1',54),
('SPK120925','PMv',172),
('SPK120924','PMd',108),
('RUS120523','M1',38),
('RUS120523','PMv',181),
('RUS120521','PMd',19),
('RUS120518','PMv',2),
('SPK120925','M1',14),
('SPK120918','PMv',205),
('SPK120925','PMd',116),
('SPK120924','M1',100),
('SPK120918','PMv',38),
('SPK120925','PMv',211),
('SPK120924','PMv',159),
('RUS120523','M1',1),
('RUS120523','PMv',148),
('SPK120918','PMv',127),
('SPK120918','PMd',111),
('SPK120925','M1',17),
('SPK120918','PMv',236),
('RUS120518','PMd',93),
('SPK120918','PMd',17),
('SPK120924','PMd',100),
('SPK120918','M1',49),
('RUS120518','M1',64),
('SPK120925','PMv',8),
('SPK120925','PMv',157),
('RUS120523','M1',87),
('RUS120518','M1',121),
('SPK120918','PMv',227),
('SPK120918','PMd',98),
('SPK120925','M1',95),
('SPK120925','M1',48),
('RUS120518','M1',18),
('RUS120518','M1',50),
('SPK120918','M1',80),
('RUS120518','PMd',69),
('SPK120918','PMv',222),
('RUS120523','PMv',165),
('SPK120924','PMv',221),
('RUS120518','PMd',110),
('RUS120521','PMv',43),
('SPK120918','PMd',119),
('RUS120523','PMv',104),
('SPK120918','PMv',22),
('SPK120925','M1',56),
('SPK120925','PMv',163),
('SPK120924','PMv',175),
('RUS120523','M1',17),
('RUS120523','PMv',49),
('SPK120925','PMd',62),
('SPK120924','PMv',51),
('RUS120518','PMd',77),
('SPK120924','PMd',7),
('SPK120925','PMd',7),
('SPK120918','PMv',53),
('RUS120518','PMd',22),
('RUS120521','PMv',166),
('RUS120518','M1',9),
('SPK120925','PMv',86),
('SPK120925','PMv',171),
('SPK120924','PMd',117),
('SPK120918','PMd',36),
('RUS120523','PMv',86),
('SPK120918','PMv',212),
('RUS120518','PMv',145),
('SPK120918','PMv',185),
('RUS120521','PMd',51),
('RUS120518','PMv',34),
('RUS120518','PMv',180),
('RUS120518','PMd',30),
('RUS120518','M1',1),
('SPK120924','PMd',19),
('SPK120918','PMd',61),
('SPK120924','PMd',112),
('SPK120925','PMv',179),
('SPK120918','PMv',84),
('SPK120924','PMd',92),
('RUS120521','M1',55),
('SPK120925','PMd',46),
('RUS120523','PMv',94),
('SPK120924','PMv',86),
('RUS120521','PMd',59),
('RUS120523','PMv',202),
('SPK120918','M1',52),
('RUS120523','PMv',103),
('SPK120925','M1',86),
('SPK120918','PMv',37),
('SPK120924','M1',16),
('SPK120924','PMv',155),
('RUS120523','M1',119),
('SPK120918','PMv',138),
('RUS120521','M1',105),
('RUS120518','PMv',18),
('SPK120924','M1',47),
('SPK120924','PMv',214),
('RUS120523','PMd',81),
('RUS120521','M1',66),
('RUS120521','PMd',26),
('RUS120518','PMv',21),
('RUS120523','PMv',41),
('SPK120918','PMv',87),
('SPK120925','PMd',22),
('RUS120521','PMv',75),
('RUS120523','PMd',118),
('SPK120924','PMv',110),
('RUS120521','PMv',34),
('SPK120918','PMd',126),
('SPK120925','PMv',69),
('SPK120924','PMv',1),
('SPK120925','M1',51),
('SPK120918','PMd',87),
('RUS120521','M1',97),
('SPK120918','M1',83),
('SPK120925','PMd',68),
('SPK120918','PMv',118),
('RUS120518','M1',66),
('RUS120518','PMd',68),
('SPK120924','PMv',79),
('RUS120523','PMv',164),
('SPK120918','PMd',12),
('RUS120518','PMv',29),
('SPK120925','PMd',30),
('SPK120918','M1',46),
('SPK120924','PMv',102),
('RUS120518','M1',122),
('SPK120924','PMd',35),
('RUS120518','M1',55),
('SPK120925','PMd',103),
('RUS120523','PMv',23),
('SPK120925','M1',59),
('SPK120925','PMv',54),
('SPK120925','PMd',57),
('SPK120925','PMv',139),
('SPK120924','PMv',193),
('RUS120518','PMv',136),
('RUS120523','PMv',172),
('SPK120918','PMd',4),
('RUS120521','M1',31),
('SPK120925','PMv',31),
('RUS120523','PMv',118),
('SPK120925','M1',105),
('SPK120925','PMv',108),
('RUS120518','PMd',21),
('RUS120521','PMv',19),
('RUS120521','PMv',161),
('RUS120521','PMd',42),
('SPK120925','PMv',85),
('SPK120924','M1',56),
('SPK120924','PMv',17),
('RUS120518','PMd',83),
('RUS120521','M1',81),
('RUS120523','PMv',139),
('SPK120924','PMd',25),
('SPK120925','PMv',62),
('SPK120925','PMd',33),
('SPK120924','M1',87),
('SPK120925','M1',8),
('RUS120523','PMd',121),
('SPK120918','PMv',184),
('RUS120521','PMd',99),
('SPK120924','PMd',48),
('RUS120521','PMd',31),
('SPK120925','PMd',14),
('SPK120924','M1',110),
('RUS120518','M1',45),
('SPK120925','M1',81),
('SPK120924','M1',25),
('SPK120925','PMv',201),
('SPK120924','PMd',87),
('SPK120925','PMv',133),
('SPK120918','PMd',70),
('RUS120518','PMd',38),
('SPK120925','PMv',178),
('RUS120518','PMd',91),
('RUS120521','M1',73),
('SPK120925','PMv',210),
('SPK120918','M1',59),
('RUS120523','PMv',32),
('RUS120521','PMv',185),
('SPK120925','PMd',41),
('RUS120518','PMv',77),
('RUS120521','PMv',136),
('SPK120918','M1',33),
('RUS120518','M1',127),
('SPK120918','PMd',97),
('RUS120521','PMd',89),
('SPK120924','PMv',14),
('RUS120523','M1',118),
('RUS120523','PMd',59),
('RUS120521','PMv',23),
('RUS120521','PMv',105),
('SPK120918','M1',90),
('SPK120924','M1',40),
('SPK120925','M1',19),
('SPK120918','PMv',226),
('SPK120924','PMd',102),
('RUS120521','PMd',25),
('RUS120518','M1',94),
('RUS120521','PMd',27),
('SPK120924','PMd',23),
('SPK120918','PMd',131),
('RUS120523','M1',81),
('SPK120918','PMv',168),
('RUS120521','PMv',45),
('SPK120918','PMd',121),
('SPK120924','PMd',69),
('RUS120523','PMd',61),
('SPK120924','PMv',6),
('RUS120523','M1',142),
('SPK120924','PMv',149),
('SPK120918','PMd',69),
('SPK120918','M1',82),
('SPK120924','PMv',177),
('RUS120518','PMv',25),
('SPK120924','PMd',65),
('SPK120925','M1',27),
('RUS120523','PMd',104),
('RUS120518','PMd',75),
('RUS120521','PMv',147),
('SPK120924','PMv',211),
('SPK120925','PMv',22),
('RUS120523','PMv',77),
('RUS120523','PMd',119),
('SPK120918','PMv',51),
('RUS120523','M1',105),
('SPK120918','PMv',160),
('SPK120925','PMd',102),
('RUS120523','PMv',22),
('SPK120925','M1',73),
('SPK120925','PMv',76),
('SPK120925','M1',58),
('RUS120523','PMd',75),
('RUS120523','PMd',96),
('RUS120521','PMv',155),
('RUS120521','M1',57),
('RUS120518','M1',110),
('SPK120924','PMv',185),
('SPK120925','PMv',30),
('RUS120521','PMv',160),
('RUS120523','PMd',23),
('SPK120924','PMv',130),
('SPK120918','PMv',12),
('SPK120924','PMv',22),
('RUS120518','PMd',82),
('SPK120924','PMv',165),
('RUS120523','M1',43),
('SPK120924','PMv',200),
('RUS120521','PMv',35),
('SPK120918','PMv',101),
('SPK120924','M1',80),
('RUS120523','M1',64),
('RUS120518','M1',102),
('SPK120918','PMd',104),
('RUS120518','PMv',169),
('SPK120925','M1',80),
('RUS120518','PMv',22),
('SPK120924','M1',26),
('RUS120518','PMv',184),
('RUS120523','PMd',34),
('SPK120924','PMd',72),
('SPK120925','PMd',86),
('SPK120924','M1',134),
('SPK120918','PMv',4),
('SPK120924','M1',49),
('SPK120924','PMv',46),
('RUS120518','PMd',90),
('SPK120924','PMd',111),
('SPK120918','PMd',62),
('RUS120518','PMv',15),
('RUS120523','PMv',47),
('RUS120523','PMv',92),
('RUS120523','PMd',112),
('SPK120918','PMd',23),
('SPK120924','PMd',41),
('SPK120925','PMd',113),
('SPK120918','PMv',27),
('SPK120924','M1',18),
('RUS120518','PMv',111),
('RUS120518','M1',106),
('SPK120918','PMd',84),
('SPK120924','PMv',146),
('SPK120918','M1',93),
('SPK120924','PMd',53),
('RUS120523','PMv',14),
('SPK120924','M1',41),
('RUS120518','PMv',161),
('RUS120521','M1',64),
('RUS120523','M1',59),
('SPK120918','PMd',54),
('SPK120925','PMd',16),
('RUS120521','PMd',74),
('RUS120518','PMv',163),
('RUS120523','M1',80),
('RUS120518','M1',118),
('SPK120924','PMv',189),
('RUS120523','PMv',109),
('SPK120918','PMv',19),
('RUS120518','PMv',200),
('SPK120924','PMd',120),
('SPK120925','PMd',70),
('SPK120918','PMv',116),
('SPK120925','PMd',51),
('SPK120924','PMv',141),
('SPK120918','PMd',14),
('SPK120924','PMv',208),
('SPK120925','PMv',106),
('RUS120523','M1',104),
('RUS120521','PMd',94),
('SPK120925','PMv',83),
('RUS120518','M1',88),
('SPK120924','PMv',31),
('RUS120523','M1',129),
('SPK120925','PMv',160),
('SPK120924','PMv',162),
('RUS120523','M1',18),
('SPK120925','PMd',78),
('SPK120918','PMv',108),
('RUS120521','PMd',65),
('RUS120521','M1',96),
('RUS120523','PMd',99),
('RUS120521','PMv',154),
('SPK120918','PMv',199),
('SPK120918','M1',2),
('RUS120523','PMv',199),
('SPK120924','M1',112),
('SPK120925','M1',107),
('RUS120521','PMv',163),
('RUS120521','M1',118),
('RUS120518','PMv',94),
('SPK120918','PMd',72),
('RUS120521','PMd',67),
('SPK120925','M1',37),
('SPK120918','M1',19),
('SPK120924','PMd',104),
('RUS120518','PMv',6),
('SPK120925','PMd',35),
('SPK120925','M1',10),
('RUS120523','M1',146),
('RUS120523','PMd',12),
('RUS120521','PMd',69),
('SPK120925','PMd',8),
('SPK120924','M1',104),
('RUS120518','PMd',94),
('SPK120925','PMv',122),
('SPK120924','PMv',116),
('RUS120523','PMd',37),
('RUS120518','M1',78),
('RUS120518','M1',30),
('RUS120523','M1',5),
('SPK120924','PMd',99),
('RUS120523','PMv',5),
('SPK120925','PMv',35),
('RUS120518','PMd',36),
('SPK120924','PMv',65),
('RUS120518','PMd',89),
('SPK120924','PMd',96),
('RUS120518','M1',99),
('RUS120518','M1',68),
('RUS120523','PMv',46),
('SPK120925','PMv',4),
('SPK120925','PMd',43),
('SPK120918','PMv',193),
('RUS120518','PMd',98),
('RUS120523','PMd',4),
('SPK120918','PMd',22),
('SPK120924','PMd',42),
('RUS120521','PMd',77),
('SPK120925','PMd',112),
('SPK120924','PMv',12),
('SPK120925','PMv',215),
('SPK120924','PMv',147),
('SPK120918','PMd',88),
('RUS120518','M1',80),
('RUS120518','PMd',44),
('SPK120924','PMv',39),
('SPK120925','M1',21),
('SPK120924','PMv',74),
('SPK120918','PMd',49),
('SPK120924','PMd',11),
('SPK120925','PMv',12),
('SPK120925','PMd',19),
('SPK120925','PMv',89),
('RUS120518','PMv',162),
('SPK120918','PMv',174),
('RUS120521','PMv',47),
('RUS120521','PMd',85),
('RUS120518','M1',42),
('RUS120523','PMd',97),
('SPK120924','M1',11),
('SPK120925','M1',52),
('SPK120918','PMv',135),
('RUS120518','PMv',100),
('RUS120523','M1',21),
('RUS120518','PMv',50),
('RUS120523','PMv',53),
('SPK120925','PMd',50),
('RUS120518','PMd',73),
('SPK120918','PMd',77),
('RUS120518','PMv',141),
('SPK120924','PMd',3),
('SPK120925','PMv',20),
('SPK120925','M1',98),
('RUS120521','PMv',186),
('RUS120521','PMv',87),
('SPK120918','PMd',115),
('RUS120521','PMd',13),
('RUS120518','PMd',59),
('RUS120523','M1',128),
('SPK120925','PMv',167),
('SPK120924','PMd',33),
('SPK120924','PMd',90),
('SPK120918','PMd',40),
('SPK120925','PMd',73),
('RUS120521','M1',102),
('SPK120924','M1',90),
('SPK120924','PMv',55),
('SPK120925','M1',5),
('SPK120918','PMv',208),
('RUS120518','PMd',113),
('SPK120924','PMv',90),
('RUS120523','M1',74),
('RUS120518','PMv',38),
('RUS120521','PMd',93),
('RUS120518','PMd',8),
('SPK120925','PMd',3),
('SPK120925','PMv',113),
('RUS120518','M1',5),
('SPK120918','M1',79),
('SPK120925','PMd',104),
('RUS120523','PMv',28),
('RUS120523','PMd',62),
('SPK120924','PMv',20),
('SPK120925','PMv',82),
('SPK120924','PMd',105),
('RUS120523','PMv',198),
('SPK120918','PMd',32),
('RUS120518','M1',75),
('RUS120523','PMv',37),
('RUS120521','PMv',10),
('SPK120924','M1',82),
('SPK120925','M1',13),
('RUS120518','PMv',46),
('SPK120925','PMd',11),
('RUS120523','PMv',123),
('RUS120521','PMv',104),
('RUS120523','M1',123),
('SPK120918','PMv',150),
('SPK120924','PMv',152),
('RUS120521','PMv',57),
('RUS120521','M1',86),
('SPK120924','PMv',108),
('SPK120924','PMd',97),
('RUS120521','M1',51),
('SPK120918','PMd',19),
('SPK120924','PMv',194),
('RUS120523','PMd',114),
('RUS120518','PMv',165),
('SPK120918','M1',21),
('RUS120521','PMv',38),
('RUS120518','PMv',115),
('SPK120925','M1',46),
('SPK120918','M1',102),
('SPK120925','PMd',115),
('RUS120523','PMv',99),
('SPK120925','M1',90),
('RUS120523','PMd',60),
('RUS120521','M1',101),
('RUS120523','M1',28),
('SPK120918','PMv',114),
('SPK120924','M1',43),
('SPK120924','PMv',36),
('SPK120925','PMv',191),
('SPK120924','PMv',140),
('RUS120523','PMv',184),
('RUS120521','PMd',30),
('RUS120518','M1',91),
('SPK120924','PMd',122),
('RUS120518','PMd',105),
('RUS120518','PMv',62),
('SPK120918','PMd',122),
('SPK120925','PMd',123),
('RUS120521','PMv',56),
('SPK120924','M1',4),
('SPK120924','PMv',5),
('RUS120521','PMv',133),
('RUS120521','M1',93),
('SPK120918','PMd',83),
('SPK120925','PMd',64),
('SPK120924','M1',35),
('SPK120918','PMv',223),
('RUS120518','PMd',72),
('RUS120521','PMv',156),
('RUS120518','M1',6),
('RUS120521','M1',35),
('SPK120924','PMd',59),
('RUS120518','PMv',101),
('SPK120918','PMv',48),
('RUS120523','M1',106),
('SPK120918','PMv',165),
('SPK120924','PMd',91),
('SPK120918','M1',118),
('RUS120523','PMd',65),
('SPK120925','M1',74),
('RUS120518','PMv',210),
('SPK120925','M1',63),
('SPK120925','PMv',166),
('SPK120924','PMd',114),
('SPK120924','PMd',80),
('SPK120918','PMd',43),
('SPK120918','M1',39),
('SPK120925','PMv',58),
('SPK120925','PMv',143),
('SPK120924','PMv',91),
('RUS120518','M1',107),
('SPK120918','PMd',109),
('SPK120925','M1',109),
('SPK120924','PMv',114),
('RUS120523','M1',98),
('SPK120918','PMv',157),
('SPK120924','PMd',83),
('SPK120918','PMd',74),
('SPK120925','PMd',107),
('SPK120918','PMv',1),
('RUS120518','PMd',34),
('SPK120924','PMv',184),
('SPK120918','PMd',35),
('RUS120523','PMv',36),
('SPK120925','PMd',45),
('SPK120924','M1',83),
('RUS120523','PMd',125),
('SPK120924','PMv',199),
('SPK120918','PMv',81),
('SPK120918','PMv',188),
('RUS120521','M1',19),
('SPK120925','M1',85),
('SPK120924','PMv',10),
('SPK120925','PMv',140),
('RUS120523','PMd',39),
('SPK120924','PMd',75),
('SPK120918','M1',70),
('RUS120521','PMv',15),
('SPK120925','PMv',26),
('RUS120523','PMd',92),
('RUS120518','PMd',95),
('RUS120523','PMv',191),
('RUS120521','PMv',24),
('RUS120523','PMv',44),
('RUS120523','PMd',11),
('RUS120523','PMd',54),
('RUS120523','PMd',117),
('SPK120925','PMv',159),
('RUS120518','PMd',96),
('SPK120924','PMd',8),
('SPK120918','PMd',94),
('RUS120521','PMd',62),
('RUS120518','PMv',49),
('SPK120918','M1',121),
('SPK120924','PMd',63),
('RUS120521','M1',91),
('SPK120924','PMv',2),
('RUS120523','M1',114),
('RUS120518','M1',20),
('RUS120521','PMv',92),
('SPK120918','PMv',50),
('SPK120925','PMv',41),
('SPK120924','M1',36),
('SPK120925','PMv',190),
('SPK120924','PMv',72),
('SPK120918','M1',55),
('RUS120518','PMv',16),
('SPK120925','PMv',18),
('RUS120523','PMv',65),
('RUS120521','M1',11),
('SPK120918','PMv',63),
('RUS120518','PMv',172),
('RUS120523','PMv',136),
('RUS120521','PMv',41),
('SPK120924','PMv',26),
('RUS120523','PMd',55),
('RUS120518','PMv',1),
('SPK120924','PMv',169),
('RUS120521','M1',37),
('SPK120918','M1',86),
('SPK120925','PMd',67),
('SPK120924','M1',92),
('SPK120924','PMv',61),
('RUS120523','PMd',108),
('RUS120518','PMd',79),
('RUS120521','M1',24),
('RUS120523','M1',76),
('SPK120918','PMd',55),
('SPK120925','PMd',34),
('RUS120521','PMv',4),
('SPK120924','PMv',123),
('RUS120523','M1',101),
('RUS120523','PMd',22),
('RUS120521','PMv',81),
('SPK120918','M1',73),
('SPK120925','PMd',98),
('SPK120925','M1',77),
('SPK120925','PMv',80),
('SPK120924','M1',61),
('SPK120924','PMv',18),
('RUS120523','M1',130),
('SPK120924','PMd',115),
('SPK120918','PMd',42),
('RUS120523','PMv',59),
('RUS120521','M1',100),
('SPK120918','PMv',187),
('RUS120521','PMd',45),
('RUS120521','M1',26),
('SPK120918','PMd',108),
('SPK120925','PMv',94),
('SPK120924','M1',30),
('SPK120925','PMv',196),
('SPK120924','PMv',134),
('SPK120925','PMd',91),
('SPK120918','M1',65),
('SPK120925','PMd',106),
('RUS120523','PMv',26),
('SPK120924','PMv',87),
('SPK120924','M1',53),
('SPK120924','PMv',42),
('RUS120523','PMd',71),
('RUS120518','PMd',86),
('SPK120924','PMd',107),
('RUS120518','M1',73),
('SPK120918','PMv',235),
('SPK120925','PMv',1),
('RUS120523','PMv',80),
('SPK120925','M1',15),
('SPK120918','PMv',206),
('RUS120523','M1',92),
('RUS120523','PMd',1),
('RUS120518','M1',98),
('SPK120918','M1',96),
('SPK120924','M1',107),
('SPK120925','M1',84),
('SPK120918','PMv',39),
('SPK120925','PMv',195),
('RUS120523','M1',6),
('RUS120523','PMv',149),
('SPK120925','PMv',52),
('SPK120925','PMv',32),
('SPK120924','PMd',16),
('SPK120925','PMv',181),
('RUS120521','PMv',29),
('RUS120521','PMv',110),
('RUS120523','M1',63),
('RUS120523','PMv',190),
('RUS120518','M1',65),
('SPK120918','PMd',82),
('SPK120925','PMv',9),
('SPK120925','PMd',20),
('RUS120518','PMv',36),
('SPK120918','PMv',198),
('RUS120518','PMd',103),
('RUS120523','PMd',36),
('SPK120924','PMv',104),
('SPK120918','PMv',171),
('SPK120918','M1',23),
('SPK120924','PMd',45),
('RUS120521','M1',111),
('RUS120523','PMv',97),
('RUS120521','M1',42),
('SPK120925','PMv',212),
('SPK120924','PMv',150),
('RUS120523','PMv',157),
('SPK120925','M1',22),
('SPK120925','PMv',189),
('SPK120918','PMd',50),
('SPK120924','PMv',220),
('SPK120924','PMv',188),
('SPK120925','PMv',17),
('RUS120518','PMv',90),
('RUS120521','PMv',191),
('SPK120925','PMv',103),
('SPK120924','M1',6),
('RUS120521','PMv',196),
('RUS120523','M1',133),
('SPK120925','PMv',220),
('SPK120924','PMd',124),
('SPK120918','PMd',45),
('SPK120924','PMv',166),
('SPK120924','M1',93),
('SPK120924','PMv',50),
('SPK120925','M1',30),
('SPK120918','PMv',221),
('RUS120521','PMv',158),
('RUS120523','PMv',174),
('RUS120523','PMd',30),
('SPK120925','PMv',25),
('SPK120918','PMv',54),
('SPK120924','PMv',120),
('SPK120918','PMv',155),
('SPK120924','PMd',38),
('RUS120521','M1',62),
('RUS120521','M1',3),
('RUS120521','PMd',98),
('SPK120925','PMd',109),
('RUS120523','PMv',17),
('SPK120925','PMv',87),
('SPK120924','M1',62),
('SPK120925','M1',33),
('RUS120523','PMd',78),
('RUS120518','M1',35),
('SPK120918','PMv',218),
('SPK120924','PMd',46),
('RUS120523','PMv',58)}
published={('RUS120521','PMd',45),
('SPK120924','M1',16),
('RUS120523','PMv',157),
('SPK120924','PMv',114),
('SPK120925','PMv',190),
('SPK120918','PMd',111),
('SPK120918','PMd',84),
('RUS120521','PMd',51),
('SPK120925','PMv',189),
('RUS120523','PMv',149),
('RUS120518','PMv',90),
('SPK120925','PMd',35),
('RUS120523','PMd',97),
('SPK120925','PMv',103),
('SPK120925','PMd',14),
('SPK120924','M1',110),
('SPK120924','PMv',42),
('SPK120918','PMd',119),
('RUS120521','PMv',56),
('SPK120924','PMv',185),
('SPK120918','PMv',116),
('SPK120924','M1',26),
('RUS120523','PMd',55),
('RUS120521','PMv',57),
('SPK120925','PMd',86),
('RUS120518','PMd',73),
('RUS120523','M1',92),
('RUS120518','PMv',141),
('SPK120924','PMd',35),
('SPK120918','PMv',205),
('SPK120924','PMd',53),
('RUS120518','PMd',36),
('SPK120924','PMd',20),
('SPK120924','PMv',102),
('SPK120925','PMd',41),
('RUS120523','PMv',22),
('RUS120523','PMv',47),
('SPK120925','PMv',133),
('RUS120518','PMv',101),
('SPK120924','PMv',193),
('RUS120521','PMv',87),
('RUS120523','PMv',46),
('RUS120523','M1',106),
('RUS120523','PMd',23),
('SPK120924','PMv',110),
('RUS120523','M1',1),
('RUS120523','PMv',148),
('SPK120918','PMd',23),
('RUS120518','M1',1),
('RUS120518','PMd',59),
('RUS120523','PMv',92),
('SPK120918','PMv',51),
('RUS120523','PMd',114),
('RUS120521','PMv',23),
('SPK120925','PMv',195),
('SPK120918','PMv',53),
('SPK120918','PMv',199),
('SPK120924','PMv',51),
('SPK120924','M1',112),
('SPK120925','PMd',20),
('SPK120925','PMv',108),
('SPK120925','PMv',30),
('SPK120918','PMv',208),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('SPK120918','PMd',88),
('RUS120523','PMv',59),
('RUS120518','PMv',38),
('SPK120918','PMd',55),
('RUS120523','PMv',58),
('SPK120925','M1',21),
('SPK120924','PMv',31),('SPK120924','PMd',81),
('RUS120521','PMd',67),
('SPK120918','PMd',72),
('RUS120523','PMv',139),
('RUS120521','M1',86),
('SPK120925','PMd',51),
('SPK120918','PMd',109),
('SPK120925','PMd',33),
('SPK120924','PMd',104),
('SPK120925','PMd',107),
('SPK120924','PMd',41),
('SPK120924','PMd',48),
('SPK120925','M1',10),
('RUS120521','PMd',13),
('SPK120924','M1',25),
('RUS120523','PMv',36),
('SPK120925','PMd',45),
('SPK120925','PMv',172),
('SPK120925','PMd',8),
('SPK120924','M1',104),
('RUS120523','M1',38),
('SPK120918','PMv',188),
('SPK120924','PMv',116),
('RUS120523','PMv',86),
('SPK120925','PMv',178),
('RUS120518','PMd',91),
('SPK120925','M1',14),
('SPK120924','PMv',188),
('SPK120925','M1',85),
('RUS120523','PMv',5),
('SPK120925','PMv',35),
('SPK120924','PMd',75),
('RUS120518','PMd',89),
('SPK120925','PMd',116),
('SPK120924','M1',100),
('SPK120918','PMv',38),
('RUS120518','M1',68),
('RUS120518','M1',66),
('RUS120523','M1',118),
('SPK120924','PMd',42),
('RUS120523','PMd',11),
('SPK120918','PMv',236),
('SPK120925','PMd',112),
('RUS120523','PMd',117),
('RUS120521','PMd',62),
('RUS120518','PMv',49),
('SPK120918','PMd',35),
('SPK120924','PMd',124),
('RUS120523','M1',114),
('RUS120521','PMv',35),
('RUS120523','M1',81),
('SPK120925','M1',48),
('SPK120918','PMv',84),
('SPK120918','M1',80),
('RUS120518','PMv',16),
('RUS120521','M1',96),
('RUS120521','PMd',85),
('SPK120924','PMv',221),
('RUS120518','PMd',75),
('SPK120925','PMd',50),
('SPK120925','PMd',67),
('SPK120924','PMd',3),
('SPK120925','M1',56),
('SPK120925','PMv',163),
('SPK120925','M1',98),
('SPK120925','M1',73),
('SPK120925','PMd',22),
('SPK120925','M1',58),
('RUS120518','PMd',77),
('SPK120925','PMv',82),
('SPK120924','PMd',7),
('RUS120521','PMv',81),
('RUS120523','PMd',12),
('SPK120925','PMd',7),
('RUS120523','PMd',96),
('RUS120521','PMd',99),
('SPK120925','PMd',98),
('RUS120518','PMd',22),
('SPK120925','M1',77),
('SPK120925','PMd',57),
('RUS120518','M1',110),
('RUS120523','M1',130),
('SPK120924','PMv',140),
('RUS120523','M1',74),
('SPK120925','PMv',86),
('SPK120924','PMd',59),
('SPK120925','PMv',171),
('SPK120925','PMd',3),
('SPK120924','M1',83),
('RUS120518','M1',106),
('SPK120918','PMd',108),
('RUS120523','M1',43),
('SPK120925','PMd',104),
('SPK120918','PMv',185),
('RUS120521','M1',24),
('SPK120924','PMd',105),
('RUS120523','PMd',62),
('RUS120518','PMd',30),
('RUS120518','M1',75),
('SPK120924','M1',53),
('SPK120925','PMd',34),
('RUS120523','PMd',71),
('RUS120521','PMd',31),
('SPK120925','M1',13),
('SPK120925','PMv',94),
('RUS120518','PMv',46),
('SPK120924','PMd',112),
('SPK120924','PMd',72),
('SPK120925','PMd',11),
('SPK120925','PMv',167),
('RUS120523','PMv',123),
('SPK120925','M1',15),
('SPK120924','PMd',16),
('SPK120925','M1',74),
('RUS120523','M1',123),
('SPK120924','M1',49),
('SPK120925','PMd',46),
('RUS120523','PMv',94),
('SPK120925','M1',84),
('SPK120918','PMd',62),
('RUS120518','PMv',15),
('RUS120521','PMd',59),
('SPK120924','PMd',97),
('SPK120925','M1',86),
('RUS120523','M1',119),
('RUS120523','PMd',56),
('SPK120925','PMd',113),
('RUS120518','M1',65),
('RUS120521','M1',66),
('SPK120924','PMd',8),
('SPK120924','PMd',45),
('SPK120924','M1',61),
('RUS120521','M1',101),
('RUS120523','M1',28),
('RUS120523','PMd',118),
('SPK120924','PMd',92),
('SPK120924','M1',35),
('RUS120518','PMv',50),
('SPK120925','PMv',191),
('RUS120518','PMd',90),
('RUS120518','PMd',21),
('SPK120924','PMv',1),
('SPK120925','M1',51),
('SPK120924','PMd',33),
('SPK120925','M1',22),
('SPK120925','PMd',68),
('SPK120918','PMd',50),
('SPK120924','PMv',220),
('SPK120918','PMd',122),
('SPK120925','PMd',123),
('SPK120918','PMd',12),
('SPK120924','M1',4),
('RUS120521','PMv',34),
('SPK120925','PMd',30),
('RUS120521','M1',93),
('RUS120523','PMv',118),
('SPK120918','PMd',83),
('RUS120521','M1',37),
('RUS120523','M1',133),
('SPK120925','PMv',220),
('RUS120521','M1',91),
('RUS120518','PMd',72),
('SPK120925','PMd',103),
('SPK120925','M1',59),
('SPK120924','M1',93),
('SPK120924','PMv',50),
('SPK120925','M1',30),
('SPK120924','PMv',87),
('SPK120925','PMv',54),
('RUS120523','M1',129),
('SPK120925','PMv',160),
('SPK120925','PMv',25),
('SPK120924','PMv',162),
('RUS120523','PMd',65),
('SPK120925','PMv',179),
('SPK120924','PMv',120),
('RUS120521','PMd',42),
('RUS120521','M1',31),
('SPK120925','PMv',52),
('SPK120925','PMv',166),
('SPK120924','PMd',114),
('RUS120521','PMd',98),
('SPK120925','M1',105),
('SPK120918','PMd',43),
('RUS120523','PMd',54),
('SPK120918','M1',2),
('SPK120924','M1',62),
('SPK120924','PMv',141),
('SPK120924','M1',56),
('SPK120924','PMv',166),
('SPK120925','PMv',85),
('RUS120518','M1',107)}
useable={('RUS120521','PMd',45),
('SPK120924','M1',16),
('RUS120523','PMv',157),
('SPK120924','PMv',114),
('SPK120925','PMv',190),
('SPK120918','PMd',111),
('SPK120918','PMd',84),
('RUS120521','PMd',51),
('SPK120925','PMv',189),
('RUS120523','PMv',149),
('RUS120518','PMv',90),
('SPK120925','PMd',35),
('RUS120523','PMd',97),
('SPK120925','PMv',103),
('SPK120925','PMd',14),
('SPK120924','M1',110),
('SPK120924','PMv',42),
('SPK120918','PMd',119),
('RUS120521','PMv',56),
('SPK120924','PMv',185),
('SPK120918','PMv',116),
('SPK120924','M1',26),
('RUS120523','PMd',55),
('RUS120521','PMv',57),
('SPK120925','PMd',86),
('RUS120518','PMd',73),
('RUS120523','M1',92),
('RUS120518','PMv',141),
('SPK120924','PMd',35),
('SPK120918','PMv',205),
('SPK120924','PMd',53),
('RUS120518','PMd',36),
('SPK120924','PMd',20),
('SPK120924','PMv',102),
('SPK120925','PMd',41),
('RUS120523','PMv',22),
('RUS120523','PMv',47),
('SPK120925','PMv',133),
('RUS120518','PMv',101),
('SPK120924','PMv',193),
('RUS120521','PMv',87),
('RUS120523','PMv',46),
('RUS120523','M1',106),
('RUS120523','PMd',23),
('SPK120924','PMv',110),
('RUS120523','M1',1),
('RUS120523','PMv',148),
('SPK120918','PMd',23),
('RUS120518','M1',1),
('RUS120518','PMd',59),
('RUS120523','PMv',92),
('SPK120918','PMv',51),
('RUS120523','PMd',114),
('RUS120521','PMv',23),
('SPK120925','PMv',195),
('SPK120918','PMv',53),
('SPK120918','PMv',199),
('SPK120924','PMv',51),
('SPK120924','M1',112),
('SPK120925','PMd',20),
('SPK120925','PMv',108),
('SPK120925','PMv',30),
('SPK120918','PMv',208),
('RUS120521','M1',102),
('RUS120521','PMd',26),
('SPK120918','PMd',88),
('RUS120523','PMv',59),
('RUS120518','PMv',38),
('SPK120918','PMd',55),
('RUS120523','PMv',58),
('SPK120925','M1',21),
('SPK120924','PMv',31),
('SPK120924','PMd',81),
('RUS120521','PMd',67),
('SPK120918','PMd',72),
('RUS120523','PMv',139),
('RUS120521','M1',86),
('SPK120925','PMd',51),
('SPK120918','PMd',109),
('SPK120925','PMd',33),
('SPK120924','PMd',104),
('SPK120925','PMd',107),
('SPK120924','PMd',41),
('SPK120924','PMd',48),
('SPK120925','M1',10),
('RUS120521','PMd',13),
('SPK120924','M1',25),
('RUS120523','PMv',36),
('SPK120925','PMd',45),
('SPK120925','PMv',172),
('SPK120925','PMd',8),
('SPK120924','M1',104),
('RUS120523','M1',38),
('SPK120918','PMv',188),
('SPK120924','PMv',116),
('RUS120523','PMv',86),
('SPK120925','PMv',178),
('RUS120518','PMd',91),
('SPK120925','M1',14),
('SPK120924','PMv',188),
('SPK120925','M1',85),
('RUS120523','PMv',5),
('SPK120925','PMv',35),
('SPK120924','PMd',75),
('RUS120518','PMd',89),
('SPK120925','PMd',116),
('SPK120924','M1',100),
('SPK120918','PMv',38),
('RUS120518','M1',68),
('RUS120518','M1',66),
('RUS120523','M1',118),
('SPK120924','PMd',42),
('RUS120523','PMd',11),
('SPK120918','PMv',236),
('SPK120925','PMd',112),
('RUS120523','PMd',117),
('RUS120521','PMd',62),
('RUS120518','PMv',49),
('SPK120918','PMd',35),
('SPK120924','PMd',124),
('RUS120523','M1',114),
('RUS120521','PMv',35),
('RUS120523','M1',81),
('SPK120925','M1',48),
('SPK120918','PMv',84),
('SPK120918','M1',80),
('RUS120518','PMv',16),
('RUS120521','M1',96),
('RUS120521','PMd',85),
('SPK120924','PMv',221),
('RUS120518','PMd',75),
('SPK120925','PMd',50),
('SPK120925','PMd',67),
('SPK120924','PMd',3),
('SPK120925','M1',56),
('SPK120925','PMv',163),
('SPK120925','M1',98),
('SPK120925','M1',73),
('SPK120925','PMd',22),
('SPK120925','M1',58),
('RUS120518','PMd',77),
('SPK120925','PMv',82),
('SPK120924','PMd',7),
('RUS120521','PMv',81),
('RUS120523','PMd',12),
('SPK120925','PMd',7),
('RUS120523','PMd',96),
('RUS120521','PMd',99),
('SPK120925','PMd',98),
('RUS120518','PMd',22),
('SPK120925','M1',77),
('SPK120925','PMd',57),
('RUS120518','M1',110),
('RUS120523','M1',130),
('SPK120924','PMv',140),
('RUS120523','M1',74),
('SPK120925','PMv',86),
('SPK120924','PMd',59),
('SPK120925','PMv',171),
('SPK120925','PMd',3),
('SPK120924','M1',83),
('RUS120518','M1',106),
('SPK120918','PMd',108),
('RUS120523','M1',43),
('SPK120925','PMd',104),
('SPK120918','PMv',185),
('RUS120521','M1',24),
('SPK120924','PMd',105),
('RUS120523','PMd',62),
('RUS120518','PMd',30),
('RUS120518','M1',75),
('SPK120924','M1',53),
('SPK120925','PMd',34),
('RUS120523','PMd',71),
('RUS120521','PMd',31),
('SPK120925','M1',13),
('SPK120925','PMv',94),
('RUS120518','PMv',46),
('SPK120924','PMd',112),
('SPK120924','PMd',72),
('SPK120925','PMd',11),
('SPK120925','PMv',167),
('RUS120523','PMv',123),
('SPK120925','M1',15),
('SPK120924','PMd',16),
('SPK120925','M1',74),
('RUS120523','M1',123),
('SPK120924','M1',49),
('SPK120925','PMd',46),
('RUS120523','PMv',94),
('SPK120925','M1',84),
('SPK120918','PMd',62),
('RUS120518','PMv',15),
('RUS120521','PMd',59),
('SPK120924','PMd',97),
('SPK120925','M1',86),
('RUS120523','M1',119),
('RUS120523','PMd',56),
('SPK120925','PMd',113),
('RUS120518','M1',65),
('RUS120521','M1',66),
('SPK120924','PMd',8),
('SPK120924','PMd',45),
('SPK120924','M1',61),
('RUS120521','M1',101),
('RUS120523','M1',28),
('RUS120523','PMd',118),
('SPK120924','PMd',92),
('SPK120924','M1',35),
('RUS120518','PMv',50),
('SPK120925','PMv',191),
('RUS120518','PMd',90),
('RUS120518','PMd',21),
('SPK120924','PMv',1),
('SPK120925','M1',51),
('SPK120924','PMd',33),
('SPK120925','M1',22),
('SPK120925','PMd',68),
('SPK120918','PMd',50),
('SPK120924','PMv',220),
('SPK120918','PMd',122),
('SPK120925','PMd',123),
('SPK120918','PMd',12),
('SPK120924','M1',4),
('RUS120521','PMv',34),
('SPK120925','PMd',30),
('RUS120521','M1',93),
('RUS120523','PMv',118),
('SPK120918','PMd',83),
('RUS120521','M1',37),
('RUS120523','M1',133),
('SPK120925','PMv',220),
('RUS120521','M1',91),
('RUS120518','PMd',72),
('SPK120925','PMd',103),
('SPK120925','M1',59),
('SPK120924','M1',93),
('SPK120924','PMv',50),
('SPK120925','M1',30),
('SPK120924','PMv',87),
('SPK120925','PMv',54),
('RUS120523','M1',129),
('SPK120925','PMv',160),
('SPK120925','PMv',25),
('SPK120924','PMv',162),
('RUS120523','PMd',65),
('SPK120925','PMv',179),
('SPK120924','PMv',120),
('RUS120521','PMd',42),
('RUS120521','M1',31),
('SPK120925','PMv',52),
('SPK120925','PMv',166),
('SPK120924','PMd',114),
('RUS120521','PMd',98),
('SPK120925','M1',105),
('SPK120918','PMd',43),
('RUS120523','PMd',54),
('SPK120918','M1',2),
('SPK120924','M1',62),
('SPK120924','PMv',141),
('SPK120924','M1',56),
('SPK120924','PMv',166),
('SPK120925','PMv',85),
('RUS120518','M1',107),
('SPK120924','PMv',74),
('SPK120924','M1',43),
('SPK120918','PMv',222),
('SPK120924','M1',54),
('RUS120523','M1',76),
('SPK120925','M1',37),
('SPK120924','PMd',90),
('SPK120924','M1',80),
('RUS120523','PMd',36),
('SPK120925','PMv',212),
('RUS120521','PMv',185),
('SPK120924','PMd',23),
('SPK120924','M1',40),
('SPK120918','PMd',70),
('SPK120924','PMv',147),
('RUS120523','PMv',199),
('RUS120523','PMd',60),
('RUS120518','PMd',95),
('RUS120521','PMv',4),
('SPK120925','PMv',210),
('SPK120918','PMv',223),
('RUS120523','PMd',39),
('SPK120924','PMd',25),
('RUS120521','PMd',27),
('SPK120918','PMv',221),
('SPK120924','PMd',115),
('SPK120918','M1',23),
('SPK120925','PMd',106)}