Comments (10)
Now checking the entire decay chain recursively for every particle (not only muons and neutrinos) by changing the code for "get_decay_products" function to:
def get_decay_products(particle, event_num):
"""
Get the decay products of a given particle.
"""
proper_decay_products = []
other_particles = []
# Check if the particle itself is a final state muon or neutrino
if particle.status == 1 and abs(particle.pid) in [13, 12, -13, -12]:
proper_decay_products.append((particle.pid, particle.id))
elif particle.status == 1:
other_particles.append((particle.pid, particle.id, event_num))
# If the particle has an end vertex, check its descendants recursively
if particle.end_vertex:
for p in particle.end_vertex.particles_out:
proper, other = get_decay_products(p, event_num)
proper_decay_products.extend(proper)
other_particles.extend(other)
return proper_decay_products, other_particles
We get this output:
Proper decay products: Counter({13: 134849, -13: 134840, -12: 134688, 12: 134679})
Other particles: Counter({22: 190298, -211: 956, 211: 929, 11: 535, -11: 533, 321: 132,
2212: 111, -2212: 111,
130: 92, -321: 87, -14: 19, -2112: 10, 2112: 10, -16: 5, 16: 5, 14: 3})
Please notice that this code overcounts particles, I'm looking into it, but I looked at one of the desired events. I already verified that such particles only appear when children particles of photons are considered.
from llp_sleptons_rpv_susy.
Currently looking at the events of those unexpected particles
from llp_sleptons_rpv_susy.
# Find the first non photon PID in other_particles, which is the first particle that is out of a photon end_vertix
first_non_22 = next((pid, id, event_num) for pid, id, event_num in other_particles if pid != 22)
Output:
First non-22 PID in 'other_particles': PID = 11, ID = 1399, Event number = 11
from llp_sleptons_rpv_susy.
from llp_sleptons_rpv_susy.
The photon created a pair of e- e+ with different energies.
from llp_sleptons_rpv_susy.
This is the code used to print the event diagram with the particle of interest and its vertex position
import pyhepmc
from pyhepmc.view import savefig
filename = "tag_1_pythia8_events.hepmc"
# pyhepmc.open can read most HepMC formats using auto-detection
with pyhepmc.open(filename) as f:
# Loop over events in the file
for i, event in enumerate(f):
# Save a plot of the event
if event.event_number == 11:
savefig(event, f"event{i}.svg")
savefig(event, f"event{i}.png")
savefig(event, f"event{i}.pdf")
for particle in event.particles:
if particle.id ==1399 and particle.pid==11:
ver=particle.production_vertex
print(ver)
However, it creates diagrams for every single event with "11" in its event number, but that is not a concern at the moment.
from llp_sleptons_rpv_susy.
I checked one of the events where we get unexpected particles, this particular one was an electron coming from an e- e+ pair production by a photon. All unexpected particles in the chain come from a vertex with a photon coming in. The process is detailed above for your reference. Thanks!
@trholmes
from llp_sleptons_rpv_susy.
I don't totally understand from looking at your diagram. Yes the one of your muons radiates two photons, and the photons create e+e- pairs, but I don't see any other particles coming from those chains. Can you clarify?
from llp_sleptons_rpv_susy.
This is an event where unexpected particles appear, the e+ e- themselves are not expected particles from our decay chain (smuon ---> mu + neutrino). It just happened that the photon in this event was creating two final status electrons, but if I pick another event, I could for example see a photon creating pi+ pi-
I hope this clarifies what I meant.
from llp_sleptons_rpv_susy.
e-e+ was not the best example. Here is another example, the most observed unexpected particle in our chain is -211 (pi-)
Here is a pi- in a photon chain (where the photon itself is from a muon coming from smuon)
PID = -211, ID = 552, Event number = 868
production vertex: (FourVector(-10.8, -21.9, 22.9, 46.4))
Event #868
Here is the decay chain:
The same picture but clearer (the chain of the photon with 58 GeV)
This one has a lot of unexpected particles compared to photon -----> e- e+
e-e+ are also not expected from our smuon decay chain and I should have mentioned that.
I hope this clarifies everything.
from llp_sleptons_rpv_susy.
Related Issues (20)
- Fraction of surviving events
- Automation HOT 4
- Modifying the UFO Model Parameters. Tuesday 09/26 HOT 1
- Getting started
- Efficiencies as weights. October 3rd
- Selection cuts and efficiency code using weights
- Turning on Pythia via TXT Problem Solved BUT New SERF Problem HOT 1
- Oct 24 Update: Pythia Issues Resolved HOT 1
- For my own reference
- Weights
- Expected Events from Weights: Nov 21st HOT 1
- Tuesday Dec 5th
- Acceptance_Weights_Expected Gives Inaccurate Output HOT 1
- Wrong Units in .HEPMC File HOT 1
- Cutflow: pass trigger and at least 2 baseline leptons selection
- Condor Submission to and from EOS HOT 4
- Acceptance, Friday March 1st HOT 2
- Generation via CMVFS/LCG Source HOT 7
- Immortal Smuons HOT 1
- Fri Apr. 26 - 10 ns
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from llp_sleptons_rpv_susy.