feat: calculate PhPerp for SIDIS hadrons and dihadrons

JeffersonLab · Oct 28, 2024 · bf0840a · bf0840a
1 parent 87e2fc3
commit bf0840a
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 0 deletions.
diff --git a/src/iguana/algorithms/physics/DihadronKinematics/Algorithm.cc b/src/iguana/algorithms/physics/DihadronKinematics/Algorithm.cc
@@ -25,6 +25,7 @@ namespace iguana::physics {
           "pdg_b/I",
           "Mh/D",
           "z/D",
+          "PhPerp/D",
           "MX/D",
           "xF/D",
           "phiH/D",
@@ -39,6 +40,7 @@ namespace iguana::physics {
     i_pdg_b    = result_schema.getEntryOrder("pdg_b");
     i_Mh       = result_schema.getEntryOrder("Mh");
     i_z        = result_schema.getEntryOrder("z");
+    i_PhPerp   = result_schema.getEntryOrder("PhPerp");
     i_MX       = result_schema.getEntryOrder("MX");
     i_xF       = result_schema.getEntryOrder("xF");
     i_phiH     = result_schema.getEntryOrder("phiH");
@@ -135,6 +137,9 @@ namespace iguana::physics {
       // calculate z
       double z = p_target.Dot(p_Ph) / p_target.Dot(p_q);
 
+      // calculate PhPerp
+      double PhPerp = tools::RejectVector(p_Ph.Vect(), p_q.Vect()).R();
+
       // calculate Mh
       double Mh = p_Ph.M();
 
@@ -190,6 +195,7 @@ namespace iguana::physics {
       result_bank.putInt(i_pdg_b,      dih_row, had_b.pdg);
       result_bank.putDouble(i_Mh,      dih_row, Mh);
       result_bank.putDouble(i_z,       dih_row, z);
+      result_bank.putDouble(i_PhPerp,  dih_row, PhPerp);
       result_bank.putDouble(i_MX,      dih_row, MX);
       result_bank.putDouble(i_xF,      dih_row, xF);
       result_bank.putDouble(i_phiH,    dih_row, phiH);

diff --git a/src/iguana/algorithms/physics/DihadronKinematics/Algorithm.h b/src/iguana/algorithms/physics/DihadronKinematics/Algorithm.h
@@ -21,6 +21,8 @@ namespace iguana::physics {
     double Mh;
     /// @brief @latex{z}: Momentum fraction of the fragmenting parton carried by the dihadron
     double z;
+    /// @brief @latex{P_h^\perp}: transverse momentum of the dihadron in the @latex{\perp}-frame (transverse to @latex{\vec{q}})
+    double PhPerp;
     /// @brief @latex{M_X(ehhX)}: Missing mass of the dihadron
     double MX;
     /// @brief @latex{x_F}: Feynman-x of the dihadron
@@ -95,6 +97,7 @@ namespace iguana::physics {
       int i_pdg_b;
       int i_Mh;
       int i_z;
+      int i_PhPerp;
       int i_MX;
       int i_xF;
       int i_phiH;