Draft: Resolve "Add Roman Pots Reconstruction Matrix"

Closes #70

benchmarks/far_forward/scripts/hits_far_forward_protons.cxx

@@ -48,9 +48,17 @@ auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
 };
 auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
   std::vector<double> result;
+  ROOT::Math::PxPyPzMVector transFourVect;
   for (size_t i = 0; i < in.size(); ++i) {
-   in[i].RotationY(-crossingAngle);
-   result.push_back(in[i].Theta()*1000);
+    Double_t s = std::sin(-crossingAngle);
+    Double_t c = std::cos(-crossingAngle);
+    Double_t zz = in[i].Pz();
+    Double_t xx = in[i].Px();
+    pz = c*zz - s*xx;
+    px = s*zz + c*xx;
+    transFourVect.SetPxPyPzE(in[i].Px(), in[i].Py(), in[i].Pz, in[i].E());
+
+    result.push_back(transFourVect.Theta()*1000);
   }
   return result;
 };