StMuDSTMaker/COMMON/StMuMcTrack.h

// $Id: StMuMcTrack.h,v 1.3 2012/05/07 14:47:06 fisyak Exp $
#ifndef __StMuMcTrack_h__
#define __StMuMcTrack_h__
#include "tables/St_g2t_track_Table.h" 
#include "StThreeVectorF.hh"
class StMuMcTrack : public TObject {
 public:
  enum EHIT {ktpc, ksvt, kssd,
             kctb, keem, kemc, kesm, kftp, kgem, khpd, kist, kigt, kfst, 
             kfgt, kfpd, kmwc, kpgc, kpmd, ksmd, kpix, ktof, kvpd, ktot};
  StMuMcTrack(const g2t_track_st &t) : TObject(), mGePid(t.ge_pid), mId(t.id), mIsShower(t.is_shower), mItrmdVertex(t.itrmd_vertex_p),
    mIdVx(t.start_vertex_p), mIdVxEnd(t.stop_vertex_p), mCharge(t.charge), mE(t.e), mEta(t. eta), mPxyz(t.p), mpT(t.pt), mPtot(t.ptot), 
    mRapidity(t.rapidity) {
    mHits[kctb] = 0xff & t.n_ctb_hit;  /* Nhits in ctb */
    mHits[keem] = 0xff & t.n_eem_hit;  /* Nhits in eem (endcap em cal) */
    mHits[kemc] = 0xff & t.n_emc_hit;  /* Nhits in emc */
    mHits[kesm] = 0xff & t.n_esm_hit;  /* Nhits in esm (endcap shower max) */
    mHits[kftp] = 0xff & t.n_ftp_hit;  /* Nhits in forward tpc */
    mHits[kgem] = 0xff & t.n_gem_hit;  /* Nhits in gem barrel */
    mHits[khpd] = 0xff & t.n_hpd_hit;  /* Nhits in hpd */
    mHits[kist] = 0xff & t.n_ist_hit;  /* Nhits in ist */
    mHits[kigt] = 0xff & t.n_igt_hit;  /* Nhits in igt */
    mHits[kfst] = 0xff & t.n_fst_hit;  /* Nhits in fst */
    mHits[kfgt] = 0xff & t.n_fgt_hit;  /* Nhits in fgt */
    mHits[kfpd] = 0xff & t.n_fpd_hit;  /* Nhits in fpd */
    mHits[kmwc] = 0xff & t.n_mwc_hit;  /* Nhits in mwc */
    mHits[kpgc] = 0xff & t.n_pgc_hit;  /* Nhits in pgc  ???  */
    mHits[kpmd] = 0xff & t.n_pmd_hit;  /* Nhits in pmd (PMD) */
    mHits[ksmd] = 0xff & t.n_smd_hit;  /* number of hits in shower max */
    mHits[kssd] = 0xff & t.n_ssd_hit;  /* Nhits in ssd */
    mHits[ksvt] = 0xff & t.n_svt_hit;  /* Nhits in svt */
    mHits[kpix] = 0xff & t.n_pix_hit;  /* Nhits in pix */
    mHits[ktof] = 0xff & t.n_tof_hit;  /* Nhits in tof */
    mHits[ktpc] = 0xff & t.n_tpc_hit;  /* Nhits in tpc */
    mHits[kvpd] = 0xff & t.n_vpd_hit;  /* Nhits in vpd */
  }
  StMuMcTrack() {}
  virtual ~StMuMcTrack() {}
  Int_t                 GePid()        const {return mGePid;} /* GEANT particle id */           
  Int_t                 Id()           const {return mId;}    /* primary key */                
  Bool_t                IsShower()     const {return mIsShower;} /* 1 if shower track, 0 if not */
  Int_t                 NoHits()       const {Int_t n = 0; for (Int_t i = ktpc; i < ktot; i++) n+= NoHits(i); return n;}
  UChar_t               NoHits(Int_t k)const {return mHits[k];}                                    
  UChar_t               No_ctb_hit()   const {return NoHits(kctb);}   /* Nhits in ctb */                           
  UChar_t               No_eem_hit()   const {return NoHits(keem);}   /* Nhits in eem (endcap em cal) */           
  UChar_t               No_emc_hit()   const {return NoHits(kemc);}   /* Nhits in emc */                           
  UChar_t               No_esm_hit()   const {return NoHits(kesm);}   /* Nhits in esm (endcap shower max) */ 
  UChar_t               No_ftp_hit()   const {return NoHits(kftp);}   /* Nhits in forward tpc */                   
  UChar_t               No_gem_hit()   const {return NoHits(kgem);}   /* Nhits in gem barrel */            
  UChar_t               No_hpd_hit()   const {return NoHits(khpd);}   /* Nhits in hpd */                           
  UChar_t               No_ist_hit()   const {return NoHits(kist);}   /* Nhits in ist */                           
  UChar_t               No_igt_hit()   const {return NoHits(kigt);}   /* Nhits in igt */                           
  UChar_t               No_fst_hit()   const {return NoHits(kfst);}   /* Nhits in fst */                           
  UChar_t               No_fgt_hit()   const {return NoHits(kfgt);}   /* Nhits in fgt */                           
  UChar_t               No_fpd_hit()   const {return NoHits(kfpd);}   /* Nhits in fpd */                           
  UChar_t               No_mwc_hit()   const {return NoHits(kmwc);}   /* Nhits in mwc */                           
  UChar_t               No_pgc_hit()   const {return NoHits(kpgc);}   /* Nhits in pgc  ???  */             
  UChar_t               No_pmd_hit()   const {return NoHits(kpmd);}   /* Nhits in pmd (PMD) */             
  UChar_t               No_smd_hit()   const {return NoHits(ksmd);}   /* number of hits in shower max */           
  UChar_t               No_ssd_hit()   const {return NoHits(kssd);}   /* Nhits in ssd */                           
  UChar_t               No_svt_hit()   const {return NoHits(ksvt);}   /* Nhits in svt */                           
  UChar_t               No_pix_hit()   const {return NoHits(kpix);}   /* Nhits in pix */                           
  UChar_t               No_tof_hit()   const {return NoHits(ktof);}   /* Nhits in tof */                           
  UChar_t               No_tpc_hit()   const {return NoHits(ktpc);}   /* Nhits in tpc */                           
  UChar_t               No_vpd_hit()   const {return NoHits(kvpd);}   /* Nhits in vpd */                     
  Int_t                 ItrmdVertex()  const {return mItrmdVertex;} /* First intermediate vertex */        
  Int_t                 IdVx       ()  const {return mIdVx;       } /* Id of start vertex of track */      
  Int_t                 IdVxEnd    ()  const {return mIdVxEnd;    } /* Id of stop vertex of this track */ 
  Char_t                Charge     ()  const {return mCharge;     } /* Charge */                           
  Float_t               E          ()  const {return mE;          } /* Energy */                           
  Float_t               Eta        ()  const {return mEta;        } /* Pseudorapidity */                  
  const StThreeVectorF& Pxyz       ()  const {return *&mPxyz;     } /* Momentum */                        
  Float_t               pT         ()  const {return mpT;         } /* Transverse momentum */              
  Float_t               Ptot       ()  const {return mPtot;       } /* Total momentum */                   
  Float_t               Rapidity   ()  const {return mRapidity;   } /* Rapidity */                        
  virtual void          Print(Option_t* option = "") const;  
  static Int_t          CorrectGePid(Int_t gePid);
  const Char_t         *GeName();
 private:
  Int_t          mGePid;       /* GEANT particle id */
  Int_t          mId;          /* primary key */
  Bool_t         mIsShower;    /* 1 if shower track, 0 if not */
  UChar_t        mHits[ktot];  /* Nhits in a detector */
  Int_t          mItrmdVertex; /* First intermediate vertex */
  Int_t          mIdVx;        /* Id of start vertex of track */
  Int_t          mIdVxEnd;     /* Id of stop vertex of this track */
  Char_t         mCharge;      /* Charge */
  Float_t        mE;           /* Energy */
  Float_t        mEta;         /* Pseudorapidity */
  StThreeVectorF mPxyz;        /* Momentum */
  Float_t        mpT;          /* Transverse momentum */
  Float_t        mPtot;        /* Total momentum */
  Float_t        mRapidity;    /* Rapidity */
  ClassDef(StMuMcTrack,1)
};
ostream&              operator<<(ostream& os, StMuMcTrack const & v);
#endif
// $Log: StMuMcTrack.h,v $
// Revision 1.3  2012/05/07 14:47:06  fisyak
// Add handles for track to fast detector matching
//
// Revision 1.2  2011/10/17 00:19:14  fisyak
// Active handing of IdTruth
//

---