lcsim/src/org/lcsim/contrib/uiowa/template
diff -u -r1.1 -r1.2
--- ChargedHadronIdentifier.java 27 Jan 2006 23:54:01 -0000 1.1
+++ ChargedHadronIdentifier.java 4 Feb 2006 01:22:24 -0000 1.2
@@ -112,6 +112,11 @@
BasicReconstructedParticle part = new BasicReconstructedParticle();
part.addTrack(tr);
part.addCluster(matchedCluster);
+ double[] trackMomentum = tr.getMomentum();
+ double trackMomentumMagSq = (trackMomentum[0]*trackMomentum[0] + trackMomentum[1]*trackMomentum[1] + trackMomentum[2]*trackMomentum[2]);
+ double mass = 0.140; // FIXME: Everything is a pion!
+ double energy = Math.sqrt(trackMomentumMagSq + mass*mass);
+ part.setEnergy(energy);
outputParticleList.add(part);
} else {
// No match -- put the track into the list of unused tracks