lcsim/src/org/lcsim/mc/fast/cluster/ronan

diff -u -r1.7 -r1.8
--- MCFastRonan.java	9 Aug 2005 18:34:45 -0000	1.7
+++ MCFastRonan.java	24 Sep 2005 01:21:16 -0000	1.8
@@ -1,17 +1,15 @@

 package org.lcsim.mc.fast.cluster.ronan;
 

-import hep.physics.particle.Particle;
-import hep.physics.particle.properties.ParticleType;

 import java.util.ArrayList;

-import java.util.Iterator;

 import java.util.List;
 import java.util.Random;

+

 import org.lcsim.conditions.ConditionsEvent;
 import org.lcsim.conditions.ConditionsListener;
 import org.lcsim.conditions.ConditionsSet;

-import org.lcsim.event.MCParticle;

 import org.lcsim.event.Cluster;
 import org.lcsim.event.EventHeader;

+import org.lcsim.event.MCParticle;

 import org.lcsim.util.Driver;
 
 

@@ -20,6 +18,7 @@

 /**
  * Fast Monte Carlo cluster simulator
  * @author M.Ronan  Oct 2000 - Added "refined" cluster simulation

+ * @version

  */
 public class MCFastRonan extends Driver implements ConditionsListener
 {

@@ -27,6 +26,12 @@

     private final static int NuMuID = 14;
     private final static int NuTauID = 16;
     private final static int PhotonID = 22;

+    private final static int Neutralino1 = 1000022;
+    private final static int Neutralino2 = 1000023;
+    private final static int Neutralino3 = 1000025;
+    private final static int Neutralino4 = 1000035;
+    
+

     private ClusterResolutionTables clusterParm;
     
     protected void process(EventHeader event)

@@ -52,10 +57,13 @@

             }
             
             int PDGID = p.getPDGID();

+            int absPDGID = Math.abs(PDGID);

             double charge = p.getCharge();
             

+            

             // filter neutrinos

-            boolean neutrino = (Math.abs(PDGID) == NuEID) || (Math.abs(PDGID) == NuMuID) || (Math.abs(PDGID) == NuTauID);

+            boolean neutrino = absPDGID == NuEID || absPDGID == NuMuID || absPDGID == NuTauID
+                    || absPDGID == Neutralino1 || absPDGID == Neutralino2 || absPDGID == Neutralino3 || absPDGID == Neutralino4;

             if (neutrino)
             {
                 continue;