PndSttHit.h

Go to the documentation of this file.

#ifndef PNDSTTHIT_H
#define PNDSTTHIT_H 1

#include "TVector3.h"
#include "FairHit.h"
#include "PndTrack.h"

class PndSttHit : public FairHit {

 public:
  PndSttHit();

  // THIS ONE!
  PndSttHit(Int_t detID, Int_t tubeID, Int_t mcindex, TVector3 &pos, TVector3 &dpos, Double_t p, Double_t isochrone, Double_t isochroneError, Double_t chDep);

  virtual ~PndSttHit();

  virtual void Print(const Option_t *opt = nullptr) const
  {
    std::cout << " opt = " << opt << std::endl;
    return;
  }

  void Clear();

  Double_t GetIsochrone() const { return fIsochrone; };
  Double_t GetIsochroneError() const { return fIsochroneError; };
  Double_t GetPulse() const { return fPulse; };
  Double_t GetDepCharge() const { return fDepCharge; };
  Double_t GetEnergyLoss() const { return fDepCharge / 1e6; };

  void SetIsochrone(Double_t isochrone) { fIsochrone = isochrone; };
  void SetIsochroneError(Double_t isochroneError) { fIsochroneError = isochroneError; };
  void SetDepCharge(Double_t depcharge) { fDepCharge = depcharge; }

  // tube ID // CHECK added
  void SetTubeID(Int_t tubeid) { fTubeID = tubeid; }
  Int_t GetTubeID() const { return fTubeID; }

  // computation of dE/dx
  Double_t ComputedEdx(PndTrack *track, Double_t tuberadius);

  virtual bool equal(FairTimeStamp *data)
  {
    PndSttHit *myDigi = dynamic_cast<PndSttHit *>(data);
    if (myDigi != nullptr) {
      if (fTubeID == myDigi->GetTubeID())
        return true;
    }
    return false;
  }

  virtual bool operator<(const PndSttHit &myDigi) const
  {
    if (fTubeID < myDigi.GetTubeID())
      return true;
    else
      return false;
  }

  friend std::ostream &operator<<(std::ostream &out, PndSttHit &digi)
  {
    out << "PndSttHit in Tube: " << digi.GetTubeID() << " Isochrone: " << digi.GetIsochrone() << " +/- " << digi.GetIsochroneError() << " Charge: " << digi.GetDepCharge()
        << " Pulse: " << digi.GetPulse() << std::endl;
    return out;
  }

 protected:
  Int_t fTubeID; // CHECK added
  Double_t fPulse;
  Double_t fIsochrone;
  Double_t fIsochroneError;
  Double_t fDepCharge;

  ClassDef(PndSttHit, 1);
};

#endif