ERGadastCsIPoint.cxx

/********************************************************************************
 *              Copyright (C) Joint Institute for Nuclear Research              *
 *                                                                              *
 *              This software is distributed under the terms of the             *
 *         GNU Lesser General Public Licence version 3 (LGPL) version 3,        *
 *                  copied verbatim in the file "LICENSE"                       *
 ********************************************************************************/

#include "ERGadastCsIPoint.h"

#include "FairLogger.h"

//--------------------------------------------------------------------------------------------------
ERGadastCsIPoint::ERGadastCsIPoint():
  FairMCPoint(),
  fX_out(0.), fY_out(0.), fZ_out(0.),
  fPx_out(0.), fPy_out(0.), fPz_out(0.) 
{
}

//--------------------------------------------------------------------------------------------------
ERGadastCsIPoint::ERGadastCsIPoint(Int_t eventID, Int_t trackID,
    Int_t mot0trackID,
    Double_t mass,
    TVector3 posIn,
    TVector3 posOut, TVector3 momIn, TVector3 momOut,
    Double_t tof, Double_t length, Double_t eLoss, Int_t pdg,
    Int_t wall, Int_t block, Int_t cell):
  FairMCPoint(trackID, -1., posIn, momIn, tof, length, eLoss),
  fEventID(eventID),
  fX_out(posOut.X()), fY_out(posOut.Y()), fZ_out(posOut.Z()),
  fPx_out(momOut.X()), fPy_out(momOut.Y()), fPz_out(momOut.Z()),
  fPDG(pdg), fWall(wall), fBlock(block), fCell(cell)
{
}

//--------------------------------------------------------------------------------------------------
ERGadastCsIPoint::ERGadastCsIPoint(const ERGadastCsIPoint& right):
  FairMCPoint(right),
  fX_out(right.fX_out), fY_out(right.fY_out), fZ_out(right.fZ_out),
  fPx_out(right.fPx_out), fPy_out(right.fPy_out), fPz_out(right.fPz_out),
  fPDG(right.fPDG), fWall(right.fWall), fBlock(right.fBlock), fCell(right.fCell)
{
}

//--------------------------------------------------------------------------------------------------
ERGadastCsIPoint::~ERGadastCsIPoint() {
}

//--------------------------------------------------------------------------------------------------
void ERGadastCsIPoint::Print(const Option_t* opt /* = 0*/) const {
  LOG(INFO) << "ERGadastCsIPoint: track " << fTrackID;
  LOG(INFO) << " mother track = " << fMot0TrackID << FairLogger::endl;
  LOG(INFO) << "    Position (" << fX << ", " << fY << ", " << fZ << ") cm" << FairLogger::endl;
  LOG(INFO) << "    Momentum (" << fPx << ", " << fPy << ", " << fPz << ") GeV" << FairLogger::endl;
  LOG(INFO) << "    Time " << fTime << " ns,  Length " << fLength << " cm" << FairLogger::endl;
  LOG(INFO) << "    Energy loss " << fELoss << " keV "<< FairLogger::endl;
}

//--------------------------------------------------------------------------------------------------
Double_t ERGadastCsIPoint::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) );
}

//--------------------------------------------------------------------------------------------------
Double_t ERGadastCsIPoint::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) );
}

//--------------------------------------------------------------------------------------------------
Bool_t ERGadastCsIPoint::IsUsable() const {
  Double_t dz = fZ_out - fZ;
  if ( TMath::Abs(dz) < 1.e-4 ) return kFALSE;
  return kTRUE;
}

//--------------------------------------------------------------------------------------------------
ClassImp(ERGadastCsIPoint)