Skip to content
Snippets Groups Projects
Commit d51e685c authored by Marshall Scott's avatar Marshall Scott
Browse files

fixed pions back to generating pi+

parent 6eec429b
No related branches found
No related tags found
1 merge request!38Fixed Pions 2
Hide whitespace changes
Inline Side-by-side
benchmarks/barrel_ecal/scripts/emcal_barrel_pions.cxx
+ 2
− 4
View file @ d51e685c
...@@ -56,15 +56,13 @@ void emcal_barrel_pions(int n_events = 1e6, double e_start = 0.0, double e_end = ...@@ -56,15 +56,13 @@ void emcal_barrel_pions(int n_events = 1e6, double e_start = 0.0, double e_end =
// type 1 is final state // type 1 is final state
// pdgid 211 - pion+ 139.570 MeV/c^2 // pdgid 211 - pion+ 139.570 MeV/c^2
// pdgid 111 - pion0 134.977 MeV/c^2 // pdgid 111 - pion0 134.977 MeV/c^2
//GenParticlePtr p3 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.139570 * 0.139570))), 211, 1); GenParticlePtr p3 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.139570 * 0.139570))), 211, 1);
GenParticlePtr p4 = std::make_shared<GenParticle>(FourVector(px, py, pz, sqrt(p * p + (0.134977 * 0.134977))), 111, 1);
GenVertexPtr v1 = std::make_shared<GenVertex>(); GenVertexPtr v1 = std::make_shared<GenVertex>();
v1->add_particle_in(p1); v1->add_particle_in(p1);
v1->add_particle_in(p2); v1->add_particle_in(p2);
//v1->add_particle_out(p3); v1->add_particle_out(p3);
v1->add_particle_out(p4);
evt.add_vertex(v1); evt.add_vertex(v1);
if (events_parsed == 0) { if (events_parsed == 0) {
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment