@@ -155,9 +155,9 @@ namespace BlipUtils {
155155 float totE = tblip.Energy + pinfo.depEnergy ;
156156 float w1 = tblip.Energy /totE;
157157 float w2 = pinfo.depEnergy /totE;
158- tblip.Position .SetXYZ ( w1*tblip.Position .X () + w2*pinfo.position .X (),
159- w1*tblip.Position .Y () + w2*pinfo.position .Y (),
160- w1*tblip.Position .Z () + w2*pinfo.position .Z ());
158+ tblip.Position .SetXYZ ( w1*tblip.Position .X () + w2*pinfo.Position .X (),
159+ w1*tblip.Position .Y () + w2*pinfo.Position .Y (),
160+ w1*tblip.Position .Z () + w2*pinfo.Position .Z ());
161161 tblip.Time = w1*tblip.Time + w2*pinfo.time ;
162162 tblip.LeadCharge = pinfo.depElectrons ;
163163 // ... if the particle isn't a match, show's over
@@ -206,11 +206,9 @@ namespace BlipUtils {
206206 float w1 = blip_i.DepElectrons /totQ;
207207 float w2 = blip_j.DepElectrons /totQ;
208208 blip_i.Energy += blip_j.Energy ;
209- geo::vect::MiddlePointAccumulator mpalg;
210- mpalg.add (&blip_i.Position , w1);
211- mpalg.add (&blip_j.Position , w2);
212- blip_i.Position = mpalg.middlePoint ();
213- mpalg.clear ();
209+ blip_i.Position .SetXYZ ( w1*blip_i.Position .X () + w2*blip_j.Position .X (),
210+ w1*blip_i.Position .Y () + w2*blip_j.Position .Y (),
211+ w1*blip_i.Position .Z () + w2*blip_j.Position .Z ());
214212 blip_i.DriftTime = w1*blip_i.DriftTime + w2*blip_j.DriftTime ;
215213 blip_i.Time = w1*blip_i.Time + w2*blip_j.Time ;
216214 blip_i.DepElectrons += blip_j.DepElectrons ;
0 commit comments