README


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
generated by cgit  (git 2.34.1) at 2025-10-13 23:35:17 +0000
 

n>) & (ev.gtr == 0) & (ev.trg_type == 0))]
        orphans_per_run[run] = len(orphans[orphans.nhit >= args.nhit_thresh])

    if args.output:
        np.savetxt(args.output,sorted([run for run, orphans in orphans_per_run.iteritems() if orphans < 100]),fmt='%i')

    fig = plt.figure()
    plt.hist(orphans_per_run.values(),bins=np.linspace(0,1000,101),histtype='step')
    plt.axvline(x=100,ls='--',color='k')
    plt.xlabel("Number of orphans with Nhit > %i" % args.nhit_thresh)
    plt.gca().set_yscale("log")
    despine(fig,trim=True)
    plt.tight_layout()
    if args.save:
        plt.savefig("orphans.pdf")
        plt.savefig("orphans.eps")
    else:
        plt.title(r"Number of High Nhit Orphans Per Run")
        plt.show()