Print

Print


Author: [log in to unmask]
Date: Wed Dec  2 19:00:08 2015
New Revision: 4004

Log:
one more cut

Modified:
    java/trunk/recon/src/main/java/org/hps/recon/filtering/V0CandidateFilter.java

Modified: java/trunk/recon/src/main/java/org/hps/recon/filtering/V0CandidateFilter.java
 =============================================================================
--- java/trunk/recon/src/main/java/org/hps/recon/filtering/V0CandidateFilter.java	(original)
+++ java/trunk/recon/src/main/java/org/hps/recon/filtering/V0CandidateFilter.java	Wed Dec  2 19:00:08 2015
@@ -32,6 +32,8 @@
     private double v0Chi2Cut = 10.0;
     private double trackChi2Cut = 20.0;
     private double trackDtCut = 5.0;
+    private double trackPMax = 0.9;
+    private double v0PMax = 1.4;
 
     private boolean _tight = false;
     private boolean _keepEpicsDataEvents = false;
@@ -69,6 +71,12 @@
                 continue;
             }
             if (electron.getTracks().get(0).getChi2() > trackChi2Cut || positron.getTracks().get(0).getChi2() > trackChi2Cut) {
+                continue;
+            }
+            if (electron.getMomentum().magnitude() > trackPMax || positron.getMomentum().magnitude() > trackPMax) {
+                continue;
+            }
+            if (v0.getMomentum().magnitude() > v0PMax) {
                 continue;
             }
             double eleTime = TrackData.getTrackTime(TrackData.getTrackData(event, electron.getTracks().get(0)));
@@ -130,6 +138,20 @@
     }
 
     /**
+     * Maximum track momentum for a V0 to be counted. A V0 is rejected if either
+     * of the final state tracks has momentum exceeding this cut.
+     *
+     * @param trackPMax units of GeV, default of 0.9
+     */
+    public void setTrackPMax(double trackPMax) {
+        this.trackPMax = trackPMax;
+    }
+
+    public void setV0PMax(double v0PMax) {
+        this.v0PMax = v0PMax;
+    }
+
+    /**
      * Maximum difference in Calorimeter Cluster Seed Hit times [ns]
      *
      * @param d