StjeTrackListToStMuTrackFourVecList.cxx

// $Id: StjeTrackListToStMuTrackFourVecList.cxx,v 1.10 2012/03/23 05:44:20 pibero Exp $
#include "StjeTrackListToStMuTrackFourVecList.h"

#include "StjTrackToTLorentzVector.h"

#include "../emulator/StMuTrackEmu.h"
#include "../emulator/StMuTrackFourVec.h"

StjeTrackListToStMuTrackFourVecList::StjeTrackListToStMuTrackFourVecList() 
  : _trackTo4p(*(new StjTrackToTLorentzVector)) { }


FourList StjeTrackListToStMuTrackFourVecList::operator()(const StjTrackList& trackList)
{
  FourList ret;

  for(StjTrackList::const_iterator track = trackList.begin(); track != trackList.end(); ++track) {

    StMuTrackEmu *trackEmu = createTrackEmu(*track);

    TLorentzVector p4 = _trackTo4p(*track);

    StMuTrackFourVec* pmu = new StMuTrackFourVec(trackEmu, 0, 0, p4, (*track).charge, (*track).trackIndex, kTpcId);
    ret.push_back(pmu);
  }

  return ret;
}

StMuTrackEmu* StjeTrackListToStMuTrackFourVecList::createTrackEmu(const StjTrack& track)
{
  StMuTrackEmu *ret = new StMuTrackEmu();

  TVector3 momentum;
  momentum.SetPtEtaPhi(track.pt, track.eta, track.phi);

  ret->_px             =  momentum.Px()               ;
  ret->_py             =  momentum.Py()               ;
  ret->_pz             =  momentum.Pz()               ;
  ret->_flag           =  track.flag                  ;
  ret->_nHits          =  track.nHits                 ;
  ret->_charge         =  track.charge                ;
  ret->_nHitsPoss      =  track.nHitsPoss             ;
  ret->_nHitsDedx      =  track.nHitsDedx             ;
  ret->_nHitsFit       =  track.nHitsFit              ;
  ret->_nSigmaPion     =  track.nSigmaPion            ;
  ret->_nSigmaKaon     =  track.nSigmaKaon            ;
  ret->_nSigmaProton   =  track.nSigmaProton          ;
  ret->_nSigmaElectron =  track.nSigmaElectron        ;
  ret->_Tdca           =  track.Tdca                  ;
  ret->_dcaX           =  track.dcaX                  ;
  ret->_dcaY           =  track.dcaY                  ;
  ret->_dcaZ           =  track.dcaZ                  ;
  ret->_dcaD           =  track.dcaD                  ;
  ret->_chi2           =  track.chi2                  ;
  ret->_chi2prob       =  track.chi2prob              ;
  ret->_BField         =  track.BField                ;
  ret->_bemcRadius     =  track.bemcRadius            ;
  ret->_etaext         =  track.exitEta               ;
  ret->_phiext         =  track.exitPhi               ;
  ret->_exitTowerId    =  track.exitTowerId           ;
  ret->_exitDetectorId =  track.exitDetectorId        ;
  ret->_dEdx           =  track.dEdx                  ;
  ret->_beta           =  track.beta                  ;
  ret->_firstPoint     =  track.firstPoint            ;
  ret->_lastPoint      =  track.lastPoint             ;
  ret->_trackIndex     =  track.trackIndex            ; 
  ret->_id             =  track.id                    ;
  ret->_idTruth        =  track.idTruth               ;
  ret->_qaTruth        =  track.qaTruth               ;
  ret->_detectorId     =  track.detectorId            ;
  ret->_btofTrayId         =  track.btofTrayId        ; 
  ret->_nSigmaTofPion      =  track.nSigmaTofPion     ;
  ret->_nSigmaTofKaon      =  track.nSigmaTofKaon     ;
  ret->_nSigmaTofProton    =  track.nSigmaTofProton   ;
  ret->_nSigmaTofElectron  =  track.nSigmaTofElectron ;

  return ret;
}