Chaos/BGO
12
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity

Compare commits

base: Chaos:9bc50667a37fb7193ff2f33db07f792673c08496
Branches Tags
Chaos:main
Chaos:temp_fix
..
compare: Chaos:aec9fdb142ff55f7536c1c0759d79667a918d731
Branches Tags
Chaos:main
Chaos:temp_fix

2 commits
9bc50667a3 ... aec9fdb142

Author SHA1 Message Date
ava
aec9fdb142 CHANGE BigBGO.py 2023-10-24 12:28:49 +02:00
ava
8c7e33b1cf Change BigBGO.py 2023-10-18 09:09:55 +02:00
1 changed files with 33 additions and 25 deletions
Download patch file Download diff file Expand all files Collapse all files

24
BigBGO.py
View file

@ -90,9 +90,8 @@ def main():
l = line.split() l = line.split()
if not 'EOF' in l: if not 'EOF' in l:
if l[0] == 'H': if l[0] == 'H':
T = adc_to_T(l[8]) T = adc_to_T(int(l[8]))
if l[0] == 'EI': if l[0] == 'EI' and not T=='unknown':
if not T=='unknown':
if (not args.xray and not args.muon) or(args.xray and isXray(l)) or (args.muon and ismuon(l)): if (not args.xray and not args.muon) or(args.xray and isXray(l)) or (args.muon and ismuon(l)):
if not args.trigger or (args.trigger and CmV(l,ch[3]) >= thr[3]): if not args.trigger or (args.trigger and CmV(l,ch[3]) >= thr[3]):
for i in range(0,len(ch)): for i in range(0,len(ch)):
@ -118,15 +117,14 @@ def main():
if not 'EOF' in l: if not 'EOF' in l:
if l[0] == 'H': if l[0] == 'H':
T = adc_to_T(int(l[8])) T = adc_to_T(int(l[8]))
if l[0] == 'EI': if l[0] == 'EI' and not T=='unknown':
if not T=='unknown':
if not args.trigger or (args.trigger and CmV(l,ch[3]) >= thr[3]): if not args.trigger or (args.trigger and CmV(l,ch[3]) >= thr[3]):
if not args.muon or (args.muon and ismuon(l)): if not args.muon or (args.muon and ismuon(l)):
if not args.two: if not args.two:
data[0].append(CmV(l,ch[0])*u[0]+CmV(l,ch[1])*u[1]+CmV(l,ch[2])*u[2]) data[0].append(CmV(l,ch[0])*u[0]+CmV(l,ch[1])*u[1]+CmV(l,ch[2])*u[2])
else: else:
#data[0].append(CmV(l,ch[1])*u[1]+CmV(l,ch[2])*u[2]) data[0].append(CmV(l,ch[1])*u[1]+CmV(l,ch[2])*u[2])
data[0].append(CmV(l,ch[0])*u[0]+CmV(l,ch[1])*u[1]) #data[0].append(CmV(l,ch[0])*u[1]+CmV(l,ch[2])*u[1])
data[1].append(CmV(l,ch[5])*u[5]+CmV(l,ch[6])*u[6]) data[1].append(CmV(l,ch[5])*u[5]+CmV(l,ch[6])*u[6])
if args.zero: if args.zero:
@ -197,9 +195,19 @@ def ismuon(l):
if not args.two and CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6] and \ if not args.two and CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6] and \
CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1] and CmV(l,A3) >= thr[2]: return True CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1] and CmV(l,A3) >= thr[2]: return True
if args.two and CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6] and \ if args.two and CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6] and \
CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1]: return True (CmV(l,A2) >= thr[1] and CmV(l,A3) >= thr[2]): return True
#CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1]: return True
else: return False else: return False
# def ismuon(l): #single BGO
# if not args.two and ( (CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6]) or \
# (CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1] and CmV(l,A3) >= thr[2])): return True
# if args.two and ( (CmV(l,B1) >= thr[5] and CmV(l,B2) >= thr[6]) or \
# #(CmV(l,A1) >= thr[0] and CmV(l,A2) >= thr[1])): return True
# (CmV(l,A2) >= thr[1] and CmV(l,A3) >= thr[2])): return True
# else: return False
def polynom(x, a,b,c,d): def polynom(x, a,b,c,d):
return np.multiply(np.multiply(np.multiply(a,x),x),x)+np.multiply(np.multiply(b,x),x)+np.multiply(c,x)+d return np.multiply(np.multiply(np.multiply(a,x),x),x)+np.multiply(np.multiply(b,x),x)+np.multiply(c,x)+d