lcsim/src/org/lcsim/mc/fast/reconstructedparticle
diff -u -r1.11 -r1.12
--- MCFastReconstructedParticle.java 18 Sep 2007 03:46:27 -0000 1.11
+++ MCFastReconstructedParticle.java 29 Nov 2007 09:50:08 -0000 1.12
@@ -45,14 +45,22 @@
private BasicHepLorentzVector p_reco = new BasicHepLorentzVector();
private BasicHepLorentzVector p_track = new BasicHepLorentzVector();
- public MCFastReconstructedParticle(Track t, ParticleType type, Particle p, Cluster assoc_c, double wtcal)
+ public MCFastReconstructedParticle(Track t, ParticleType type, Particle p, Cluster assoc_c, double wtcal, boolean refPoint000)
{
MCFast.log.info(" PDGID= "+type.getPDGID()+" t.getPX,...= "+t.getPX()+" "+t.getPY()+" "+t.getPZ());
_mass = type.getMass();
addTrack(t);
_charge = t.getCharge();
- // Use true point of origin for reference point for now.
+ // Use (0,0,0) for reference point if refPoint000=true
+ // Use true point of origin for reference point if refPoint000=false
+ if(refPoint000)
+ {
+ _referencePoint = new BasicHep3Vector(0, 0, 0);
+ }
+ else
+ {
_referencePoint = p.getOrigin();
+ }
e_track = sqrt(((ReconTrack)t).getDocaMomentumVec(_referencePoint).magnitudeSquared()+ _mass*_mass);
p_track.setV3(e_track, ((ReconTrack)t).getDocaMomentumVec(_referencePoint));
p3_track = p_track.v3();