Skip to content
Snippets Groups Projects
Commit aa8ad68b authored by Ziyue Zhang's avatar Ziyue Zhang
Browse files

Test on t

parent d35832cf
No related branches found
No related tags found
1 merge request!25Add t in dvmp; add a place holder for InvQuant in REC
...@@ -28,7 +28,7 @@ namespace util { ...@@ -28,7 +28,7 @@ namespace util {
// for structure functions // for structure functions
struct inv_quant { // add more when needed struct inv_quant { // add more when needed
double nu, Q2, x; double nu, Q2, x, t;
}; };
// for simu // for simu
...@@ -36,10 +36,12 @@ namespace util { ...@@ -36,10 +36,12 @@ namespace util {
{ {
ROOT::Math::PxPyPzMVector q(parts[0] - parts[2]); ROOT::Math::PxPyPzMVector q(parts[0] - parts[2]);
ROOT::Math::PxPyPzMVector P(parts[3]); ROOT::Math::PxPyPzMVector P(parts[3]);
ROOT::Math::PxPyPzMVector Delta(parts[6] - parts[3]);
double nu = q.Dot(P) / P.mass(); double nu = q.Dot(P) / P.mass();
double Q2 = -q.Dot(q); double Q2 = -q.Dot(q);
inv_quant quantities = {nu, Q2, Q2 / 2. / P.mass() / nu}; double t = Delta.Dot(Delta);
inv_quant quantities = {nu, Q2, Q2 / 2. / P.mass() / nu, t};
return quantities; return quantities;
} }
......
...@@ -179,6 +179,8 @@ int vm_invar(const std::string& config_name) ...@@ -179,6 +179,8 @@ int vm_invar(const std::string& config_name)
t23->Draw(); t23->Draw();
c.Print(fmt::format("{}InvariantQuantities.png", output_prefix).c_str()); c.Print(fmt::format("{}InvariantQuantities.png", output_prefix).c_str());
} }
// TODO we're not actually getting the resolutions yet // TODO we're not actually getting the resolutions yet
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Please register or to comment