add DCA, mc info and vtx index info to StMuFwdTrack (#733)

to mirror what is in StEvent structure (StFwdTrack)

Author: jdbrice

StRoot/StMuDSTMaker/COMMON/StMuFwdTrack.cxx

@@ -39,6 +39,12 @@ void StMuFwdTrack::set( StFwdTrack * evTrack) {
             StMuFwdTrackSeedPoint( TVector3( sp.mXYZ.x(), sp.mXYZ.y(), sp.mXYZ.z() ), sp.mSector, sp.mTrackId, sp.mCov )
         );
     }
+
+    setDCA(evTrack->dca().x(), evTrack->dca().y(), evTrack->dca().z());
+    mIdTruth = evTrack->idTruth();
+    mQATruth = evTrack->qaTruth();
+    mVtxIndex = evTrack->vertexIndex();
+
 }
 
 

StRoot/StMuDSTMaker/COMMON/StMuFwdTrack.h

@@ -15,6 +15,7 @@
 #include "TRefArray.h"
 
 #include "StMuFcsCluster.h"
+#include <climits>
 
 
 class StFwdTrack;
@@ -152,6 +153,18 @@ class StMuFwdTrack : public TObject {
 
     TRefArray mEcalClusters;
     TRefArray mHcalClusters;
+
+    void setMc( UShort_t idt, UShort_t qual ) { mIdTruth = idt; mQATruth = qual; }
+    void setDCA( TVector3 dca ) { mDCA[0] = dca.X(); mDCA[1] = dca.Y(); mDCA[2] = dca.Z(); }
+    void setDCA( float dcaX, float dcaY, float dcaZ ) { mDCA[0] = dcaX; mDCA[1] = dcaY; mDCA[2] = dcaZ; }
+    void setVtxIndex( UChar_t vtxIndex ) { mVtxIndex = vtxIndex; }
+
+
+    UShort_t idTruth() const { return mIdTruth; }
+    UShort_t qaTruth() const { return mQATruth; }
+    TVector3 dca() const { return TVector3( mDCA[0], mDCA[1], mDCA[2] ); }
+    UChar_t vertexIndex() const { return mVtxIndex; }
+    bool isPrimary() const { return mVtxIndex != UCHAR_MAX; }
 
 protected:
 
@@ -166,12 +179,14 @@ class StMuFwdTrack : public TObject {
     float mPval;
     short mCharge;
     TVector3 mPrimaryMomentum;
-    
 
-    // StPtrVecFcsCluster mEcalClusters;
-    // StPtrVecFcsCluster mHcalClusters;
+    UShort_t mIdTruth; // MC track id
+    UShort_t mQATruth; // MC track quality (percentage of hits coming from corresponding MC track)
+
+    float mDCA[3]; // DCA to the primary vertex
+    UChar_t mVtxIndex;
 
-    ClassDef(StMuFwdTrack,2)
+    ClassDef(StMuFwdTrack,3)
 
 };