Digit.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
#ifndef ALICEO2_PHOS_DIGIT_H_
#define ALICEO2_PHOS_DIGIT_H_
 
#include <cmath>
#include "CommonDataFormat/TimeStamp.h"
 
namespace o2
{
 
namespace phos
{
class Hit;
 
using DigitBase = o2::dataformats::TimeStamp<double>;
class Digit : public DigitBase
{
 
 public:
  static constexpr int kTimeGate = 25; // Time in ns between digits to be added as one signal.
                                       // Should it be readout time (6000 ns???): to be tested
 
  Digit() = default;
 
  Digit(short cell, float amplitude, float time, int label);
 
  Digit(short cell, float amplitude, float time, bool isTrigger2x2, int label);
 
  Digit(const Hit& hit, int label);
 
  ~Digit() = default; // override
 
  void fillFromHit(const Hit& hit);
 
  inline bool operator<(const Digit& other) const
  {
    if (fabs(getTimeStamp() - other.getTimeStamp()) < kTimeGate) {
      return getAbsId() < other.getAbsId();
    } else {
      return getTimeStamp() < other.getTimeStamp();
    }
  }
 
  inline bool operator>(const Digit& other) const
  {
    if (fabs(getTimeStamp() - other.getTimeStamp()) <= kTimeGate) {
      return getAbsId() > other.getAbsId();
    } else {
      return getTimeStamp() > other.getTimeStamp();
    }
  }
 
  inline bool operator==(const Digit& other) const
  {
    return ((fabs(getTimeStamp() - other.getTimeStamp()) <= kTimeGate) &&
            getAbsId() == other.getAbsId());
  }
 
  bool canAdd(const Digit other) const;
  Digit& operator+=(const Digit& other); //
 
  void addEnergyTime(float energy, float time);
 
  // true if tru and not readount digit
  bool isTRU() const { return mAbsId >= NREADOUTCHANNELS; }
 
  short getAbsId() const { return mAbsId; }
  void setAbsId(short cellId) { mAbsId = cellId; }
 
  short getTRUId() const { return mAbsId - NREADOUTCHANNELS; }
  void setTRUId(short cellId) { mAbsId = cellId + NREADOUTCHANNELS; }
 
  float getAmplitude() const { return mAmplitude; }
  void setAmplitude(float amplitude) { mAmplitude = amplitude; }
 
  float getTime() const { return mTime; }
  void setTime(float time) { mTime = time; }
 
  bool isHighGain() const { return mIsHighGain; }
  void setHighGain(Bool_t isHG) { mIsHighGain = isHG; }
 
  bool is2x2Tile() const { return isTRU() && isHighGain(); }
 
  int getLabel() const { return mLabel; }
  void setLabel(int l) { mLabel = l; }
 
  void reset()
  {
    mIsHighGain = true;
    mAbsId = 0;
    mLabel = -1;
    mAmplitude = 0;
    mTime = 0;
  }
 
  void PrintStream(std::ostream& stream) const;
 
 private:
  static constexpr short NREADOUTCHANNELS = 14337; 
 
  bool mIsHighGain = true; 
  short mAbsId = 0;        
  int mLabel = -1;         
  float mAmplitude = 0;    
  float mTime = 0.;        
 
  ClassDefNV(Digit, 1);
};
 
std::ostream& operator<<(std::ostream& stream, const Digit& dig);
} // namespace phos
} // namespace o2
#endif