Add Z reco-truth histograms for identified Primary Vertex

InnerDetector/InDetValidation/InDetPhysValMonitoring/share/InDetPVMPlotDefCommon.xml

@@ -756,6 +756,14 @@
   <x title="Vertex Truth-Matched type" n="5" lo="-0.5" hi="4.5"/>
   <y title="Entries"/>
 </h>
+<h id="vx_z_diff" type="TH1F" title="PV selected-truth Z difference">
+  <x title="Vertex Truth-Matched type" n="100" lo="-0.5" hi="0.5"/>
+  <y title="Entries"/>
+</h>
+<h id="vx_z_diff_pull" type="TH1F" title="PV selected-truth Z difference pull">
+  <x title="Vertex Truth-Matched type" n="100" lo="-5.0" hi="5.0"/>
+  <y title="Entries"/>
+</h>
 <h id="vx_hs_classification" type="TH1F" title="Event hardscatter vertex classification">
   <x title="Hardscatter classification" n="5" lo="-0.5" hi="4.5"/>
   <y title="Entries"/>

InnerDetector/InDetValidation/InDetPhysValMonitoring/src/InDetPerfPlot_VertexTruthMatching.cxx

@@ -81,6 +81,8 @@ InDetPerfPlot_VertexTruthMatching::InDetPerfPlot_VertexTruthMatching(InDetPlotBa
 void InDetPerfPlot_VertexTruthMatching::initializePlots() {
 
     book(m_vx_type_truth,"vx_type_truth");
+    book(m_vx_z_diff,"vx_z_diff");
+    book(m_vx_z_diff_pull,"vx_z_diff_pull");
     if (m_iDetailLevel >= 200) {
         book(m_vx_hs_classification,"vx_hs_classification");
         book(m_vx_nReco_vs_nTruth_inclusive,"vx_nReco_vs_nTruth_inclusive");
@@ -314,7 +316,15 @@ const xAOD::TruthVertex* InDetPerfPlot_VertexTruthMatching::getTruthVertex(const
     return truthVtx;
 }
 
-void InDetPerfPlot_VertexTruthMatching::fill(const xAOD::Vertex& vertex) {
+void InDetPerfPlot_VertexTruthMatching::fill(const xAOD::Vertex& vertex, const xAOD::TruthVertex * tvrt) {
+
+    if(tvrt){
+      float diff_z=vertex.z()-tvrt->z();
+      const AmgSymMatrix(3)& covariance = vertex.covariancePosition();
+      float err_z = fabs(Amg::error(covariance, 2)) > 1e-7 ? Amg::error(covariance, 2) : 1000.;
+      fillHisto(m_vx_z_diff,diff_z);
+      fillHisto(m_vx_z_diff_pull,diff_z/err_z);
+    }
 
     // Get the match type info for each vertex:
     const static xAOD::Vertex::Decorator<InDetVertexTruthMatchUtils::VertexMatchType> recoVtxMatchTypeInfo("VertexMatchType");

InnerDetector/InDetValidation/InDetPhysValMonitoring/src/InDetPerfPlot_VertexTruthMatching.h

@@ -41,13 +41,15 @@
 class InDetPerfPlot_VertexTruthMatching: public InDetPlotBase {
 public:
     InDetPerfPlot_VertexTruthMatching(InDetPlotBase* pParent, const std::string& dirName, const int iDetailLevel = 10);
-    void fill(const xAOD::Vertex& vertex);
+    void fill(const xAOD::Vertex& vertex, const xAOD::TruthVertex * tvrt =0);
     void fill(const xAOD::VertexContainer& vertexContainer, const std::vector<const xAOD::TruthVertex*>& truthHSVertices, const std::vector<const xAOD::TruthVertex*>& truthPUVertices);
 private:
     int m_iDetailLevel;
     float m_cutMinTruthRecoRadialDiff = 0.1;
     ///truth type
     TH1* m_vx_type_truth;
+    TH1* m_vx_z_diff;
+    TH1* m_vx_z_diff_pull;
     ///hardscatter classification
     TH1* m_vx_hs_classification;
     ///vertex reco efficiency

InnerDetector/InDetValidation/InDetPhysValMonitoring/src/InDetRttPlots.cxx

@@ -316,7 +316,8 @@ InDetRttPlots::fill(const xAOD::VertexContainer& vertexContainer, const std::vec
     ATH_MSG_DEBUG("IN InDetRttPlots::fill, filling for all vertices");
     if (vtx->vertexType() == xAOD::VxType::PriVtx) {
       m_hardScatterVertexPlots.fill(*vtx);
-      m_hardScatterVertexTruthMatchingPlots.fill(*vtx);
+      if(truthHSVertices.size()>0)m_hardScatterVertexTruthMatchingPlots.fill(*vtx,truthHSVertices[0]);
+      else m_hardScatterVertexTruthMatchingPlots.fill(*vtx);
       ATH_MSG_DEBUG("IN InDetRttPlots::fill, filling for all HS vertex");
     }
   }