Skip to content
Snippets Groups Projects

Add J/psi resolution; re-write RC with Dummy instead of rc tracks

Merged Ziyue Zhang requested to merge ziyue_work_branch into master
2 files
+ 8
9
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -84,10 +84,10 @@ int vm_mass(const std::string& config_name)
@@ -84,10 +84,10 @@ int vm_mass(const std::string& config_name)
// util::PrintGeant4(mcparticles2);
// util::PrintGeant4(mcparticles2);
// Define analysis flow
// Define analysis flow
auto d_im = d.Define("p_rec", momenta_from_tracking, {"outputTrackParameters"})
auto d_im = d.Define("p_rec", momenta_from_tracking, {"outputTrackParameters"})
.Define("p_rec_test", util::p_dummy, {"DummyReconstructedParticles"})
.Define("p_rec_test", util::momenta_RC, {"DummyReconstructedParticles"})
.Define("N", "p_rec.size()")
.Define("N", "p_rec.size()")
.Define("p_sim", util::momenta_from_simulation, {"mcparticles2"})
.Define("p_sim", util::momenta_from_simulation, {"mcparticles2"})
.Define("decay_pair_rec", find_decay_pair, {"p_rec"})
.Define("decay_pair_rec", find_decay_pair, {"p_rec_test"})
.Define("decay_pair_sim", find_decay_pair, {"p_sim"})
.Define("decay_pair_sim", find_decay_pair, {"p_sim"})
.Define("p_vm_rec", "decay_pair_rec.first + decay_pair_rec.second")
.Define("p_vm_rec", "decay_pair_rec.first + decay_pair_rec.second")
.Define("p_vm_sim", "decay_pair_sim.first + decay_pair_sim.second")
.Define("p_vm_sim", "decay_pair_sim.first + decay_pair_sim.second")
Loading