Skip to content
Snippets Groups Projects

Add t in dvmp; add a place holder for InvQuant in REC

Merged Ziyue Zhang requested to merge ziyue_work_branch into master
Compare and
2 files
+ 67
35
Compare changes
  • Side-by-side
  • Inline
Files
2
@@ -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,16 +36,19 @@ namespace util {
@@ -36,16 +36,19 @@ 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;
}
}
inline double get_nu_simu(inv_quant quantities) { return quantities.nu / 1000.; }
inline double get_nu(inv_quant quantities) { return quantities.nu / 1000.; }
inline double get_Q2_simu(inv_quant quantities) { return quantities.Q2; }
inline double get_Q2(inv_quant quantities) { return quantities.Q2; }
inline double get_x_simu(inv_quant quantities) { return quantities.x; }
inline double get_x(inv_quant quantities) { return quantities.x; }
 
inline double get_t(inv_quant quantities) { return quantities.t; }
// for tracking, add later
// for tracking, add later
Loading