Draft: Resolve "Add Roman Pots Reconstruction Matrix"

Closes #70

benchmarks/far_forward/scripts/gen_far_forward_protons.cxx

@@ -21,8 +21,8 @@ void gen_far_forward_protons(int n_events = 100,
   double minMomentum = 270.0; //GeV
   double maxMomentum = 275.0; //GeV
   
-  double cos_theta_min = std::cos(0.0*(M_PI/180.0)); //we need MRAD for the units - syntax is degrees*Pi/180
-  double cos_theta_max = std::cos(0.2864789*(M_PI/180.0)); //1.15 degrees = 20mrad FF region is 0.0 to 20.0 mrad
+  double cos_theta_min = std::cos(0.0); //we need MRAD for the units - syntax is degrees*Pi/180
+  double cos_theta_max = std::cos(0.005); //1.15 degrees = 20mrad FF region is 0.0 to 20.0 mrad
 
   WriterAscii hepmc_output(out_fname);
   int events_parsed = 0;