Updating for position reconstruction

libPISTA/inc/Pista.h

@@ -39,6 +39,9 @@ public:
     inline void SetTS(ULong64_t& TS){Pista_TS = TS;};
     inline ULong64_t GetTS(void){return Pista_TS;};
 
+    void ReadGeometry(string filename);
+    void CreateHistogramsCal2D(TDirectory* Dir); 
+
 private:
     UShort_t NStrips_DE;
     UShort_t NStrips_E;
@@ -56,6 +59,7 @@ private:
     // TimeStamps of Mesytec Crate
     ULong64_t Pista_TS;
 
+   
 
 #ifdef WITH_NPTOOL
     TPISTAData* fPISTAData;

libPISTA/inc/PistaTelescope.h

@@ -3,6 +3,8 @@
 
 #include "BaseDetector.hh"
 #include "MesytecParameters.hh"
+#include "TVector3.h"
+
 
 class PistaTelescope : public BaseDetector
 {
@@ -24,6 +26,21 @@ public:
 
     //! SetParameters
     void SetParameters(Parameters* Par,Map* Map);
+    void AddGeometry(TVector3 A, TVector3 B, TVector3 C, TVector3 D);
+    
+    double GetStripPositionX(int X, int Y){
+      return m_StripPositionX[X][Y];
+    };
+    double GetStripPositionY(int X, int Y){
+      return m_StripPositionY[X][Y];
+    };
+    double GetStripPositionZ(int X, int Y){
+      return m_StripPositionZ[X][Y];
+    };
+
+    TVector3 CalculatePositionOfInteraction(int stripE, int stripDE);
+    TVector3 GetDetectorNormal();
+    vector<TVector3> GetPositionOfInteraction() {return m_PositionOfInteraction;};
 
 #ifdef WITH_ROOT
     void SetOpt(TTree *OutTTree, TTree *InTTree);
@@ -32,6 +49,16 @@ public:
 private:
     UShort_t NStrips_DE;
     UShort_t NStrips_E;
+    
+     // for position reconstruction
+    TVector3 m_A; //!
+    TVector3 m_B; //!
+    TVector3 m_C; //!
+    TVector3 m_D; //!
+    vector<vector<double>> m_StripPositionX; //!
+    vector<vector<double>> m_StripPositionY; //!
+    vector<vector<double>> m_StripPositionZ; //!
+    vector<TVector3> m_PositionOfInteraction; //!
 };
 
 #endif // PISTATELESCOPE_H