CbmMvdDigi.cxx

Go to the documentation of this file.

/* Copyright (C) 2008-2020 GSI Helmholtzzentrum fuer Schwerionenforschung, Darmstadt
   SPDX-License-Identifier: GPL-3.0-only
   Authors: Christina Dritsa [committer], Florian Uhlig, Philipp Sitzmann */
 
// -----------------------------------------------------------------------
// -----               CbmMvdDigi source file                        -----
// -----              Created 17/04/08  by C. Dritsa                 -----
// -----------------------------------------------------------------------
 
#include "CbmMvdDigi.h"
 
#include <Logger.h>  // for LOG, Logger
 
#include <TObject.h>  // for TObject
 
#include <cmath>
#include <limits>
 
// -----   Default constructor   -------------------------------------------
CbmMvdDigi::CbmMvdDigi()
  : TObject()
  , fDigiTime(0.)
  , fCharge(0.)
  , fPixelSizeX(0.)
  , fPixelSizeY(0.)
  , fChannelNr(0)
  , fFrameNumber(0)
  , fRefId(-1)
  , fAddress(0)
  , fDigiFlag(-1)
{
}
// -------------------------------------------------------------------------
 
 
/*
// -----   Constructor with parameters   -----------------------------------
CbmMvdDigi::CbmMvdDigi(int32_t iStation, int32_t iChannelNrX, int32_t iChannelNrY, float charge,
                       float pixelSizeX, float pixelSizeY)
  : CbmDigi(kMVD, 0),
    CbmMvdDetectorId(),
    fCharge(charge),
    fChannelNrX(iChannelNrX),
    fChannelNrY(iChannelNrY),
    fTrackID(-1),
    fPointID(0),
    fPixelSizeX(pixelSizeX),
    fPixelSizeY(pixelSizeY),
    fDetectorId(0),
    fChannelNr(0),
    fDigiTime(0.),
    fFrameNumber(0),
    fRefId(-1),
    fDigiFlag(-1)
{
    // Check range for station
    if ( ! ( iStation >= 0 && iStation <= 255 ) ) {
        LOG(fatal) << "Illegal station number " << iStation;
    }
 
    fDetectorId = DetectorId(iStation);
 
    fChannelNrY = iChannelNrY;
    fCharge  = charge;
    fChannelNrX=iChannelNrX;
    fChannelNrY=iChannelNrY;
    fPixelSizeX=pixelSizeX;
    fPixelSizeY=pixelSizeY;
    fDigiFlag=-1;
 
}
// -------------------------------------------------------------------------
*/
 
// -----   Constructor with parameters  --> used only due to error TODO include correct version -----------------------------------
CbmMvdDigi::CbmMvdDigi(int32_t sensorNr, int32_t iChannelNrX, int32_t iChannelNrY, float charge, float pixelSizeX,
                       float pixelSizeY, double time, int32_t frame)
  : TObject()
  , fDigiTime(time)
  , fCharge(charge)
  , fPixelSizeX(pixelSizeX)
  , fPixelSizeY(pixelSizeY)
  , fChannelNr(0)
  , fFrameNumber(frame)
  , fRefId(-1)
  , fAddress(CbmMvdAddress::GetAddressFromSensorNrAndPixelXY(sensorNr, iChannelNrY, iChannelNrX))
  , fDigiFlag(-1)
{
  // Check range for sensorNr
  if (!(sensorNr >= 0 && sensorNr <= 600)) {
    LOG(fatal) << "Illegal sensor number " << sensorNr;
  }
}
// -------------------------------------------------------------------------
 
// -----   Conversion Constructor from  CbmMvdRawDigi  ---------------------
CbmMvdDigi::CbmMvdDigi(const CbmMvdRawDigi& rhs, float pixelSizeX, float pixelSizeY)
  : TObject()
  , fCharge(10000)  // Based on email by M. Deveaux, value high enough to be above threshold
  , fPixelSizeX(pixelSizeX)
  , fPixelSizeY(pixelSizeY)
  , fChannelNr(0)  // Deprecated
  , fRefId(-1)     // set only by consumers (ClusterFinder)
  , fDigiFlag(-1)
{
  fAddress     = rhs.GetAddress();
  fFrameNumber = rhs.GetFrame();
 
  if (std::numeric_limits<int32_t>::max() < rhs.GetFrame()) {
    LOG(fatal) << "Frame number in Raw Digi (uint32) too large for MvdDigi field (int32): " << rhs.GetFrame()
               << " (max allowed: " << std::numeric_limits<int32_t>::max() << ")";
  }
  fDigiTime = rhs.GetTime();  // !!!! Conversion from uint32_t to double, range/precision to check? !!!!
}
// -------------------------------------------------------------------------
 
// -----   Conversion Constructor to  CbmMvdRawDigi  ---------------------
CbmMvdRawDigi CbmMvdDigi::ConvertToRawDigi(const CbmMvdDigi& rhs)
{
  if (rhs.GetTime() < 0) {
    LOG(fatal) << " Cannot convert CbmMvdDigi to CbmMvdRawDigi as time negative (conversion to uint64_t "
               << rhs.GetTime();
  }

  return CbmMvdRawDigi(rhs.GetAddress(),
                       static_cast<uint32_t>(rhs.GetTime()),  // Convert double in ns to uint32_t rounded to the ns
                       rhs.GetFrameNumber());
}
// -------------------------------------------------------------------------
 
int32_t CbmMvdDigi::GetAdcCharge(int32_t adcDynamic, int32_t adcOffset, int32_t adcBits) const
{
 
  int32_t adcCharge;
 
  if (fCharge < adcOffset) { return 0; };
 
 
  double stepSize;
  //    int32_t adcMax = adcOffset + adcDynamic;
 
  stepSize  = adcDynamic / pow(2, adcBits);
  adcCharge = int((fCharge - adcOffset) / stepSize);
 
 
  if (adcCharge > int(pow(2, adcBits) - 1)) { adcCharge = (int) pow(2, adcBits) - 1; }
 
  if (gDebug > 0) { LOG(debug) << "CbmMvdDigi::GetAdcCharge() " << adcCharge; }
 
  return adcCharge;
}
 
 
// -------------------------------------------------------------------------
 
ClassImp(CbmMvdDigi)