Print

---

Author: [log in to unmask]
Date: Thu Jul 30 11:41:30 2015
New Revision: 3309

Log:
Simple Driver to strip events containing a single Full Energy Electron candidate.

Added:
    java/trunk/analysis/src/main/java/org/hps/analysis/examples/StripSingleFeeDriver.java

Added: java/trunk/analysis/src/main/java/org/hps/analysis/examples/StripSingleFeeDriver.java
 =============================================================================
--- java/trunk/analysis/src/main/java/org/hps/analysis/examples/StripSingleFeeDriver.java	(added)
+++ java/trunk/analysis/src/main/java/org/hps/analysis/examples/StripSingleFeeDriver.java	Thu Jul 30 11:41:30 2015
@@ -0,0 +1,67 @@
+package org.hps.analysis.examples;
+
+import java.util.List;
+import org.lcsim.event.EventHeader;
+import org.lcsim.event.ReconstructedParticle;
+import org.lcsim.util.Driver;
+
+/**
+ * Simple Driver to select events containing only one Full Energy Electron
+ * (FEE). This is defined as one ReconstructedParticle which has been identified
+ * as an electron and has an energy equal to or greater than the energyCut
+ * (default is 0.85) and a track which has greater than or equal the
+ * numberOfHitsOnTrack (default is 6).
+ *
+ * @author Norman A Graf
+ *
+ * @version $Id:
+ */
+public class StripSingleFeeDriver extends Driver
+{
+
+    private double _energyCut = 0.85;
+    private int _nHitsOnTrack = 6;
+
+    /**
+     *
+     * @param event
+     */
+    @Override
+    protected void process(EventHeader event)
+    {
+        boolean skipEvent = true;
+        List<ReconstructedParticle> rps = event.get(ReconstructedParticle.class, "FinalStateParticles");
+        if (rps.size() == 1) {
+            ReconstructedParticle rp = rps.get(0);
+            double energy = rp.getEnergy();
+            int pdgId = rp.getParticleIDUsed().getPDG();
+            if (pdgId == 11) {
+                if (energy >= _energyCut && rp.getTracks().get(0).getTrackerHits().size() >= _nHitsOnTrack) // electron
+                {
+                    skipEvent = false;
+                }
+            }
+        }
+    }
+
+    /**
+     * Electrons having energy below the cut will be rejected.
+     *
+     * @param cut
+     */
+    public void setEnergyCut(double cut)
+    {
+        _energyCut = cut;
+    }
+
+    /**
+     * Tracks having fewer than the number of hits will be rejected.
+     *
+     * @param cut
+     */
+    public void setNumberOfHitsOnTrack(int cut)
+    {
+        _nHitsOnTrack = cut;
+    }
+
+}