Skip to content
Snippets Groups Projects

Draft: Resolve "Add Roman Pots Reconstruction Matrix"

Open Alex Jentsch requested to merge 70-add-roman-pots-reconstruction-matrix into master
1 file
+ 4
3
Compare changes
  • Side-by-side
  • Inline
@@ -20,7 +20,7 @@ using ROOT::RDataFrame;
using namespace ROOT::VecOps;
double local_x_offset = -833.3878326;
double crossingAngle = -0.025;
auto p_track = [](std::vector<eic::TrackParametersData> const& in) {
std::vector<double> result;
@@ -49,7 +49,8 @@ auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
std::vector<double> result;
for (size_t i = 0; i < in.size(); ++i) {
result.push_back(in[i].Theta()*1000);
double unRotateTheta = in[i].Theta() - crossingAngle;
result.push_back(unRotateTheta*1000);
}
return result;
};
@@ -148,7 +149,7 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root")
.Define("y_pos", y_pos, {"xy_hit_pos"})
;
auto h_local_pos = df0.Histo2D({"h_local_pos", ";x [mm]; y [mm] ", 100, -1200.0, -500.0, 100, -100.0, 100.0}, "x_pos", "y_pos");
auto h_local_pos = df0.Histo2D({"h_local_pos", ";x [mm]; y [mm] ", 100, -1000.0, -650.0, 100, -100.0, 100.0}, "x_pos", "y_pos");
Loading