Skip to content
Snippets Groups Projects
Commit 75ec438c authored by Ziyue Zhang's avatar Ziyue Zhang
Browse files

WIP: Test on PID instead of pmass

parent cb2c25ec
Branches
Tags
1 merge request!39Replacing nu with y
This commit is part of merge request !39. Comments created here will be created in the context of that merge request.
...@@ -40,7 +40,7 @@ namespace util { ...@@ -40,7 +40,7 @@ namespace util {
} }
//========================================import Geant4 and set the wanted particles in the intended order======================== //========================================import Geant4 and set the wanted particles in the intended order========================
//0:e0 1:p0 2:e1 3:p1 4:recoil system (without p1) 5:l1 from 4 6:l2 from 4 //0:e0 1:p0 2:e1 3:p1 4:recoil system (without p1) 5:l1 from 4 6:l2 from 4
/*inline auto momenta_ordered_sim(const std::vector<dd4pod::Geant4ParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events inline auto momenta_ordered_sim(const std::vector<dd4pod::Geant4ParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
std::vector<ROOT::Math::PxPyPzMVector> momenta{7}; std::vector<ROOT::Math::PxPyPzMVector> momenta{7};
int order_map[7] = {0, 3, 2, 6, 5, 7, 8}; int order_map[7] = {0, 3, 2, 6, 5, 7, 8};
for(int i = 0 ; i < 7 ; i++){ for(int i = 0 ; i < 7 ; i++){
...@@ -52,7 +52,7 @@ namespace util { ...@@ -52,7 +52,7 @@ namespace util {
momenta[i].SetPxPyPzE(px, py, pz, e); momenta[i].SetPxPyPzE(px, py, pz, e);
} }
return momenta; return momenta;
}*/ }
//import Reconstructed particles and set the wanted particles in the intended order======================== //import Reconstructed particles and set the wanted particles in the intended order========================
/*inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events /*inline auto momenta_ordered_rec(const std::vector<eic::ReconstructedParticleData>& parts, std::string_view mother, std::string_view daughter){//mother and daughter are not used yet; will be useful when generater is different and/or when the mcparticles doesn't follow the same order in all events
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment