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;
}
// -------------------------------------------------------------------------