@@ -264,6 +264,20 @@ DECLARE_SOA_DYNAMIC_COLUMN(OPAngle, opAngle,
264264 return v1.Angle (v2);
265265 });
266266
267+ DECLARE_SOA_DYNAMIC_COLUMN (OPAngleKShortKStar, opAngleKShortKStar,
268+ [](float photonPx, float photonPy, float photonPz, float kshortPx, float kshortPy, float kshortPz) {
269+ TVector3 v1 (photonPx + kshortPx, photonPy + kshortPy, photonPz + kshortPz);
270+ TVector3 v2 (kshortPx, kshortPy, kshortPz);
271+ return v1.Angle (v2);
272+ });
273+
274+ DECLARE_SOA_DYNAMIC_COLUMN (OPAnglePhotonKStar, opAnglePhotonKStar,
275+ [](float photonPx, float photonPy, float photonPz, float kshortPx, float kshortPy, float kshortPz) {
276+ TVector3 v1 (photonPx, photonPy, photonPz);
277+ TVector3 v2 (kshortPx + photonPx, kshortPy + photonPy, kshortPz + photonPz);
278+ return v1.Angle (v2);
279+ });
280+
267281// Photon
268282DECLARE_SOA_DYNAMIC_COLUMN (PhotonPt, photonPt, // ! Transverse momentum in GeV/c
269283 [](float photonPx, float photonPy) -> float {
@@ -332,6 +346,8 @@ DECLARE_SOA_TABLE(KStarCores, "AOD", "KSTARCORES",
332346 kstarCore::Eta<kstarCore::PhotonPx, kstarCore::PhotonPy, kstarCore::PhotonPz, kstarCore::KShortPx, kstarCore::KShortPy, kstarCore::KShortPz>,
333347 kstarCore::Radius<kstarCore::X, kstarCore::Y>,
334348 kstarCore::OPAngle<kstarCore::PhotonPx, kstarCore::PhotonPy, kstarCore::PhotonPz, kstarCore::KShortPx, kstarCore::KShortPy, kstarCore::KShortPz>,
349+ kstarCore::OPAngleKShortKStar<kstarCore::PhotonPx, kstarCore::PhotonPy, kstarCore::PhotonPz, kstarCore::KShortPx, kstarCore::KShortPy, kstarCore::KShortPz>,
350+ kstarCore::OPAnglePhotonKStar<kstarCore::PhotonPx, kstarCore::PhotonPy, kstarCore::PhotonPz, kstarCore::KShortPx, kstarCore::KShortPy, kstarCore::KShortPz>,
335351
336352 kstarCore::PhotonPt<kstarCore::PhotonPx, kstarCore::PhotonPy>,
337353 kstarCore::PhotonP<kstarCore::PhotonPx, kstarCore::PhotonPy, kstarCore::PhotonPz>,
0 commit comments