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

diff -u -r1.12 -r1.13
--- MCFastRonan.java	16 Feb 2007 01:34:34 -0000	1.12
+++ MCFastRonan.java	30 Aug 2007 19:40:22 -0000	1.13
@@ -9,9 +9,11 @@

 import org.lcsim.conditions.ConditionsListener;
 import org.lcsim.conditions.ConditionsSet;
 import org.lcsim.event.Cluster;

+import org.lcsim.event.LCRelation;

 import org.lcsim.event.EventHeader;
 import org.lcsim.event.MCParticle;
 import org.lcsim.util.Driver;

+import org.lcsim.digisim.MyLCRelation;

@@ -46,6 +48,7 @@

             clusterParm = new ClusterResolutionTables(conditions);
         }
         List<Cluster> cl = new ArrayList<Cluster>();

+	List<LCRelation> lcrelationList = new ArrayList<LCRelation>();

 
         boolean hist = getHistogramLevel() > 0;
         

@@ -99,6 +102,7 @@

                 }
                 
                 cl.add(new ReconEMCluster(clusterParm, rand, p, hist));

+		lcrelationList.add(new MyLCRelation(cl.get(cl.size()-1),(MCParticle) p));

                 
             }
             

@@ -122,6 +126,7 @@

                 }
                 
                 cl.add(new ReconHADCluster(clusterParm, rand, p, hist));

+		lcrelationList.add(new MyLCRelation(cl.get(cl.size()-1),(MCParticle) p));

             }
         }
 	double neg_energy_total = 0.;

@@ -136,6 +141,7 @@

 	if( neg_energy_total < -Double.MIN_VALUE ) for (Cluster rcl :  cl ) if( Math.abs(((ReconCluster)rcl).getMCParticle().getCharge()) < Double.MIN_VALUE && ((ReconCluster)rcl).getNegEnergy() >= 0. ) 
                                                                                ((ReconCluster)rcl).adjustEnergy(neg_energy_total,pos_energy_weight_total);
         event.put(EventHeader.CLUSTERS, cl, Cluster.class, 0);

+	event.put("ClustersToMCP", lcrelationList, LCRelation.class, 0);

     }