// ------------------------------------------------------------------------- // ----- CbmStsPoint source file ----- // ----- Created 26/07/04 by V. Friese ----- // ------------------------------------------------------------------------- #include "CbmStsPoint.h" #include using std::cout; using std::endl; using std::flush; // ----- Default constructor ------------------------------------------- CbmStsPoint::CbmStsPoint() : FairMCPoint() { fX_out = fY_out = fZ_out = 0.; fPx_out = fPy_out = fPz_out = 0.; } // ------------------------------------------------------------------------- // ----- Standard constructor ------------------------------------------ CbmStsPoint::CbmStsPoint(Int_t trackID, Int_t detID, TVector3 posIn, TVector3 posOut, TVector3 momIn, TVector3 momOut, Double_t tof, Double_t length, Double_t eLoss) : FairMCPoint(trackID, detID, posIn, momIn, tof, length, eLoss) { fX_out = posOut.X(); fY_out = posOut.Y(); fZ_out = posOut.Z(); fPx_out = momOut.Px(); fPy_out = momOut.Py(); fPz_out = momOut.Pz(); } // ------------------------------------------------------------------------- // ----- Destructor ---------------------------------------------------- CbmStsPoint::~CbmStsPoint() { } // ------------------------------------------------------------------------- // ----- Public method Print ------------------------------------------- void CbmStsPoint::Print(const Option_t* opt) const { cout << "-I- CbmStsPoint: STS Point for track " << fTrackID << " in detector " << fDetectorID << endl; cout << " Position (" << fX << ", " << fY << ", " << fZ << ") cm" << endl; cout << " Momentum (" << fPx << ", " << fPy << ", " << fPz << ") GeV" << endl; cout << " Time " << fTime << " ns, Length " << fLength << " cm, Energy loss " << fELoss*1.0e06 << " keV" << endl; } // ------------------------------------------------------------------------- // ----- Point x coordinate from linear extrapolation ------------------ Double_t CbmStsPoint::GetX(Double_t z) const { // cout << fZ << " " << z << " " << fZ_out << endl; if ( (fZ_out-z)*(fZ-z) >= 0. ) return (fX_out+fX)/2.; Double_t dz = fZ_out - fZ; return ( fX + (z-fZ) / dz * (fX_out-fX) ); } // ------------------------------------------------------------------------- // ----- Point y coordinate from linear extrapolation ------------------ Double_t CbmStsPoint::GetY(Double_t z) const { if ( (fZ_out-z)*(fZ-z) >= 0. ) return (fY_out+fY)/2.; Double_t dz = fZ_out - fZ; // if ( TMath::Abs(dz) < 1.e-3 ) return (fY_out+fY)/2.; return ( fY + (z-fZ) / dz * (fY_out-fY) ); } // ------------------------------------------------------------------------- // ----- Public method IsUsable ---------------------------------------- Bool_t CbmStsPoint::IsUsable() const { Double_t dz = fZ_out - fZ; if ( TMath::Abs(dz) < 1.e-4 ) return kFALSE; return kTRUE; } // ------------------------------------------------------------------------- ClassImp(CbmStsPoint)