diff --git a/benchmarks/u_omega/analysis/demo.cxx b/benchmarks/u_omega/analysis/demo.cxx index 82432a671a24766973b7ebaa348c6161824bf11b..46f7b472b536615cfa89200d365a808c8a9e39ba 100644 --- a/benchmarks/u_omega/analysis/demo.cxx +++ b/benchmarks/u_omega/analysis/demo.cxx @@ -32,7 +32,6 @@ auto p_track = [](std::vector<edm4eic::TrackParametersData> const& in) { return result; }; - auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { std::vector<double> result; for (size_t i = 0; i < in.size(); ++i) { @@ -40,6 +39,7 @@ auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { } return result; }; + auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { std::vector<double> result; for (size_t i = 0; i < in.size(); ++i) { @@ -47,15 +47,7 @@ auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { } return result; }; -auto fourvec = [](ROOT::VecOps::RVec<dd4pod::Geant4ParticleData> const& in) { - std::vector<ROOT::Math::PxPyPzMVector> result; - ROOT::Math::PxPyPzMVector lv; - for (size_t i = 0; i < in.size(); ++i) { - lv.SetCoordinates(in[i].ps.x, in[i].ps.y, in[i].ps.z, in[i].mass); - result.push_back(lv); - } - return result; -}; + auto recfourvec = [](ROOT::VecOps::RVec<edm4eic::ReconstructedParticleData> const& in) { std::vector<ROOT::Math::PxPyPzMVector> result; ROOT::Math::PxPyPzMVector lv; @@ -89,14 +81,6 @@ void demo(const char* fname = "rec_dvcs.root"){ PxPyPzMVector p_ebeam = {0,0,-10, 0.000511}; PxPyPzMVector p_pbeam = {0,0,275, 0.938 }; - auto eprime = [](ROOT::VecOps::RVec<dd4pod::Geant4ParticleData> const& in) { - for(const auto& p : in){ - if(p.pdgID == 11 ) { - return PxPyPzMVector(p.ps.x,p.ps.y,p.ps.z,p.mass); - } - } - return PxPyPzMVector(0,0,0,0); - }; auto q_vec = [=](PxPyPzMVector const& p) { return p_ebeam - p; };