Draft: electron_incl_kin DEBUG

assigned to @wdconinc

electron_incl_kin   DEBUG pi = -3.66374e-15 6.55032e-15 274.999 275 
electron_incl_kin   DEBUG ei = 2.13163e-14 1.42109e-14 -18 18 
electron_incl_kin   DEBUG ef = [(-1.44807 1.84416 -18.0475 18.1992 , 6 1233367242 )]
electron_incl_kin   DEBUG q = 1.44807 -1.84416 0.0474999 -0.199178 
electron_incl_kin   DEBUG x,y,Q2,W,nu = -nan,-0.00685218,-nan,-nan,-72.2994

      const auto q = ei.subtract(ef.front().first);
      kin.Q2(q.mass());
      kin.y(q.dot(pi) / ei.dot(pi));
      kin.nu(q.dot(pi) / m_proton);
      kin.x(kin.Q2() / (2. * q.dot(pi)));
      kin.W(sqrt(pi.add(q).mass()));

mentioned in merge request EIC/juggler!315 (merged)

Updated to

      const auto [ef, scatID] = electrons.front();
      const auto q = ei.subtract(ef);
      const auto q_dot_pi = q.dot(pi);
      kin.Q2(-q.dot(q));
      kin.y(q_dot_pi / ei.dot(pi));
      kin.nu(q_dot_pi / m_proton);
      kin.x(kin.Q2() / (2. * q_dot_pi));
      kin.W((pi.add(q)).mass());

electron_incl_kin   DEBUG pi = -3.66374e-15 6.55032e-15 274.999 275 
electron_incl_kin   DEBUG ei = 2.13163e-14 1.42109e-14 -18 18 
electron_incl_kin   DEBUG ef = -1.42251 1.80934 -17.7516 17.9002 
electron_incl_kin   DEBUG q = 1.42251 -1.80934 -0.248393 0.0998085 
electron_incl_kin   DEBUG x,y,Q2,W,nu = 0.0279306,0.00967225,5.34898,187.041,102.055

closed