MuonPRDTest - Add TGC measurement uncertainty to the tree

MuonPRDTest - Add TGC measurement uncertainty to the tree

MuonSpectrometer/MuonValidation/MuonPRDTest/MuonPRDTest/TGCPRDVariables.h

@@ -28,6 +28,7 @@ namespace MuonPRDTest{
         TwoVectorBranch m_TGC_PRD_localPos{parent(), "PRD_TGC_localPos"};
         TgcIdentifierBranch m_TGC_PRD_id{parent(), "PRD_TGC"};
         VectorBranch<uint8_t>& m_TGC_PRD_bcId{parent().newVector<uint8_t>("PRD_TGC_bcId")};
+        VectorBranch<float>& m_TGC_PRD_cov{parent().newVector<float>("PRD_TGC_cov")};
     };
 };
 

MuonSpectrometer/MuonValidation/MuonPRDTest/Root/TGCPRDVariables.cxx

 /*
-  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
+  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
 */
 
 #include "MuonPRDTest/TGCPRDVariables.h"
 
 #include "MuonReadoutGeometry/TgcReadoutElement.h"
-
+#include "EventPrimitives/EventPrimitivesHelpers.h"
 namespace MuonPRDTest {
     TGCPRDVariables::TGCPRDVariables(MuonTesterTree& tree, const std::string& container_name, MSG::Level msglvl) :
         PrdTesterModule(tree, "PRD_TGC", true, msglvl), m_key{container_name} {}
@@ -41,6 +41,7 @@ namespace MuonPRDTest {
                 m_TGC_PRD_globalPos.push_back(pos);
                 m_TGC_PRD_localPos.push_back(loc_pos);                
                 m_TGC_PRD_bcId.push_back(prd->getBcBitMap());
+                m_TGC_PRD_cov.push_back(Amg::error(prd->localCovariance(), Trk::locX));
                 ++n_PRD;
             }
         }

MuonSpectrometer/MuonValidation/MuonPRDTest/Root/TGCSimHitVariables.cxx

@@ -55,7 +55,7 @@ namespace MuonPRDTest {
 
             const Amg::Vector3D& localPosition = hit.localPosition();
             m_TGC_hitLocalPosition.push_back(localPosition);
-            m_TGC_hitGlobalPosition.push_back(tgcdet->localToGlobalCoords(localPosition, offid));
+            m_TGC_hitGlobalPosition.push_back(tgcdet->localToGlobalTransf(offid) *localPosition);
             m_TGC_detector_globalPosition.push_back(tgcdet->globalPosition());
 
             m_TGC_particleEncoding.push_back(hit.particleEncoding());