def reconstruct_track(self, track, clusters=None): vertex = track.path.points['vertex'] pdg_id = 211 * track.charge mass, charge = particle_data[pdg_id] p4 = TLorentzVector() p4.SetVectM(track.p3, mass) particle = Particle(p4, vertex, charge, pdg_id) particle.set_path(track.path) particle.clusters = clusters track.locked = True return particle