-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathJPAS-diagram-4f.py
More file actions
1440 lines (1368 loc) · 85.6 KB
/
JPAS-diagram-4f.py
File metadata and controls
1440 lines (1368 loc) · 85.6 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
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
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
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
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
'''
Make color-color diagram for JPAS
'''
from __future__ import print_function
import numpy as np
import glob
import json
import matplotlib.pyplot as plt
import seaborn as sns
import sys
pattern = "*-spectros/*-JPAS17-magnitude.json"
file_list = glob.glob(pattern)
def filter_mag(e, s, f1, f2, f3, f4):
col, col0 = [], []
if data['id'].endswith(e):
if data['id'].startswith(str(s)):
filter1 = data[f1]
filter2 = data[f2]
filter3 = data[f3]
filter4 = data[f4]
diff = filter1 - filter2
diff0 = filter3 - filter4
col.append(diff)
col0.append(diff0)
return col, col0
def plot_mag(f1, f2, f3, f4):
x, y = filter_mag("HPNe", "", f1, f2, f3, f4)
x1, y1 = filter_mag("CV", "", f1, f2, f3, f4)
x2, y2 = filter_mag("E00", "DdDm1_L2", f1, f2, f3, f4)
x3, y3 = filter_mag("E00", "DdDm1_L3", f1, f2, f3, f4)
x4, y4 = filter_mag("E00", "DdDm1_L4", f1, f2, f3, f4)
x5, y5 = filter_mag("E00", "DdDm1_L5", f1, f2, f3, f4)
x6, y6 = filter_mag("E00", "N2242_L2", f1, f2, f3, f4)
x7, y7 = filter_mag("E00", "N2242_L3", f1, f2, f3, f4)
x8, y8 = filter_mag("E00", "N2242_L4", f1, f2, f3, f4)
x9, y9 = filter_mag("E00", "N2242_L5", f1, f2, f3, f4)
x10, y10 = filter_mag("E00", "K648_L2", f1, f2, f3, f4)
x11, y11 = filter_mag("E00", "K648_L3", f1, f2, f3, f4)
x12, y12 = filter_mag("E00", "K648_L4", f1, f2, f3, f4)
x13, y13 = filter_mag("E00", "K648_L5", f1, f2, f3, f4)
x14, y14 = filter_mag("E00", "BB1_L2", f1, f2, f3, f4)
x15, y15 = filter_mag("E00", "BB1_L3", f1, f2, f3, f4)
x16, y16 = filter_mag("E00", "BB1_L4", f1, f2, f3, f4)
x17, y17 = filter_mag("E00", "BB1_L5", f1, f2, f3, f4)
x18, y18 = filter_mag("E00", "Typ_L2", f1, f2, f3, f4)
x19, y19 = filter_mag("E00", "Typ_L3", f1, f2, f3, f4)
x20, y20 = filter_mag("E00", "Typ_L4", f1, f2, f3, f4)
x21, y21 = filter_mag("E00", "Typ_L5", f1, f2, f3, f4)
x22, y22 = filter_mag("E01", "DdDm1_L2", f1, f2, f3, f4)
x23, y23 = filter_mag("E01", "DdDm1_L3", f1, f2, f3, f4)
x24, y24 = filter_mag("E01", "DdDm1_L4", f1, f2, f3, f4)
x25, y25 = filter_mag("E01", "DdDm1_L5", f1, f2, f3, f4)
x26, y26 = filter_mag("E01", "N2242_L2", f1, f2, f3, f4)
x27, y27 = filter_mag("E01", "N2242_L3", f1, f2, f3, f4)
x28, y28 = filter_mag("E01", "N2242_L4", f1, f2, f3, f4)
x29, y29 = filter_mag("E01", "N2242_L5", f1, f2, f3, f4)
x30, y30 = filter_mag("E01", "K648_L2", f1, f2, f3, f4)
x31, y31 = filter_mag("E01", "K648_L3", f1, f2, f3, f4)
x32, y32 = filter_mag("E01", "K648_L4", f1, f2, f3, f4)
x33, y33 = filter_mag("E01", "K648_L5", f1, f2, f3, f4)
x34, y34 = filter_mag("E01", "BB1_L2", f1, f2, f3, f4)
x35, y35 = filter_mag("E01", "BB1_L3", f1, f2, f3, f4)
x36, y36 = filter_mag("E01", "BB1_L4", f1, f2, f3, f4)
x37, y37 = filter_mag("E01", "BB1_L5", f1, f2, f3, f4)
x38, y38 = filter_mag("E01", "Typ_L2", f1, f2, f3, f4)
x39, y39 = filter_mag("E01", "Typ_L3", f1, f2, f3, f4)
x40, y40 = filter_mag("E01", "Typ_L4", f1, f2, f3, f4)
x41, y41 = filter_mag("E01", "Typ_L5", f1, f2, f3, f4)
x42, y42 = filter_mag("E02", "DdDm1_L2", f1, f2, f3, f4)
x43, y43 = filter_mag("E02", "DdDm1_L3", f1, f2, f3, f4)
x44, y44 = filter_mag("E02", "DdDm1_L4", f1, f2, f3, f4)
x45, y45 = filter_mag("E02", "DdDm1_L5", f1, f2, f3, f4)
x46, y46 = filter_mag("E02", "N2242_L2", f1, f2, f3, f4)
x47, y47 = filter_mag("E02", "N2242_L3", f1, f2, f3, f4)
x48, y48 = filter_mag("E02", "N2242_L4", f1, f2, f3, f4)
x49, y49 = filter_mag("E02", "N2242_L5", f1, f2, f3, f4)
x50, y50 = filter_mag("E02", "K648_L2", f1, f2, f3, f4)
x51, y51 = filter_mag("E02", "K648_L3", f1, f2, f3, f4)
x52, y52 = filter_mag("E02", "K648_L4", f1, f2, f3, f4)
x53, y53 = filter_mag("E02", "K648_L5", f1, f2, f3, f4)
x54, y54 = filter_mag("E02", "BB1_L2", f1, f2, f3, f4)
x55, y55 = filter_mag("E02", "BB1_L3", f1, f2, f3, f4)
x56, y56 = filter_mag("E02", "BB1_L4", f1, f2, f3, f4)
x57, y57 = filter_mag("E02", "BB1_L5", f1, f2, f3, f4)
x58, y58 = filter_mag("E02", "Typ_L2", f1, f2, f3, f4)
x59, y59 = filter_mag("E02", "Typ_L3", f1, f2, f3, f4)
x60, y60 = filter_mag("E02", "Typ_L4", f1, f2, f3, f4)
x61, y61 = filter_mag("E02", "Typ_L5", f1, f2, f3, f4)
x62, y62 = filter_mag("E00_900", "DdDm1_L2", f1, f2, f3, f4)
x63, y63 = filter_mag("E00_900", "DdDm1_L3", f1, f2, f3, f4)
x64, y64 = filter_mag("E00_900", "DdDm1_L4", f1, f2, f3, f4)
x65, y65 = filter_mag("E00_900", "DdDm1_L5", f1, f2, f3, f4)
x66, y66 = filter_mag("E00_900", "N2242_L2", f1, f2, f3, f4)
x67, y67 = filter_mag("E00_900", "N2242_L3", f1, f2, f3, f4)
x68, y68 = filter_mag("E00_900", "N2242_L4", f1, f2, f3, f4)
x69, y69 = filter_mag("E00_900", "N2242_L5", f1, f2, f3, f4)
x70, y70 = filter_mag("E0_900", "K648_L2", f1, f2, f3, f4)
x71, y71 = filter_mag("E00_900", "K648_L3", f1, f2, f3, f4)
x72, y72 = filter_mag("E00_900", "K648_L4", f1, f2, f3, f4)
x73, y73 = filter_mag("E00_900", "K648_L5", f1, f2, f3, f4)
x74, y74 = filter_mag("E00_900", "BB1_L2", f1, f2, f3, f4)
x75, y75 = filter_mag("E00_900", "BB1_L3", f1, f2, f3, f4)
x76, y76 = filter_mag("E00_900", "BB1_L4", f1, f2, f3, f4)
x77, y77 = filter_mag("E00_900", "BB1_L5", f1, f2, f3, f4)
x78, y78 = filter_mag("E00_900", "Typ_L2", f1, f2, f3, f4)
x79, y79 = filter_mag("E00_900", "Typ_L3", f1, f2, f3, f4)
x80, y80 = filter_mag("E00_900", "Typ_L4", f1, f2, f3, f4)
x81, y81 = filter_mag("E00_900", "Typ_L5", f1, f2, f3, f4)
x82, y82 = filter_mag("E01_900", "DdDm1_L2", f1, f2, f3, f4)
x83, y83 = filter_mag("E01_900", "DdDm1_L3", f1, f2, f3, f4)
x84, y84 = filter_mag("E01_900", "DdDm1_L4", f1, f2, f3, f4)
x85, y85 = filter_mag("E01_900", "DdDm1_L5", f1, f2, f3, f4)
x86, y86 = filter_mag("E01_900", "N2242_L2", f1, f2, f3, f4)
x87, y87 = filter_mag("E01_900", "N2242_L3", f1, f2, f3, f4)
x88, y88 = filter_mag("E01_900", "N2242_L4", f1, f2, f3, f4)
x89, y89 = filter_mag("E01_900", "N2242_L5", f1, f2, f3, f4)
x90, y90 = filter_mag("E01_900", "K648_L2", f1, f2, f3, f4)
x91, y91 = filter_mag("E01_900", "K648_L3", f1, f2, f3, f4)
x92, y92 = filter_mag("E01_900", "K648_L4", f1, f2, f3, f4)
x93, y93 = filter_mag("E01_900", "K648_L5", f1, f2, f3, f4)
x94, y94 = filter_mag("E01_900", "BB1_L2", f1, f2, f3, f4)
x95, y95 = filter_mag("E01_900", "BB1_L3", f1, f2, f3, f4)
x96, y96 = filter_mag("E01_900", "BB1_L4", f1, f2, f3, f4)
x97, y97 = filter_mag("E01_900", "BB1_L5", f1, f2, f3, f4)
x98, y98 = filter_mag("E01_900", "Typ_L2", f1, f2, f3, f4)
x99, y99 = filter_mag("E01_900", "Typ_L3", f1, f2, f3, f4)
x100, y100 = filter_mag("E01_900", "Typ_L4", f1, f2, f3, f4)
x101, y101 = filter_mag("E01_900", "Typ_L5", f1, f2, f3, f4)
x102, y102 = filter_mag("E02_900", "DdDm1_L2", f1, f2, f3, f4)
x103, y103 = filter_mag("E02_900", "DdDm1_L3", f1, f2, f3, f4)
x104, y104 = filter_mag("E02_900", "DdDm1_L4", f1, f2, f3, f4)
x105, y105 = filter_mag("E02_900", "DdDm1_L5", f1, f2, f3, f4)
x106, y106 = filter_mag("E02_900", "N2242_L2", f1, f2, f3, f4)
x107, y107 = filter_mag("E02_900", "N2242_L3", f1, f2, f3, f4)
x108, y108 = filter_mag("E02_900", "N2242_L4", f1, f2, f3, f4)
x109, y109 = filter_mag("E02_900", "N2242_L5", f1, f2, f3, f4)
x110, y110 = filter_mag("E02_900", "K648_L2", f1, f2, f3, f4)
x111, y111 = filter_mag("E02_900", "K648_L3", f1, f2, f3, f4)
x112, y112 = filter_mag("E02_900", "K648_L4", f1, f2, f3, f4)
x113, y113 = filter_mag("E02_900", "K648_L5", f1, f2, f3, f4)
x114, y114 = filter_mag("E02_900", "BB1_L2", f1, f2, f3, f4)
x115, y115 = filter_mag("E02_900", "BB1_L3", f1, f2, f3, f4)
x116, y116 = filter_mag("E02_900", "BB1_L4", f1, f2, f3, f4)
x117, y117 = filter_mag("E02_900", "BB1_L5", f1, f2, f3, f4)
x118, y118 = filter_mag("E02_900", "Typ_L2", f1, f2, f3, f4)
x119, y119 = filter_mag("E02_900", "Typ_L3", f1, f2, f3, f4)
x120, y120 = filter_mag("E02_900", "Typ_L4", f1, f2, f3, f4)
x121, y121 = filter_mag("E02_900", "Typ_L5", f1, f2, f3, f4)
x122, y122= filter_mag("-DPNe", "", f1, f2, f3, f4)
x123, y123= filter_mag("QSOs-13", "", f1, f2, f3, f4)
x124, y124 = filter_mag("QSOs-24", "", f1, f2, f3, f4)
x125, y125 = filter_mag("QSOs-32", "", f1, f2, f3, f4)
x126, y126 = filter_mag("-SFGs", "", f1, f2, f3, f4)
x127, y127 = filter_mag("-sys", "", f1, f2, f3, f4)
x128, y128 = filter_mag("-sys-IPHAS", "", f1, f2, f3, f4)
x129, y129 = filter_mag("-ExtHII", "", f1, f2, f3, f4)
x130, y130 = filter_mag("-sys-Ext", '', f1, f2, f3, f4)
x131, y131 = filter_mag("-survey", '', f1, f2, f3, f4)
x132, y132 = filter_mag("-SNR", '', f1, f2, f3, f4)
x133, y133 = filter_mag("extr-SySt-raman", '', f1, f2, f3, f4)
x134, y134 = filter_mag("-extr-SySt", '', f1, f2, f3, f4)
x135, y135 = filter_mag("-sys-raman", "", f1, f2, f3, f4)
x136, y136 = filter_mag("-sys-IPHAS-raman", "", f1, f2, f3, f4)
x137, y137 = filter_mag("ngc185-raman", "", f1, f2, f3, f4)
x138, y138 = filter_mag("SySt-ic10", "", f1, f2, f3, f4)
x139, y139 = filter_mag("LOAN-HPNe-", "", f1, f2, f3, f4)
x140, y140 = filter_mag("1359559-HPNe-", "", f1, f2, f3, f4)
x141, y141 = filter_mag("DdDm-1-HPNe-", "", f1, f2, f3, f4)
x142, y142 = filter_mag("QSOs-30", "", f1, f2, f3, f4)
for a, b in zip(x, y):
A1[0].append(a)
B1[0].append(b)
for a, b in zip(x1, y1):
A1[1].append(a)
B1[1].append(b)
for a, b in zip(x2, y2):
A1[2].append(a)
B1[2].append(b)
for a, b in zip(x3, y3):
A1[3].append(a)
B1[3].append(b)
for a, b in zip(x4, y4):
A1[4].append(a)
B1[4].append(b)
for a, b in zip(x5, y5):
A1[5].append(a)
B1[5].append(b)
for a, b in zip(x6, y6):
A1[6].append(a)
B1[6].append(b)
for a, b in zip(x7, y7):
A1[7].append(a)
B1[7].append(b)
for a, b in zip(x8, y8):
A1[8].append(a)
B1[8].append(b)
for a, b in zip(x9, y9):
A1[9].append(a)
B1[9].append(b)
for a, b in zip(x10, y10):
A1[10].append(a)
B1[10].append(b)
for a, b in zip(x11, y11):
A1[11].append(a)
B1[11].append(b)
for a, b in zip(x12, y12):
A1[12].append(a)
B1[12].append(b)
for a, b in zip(x13, y13):
A1[13].append(a)
B1[13].append(b)
for a, b in zip(x14, y14):
A1[14].append(a)
B1[14].append(b)
for a, b in zip(x15, y15):
A1[15].append(a)
B1[15].append(b)
for a, b in zip(x16, y16):
A1[16].append(a)
B1[16].append(b)
for a, b in zip(x17, y17):
A1[17].append(a)
B1[17].append(b)
for a, b in zip(x18, y18):
A1[18].append(a)
B1[18].append(b)
for a, b in zip(x19, y19):
A1[19].append(a)
B1[19].append(b)
for a, b in zip(x20, y20):
A1[20].append(a)
B1[20].append(b)
for a, b in zip(x21, y21):
A1[21].append(a)
B1[21].append(b)
for a, b in zip(x22, y22):
A1[22].append(a)
B1[22].append(b)
for a, b in zip(x23, y23):
A1[23].append(a)
B1[23].append(b)
for a, b in zip(x24, y24):
A1[24].append(a)
B1[24].append(b)
for a, b in zip(x25, y25):
A1[25].append(a)
B1[25].append(b)
for a, b in zip(x26, y26):
A1[26].append(a)
B1[26].append(b)
for a, b in zip(x27, y27):
A1[27].append(a)
B1[27].append(b)
for a, b in zip(x28, y28):
A1[28].append(a)
B1[28].append(b)
for a, b in zip(x29, y29):
A1[29].append(a)
B1[29].append(b)
for a, b in zip(x30, y30):
A1[30].append(a)
B1[30].append(b)
for a, b in zip(x31, y31):
A1[31].append(a)
B1[31].append(b)
for a, b in zip(x32, y32):
A1[32].append(a)
B1[32].append(b)
for a, b in zip(x33, y33):
A1[33].append(a)
B1[33].append(b)
for a, b in zip(x34, y34):
A1[34].append(a)
B1[34].append(b)
for a, b in zip(x35, y35):
A1[35].append(a)
B1[35].append(b)
for a, b in zip(x36, y36):
A1[36].append(a)
B1[36].append(b)
for a, b in zip(x37, y37):
A1[37].append(a)
B1[37].append(b)
for a, b in zip(x38, y38):
A1[38].append(a)
B1[38].append(b)
for a, b in zip(x39, y39):
A1[39].append(a)
B1[39].append(b)
for a, b in zip(x40, y40):
A1[40].append(a)
B1[40].append(b)
for a, b in zip(x41, y41):
A1[41].append(a)
B1[41].append(b)
for a, b in zip(x42, y42):
A1[42].append(a)
B1[42].append(b)
for a, b in zip(x43, y43):
A1[43].append(a)
B1[43].append(b)
for a, b in zip(x44, y44):
A1[44].append(a)
B1[44].append(b)
for a, b in zip(x45, y45):
A1[45].append(a)
B1[45].append(b)
for a, b in zip(x46, y46):
A1[46].append(a)
B1[46].append(b)
for a, b in zip(x47, y47):
A1[47].append(a)
B1[47].append(b)
for a, b in zip(x48, y48):
A1[48].append(a)
B1[48].append(b)
for a, b in zip(x49, y49):
A1[49].append(a)
B1[49].append(b)
for a, b in zip(x50, y50):
A1[50].append(a)
B1[50].append(b)
for a, b in zip(x51, y51):
A1[51].append(a)
B1[51].append(b)
for a, b in zip(x52, y52):
A1[52].append(a)
B1[52].append(b)
for a, b in zip(x53, y53):
A1[53].append(a)
B1[53].append(b)
for a, b in zip(x54, y54):
A1[54].append(a)
B1[54].append(b)
for a, b in zip(x55, y55):
A1[55].append(a)
B1[55].append(b)
for a, b in zip(x56, y56):
A1[56].append(a)
B1[56].append(b)
for a, b in zip(x57, y57):
A1[57].append(a)
B1[57].append(b)
for a, b in zip(x58, y58):
A1[58].append(a)
B1[58].append(b)
for a, b in zip(x59, y59):
A1[59].append(a)
B1[59].append(b)
for a, b in zip(x60, y60):
A1[60].append(a)
B1[60].append(b)
for a, b in zip(x61, y61):
A1[61].append(a)
B1[61].append(b)
for a, b in zip(x62, y62):
A1[62].append(a)
B1[62].append(b)
for a, b in zip(x63, y63):
A1[63].append(a)
B1[63].append(b)
for a, b in zip(x64, y64):
A1[64].append(a)
B1[64].append(b)
for a, b in zip(x65, y65):
A1[65].append(a)
B1[65].append(b)
for a, b in zip(x66, y66):
A1[66].append(a)
B1[66].append(b)
for a, b in zip(x67, y67):
A1[67].append(a)
B1[67].append(b)
for a, b in zip(x68, y68):
A1[68].append(a)
B1[68].append(b)
for a, b in zip(x69, y69):
A1[69].append(a)
B1[69].append(b)
for a, b in zip(x70, y70):
A1[70].append(a)
B1[70].append(b)
for a, b in zip(x72, y72):
A1[72].append(a)
B1[72].append(b)
for a, b in zip(x73, y73):
A1[73].append(a)
B1[73].append(b)
for a, b in zip(x74, y74):
A1[74].append(a)
B1[74].append(b)
for a, b in zip(x75, y75):
A1[75].append(a)
B1[75].append(b)
for a, b in zip(x76, y76):
A1[76].append(a)
B1[76].append(b)
for a, b in zip(x77, y77):
A1[77].append(a)
B1[77].append(b)
for a, b in zip(x78, y78):
A1[78].append(a)
B1[78].append(b)
for a, b in zip(x79, y79):
A1[79].append(a)
B1[79].append(b)
for a, b in zip(x80, y80):
A1[80].append(a)
B1[80].append(b)
for a, b in zip(x81, y81):
A1[81].append(a)
B1[81].append(b)
for a, b in zip(x82, y82):
A1[82].append(a)
B1[82].append(b)
for a, b in zip(x83, y83):
A1[83].append(a)
B1[83].append(b)
for a, b in zip(x84, y84):
A1[84].append(a)
B1[84].append(b)
for a, b in zip(x85, y85):
A1[85].append(a)
B1[85].append(b)
for a, b in zip(x86, y86):
A1[86].append(a)
B1[86].append(b)
for a, b in zip(x87, y87):
A1[87].append(a)
B1[87].append(b)
for a, b in zip(x88, y88):
A1[88].append(a)
B1[88].append(b)
for a, b in zip(x89, y89):
A1[89].append(a)
B1[89].append(b)
for a, b in zip(x90, y90):
A1[90].append(a)
B1[90].append(b)
for a, b in zip(x91, y91):
A1[91].append(a)
B1[91].append(b)
for a, b in zip(x92, y92):
A1[92].append(a)
B1[92].append(b)
for a, b in zip(x93, y93):
A1[93].append(a)
B1[93].append(b)
for a, b in zip(x94, y94):
A1[94].append(a)
B1[94].append(b)
for a, b in zip(x95, y95):
A1[95].append(a)
B1[95].append(b)
for a, b in zip(x96, y96):
A1[96].append(a)
B1[96].append(b)
for a, b in zip(x97, y97):
A1[97].append(a)
B1[97].append(b)
for a, b in zip(x98, y98):
A1[98].append(a)
B1[98].append(b)
for a, b in zip(x99, y99):
A1[99].append(a)
B1[99].append(b)
for a, b in zip(x100, y100):
A1[100].append(a)
B1[100].append(b)
for a, b in zip(x101, y101):
A1[101].append(a)
B1[101].append(b)
for a, b in zip(x102, y102):
A1[102].append(a)
B1[102].append(b)
for a, b in zip(x103, y103):
A1[103].append(a)
B1[103].append(b)
for a, b in zip(x104, y104):
A1[104].append(a)
B1[104].append(b)
for a, b in zip(x105, y105):
A1[105].append(a)
B1[105].append(b)
for a, b in zip(x106, y106):
A1[106].append(a)
B1[106].append(b)
for a, b in zip(x107, y107):
A1[107].append(a)
B1[107].append(b)
for a, b in zip(x108, y108):
A1[108].append(a)
B1[108].append(b)
for a, b in zip(x109, y109):
A1[109].append(a)
B1[109].append(b)
for a, b in zip(x110, y110):
A1[110].append(a)
B1[110].append(b)
for a, b in zip(x111, y111):
A1[111].append(a)
B1[111].append(b)
for a, b in zip(x112, y112):
A1[112].append(a)
B1[112].append(b)
for a, b in zip(x113, y113):
A1[113].append(a)
B1[113].append(b)
for a, b in zip(x114, y114):
A1[114].append(a)
B1[114].append(b)
for a, b in zip(x115, y115):
A1[115].append(a)
B1[115].append(b)
for a, b in zip(x116, y116):
A1[116].append(a)
B1[116].append(b)
for a, b in zip(x117, y117):
A1[117].append(a)
B1[117].append(b)
for a, b in zip(x118, y118):
A1[118].append(a)
B1[118].append(b)
for a, b in zip(x119, y119):
A1[119].append(a)
B1[119].append(b)
for a, b in zip(x120, y120):
A1[120].append(a)
B1[120].append(b)
for a, b in zip(x121, y121):
A1[121].append(a)
B1[121].append(b)
for a, b in zip(x122, y122):
A1[122].append(a)
B1[122].append(b)
for a, b in zip(x123, y123):
A1[123].append(a)
B1[123].append(b)
for a, b in zip(x124, y124):
A1[124].append(a)
B1[124].append(b)
for a, b in zip(x125, y125):
A1[125].append(a)
B1[125].append(b)
for a, b in zip(x126, y126):
A1[126].append(a)
B1[126].append(b)
for a, b in zip(x127, y127):
A1[127].append(a)
B1[127].append(b)
for a, b in zip(x128, y128):
A1[128].append(a)
B1[128].append(b)
for a, b in zip(x129, y129):
A1[129].append(a)
B1[129].append(b)
for a, b in zip(x130, y130):
A1[130].append(a)
B1[130].append(b)
for a, b in zip(x131, y131):
A1[131].append(a)
B1[131].append(b)
for a, b in zip(x132, y132):
A1[132].append(a)
B1[132].append(b)
for a, b in zip(x133, y133):
A1[133].append(a)
B1[133].append(b)
for a, b in zip(x134, y134):
A1[134].append(a)
B1[134].append(b)
for a, b in zip(x135, y135):
A1[135].append(a)
B1[135].append(b)
for a, b in zip(x136, y136):
A1[136].append(a)
B1[136].append(b)
for a, b in zip(x137, y137):
A1[137].append(a)
B1[137].append(b)
for a, b in zip(x138, y138):
A1[138].append(a)
B1[138].append(b)
for a, b in zip(x139, y139):
A1[139].append(a)
B1[139].append(b)
for a, b in zip(x140, y140):
A1[140].append(a)
B1[140].append(b)
for a, b in zip(x141, y141):
A1[141].append(a)
B1[141].append(b)
for a, b in zip(x142, y142):
A1[142].append(a)
B1[142].append(b)
label = []
label1 = ["H4-1", "PNG 135.9+55.9"]
n = 143
A1, B1 = [[] for _ in range(n)], [[] for _ in range(n)]
d_644_jplus, d_768_jplus = [], []
d_644_jplus1, d_768_jplus1 = [], []
for file_name in file_list:
with open(file_name) as f:
data = json.load(f)
# if data['id'].endswith("1-HPNe"):
# label.append("DdDm-1")
# elif data['id'].endswith("SLOAN-HPNe-"):
# label.append("H4-1")
# elif data['id'].endswith("1359559-HPNe"):
# label.append("PNG 135.9+55.9")
if data['id'].startswith("ngc"):
label.append("NGC 2242")
elif data['id'].startswith("mwc"):
label.append("MWC 574")
plot_mag("Jv_0915_g_4803", "Jv0915_5101", "Jv0915_8600", "Jv_0915_u_3485")
lgd_kws = {'frameon': True, 'fancybox': True, 'shadow': True}
#sns.set(style="dark")#, context="talk")
sns.set_style('ticks')
fig = plt.figure(figsize=(7, 6))
ax1 = fig.add_subplot(111)
ax1.set_xlim(xmin=-6.5,xmax=2.2)
# ax1.set_ylim(ymin=-2.5,ymax=2.0)
#ax1.set_xlim(xmin=-2.5,xmax=2.0)
plt.tick_params(axis='x', labelsize=18)
plt.tick_params(axis='y', labelsize=18)
plt.xlabel(r'$J8600 - uSDSS$', fontsize= 18)
plt.ylabel(r'$gSDSS - J5101$', fontsize= 18)
ax1.scatter(B1[62], A1[62], c= "orange", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral )
ax1.scatter(B1[63], A1[63], c= "orange", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[64], A1[64], c= "orange", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[65], A1[65], c= "orange", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[66], A1[66], c= "green", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[67], A1[67], c= "green", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[68], A1[68], c= "green", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[69], A1[69], c= "green", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[70], A1[70], c= "brown", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[71], A1[71], c= "brown", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[72], A1[72], c= "brown", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[73], A1[73], c= "brown", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[74], A1[74], c= "cyan", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[75], A1[75], c= "cyan", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[76], A1[76], c= "cyan", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[77], A1[77], c= "cyan", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[78], A1[78], c= "magenta", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[79], A1[79], c= "magenta", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[80], A1[80], c= "magenta", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[81], A1[81], c= "magenta", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[82], A1[82], c= "orange", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[83], A1[83], c= "orange", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[84], A1[84], c= "orange", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[85], A1[85], c= "orange", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[86], A1[86], c= "green", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[87], A1[87], c= "green", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[88], A1[88], c= "green", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[89], A1[89], c= "green", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[90], A1[90], c= "brown", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[91], A1[91], c= "brown", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[92], A1[92], c= "brown", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[93], A1[93], c= "brown", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[94], A1[94], c= "cyan", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[95], A1[95], c= "cyan", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[96], A1[96], c= "cyan", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[97], A1[97], c= "cyan", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[98], A1[98], c= "magenta", alpha=.4, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[99], A1[99], c= "magenta", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[100], A1[100], c= "magenta", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[101], A1[101], c= "magenta", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax1.scatter(B1[102], A1[102], c= "orange", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L2') #
ax1.scatter(B1[103], A1[103], c= "orange", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L3')
ax1.scatter(B1[104], A1[104], c= "orange", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L4')
ax1.scatter(B1[105], A1[105], c= "orange", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L5')
ax1.scatter(B1[106], A1[106], c= "green", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L2')
ax1.scatter(B1[107], A1[107], c= "green", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L3')
ax1.scatter(B1[108], A1[108], c= "green", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L4')
ax1.scatter(B1[109], A1[109], c= "green", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L5')
ax1.scatter(B1[110], A1[110], c= "brown", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L2')
ax1.scatter(B1[111], A1[111], c= "brown", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L3')
ax1.scatter(B1[112], A1[112], c= "brown", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L4')
ax1.scatter(B1[113], A1[113], c= "brown", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L5')
ax1.scatter(B1[114], A1[114], c= "cyan", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L2')
ax1.scatter(B1[115], A1[115], c= "cyan", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L3')
ax1.scatter(B1[116], A1[116], c= "cyan", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L4')
ax1.scatter(B1[117], A1[117], c= "cyan", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L5')
ax1.scatter(B1[118], A1[118], c= "magenta", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L2')
ax1.scatter(B1[119], A1[119], c= "magenta", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L3')
ax1.scatter(B1[120], A1[120], c= "magenta", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L4')
ax1.scatter(B1[121], A1[121], c= "magenta", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L5')
ax1.scatter(B1[2], A1[2], facecolors='none', edgecolors='orange', alpha=0.8, marker='s', s=7) # c= "orange", alpha=0.8, marker='s', s=7 )
ax1.scatter(B1[3], A1[3], facecolors='none', edgecolors='orange', alpha=0.8, marker='D', s=7)
ax1.scatter(B1[4], A1[4], facecolors='none', edgecolors='orange', alpha=0.8, marker='^', s=17)
ax1.scatter(B1[5], A1[5], facecolors='none', edgecolors='orange', alpha=0.8, marker='*', s=17)
ax1.scatter(B1[6], A1[6], facecolors='none', edgecolors='green', alpha=0.8, marker='s', s=7) # c= "green", alpha=0.8, marker='s', s=7)
ax1.scatter(B1[7], A1[7], facecolors='none', edgecolors='green', alpha=0.8, marker='D', s=7)
ax1.scatter(B1[8], A1[8], facecolors='none', edgecolors='green', alpha=0.8, marker='^', s=17)
ax1.scatter(B1[9], A1[9], facecolors='none', edgecolors='green', alpha=0.8, marker='*', s=17)
ax1.scatter(B1[10], A1[10], facecolors='none', edgecolors='brown', alpha=0.8, marker='s', s=7) # c= "brown", alpha=0.8, marker='s', s=7
ax1.scatter(B1[11], A1[11], facecolors='none', edgecolors='brown', alpha=0.8, marker='D', s=7)
ax1.scatter(B1[12], A1[12], facecolors='none', edgecolors='brown', alpha=0.8, marker='^', s=17)
ax1.scatter(B1[13], A1[13], facecolors='none', edgecolors='brown', alpha=0.8, marker='*', s=17)
ax1.scatter(B1[14], A1[14], facecolors='none', edgecolors='cyan', alpha=0.8, marker='s', s=7) #c= "cyan", alpha=0.8, marker='s', s=7
ax1.scatter(B1[15], A1[15], facecolors='none', edgecolors='cyan', alpha=0.8, marker='D', s=7)
ax1.scatter(B1[16], A1[16], facecolors='none', edgecolors='cyan', alpha=0.8, marker='^', s=17)
ax1.scatter(B1[17], A1[17], facecolors='none', edgecolors='cyan', alpha=0.8, marker='*', s=17)
ax1.scatter(B1[18], A1[18], facecolors='none', edgecolors='magenta', alpha=0.8, marker='s', s=7)#c= "magenta",alpha=0.8, marker='s', s=7)
ax1.scatter(B1[19], A1[19], facecolors='none', edgecolors='magenta', alpha=0.8, marker='D', s=7)
ax1.scatter(B1[20], A1[20], facecolors='none', edgecolors='magenta', alpha=0.8, marker='^', s=17)
ax1.scatter(B1[21], A1[21], facecolors='none', edgecolors='magenta', alpha=0.8, marker='*', s=17)
ax1.scatter(B1[22], A1[22], facecolors='none', edgecolors='orange', alpha=0.8, s=13, marker='s')#c= "orange",alpha=0.8, s=13, marker='s')
ax1.scatter(B1[23], A1[23], facecolors='none', edgecolors='orange', alpha=0.8, s=13, marker='D')
ax1.scatter(B1[24], A1[24], facecolors='none', edgecolors='orange', alpha=0.8, s=23, marker='^')
ax1.scatter(B1[25], A1[25], facecolors='none', edgecolors='orange', alpha=0.8, s=23, marker='*')
ax1.scatter(B1[26], A1[26], facecolors='none', edgecolors='green', alpha=0.8, s=13, marker='s') # c= "green",alpha=0.8, s=13, marker='s'
ax1.scatter(B1[27], A1[27], facecolors='none', edgecolors='green', alpha=0.8, s=13, marker='D')
ax1.scatter(B1[28], A1[28], facecolors='none', edgecolors='green', alpha=0.8, s=23, marker='^')
ax1.scatter(B1[29], A1[29], facecolors='none', edgecolors='green', alpha=0.8, s=23, marker='*')
ax1.scatter(B1[30], A1[30], facecolors='none', edgecolors='brown', alpha=0.8, s=13, marker='s') #c= "brown", alpha=0.8, s=13, marker='s')
ax1.scatter(B1[31], A1[31], facecolors='none', edgecolors='brown', alpha=0.8, s=13, marker='D')
ax1.scatter(B1[32], A1[32], facecolors='none', edgecolors='brown', alpha=0.8, s=23, marker='^')
ax1.scatter(B1[33], A1[33], facecolors='none', edgecolors='brown', alpha=0.8, s=23, marker='*')
ax1.scatter(B1[34], A1[34], facecolors='none', edgecolors='cyan', alpha=0.8, s=13, marker='s') #c= "cyan", alpha=0.8, s=13, marker='s')
ax1.scatter(B1[35], A1[35], facecolors='none', edgecolors='cyan', alpha=0.8, s=13, marker='D')
ax1.scatter(B1[36], A1[36], facecolors='none', edgecolors='cyan', alpha=0.8, s=23, marker='^')
ax1.scatter(B1[37], A1[37], facecolors='none', edgecolors='cyan', alpha=0.8, s=23, marker='*')
ax1.scatter(B1[38], A1[38], facecolors='none', edgecolors='magenta', alpha=0.8, s=13, marker='s') #c= "magenta",alpha=0.8, s=13,marker='s')
ax1.scatter(B1[39], A1[39], facecolors='none', edgecolors='magenta', alpha=0.8, s=13, marker='D')
ax1.scatter(B1[40], A1[40], facecolors='none', edgecolors='magenta', alpha=0.8, s=23, marker='^')
ax1.scatter(B1[41], A1[41], facecolors='none', edgecolors='magenta', alpha=0.8, s=23, marker='*')
ax1.scatter(B1[42], A1[42], facecolors='none', edgecolors='orange', alpha=0.8, s=30, marker='s') #c= "orange",alpha=0.8, s=30, marker='s',
ax1.scatter(B1[43], A1[43], facecolors='none', edgecolors='orange', alpha=0.8, s=30, marker='D')
ax1.scatter(B1[44], A1[44], facecolors='none', edgecolors='orange', alpha=0.8, s=40, marker='^')
ax1.scatter(B1[45], A1[45], facecolors='none', edgecolors='orange', alpha=0.8, s=40, marker='*')
ax1.scatter(B1[46], A1[46], facecolors='none', edgecolors='green', alpha=0.8, s=30, marker='s') #c= "green", alpha=0.8, s=30, marker='s',
ax1.scatter(B1[47], A1[47], facecolors='none', edgecolors='green', alpha=0.8, s=30, marker='D')
ax1.scatter(B1[48], A1[48], facecolors='none', edgecolors='green', alpha=0.8, s=40, marker='^')
ax1.scatter(B1[49], A1[49], facecolors='none', edgecolors='green', alpha=0.8, s=40, marker='*')
ax1.scatter(B1[50], A1[50], facecolors='none', edgecolors='brown', alpha=0.8, s=30, marker='s') #c= "brown", alpha=0.8, s=30, marker='s',
ax1.scatter(B1[51], A1[51], facecolors='none', edgecolors='brown', alpha=0.8, s=30, marker='D')
ax1.scatter(B1[52], A1[52], facecolors='none', edgecolors='brown', alpha=0.8, s=40, marker='^')
ax1.scatter(B1[53], A1[53], facecolors='none', edgecolors='brown', alpha=0.8, s=40, marker='*')
ax1.scatter(B1[54], A1[54], facecolors='none', edgecolors='cyan', alpha=0.8, s=30, marker='s') #c= "cyan", alpha=0.8, s=30, marker='s',
ax1.scatter(B1[55], A1[55], facecolors='none', edgecolors='cyan', alpha=0.8, s=30, marker='D')
ax1.scatter(B1[56], A1[56], facecolors='none', edgecolors='cyan', alpha=0.8, s=40, marker='^')
ax1.scatter(B1[57], A1[57], facecolors='none', edgecolors='cyan', alpha=0.8, s=40, marker='*')
ax1.scatter(B1[58], A1[58], facecolors='none', edgecolors='magenta', alpha=0.8, s=30, marker='s') #c= "magenta", alpha=0.8, s=30, marker='s'
ax1.scatter(B1[59], A1[59], facecolors='none', edgecolors='magenta', alpha=0.8, s=30, marker='D')
ax1.scatter(B1[60], A1[60], facecolors='none', edgecolors='magenta', alpha=0.8, s=40, marker='^')
ax1.scatter(B1[61], A1[61], facecolors='none', edgecolors='magenta', alpha=0.8, s=40, marker='*')
ax1.scatter(B1[0], A1[0], c='black', alpha=0.8, s=40, label='Obs. halo PNe')
(_, caps, _) = ax1.errorbar(B1[0], A1[0], xerr=0.0, yerr=0.0, marker='.', fmt='k.', elinewidth=0.8, markeredgewidth=0.8, markersize=8,)
for cap in caps:
cap.set_markeredgewidth(1)
ax1.scatter(B1[139], A1[139], c='black', alpha=0.8, s=40)#, label='Obs. halo PNe')
ax1.errorbar(B1[139], A1[139], yerr=0.2, uplims= True, marker='.', fmt='k.', elinewidth=.8, markeredgewidth=.8, markersize=8,)
#ax1.errorbar(B1[139], A1[139], xerr=0.22, xlolims= True, marker='.', fmt='k.', elinewidth=.8, markeredgewidth=.8, markersize=8,)
ax1.scatter(B1[140], A1[140], c='black', alpha=0.8, s=40)
(_, caps, _) = ax1.errorbar(B1[140], A1[140], xerr=0.0, yerr=0.0, marker='.', fmt='k.', elinewidth=0.8, markeredgewidth=0.8, markersize=8,)
for cap in caps:
cap.set_markeredgewidth(1)
ax1.scatter(B1[141], A1[141], c='black', alpha=0.8, s=40)
(_, caps, _) = ax1.errorbar(B1[141], A1[141], xerr=0.0, yerr=0.0, marker='.', fmt='k.', elinewidth=0.8, markeredgewidth=0.8, markersize=8,)
for cap in caps:
cap.set_markeredgewidth(1)
#ax1.errorbar(B1[141], A1[141], yerr=0.1, lolims= True)
# (_, caps, _) = ax1.errorbar(B1[0], A1[0], xerr=errx_con, yerr=erry_con, marker='o', fmt='ko', elinewidth=0.8, markeredgewidth=0.8, markersize=8,)
# for cap in caps:
# cap.set_markeredgewidth(1)
#ax1.scatter(B1[62],A1[62], c= "yellow", alpha=0.8, marker='o', label='Obs. disk PN')
ax1.scatter(B1[1], A1[1], c='purple', alpha=0.8, label='SDSS CVs')
ax1.scatter(B1[123], A1[123], c= "mediumaquamarine" , alpha=0.8, marker='s', label='SDSS QSOs (1.3<z<1.4)')
ax1.scatter(B1[124], A1[124], c= "royalblue", alpha=0.8, marker='D', label='SDSS QSOs (2.4<z<2.6)')
ax1.scatter(B1[142], A1[142], c= "salmon", alpha=0.8, s=38, marker='*', label='SDSS QSOs (3.0<z<3.2)')
ax1.scatter(B1[125], A1[125], c= "goldenrod", alpha=0.8, s=38, marker='^', label='SDSS QSOs (3.2<z<3.4)')
ax1.scatter(B1[126], A1[126], c= "white", alpha=0.3, s=38, marker='^', label='SDSS SFGs ')
ax1.scatter(B1[127], A1[127], c='red', alpha=0.8, s=38, marker='s')#, label='Obs. SySts ')
ax1.scatter(B1[135], A1[135], c= "red", alpha=0.8, s=38, marker='s', label='Obs. SySts ')
#ax1.scatter(B1[72], A1[72], c= "red", alpha=0.8, marker='D', label='Symbiotics in NGC 55')
ax1.scatter(B1[133], A1[133], c= "red", alpha=0.8, s=38, marker='D', label=' Obs. SySts in NGC 205')
ax1.scatter(B1[134], A1[134], c='red', alpha=0.8, s=38, marker='D') #label=' SySts in NGC 205')
ax1.scatter(B1[128], A1[128], c='red', alpha=0.8, s=38, marker='^')#, label='IPHAS SySts')
ax1.scatter(B1[136], A1[136], c= "red", alpha=0.8, s=38, marker='^', label='IPHAS SySts')
#ax1.scatter(B1[138], A1[138], c='red', alpha=0.8, s=38, marker='o', label='Obs. SySts in IC10 ')
#ax1.scatter(B1[137], A1[137], c= "red", alpha=0.8, s=38, marker='v', label=' Obs. SySts in NGC 185')
#ax1.scatter(B1[73], A1[73], c= "red", alpha=0.8, marker='o', label='C. Buil Symbiotics')
ax1.scatter(B1[129], A1[129], c= "gray", alpha=0.8, marker='D', label='Obs. HII region in NGC 55')
#ax1.scatter(B1[74], A1[74], c= "black", alpha=0.8, marker='.', label='SN Remanents')
# ax1.scatter(d_768_jplus[0], d_644_jplus[0], c= "b", alpha=0.8, marker='o', s=40, label='J-PLUS H4-1')
# (_, caps, _) = ax1.errorbar(d_768_jplus[0], d_644_jplus[0], xerr=err_Val_h4x, yerr=err_Val_h4y, marker='.', fmt='b.', elinewidth=.8, markeredgewidth=.8, markersize=8,)
# for cap in caps:
# cap.set_markeredgewidth(1)
# ax1.scatter(d_768_jplus[1], d_644_jplus[1], c= "g", alpha=0.8, marker='o', s=40, label='J-PLUS PNG135.9+55.9 ')
# (_, caps, _) = ax1.errorbar(d_768_jplus[1], d_644_jplus[1], xerr=err_Val_pn135x, yerr=err_Val_pn135y, marker='.', fmt='g.', elinewidth=.8, markeredgewidth=.8, markersize=8,) #fmt='yo',
# for cap in caps:
# cap.set_markeredgewidth(1)
# ax1.scatter(d_768_jplus[2], d_644_jplus[2], c= "y", alpha=0.8, marker='s', label='J-PLUS CI Cyg')
# (_, caps, _) = ax1.errorbar(d_768_jplus[2], d_644_jplus[2], xerr=errx[2], yerr=erry[2], marker='s', fmt='ys', elinewidth=.8, markeredgewidth=.8, markersize=8,)
# for cap in caps:
# cap.set_markeredgewidth(1)
# ax1.scatter(d_768_jplus[3], d_644_jplus[3], c= "m", alpha=0.8, marker='s', label='J-PLUS TX CVn')
# (_, caps, _) = ax1.errorbar(d_768_jplus[3], d_644_jplus[3], xerr=errx[3], yerr=erry[3], marker='s', fmt='ms', elinewidth=.8, markeredgewidth=.8, markersize=8,) #fmt='yo',
# for cap in caps:
# cap.set_markeredgewidth(1)
for label_, x, y in zip(label, B1[0], A1[0]):
ax1.annotate(label_, (x, y), alpha=10, size=8,
xytext=(40.0, -8.0), textcoords='offset points', ha='right', va='bottom',)
ax1.annotate("H4-1", (np.array(B1[139]), np.array(A1[139])), alpha=10, size=8,
xytext=(-5, 2), textcoords='offset points', ha='right', va='bottom',)
ax1.annotate("PNG 135.9+55.9", (np.array(B1[140]), np.array(A1[140])), alpha=10, size=8,
xytext=(60, -10), textcoords='offset points', ha='right', va='bottom',)
ax1.annotate("DdDm-1", (np.array(B1[141]), np.array(A1[141])), alpha=10, size=8,
xytext=(30, -10), textcoords='offset points', ha='right', va='bottom',)
#Obersevado porJPLUS
#for label_, x, y in zip(label1, d_768_jplus[0], d_644_jplus[0]):
# ax1.annotate("H4-1", (d_768_jplus[0], d_644_jplus[0]), alpha=8, size=8,
# xytext=(-5, 3), textcoords='offset points', ha='right', va='bottom', color='blue',)
# ax1.annotate("PNG 135.9+55.9", (d_768_jplus[1], d_644_jplus[1]), alpha=8, size=8,
# xytext=(68, -10), textcoords='offset points', ha='right', va='bottom', color='green',)
# ax1.annotate("CI Cyg", (d_768_jplus[2], d_644_jplus[2]), alpha=20, size=8,
# xytext=(-5, 3), textcoords='offset points', ha='right', va='bottom', color='yellow',)
# ax1.annotate("TX CVn", (d_768_jplus[3], d_644_jplus[3]), alpha=20, size=8,
# xytext=(18, -13), textcoords='offset points', ha='right', va='bottom', color='m',)
#for label_, x, y in zip(can_alh, d_768_cAlh, d_644_cAlh):
#ax1.annotate(label_, (x, y), alpha=0.9, size=8,
#xytext=(3,-10), textcoords='offset points', ha='left', va='bottom',)
# plt.annotate(
# '', xy=(B1[2][0], A1[2][0]+0.35), xycoords='data', #
# xytext=(B1[42][0]+0.4, A1[42][0]+0.4), textcoords='data', fontsize = 7,# vector extinction
# arrowprops=dict(edgecolor='black',arrowstyle= '<-')) #
# plt.annotate(
# '', xy=(B1[2][0]+0.35, A1[2][0]+0.35), xycoords='data',
# xytext=(5, 0), textcoords='offset points', fontsize='x-small')
#for Z, x, y in zip(z, d_768_Qz, d_644_Qz):
#ax1.annotate("{:.3f}".format(Z), (x, y), fontsize='x-small',
#xytext=(5,-5), textcoords='offset points', ha='left', bbox={"boxstyle": "round", "fc": "white", "ec": "none", "alpha": 0.5}, alpha=0.7)
#ax1.set_title(" ".join([cmd_args.source]))
#ax1.grid(True)
#ax1.annotate('Higher z(3.288)', xy=(0.08749580383300781, 0.181182861328125), xytext=(-0.5, -0.58),
#arrowprops=dict(arrowstyle="->", connectionstyle="arc3"))
#ax1.annotate('Lower z(3.065)', xy=(0.3957328796386719, 0.1367034912109375), xytext=(0.5, -0.58),
#arrowprops=dict(arrowstyle="->", connectionstyle="arc3"))
ax1.minorticks_on()
#ax1.grid(which='minor')#, lw=0.3)
ax1.legend(scatterpoints=1, ncol=1, fontsize=7.0, loc='upper left', **lgd_kws)
ax1.grid()
#lgd = ax1.legend(loc='center right', bbox_to_anchor=(1.27, 0.5), fontsize=7.5, **lgd_kws)
#ax1.grid(which='minor', lw=0.5)
#sns.despine(bottom=True)
plt.tight_layout()
#plt.savefig('luis-JPLUS-Viironen.pdf')#, bbox_extra_artists=(lgd,), bbox_inches='tight')
plt.savefig('Fig5-JPAS-g-u-v1.pdf')
plt.clf()
######################################################################################################
# Fig1 Halpha and [O III]
#######################
label = []
label1 = ["H4-1", "PNG 135.9+55.9"]
n = 143
A1, B1 = [[] for _ in range(n)], [[] for _ in range(n)]
d_644_jplus, d_768_jplus = [], []
d_644_jplus1, d_768_jplus1 = [], []
for file_name in file_list:
with open(file_name) as f:
data = json.load(f)
# if data['id'].endswith("1-HPNe"):
# label.append("DdDm-1")
# elif data['id'].endswith("SLOAN-HPNe-"):
# label.append("H4-1")
# elif data['id'].endswith("1359559-HPNe"):
# label.append("PNG 135.9+55.9")
if data['id'].startswith("ngc"):
label.append("NGC 2242")
elif data['id'].startswith("mwc"):
label.append("MWC 574")
plot_mag("Jv0915_6600", "Jv0915_6098", "Jv0915_5001", "Jv0915_4600")
lgd_kws = {'frameon': True, 'fancybox': True, 'shadow': True}
#sns.set(style="dark")#, context="talk")
sns.set_style('ticks')
fig = plt.figure(figsize=(7, 6))
ax2 = fig.add_subplot(111)
ax2.set_xlim(xmin=-5.8,xmax=1.2)
ax2.set_ylim(ymin=-5.5,ymax=1.0)
#ax2.set_xlim(xmin=-2.5,xmax=2.0)
plt.tick_params(axis='x', labelsize=18)
plt.tick_params(axis='y', labelsize=18)
plt.xlabel(r'$J5001 - J4600$', fontsize= 18)
plt.ylabel(r'$J6600 - J6098$', fontsize= 18)
ax2.scatter(B1[62], A1[62], c= "orange", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral )
ax2.scatter(B1[63], A1[63], c= "orange", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[64], A1[64], c= "orange", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[65], A1[65], c= "orange", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[66], A1[66], c= "green", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[67], A1[67], c= "green", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[68], A1[68], c= "green", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[69], A1[69], c= "green", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[70], A1[70], c= "brown", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[71], A1[71], c= "brown", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[72], A1[72], c= "brown", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[73], A1[73], c= "brown", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[74], A1[74], c= "cyan", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[75], A1[75], c= "cyan", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[76], A1[76], c= "cyan", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[77], A1[77], c= "cyan", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[78], A1[78], c= "magenta", alpha=.8, marker='s', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[79], A1[79], c= "magenta", alpha=.8, marker='D', s=7, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[80], A1[80], c= "magenta", alpha=.8, marker='^', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[81], A1[81], c= "magenta", alpha=.8, marker='*', s=17, lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[82], A1[82], c= "orange", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[83], A1[83], c= "orange", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[84], A1[84], c= "orange", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[85], A1[85], c= "orange", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[86], A1[86], c= "green", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[87], A1[87], c= "green", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[88], A1[88], c= "green", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[89], A1[89], c= "green", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[90], A1[90], c= "brown", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[91], A1[91], c= "brown", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[92], A1[92], c= "brown", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[93], A1[93], c= "brown", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[94], A1[94], c= "cyan", alpha=.8, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[95], A1[95], c= "cyan", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[96], A1[96], c= "cyan", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[97], A1[97], c= "cyan", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[98], A1[98], c= "magenta", alpha=.4, s=13, marker='s', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[99], A1[99], c= "magenta", alpha=.8, s=13, marker='D', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[100], A1[100], c= "magenta", alpha=.8, s=23, marker='^', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[101], A1[101], c= "magenta", alpha=.8, s=23, marker='*', lw=0, cmap=plt.cm.spectral)
ax2.scatter(B1[102], A1[102], c= "orange", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L2') #
ax2.scatter(B1[103], A1[103], c= "orange", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L3')
ax2.scatter(B1[104], A1[104], c= "orange", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L4')
ax2.scatter(B1[105], A1[105], c= "orange", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM DdDm1 L5')
ax2.scatter(B1[106], A1[106], c= "green", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L2')
ax2.scatter(B1[107], A1[107], c= "green", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L3')
ax2.scatter(B1[108], A1[108], c= "green", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L4')
ax2.scatter(B1[109], A1[109], c= "green", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM N2242 L5')
ax2.scatter(B1[110], A1[110], c= "brown", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L2')
ax2.scatter(B1[111], A1[111], c= "brown", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L3')
ax2.scatter(B1[112], A1[112], c= "brown", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L4')
ax2.scatter(B1[113], A1[113], c= "brown", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM K648 L5')
ax2.scatter(B1[114], A1[114], c= "cyan", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L2')
ax2.scatter(B1[115], A1[115], c= "cyan", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L3')
ax2.scatter(B1[116], A1[116], c= "cyan", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L4')
ax2.scatter(B1[117], A1[117], c= "cyan", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM BB1 L5')
ax2.scatter(B1[118], A1[118], c= "magenta", alpha=.8, s=30, marker='s', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L2')
ax2.scatter(B1[119], A1[119], c= "magenta", alpha=.8, s=30, marker='D', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L3')
ax2.scatter(B1[120], A1[120], c= "magenta", alpha=.8, s=40, marker='^', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L4')
ax2.scatter(B1[121], A1[121], c= "magenta", alpha=.8, s=40, marker='*', lw=0, cmap=plt.cm.spectral)#, label='BM Typ L5')
ax2.scatter(B1[2], A1[2], facecolors='none', edgecolors='orange', alpha=0.8, marker='s', s=7) # c= "orange", alpha=0.8, marker='s', s=7 )
ax2.scatter(B1[3], A1[3], facecolors='none', edgecolors='orange', alpha=0.8, marker='D', s=7)
ax2.scatter(B1[4], A1[4], facecolors='none', edgecolors='orange', alpha=0.8, marker='^', s=17)
ax2.scatter(B1[5], A1[5], facecolors='none', edgecolors='orange', alpha=0.8, marker='*', s=17)
ax2.scatter(B1[6], A1[6], facecolors='none', edgecolors='green', alpha=0.8, marker='s', s=7) # c= "green", alpha=0.8, marker='s', s=7)
ax2.scatter(B1[7], A1[7], facecolors='none', edgecolors='green', alpha=0.8, marker='D', s=7)
ax2.scatter(B1[8], A1[8], facecolors='none', edgecolors='green', alpha=0.8, marker='^', s=17)
ax2.scatter(B1[9], A1[9], facecolors='none', edgecolors='green', alpha=0.8, marker='*', s=17)
ax2.scatter(B1[10], A1[10], facecolors='none', edgecolors='brown', alpha=0.8, marker='s', s=7) # c= "brown", alpha=0.8, marker='s', s=7
ax2.scatter(B1[11], A1[11], facecolors='none', edgecolors='brown', alpha=0.8, marker='D', s=7)
ax2.scatter(B1[12], A1[12], facecolors='none', edgecolors='brown', alpha=0.8, marker='^', s=17)
ax2.scatter(B1[13], A1[13], facecolors='none', edgecolors='brown', alpha=0.8, marker='*', s=17)
ax2.scatter(B1[14], A1[14], facecolors='none', edgecolors='cyan', alpha=0.8, marker='s', s=7) #c= "cyan", alpha=0.8, marker='s', s=7
ax2.scatter(B1[15], A1[15], facecolors='none', edgecolors='cyan', alpha=0.8, marker='D', s=7)
ax2.scatter(B1[16], A1[16], facecolors='none', edgecolors='cyan', alpha=0.8, marker='^', s=17)
ax2.scatter(B1[17], A1[17], facecolors='none', edgecolors='cyan', alpha=0.8, marker='*', s=17)
ax2.scatter(B1[18], A1[18], facecolors='none', edgecolors='magenta', alpha=0.8, marker='s', s=7)#c= "magenta",alpha=0.8, marker='s', s=7)
ax2.scatter(B1[19], A1[19], facecolors='none', edgecolors='magenta', alpha=0.8, marker='D', s=7)
ax2.scatter(B1[20], A1[20], facecolors='none', edgecolors='magenta', alpha=0.8, marker='^', s=17)
ax2.scatter(B1[21], A1[21], facecolors='none', edgecolors='magenta', alpha=0.8, marker='*', s=17)
ax2.scatter(B1[22], A1[22], facecolors='none', edgecolors='orange', alpha=0.8, s=13, marker='s')#c= "orange",alpha=0.8, s=13, marker='s')
ax2.scatter(B1[23], A1[23], facecolors='none', edgecolors='orange', alpha=0.8, s=13, marker='D')
ax2.scatter(B1[24], A1[24], facecolors='none', edgecolors='orange', alpha=0.8, s=23, marker='^')
ax2.scatter(B1[25], A1[25], facecolors='none', edgecolors='orange', alpha=0.8, s=23, marker='*')