Ongoing dis_electron Q2 analysis. An issue with Q2 being basically zero for the electrons.
All threads resolved!
So I have made the branch and got the functions necessary, but there is an issue that the electrons Q2 is essentially zero. So the sorted particles within each event with the highest momenta are basically electrons that "missed" the target. Maybe I should sort by highest Q2 and rerun.
Merge request reports
Activity
17 #include <eicd/ReconstructedParticleData.h> 18 #include <TH1D.h> 19 #include <TFitResult.h> 20 #include <TRandom3.h> 21 #include <algorithm> 22 #include <utility> 23 24 25 //Reconstuction functions 26 bool mom_sort_recon(eic::ReconstructedParticleData &part1, eic::ReconstructedParticleData &part2) 27 { 28 return (part1.p.x*part1.p.x + part1.p.y*part1.p.y + part1.p.z*part1.p.z > 29 part2.p.x*part2.p.x + part2.p.y*part2.p.y + part2.p.z*part2.p.z); 30 } 31 32 inline auto momenta_from_recon(const std::vector<eic::ReconstructedParticleData>& parts) changed this line in version 2 of the diff
14 14 #include <nlohmann/json.hpp> 15 15 #include <string> 16 16 #include <vector> 17 #include <eicd/ReconstructedParticleData.h> 18 #include <TH1D.h> 19 #include <TFitResult.h> 20 #include <TRandom3.h> 21 #include <algorithm> 22 #include <utility> 23 24 25 //Reconstuction functions 26 bool mom_sort_recon(eic::ReconstructedParticleData &part1, eic::ReconstructedParticleData &part2) changed this line in version 2 of the diff
24 25 //Reconstuction functions 26 bool mom_sort_recon(eic::ReconstructedParticleData &part1, eic::ReconstructedParticleData &part2) 27 { 28 return (part1.p.x*part1.p.x + part1.p.y*part1.p.y + part1.p.z*part1.p.z > 29 part2.p.x*part2.p.x + part2.p.y*part2.p.y + part2.p.z*part2.p.z); 30 } 31 32 inline auto momenta_from_recon(const std::vector<eic::ReconstructedParticleData>& parts) 33 { 34 std::vector <eic::ReconstructedParticleData> sort_parts = parts; 35 sort(sort_parts.begin(), sort_parts.end(), mom_sort_recon); 36 return sort_parts; 37 } 38 39 inline auto Q2_from_recon(const std::vector<eic::ReconstructedParticleData>& parts) changed this line in version 2 of the diff
added 1 commit
- 176163ae - Updated functions and changed the ebeam sign
added 1 commit
- 7db1455f - Wrote plotting function, fairly sure that the benchmark is more or less complete.
added 12 commits
-
7db1455f...33a3f85f - 6 commits from branch
master - ec7a0114 - Ongoing dis_electron Q2 analysis. An issue with Q2 being basically zero for the electrons.
- 9ed34199 - Changed the sorting function to report minimum momemta
- d9ab8dee - Updated functions. Still get negative Q2 that is a delta function near zero
- 6155842f - Updated functions and changed the ebeam sign
- 533d82c2 - Added plotting and test results
- 8e906196 - Wrote plotting function, fairly sure that the benchmark is more or less complete.
Toggle commit list-
7db1455f...33a3f85f - 6 commits from branch
added 22 commits
-
8e906196...776434d9 - 16 commits from branch
master - 3d618a33 - Ongoing dis_electron Q2 analysis. An issue with Q2 being basically zero for the electrons.
- 2a311cec - Changed the sorting function to report minimum momemta
- 0dcb7572 - Updated functions. Still get negative Q2 that is a delta function near zero
- dae9dbbc - Updated functions and changed the ebeam sign
- 4afb2a16 - Added plotting and test results
- f1c0c4b0 - Wrote plotting function, fairly sure that the benchmark is more or less complete.
Toggle commit list-
8e906196...776434d9 - 16 commits from branch
Please register or sign in to reply