CbmTofPoint.cxx

Go to the documentation of this file.

/* Copyright (C) 2006-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Volker Friese, Denis Bertini [committer] */
 
#include "CbmTofPoint.h"
 
#include <FairMCPoint.h>  // for FairMCPoint
 
#include <TVector3.h>  // for TVector3
 
#include <bitset>   // for bitset
#include <cassert>  // for assert
#include <limits>   // for numeric_limits, numeric_limits<>::digits
#include <sstream>  // for operator<<, basic_ostream, stringstream
#include <string>   // for char_traits
 
using std::endl;
using std::string;
using std::stringstream;
 
 
// -----   Default constructor   -------------------------------------------
CbmTofPoint::CbmTofPoint() : FairMCPoint(), fNofCells(0), fGapMask(0) {}
// -------------------------------------------------------------------------
 
 
// -----   Standard constructor   ------------------------------------------
CbmTofPoint::CbmTofPoint(int32_t trackID, int32_t detID, TVector3 pos, TVector3 mom, double tof, double length,
                         double eLoss)
  : FairMCPoint(trackID, detID, pos, mom, tof, length, eLoss)
  , fNofCells(0)
  , fGapMask(0)
{
}
// -------------------------------------------------------------------------
 
 
// -----   Destructor   ----------------------------------------------------
CbmTofPoint::~CbmTofPoint() {}
// -------------------------------------------------------------------------
 
 
// -----   Get the number of gaps   ----------------------------------------
int32_t CbmTofPoint::GetNGaps() const
{
  int32_t iNGaps(0);
 
  for (int32_t iGapBit = 0; iGapBit < std::numeric_limits<uint16_t>::digits; iGapBit++) {
    if (fGapMask & (0x1 << iGapBit)) { iNGaps++; }
  }
 
  return iNGaps;
}
// -------------------------------------------------------------------------
 
 
// -----   Get the index of the first gap   --------------------------------
int32_t CbmTofPoint::GetFirstGap() const
{
  for (int32_t iGapBit = 0; iGapBit < std::numeric_limits<uint16_t>::digits; iGapBit++) {
    if (fGapMask & (0x1 << iGapBit)) { return iGapBit; }
  }
 
  return -1;
}
// -------------------------------------------------------------------------
 
 
// -----   Get the index of the last gap   ---------------------------------
int32_t CbmTofPoint::GetLastGap() const
{
  int32_t iLastGap(-1);
 
  for (int32_t iGapBit = 0; iGapBit < std::numeric_limits<uint16_t>::digits; iGapBit++) {
    if (fGapMask & (0x1 << iGapBit)) { iLastGap = iGapBit; }
  }
 
  return iLastGap;
}
// -------------------------------------------------------------------------
 
 
// -----   Add one gap to the gap mask   -----------------------------------
void CbmTofPoint::SetGap(int32_t iGap)
{
  assert(0 <= iGap && std::numeric_limits<uint16_t>::digits > iGap);
  fGapMask |= 0x1 << iGap;
}
// -------------------------------------------------------------------------
 
 
// -----   String output   -------------------------------------------------
string CbmTofPoint::ToString() const
{
  stringstream ss;
  ss << "STofPoint: track ID " << fTrackID << ", detector ID " << fDetectorID << "\n";
  ss << "    Position (" << fX << ", " << fY << ", " << fZ << ") cm \n";
  ss << "    Momentum (" << fPx << ", " << fPy << ", " << fPz << ") GeV \n";
  ss << "    Time " << fTime << " ns,  Length " << fLength << " cm,  Energy loss " << fELoss * 1.0e06 << " keV \n";
  ss << "    Number of cells " << fNofCells << ", gap mask "
     << std::bitset<std::numeric_limits<uint16_t>::digits>(fGapMask) << endl;
  return ss.str();
}
// -------------------------------------------------------------------------
 
ClassImp(CbmTofPoint)