Author: [log in to unmask]
Date: Thu Oct 22 21:33:15 2015
New Revision: 3885

Log:
Added new cuts to pair process analysis driver to try and reduce background.

Modified:
    java/trunk/users/src/main/java/org/hps/users/kmccarty/TriggerProcessAnalysisDriver.java

Modified: java/trunk/users/src/main/java/org/hps/users/kmccarty/TriggerProcessAnalysisDriver.java
 =============================================================================
--- java/trunk/users/src/main/java/org/hps/users/kmccarty/TriggerProcessAnalysisDriver.java	(original)
+++ java/trunk/users/src/main/java/org/hps/users/kmccarty/TriggerProcessAnalysisDriver.java	Thu Oct 22 21:33:15 2015
@@ -21,14 +21,18 @@
 import org.lcsim.util.aida.AIDA;
 
 public class TriggerProcessAnalysisDriver extends Driver {
-	private boolean checkSVT = false;
 	private int eventsProcessed = 0;
 	private int møllersProcessed = 0;
+	private boolean checkSVT = false;
 	private int tridentsProcessed = 0;
 	private int gblMøllersProcessed = 0;
 	private int gblTridentsProcessed = 0;
 	private double timeCoincidence = 2.5;
+	private double elasticThreshold = 0.800;
+	private double møllerLowerRange = 0.900;
+	private double møllerUpperRange = 1.200;
 	private AIDA aida = AIDA.defaultInstance();
+	private boolean checkTriggerTimeWindow = false;
 	private String clusterCollectionName = "EcalClustersCorr";
 	private String particleCollectionName = "FinalStateParticles";
 	
@@ -281,6 +285,10 @@
 	
 	public void setCheckSVT(boolean state) {
 		checkSVT = state;
+	}
+	
+	public void setCheckTriggerTimeWindow(boolean state) {
+		checkTriggerTimeWindow = state;
 	}
 	
 	/**
@@ -448,10 +456,10 @@
 				continue tridentTrackLoop;
 			}
 			
-			// Require that the negative track have less than 900 MeV
-			// momentum to exclude elastic electrons.
-			if(pair[0].getCharge() < 0 && pair[0].getMomentum().magnitude() > 0.900
-					|| pair[1].getCharge() < 0 && pair[1].getMomentum().magnitude() > 0.900) {
+			// Require that the negative track have less than the
+			// elastic threshold momentum to exclude elastic electrons.
+			if(pair[0].getCharge() < 0 && pair[0].getMomentum().magnitude() > elasticThreshold
+					|| pair[1].getCharge() < 0 && pair[1].getMomentum().magnitude() > elasticThreshold) {
 				continue tridentTrackLoop;
 			}
 			
@@ -505,15 +513,32 @@
 			}
 			
 			// The clusters must within a limited time window.
+			/*
 			Cluster[] trackClusters = { pair[0].getClusters().get(0), pair[1].getClusters().get(0) };
 			if(TriggerModule.getValueTimeCoincidence(trackClusters) > timeCoincidence) {
 				continue tridentLoop;
 			}
+			*/
+			
+			// The clusters must be coincidental within an energy
+			// dependent coincidence window.
+			Cluster[] trackClusters = { pair[0].getClusters().get(0), pair[1].getClusters().get(0) };
+			if(!isCoincidental(trackClusters)) {
+				continue tridentLoop;
+			}
 			
 			// Require that the electron in the pair have an energy
-			// below 900 MeV to exclude elastic electrons.
-			if(electron.getMomentum().magnitude() >= 0.900) {
+			// below the elastic threshold to exclude elastic electrons.
+			if(electron.getMomentum().magnitude() >= elasticThreshold) {
 				continue tridentLoop;
+			}
+			
+			// Require that all clusters occur within the trigger time
+			// window to exclude accidentals.
+			if(checkTriggerTimeWindow) {
+				if(!inTriggerWindow(trackClusters[0]) || !inTriggerWindow(trackClusters[1])) {
+					continue tridentLoop;
+				}
 			}
 			
 			// If all the above conditions are met, the pair is to be
@@ -553,15 +578,15 @@
 			}
 			
 			// Require that the electrons in the pair have energies
-			// below 900 MeV to exclude elastic electrons.
-			if(pair[0].getMomentum().magnitude() > 0.900 || pair[1].getMomentum().magnitude() > 0.900) {
+			// below the elastic threshold to exclude said electrons.
+			if(pair[0].getMomentum().magnitude() > elasticThreshold || pair[1].getMomentum().magnitude() > elasticThreshold) {
 				continue møllerLoop;
 			}
 			
 			// Require that the energy of the pair be within a range
 			// that is sufficiently "Møller-like."
 			double momentumSum = VecOp.add(pair[0].getMomentum(), pair[1].getMomentum()).magnitude();
-			if(momentumSum < 0.800 || momentumSum > 1.500) {
+			if(momentumSum < møllerLowerRange || momentumSum > møllerUpperRange) {
 				continue møllerLoop;
 			}
 			
@@ -605,22 +630,39 @@
 			}
 			
 			// The clusters must within a limited time window.
+			/*
 			Cluster[] trackClusters = { pair[0].getClusters().get(0), pair[1].getClusters().get(0) };
 			if(TriggerModule.getValueTimeCoincidence(trackClusters) > timeCoincidence) {
 				continue møllerLoop;
 			}
+			*/
+			
+			// The clusters must be coincidental within an energy
+			// dependent coincidence window.
+			Cluster[] trackClusters = { pair[0].getClusters().get(0), pair[1].getClusters().get(0) };
+			if(!isCoincidental(trackClusters)) {
+				continue møllerLoop;
+			}
 			
 			// Require that the electrons in the pair have energies
-			// below 900 MeV to exclude elastic electrons.
-			if(pair[0].getMomentum().magnitude() > 0.900 || pair[1].getMomentum().magnitude() > 0.900) {
+			// below the elastic threshold to exclude said electrons.
+			if(pair[0].getMomentum().magnitude() > elasticThreshold || pair[1].getMomentum().magnitude() > elasticThreshold) {
 				continue møllerLoop;
 			}
 			
 			// Require that the energy of the pair be within a range
 			// that is sufficiently "Møller-like."
 			double momentumSum = VecOp.add(pair[0].getMomentum(), pair[1].getMomentum()).magnitude();
-			if(momentumSum < 0.800 || momentumSum > 1.500) {
-				continue møllerLoop;
+			if(momentumSum < møllerLowerRange || momentumSum > møllerUpperRange) {
+				continue møllerLoop;
+			}
+			
+			// Require that all clusters occur within the trigger time
+			// window to exclude accidentals.
+			if(checkTriggerTimeWindow) {
+				if(!inTriggerWindow(trackClusters[0]) || !inTriggerWindow(trackClusters[1])) {
+					continue møllerLoop;
+				}
 			}
 			
 			// If all the above conditions are met, the pair is to be
@@ -704,5 +746,43 @@
 	private static final double getCalculatedCoplanarity(Track[] pair) {
 		return getCalculatedCoplanarity(TrackUtils.getTrackPositionAtEcal(pair[0]).v(), TrackUtils.getTrackPositionAtEcal(pair[1]).v());
 	}
-}
-// [log in to unmask]
+	
+	private static final boolean inTriggerWindow(Cluster cluster) {
+		// Get the cluster time.
+		double clusterTime = TriggerModule.getClusterTime(cluster);
+		
+		// Check that it is within the allowed bounds.
+		return (35 <= clusterTime && clusterTime <= 50);
+	}
+	
+	private static final boolean isCoincidental(Cluster[] pair) {
+		// Get the energy sum and the time coincidence.
+		double energySum = pair[0].getEnergy() + pair[1].getEnergy();
+		double timeCoincidence = TriggerModule.getValueTimeCoincidence(pair);
+		
+		// Get the upper and lower bounds of the allowed range.
+		double mean = getTimeDependenceMean(energySum);
+		double threeSigma = 3.0 * getTimeDependenceSigma(energySum);
+		double lowerBound = mean - threeSigma;
+		double upperBound = mean + threeSigma;
+		
+		// Perform the time coincidence check.
+		return (lowerBound <= timeCoincidence && timeCoincidence <= upperBound);
+	}
+	
+	private static final double getTimeDependenceMean(double energySum) {
+		// Define the fit parameters.
+		double[] param = { 0.289337, -2.81998, 9.03475, -12.93, 8.71476, -2.26969 };
+		
+		// Calculate the mean.
+		return param[0] + energySum * (param[1] + energySum * (param[2] + energySum * (param[3] + energySum * (param[4] + energySum * (param[5])))));
+	}
+	
+	private static final double getTimeDependenceSigma(double energySum) {
+		// Define the fit parameters.
+		double[] param = { 4.3987, -24.2371, 68.9567, -98.2586, 67.562, -17.8987 };
+		
+		// Calculate the standard deviation.
+		return param[0] + energySum * (param[1] + energySum * (param[2] + energySum * (param[3] + energySum * (param[4] + energySum * (param[5])))));
+	}
+}