PWGEM/Dilepton: update muon QA

Author: dsekihat

PWGEM/Dilepton/Core/SingleTrackQC.h

@@ -266,8 +266,8 @@ struct SingleTrackQC {
       fRegistry.addClone("Track/positive/", "Track/negative/");
     } else if constexpr (pairtype == o2::aod::pwgem::dilepton::utils::pairutil::DileptonPairType::kDimuon) {
       const AxisSpec axis_pt{ConfPtlBins, "p_{T,#mu} (GeV/c)"};
-      const AxisSpec axis_eta{40, -5, -1, "#eta_{#mu}"};
-      const AxisSpec axis_phi{90, 0, 2 * M_PI, "#varphi_{#mu} (rad.)"};
+      const AxisSpec axis_eta{100, -6, -1, "#eta_{#mu}"};
+      const AxisSpec axis_phi{180, 0, 2 * M_PI, "#varphi_{#mu} (rad.)"};
       const AxisSpec axis_dca{ConfDCABins, "DCA_{#mu}^{XY} (#sigma)"};
 
       // track info

PWGEM/Dilepton/Core/SingleTrackQCMC.h

@@ -310,8 +310,8 @@ struct SingleTrackQCMC {
       }
     } else if constexpr (pairtype == o2::aod::pwgem::dilepton::utils::pairutil::DileptonPairType::kDimuon) {
       const AxisSpec axis_pt{ConfPtlBins, "p_{T,#mu} (GeV/c)"};
-      const AxisSpec axis_eta{25, -4.5, -2.0, "#eta_{#mu}"};
-      const AxisSpec axis_phi{36, 0.0, 2 * M_PI, "#varphi_{#mu} (rad.)"};
+      const AxisSpec axis_eta{100, -6, -1, "#eta_{#mu}"};
+      const AxisSpec axis_phi{180, 0.0, 2 * M_PI, "#varphi_{#mu} (rad.)"};
       const AxisSpec axis_dca{ConfDCABins, "DCA_{#mu}^{XY} (#sigma)"};
       const AxisSpec axis_charge_gen{3, -1.5, +1.5, "true charge"};
 

PWGEM/Dilepton/Utils/EMTrackUtilities.h

@@ -59,13 +59,18 @@ float fwdDcaXYinSigma(T const& track)
   float cXY = track.cXYatDCA();
   float dcaX = track.fwdDcaX(); // in cm
   float dcaY = track.fwdDcaY(); // in cm
+  float dcaXY = std::sqrt(dcaX * dcaX + dcaY * dcaY);
+  float dFdx = 2.f * dcaX / dcaXY;
+  float dFdy = 2.f * dcaY / dcaXY;
+  float sigma_dcaXY = std::sqrt(cXX * dFdx * dFdx + cYY * dFdy * dFdy + 2.f * cXY * dFdx * dFdy);
+  return dcaXY / sigma_dcaXY;
 
-  float det = cXX * cYY - cXY * cXY; // determinant
-  if (det < 0) {
-    return 999.f;
-  } else {
-    return std::sqrt(std::fabs((dcaX * dcaX * cYY + dcaY * dcaY * cXX - 2. * dcaX * dcaY * cXY) / det / 2.)); // dca xy in sigma
-  }
+  // float det = cXX * cYY - cXY * cXY; // determinant
+  // if (det < 0) {
+  //   return 999.f;
+  // } else {
+  //   return std::sqrt(std::fabs((dcaX * dcaX * cYY + dcaY * dcaY * cXX - 2. * dcaX * dcaY * cXY) / det / 2.)); // dca xy in sigma
+  // }
 }
 //_______________________________________________________________________
 template <typename T>