Compare commits

..

4 commits

3 changed files with 23 additions and 8 deletions

View file

@ -42,9 +42,17 @@ MinE = ifN ADDI C$2E
BITS C 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) # 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 # 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 CMP A$1E < A$2E
ANTI = ifC ADDI 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 # Test if the anticoicidence is above noise
# and exceed the crosstalk limit vs the hit central pixel # 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 ifC CMP MaxE < ANTI+cut_ept_ac_ratio
# C=1 S=0 for invalid event PHA # C=1 S=0 for invalid event PHA

View file

@ -161,6 +161,7 @@ het_gcr_C_min = $LOG(het_gcr_e0*1e6/eV_BGO)
cut_het_b3_max = $LOG(40e3/eV_Si) 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_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_AB_leak_thres = $LOG(800e3/eV_Si)
cut_C_ABC_electron_maxB = $LOG(500e3/eV_Si)
cut_C_thres = $LOG(1000e3/eV_BGO) cut_C_thres = $LOG(1000e3/eV_BGO)
het_path_length_IO = 0.991 # 2016-12-02 Robert E. geant4 het_path_length_IO = 0.991 # 2016-12-02 Robert E. geant4
het_path_length_OO = 0.983 # 2016-12-02 Robert E. geant4 het_path_length_OO = 0.983 # 2016-12-02 Robert E. geant4

View file

@ -75,13 +75,7 @@ het_stop_ABnC=.
.print het_stop_ABnC address $. .print het_stop_ABnC address $.
# C must be less than cut_C_thres # 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 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 ifC CMP CS < Z+cut_C_thres
ifN GOTO reject_stop ifN GOTO reject_stop
@ -244,6 +238,7 @@ Y = ADDI Y - cut_HET_ABC_e_min
# electrons: histogram 128 x 1 using total energy # electrons: histogram 128 x 1 using total energy
CMP Y < Z + cut_HET_ABC_e_max - cut_HET_ABC_e_min 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 TOTE = ifC TRIM TOT, cut_HET_ABC_TOT_e_min, 127
H = ifC ADDI hist_base + hist_HET_ABC_e H = ifC ADDI hist_base + hist_HET_ABC_e