CbmStsPoint.cxx

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                      CbmStsPoint source file                  -----
// -----                  Created 26/07/04  by V. Friese               -----
// -------------------------------------------------------------------------
 
#include "CbmStsPoint.h"
 
#include "BmnDetectorList.h"
 
#include <iostream>
 
using std::cout;
using std::endl;
using std::flush;
 
// -----   Default constructor   -------------------------------------------
CbmStsPoint::CbmStsPoint()
    : FairMCPoint()
    , fX_out(0.)
    , fY_out(0.)
    , fZ_out(0.)
    , fPx_out(0.)
    , fPy_out(0.)
    , fPz_out(0.)
    , fStation(-1)
    , fModule(-1)
{}
// -------------------------------------------------------------------------
 
// -----   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,
                         Int_t eventId)
    : FairMCPoint(trackID, detID, posIn, momIn, tof, length, eLoss, eventId)
    , fX_out(posOut.X())
    , fY_out(posOut.Y())
    , fZ_out(posOut.Z())
    , fPx_out(momOut.Px())
    , fPy_out(momOut.Py())
    , fPz_out(momOut.Pz())
    , fStation(-1)
    , fModule(-1)
{
    SetLink(FairLink(kMCTrack, trackID));
}
// -------------------------------------------------------------------------
 
// -----   Destructor   ----------------------------------------------------
CbmStsPoint::~CbmStsPoint() {}
// -------------------------------------------------------------------------
 
// -----   Copy constructor with event and epoch time   --------------------
CbmStsPoint::CbmStsPoint(const CbmStsPoint& point, Int_t eventId, Double_t eventTime, Double_t epochTime)
    : FairMCPoint(point)
    , fX_out(point.fX_out)
    , fY_out(point.fY_out)
    , fZ_out(point.fZ_out)
    , fPx_out(point.fPx_out)
    , fPy_out(point.fPy_out)
    , fPz_out(point.fPz_out)
    , fStation(point.fStation)
    , fModule(point.fModule)
{
    //  *this = point;
    if (eventId > 0)
        fEventId = eventId;
    fTime = point.GetTime() + eventTime - epochTime;
}
// -------------------------------------------------------------------------
 
// -----   Public method Print   -------------------------------------------
void CbmStsPoint::Info() const
{
    LOG(debug2) << "StsPoint: track ID " << fTrackID << ", detector ID " << fDetectorID;
    LOG(debug2) << "          IN  Position (" << fX << ", " << fY << ", " << fZ << ") cm";
    LOG(debug2) << "          OUT Position (" << fX_out << ", " << fY_out << ", " << fZ_out << ") cm";
    LOG(debug2) << "    Momentum (" << fPx << ", " << fPy << ", " << fPz << ") GeV";
    LOG(debug2) << "    Time " << fTime << " ns,  Length " << fLength << " cm,  Energy loss " << fELoss * 1.0e06
                << " keV";
    LOG(debug2) << "    StationNum " << fStation << ", ModuleNum " << fModule;
}
// -------------------------------------------------------------------------
 
// -----   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;
}
// -------------------------------------------------------------------------