Author: [log in to unmask]
Date: Mon Sep 21 17:19:52 2015
New Revision: 3661

Log:
update user dir

Added:
    java/trunk/users/src/main/java/org/hps/users/meeg/V0Filter.java
      - copied, changed from r3631, java/trunk/recon/src/main/java/org/hps/recon/filtering/MollerCandidateFilter.java

Copied: java/trunk/users/src/main/java/org/hps/users/meeg/V0Filter.java (from r3631, java/trunk/recon/src/main/java/org/hps/recon/filtering/MollerCandidateFilter.java)
 =============================================================================
--- java/trunk/recon/src/main/java/org/hps/recon/filtering/MollerCandidateFilter.java	(original)
+++ java/trunk/users/src/main/java/org/hps/users/meeg/V0Filter.java	Mon Sep 21 17:19:52 2015
@@ -1,149 +1,38 @@
-package org.hps.recon.filtering;
+package org.hps.users.meeg;
 
-import static java.lang.Math.abs;
 import java.util.List;
+import org.hps.recon.filtering.EventReconFilter;
 import org.lcsim.event.EventHeader;
 import org.lcsim.event.ReconstructedParticle;
 
-/**
- * Class to strip off Moller candidates. Currently defined as: e- e- events with
- * tracks matched to clusters. Neither electron can be a full-energy candidate
- * (momentum less than _fullEnergyCut [0.85GeV]) but the momentum sum must be
- * consistent with the beam energy (greater than _mollerMomentumSumMin and less
- * than _mollerMomentumSumMax). The Ecal cluster times must be within _timingCut
- * [2.5ns] of each other.
- *
- * @author Norman A Graf
- *
- * @version $Id:
- */
-public class MollerCandidateFilter extends EventReconFilter {
+public class V0Filter extends EventReconFilter {
 
-    private String _mollerCandidateCollectionName = "TargetConstrainedMollerCandidates";
-    private double _mollerMomentumSumMin = 0.85;
-    private double _mollerMomentumSumMax = 1.3;
-    private double _fullEnergyCut = 0.85;
-    private double _clusterTimingCut = 2.5;
-
-    private boolean _tight = false;
+    private String v0CollectionName = "UnconstrainedV0Candidates";
+    private double massMin = 0.025;
+    private double massMax = 0.030;
+    private double minPTot = 0.7 * 1.05;
 
     @Override
     protected void process(EventHeader event) {
         incrementEventProcessed();
-        if (!event.hasCollection(ReconstructedParticle.class, _mollerCandidateCollectionName)) {
+        if (!event.hasCollection(ReconstructedParticle.class, v0CollectionName)) {
             skipEvent();
         }
-        List<ReconstructedParticle> mollerCandidates = event.get(ReconstructedParticle.class, _mollerCandidateCollectionName);
-        if (mollerCandidates.size() == 0) {
+        List<ReconstructedParticle> v0List = event.get(ReconstructedParticle.class, v0CollectionName);
+        if (v0List.isEmpty()) {
             skipEvent();
         }
-
-        // tight requires ONLY ONE real vertex fit 
-        if (_tight) {
-            if (mollerCandidates.size() != 2) {
-                skipEvent();
+        boolean hasGoodVertex = false;
+        for (ReconstructedParticle rp : v0List) {
+            double mass = rp.getMass();
+            if (mass > massMin && mass < massMax && rp.getMomentum().magnitude() > minPTot) {
+                hasGoodVertex = true;
+                break;
             }
         }
-
-        for (ReconstructedParticle rp : mollerCandidates) {
-
-            ReconstructedParticle e1 = null;
-            ReconstructedParticle e2 = null;
-
-            List<ReconstructedParticle> electrons = rp.getParticles();
-            if (electrons.size() != 2) {
-                skipEvent();
-            }
-            // require both electrons to be associated with an ECal cluster
-            e1 = electrons.get(0);
-            if (e1.getClusters().size() == 0) {
-                skipEvent();
-            }
-            e2 = electrons.get(1);
-            if (e2.getClusters().size() == 0) {
-                skipEvent();
-            }
-            // remove full energy electrons
-            double p1 = e1.getMomentum().magnitude();
-            if (p1 > _fullEnergyCut) {
-                skipEvent();
-            }
-            double p2 = e2.getMomentum().magnitude();
-            if (p2 > _fullEnergyCut) {
-                skipEvent();
-            }
-
-            // require momentum sum to be approximately the beam energy
-            double pSum = p1 + p2;
-            if (pSum < _mollerMomentumSumMin || pSum > _mollerMomentumSumMax) {
-                skipEvent();
-            }
-
-            // calorimeter cluster timing cut
-            // first CalorimeterHit in the list is the seed crystal
-            double t1 = e1.getClusters().get(0).getCalorimeterHits().get(0).getTime();
-            double t2 = e2.getClusters().get(0).getCalorimeterHits().get(0).getTime();
-
-            if (abs(t1 - t2) > _clusterTimingCut) {
-                skipEvent();
-            }
+        if (!hasGoodVertex) {
+            skipEvent();
         }
         incrementEventPassed();
     }
-
-    /**
-     * Maximum difference in Calorimeter Cluster Seed Hit times [ns]
-     *
-     * @param d
-     */
-    public void setClusterTimingCut(double d) {
-        _clusterTimingCut = d;
-    }
-
-    /**
-     * Name of Moller Candidate ReconstructedParticle Collection Name
-     *
-     * @param s
-     */
-    public void setMollerCandidateCollectionName(String s) {
-        _mollerCandidateCollectionName = s;
-    }
-
-    /**
-     * Minimum value for the sum of the two electron momenta [GeV]
-     *
-     * @param d
-     */
-    public void setMollerMomentumSumMin(double d) {
-        _mollerMomentumSumMin = d;
-    }
-
-    /**
-     * Maximum value for the sum of the two electron momenta [GeV]
-     *
-     * @param d
-     */
-    public void setMollerMomentumSumMax(double d) {
-        _mollerMomentumSumMax = d;
-    }
-
-    /**
-     * Maximum value for each of two electron momenta (removes full energy
-     * electrons) [GeV]
-     *
-     * @param d
-     */
-    public void setMollerMomentumMax(double d) {
-        _fullEnergyCut = d;
-    }
-
-    /**
-     * Setting a tight constraint requires one and only one candidate in the
-     * event
-     *
-     * @param b
-     */
-    public void setTightConstraint(boolean b) {
-        _tight = b;
-    }
 }