Compare commits
4 commits
e5e1bd6d91
...
0d05073229
| Author | SHA1 | Date | |
|---|---|---|---|
| 0d05073229 | |||
| 7a5199ed67 | |||
| d6f1aee1dd | |||
| ac084f6b67 |
3 changed files with 23 additions and 8 deletions
|
|
@ -42,9 +42,17 @@ MinE = ifN ADDI C$2E
|
|||
|
||||
BITS C
|
||||
|
||||
# Test if the other segment has significant negative signal.
|
||||
|
||||
SUB Z - C$1K
|
||||
mMinK = ifN SUB Z - C$2K
|
||||
mMinE = LOG mMinK
|
||||
CMP mMinE > Z+cut_ept_ac_thres
|
||||
ifC CMP MaxE < mMinE + cut_ept_ac_ratio
|
||||
|
||||
# Check if we are above the minimum energy limit, otherwise abort (S=0 for PHA buffer selection)
|
||||
|
||||
CMP MaxE < Z+ept_hist_offset
|
||||
ifN CMP MaxE < Z+ept_hist_offset
|
||||
|
||||
# test if the central channel has seen a peak, T>1
|
||||
|
||||
|
|
@ -66,10 +74,21 @@ A$2E = LOG A$2K
|
|||
CMP A$1E < A$2E
|
||||
ANTI = ifC ADDI A$2E
|
||||
|
||||
# Test the anticoincidence channels for signifivant negative signal
|
||||
|
||||
mA$1K = SUB Z - A$1K
|
||||
mA$1E = LOG mA$1K
|
||||
mA$2K = SUB Z - A$2K
|
||||
mA$2E = LOG mA$2K
|
||||
CMP mA$1E < mA$2E
|
||||
mANTI = ifC ADDI mA$2E
|
||||
CMP mANTI > Z+cut_ept_ac_thres
|
||||
ifC CMP MaxE < mANTI + cut_ept_ac_ratio
|
||||
|
||||
# Test if the anticoicidence is above noise
|
||||
# and exceed the crosstalk limit vs the hit central pixel
|
||||
|
||||
CMP ANTI > Z+cut_ept_ac_thres
|
||||
ifN CMP ANTI > Z+cut_ept_ac_thres
|
||||
ifC CMP MaxE < ANTI+cut_ept_ac_ratio
|
||||
|
||||
# C=1 S=0 for invalid event PHA
|
||||
|
|
|
|||
|
|
@ -161,6 +161,7 @@ het_gcr_C_min = $LOG(het_gcr_e0*1e6/eV_BGO)
|
|||
cut_het_b3_max = $LOG(40e3/eV_Si)
|
||||
cut_C_AB_ratio = $LOGR(10.*BGO_Si) #TESTING! -> this is original $LOGR(10./BGO_Si) # C [125eV], AB [31.25eV]
|
||||
cut_AB_leak_thres = $LOG(800e3/eV_Si)
|
||||
cut_C_ABC_electron_maxB = $LOG(500e3/eV_Si)
|
||||
cut_C_thres = $LOG(1000e3/eV_BGO)
|
||||
het_path_length_IO = 0.991 # 2016-12-02 Robert E. geant4
|
||||
het_path_length_OO = 0.983 # 2016-12-02 Robert E. geant4
|
||||
|
|
|
|||
|
|
@ -75,13 +75,7 @@ het_stop_ABnC=.
|
|||
.print het_stop_ABnC address $.
|
||||
|
||||
# C must be less than cut_C_thres
|
||||
# and when AB has little energy, either C age must be invalid
|
||||
|
||||
CMP AB > Z+cut_AB_leak_thres
|
||||
C1T = ifN BRNG C1U{3:0}
|
||||
ifN CMP C1T <= Z+1
|
||||
C2T = ifN BRNG C2U{3:0}
|
||||
ifN CMP C2T <= Z+1
|
||||
ifC CMP CS < Z+cut_C_thres
|
||||
ifN GOTO reject_stop
|
||||
|
||||
|
|
@ -244,6 +238,7 @@ Y = ADDI Y - cut_HET_ABC_e_min
|
|||
|
||||
# electrons: histogram 128 x 1 using total energy
|
||||
CMP Y < Z + cut_HET_ABC_e_max - cut_HET_ABC_e_min
|
||||
ifC CMP B < Z + cut_C_ABC_electron_maxB
|
||||
|
||||
TOTE = ifC TRIM TOT, cut_HET_ABC_TOT_e_min, 127
|
||||
H = ifC ADDI hist_base + hist_HET_ABC_e
|
||||
|
|
|
|||
Loading…
Add table
Add a link
Reference in a new issue