Skip to content
Snippets Groups Projects
Commit 65972e64 authored by Alex Jentsch's avatar Alex Jentsch
Browse files

Hits look good now. Clean up histogram ranges.

parent 71d27e80
No related branches found
No related tags found
No related merge requests found
...@@ -20,7 +20,7 @@ using ROOT::RDataFrame; ...@@ -20,7 +20,7 @@ using ROOT::RDataFrame;
using namespace ROOT::VecOps; using namespace ROOT::VecOps;
double local_x_offset = -833.3878326; double local_x_offset = -833.3878326;
double crossingAngle = -0.025;
auto p_track = [](std::vector<eic::TrackParametersData> const& in) { auto p_track = [](std::vector<eic::TrackParametersData> const& in) {
std::vector<double> result; std::vector<double> result;
...@@ -49,7 +49,8 @@ auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { ...@@ -49,7 +49,8 @@ auto momentum = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) { auto theta = [](std::vector<ROOT::Math::PxPyPzMVector> const& in) {
std::vector<double> result; std::vector<double> result;
for (size_t i = 0; i < in.size(); ++i) { 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; return result;
}; };
...@@ -148,7 +149,7 @@ int hits_far_forward_protons(const char* fname = "sim_far_forward_protons.root") ...@@ -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"}) .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");
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment