solo/hetept/l3/ept.l3

@@ -42,17 +42,9 @@ 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)
 
-	ifN	CMP MaxE < Z+ept_hist_offset
+		CMP MaxE < Z+ept_hist_offset
 
 # test if the central channel has seen a peak, T>1 
 
@@ -74,21 +66,10 @@ 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
 
-	ifN	CMP  ANTI > Z+cut_ept_ac_thres
+		CMP  ANTI > Z+cut_ept_ac_thres
 	ifC	CMP  MaxE < ANTI+cut_ept_ac_ratio
 
 # C=1 S=0 for invalid event PHA

solo/hetept/l3/het_calib.l3

@@ -161,7 +161,6 @@ 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

solo/hetept/l3/het_stopping.l3

@@ -75,7 +75,13 @@ 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
 
@@ -238,7 +244,6 @@ 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