diff --git a/NPLib/Detectors/PISTA/TFPMWPhysics.cxx b/NPLib/Detectors/PISTA/TFPMWPhysics.cxx
index f088f42dc21ca7d78efdb2657b35e94107889ec7..1e707a51227d1d1a2aab6244b6456b51f9e2abf1 100644
--- a/NPLib/Detectors/PISTA/TFPMWPhysics.cxx
+++ b/NPLib/Detectors/PISTA/TFPMWPhysics.cxx
@@ -671,6 +671,13 @@ void TFPMWPhysics::ReadAnalysisConfig() {
 }
 
 
+///////////////////////////////////////////////////////////////////////////
+TVector3 TFPMWPhysics::GetHitPosition(int Det){
+  TVector3 Hit = TVector3(PositionX[Det], PositionY[Det], GetDetectorPositionZ(Det));
+
+  return Hit;
+}
+
 ///////////////////////////////////////////////////////////////////////////
 void TFPMWPhysics::Clear() {
   DetectorNbr.clear();
diff --git a/NPLib/Detectors/PISTA/TFPMWPhysics.h b/NPLib/Detectors/PISTA/TFPMWPhysics.h
index 959ebba21235d0097d7d337299042caa34771f07..ea0a43663bca222e75bdfe4ee159c57ec6a9f396 100644
--- a/NPLib/Detectors/PISTA/TFPMWPhysics.h
+++ b/NPLib/Detectors/PISTA/TFPMWPhysics.h
@@ -170,6 +170,7 @@ class TFPMWPhysics : public TObject, public NPL::VDetector {
     double GetDetectorPositionY(int Det) {return DetPosY[Det];}
     double GetDetectorPositionZ(int Det) {return DetPosZ[Det];}
 
+    TVector3 GetHitPosition(int Det);
 
     //////////////////////////////////////////////////////////////
     // specific methods to FPMW array