Skip to content
Snippets Groups Projects

Fix dvmp

Merged Whitney Armstrong requested to merge dvmp_fix into master
1 file
+ 4
4
Compare changes
  • Side-by-side
  • Inline
@@ -91,14 +91,14 @@ int vm_mass(const std::string& config_name)
@@ -91,14 +91,14 @@ int vm_mass(const std::string& config_name)
// types
// types
auto find_decay_pair = [vm_mass, decay_mass](const std::vector<ROOT::Math::PxPyPzMVector>& parts) {
auto find_decay_pair = [vm_mass, decay_mass](const std::vector<ROOT::Math::PxPyPzMVector>& parts) {
return util::find_decay_pair(parts, vm_mass, decay_mass);
return common_bench::find_decay_pair(parts, vm_mass, decay_mass);
};
};
// util::PrintGeant4(mcparticles2);
// common_bench::PrintGeant4(mcparticles2);
// Define analysis flow
// Define analysis flow
auto d_im = d.Define("p_rec", util::momenta_RC, {"DummyReconstructedParticles"}) //using dummy rc
auto d_im = d.Define("p_rec", common_bench::momenta_RC, {"DummyReconstructedParticles"}) //using dummy rc
.Define("N", "p_rec.size()")
.Define("N", "p_rec.size()")
.Define("p_sim", util::momenta_from_simulation, {"mcparticles2"})
.Define("p_sim", common_bench::momenta_from_simulation, {"mcparticles2"})
.Define("decay_pair_rec", find_decay_pair, {"p_rec"})
.Define("decay_pair_rec", find_decay_pair, {"p_rec"})
.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")
Loading