High Pt Muons and Final Status Muons about llp_sleptons_rpv_susy HOT 4 CLOSED

To solve this issue, I added this line to the functionand abs(particle.pid) == 13 to make it only go down the muon decay line:

def get_final_muon_descendant(particle):

    """
    Get the final state muon descendant of a given particle.
    """

    # Check if the particle itself is a final state muon

    if particle.status == 1 and abs(particle.pid) == 13:
        return particle

    # If the particle is a muon and has an end vertex, check its descendants recursively

    elif particle.end_vertex and abs(particle.pid) == 13:
        for p in particle.end_vertex.particles_out:
            final_muon = get_final_muon_descendant(p)
            if final_muon is not None:
                return final_muon

    # If the particle is not a final state muon and does not have an end vertex, return None

    return None

from llp_sleptons_rpv_susy.

This is the output after updating the code:

We are still getting the same muons with high Pt. I looked into the event with muon_Pt ~ 500 GeV

from llp_sleptons_rpv_susy.

Here is the event:

And here is the muon with Pt ~ 0.49 TeV:

As strange as it seems, it is pair produced by a photon!
The actual muon from the smuon decay is marked in green, the suspicious 0.49 TeV muon is circled in red

from llp_sleptons_rpv_susy.

@trholmes
Could you please take a look at this issue? This is the weird muon with pt ~0.5 TeV

I will look at the events with the other 7 muons that were not produced directly by a smuon decay and have pt> 65 GeV

from llp_sleptons_rpv_susy.