PndEmcApdHit.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                  PndEmcApdHit header file                     -----
// -----               Created 07/04/08  by S.Spataro                  -----
// -------------------------------------------------------------------------

//#pragma once
#ifndef PNDEMCAPDHIT_H
#define PNDEMCAPDHIT_H

#include "FairHit.h"
#include "TMath.h"
#include "TArrayI.h"
#include <vector>

class PndEmcApdPoint;

class PndEmcApdHit : public FairHit {
 public:
  PndEmcApdHit();

  PndEmcApdHit(Int_t trackid, Int_t id, Float_t energy, Float_t time, Float_t X, Float_t Y, Float_t Z);

  PndEmcApdHit(Int_t trackid, Int_t id, Float_t energy, Float_t time, Float_t X, Float_t Y, Float_t Z, Int_t npoint, Int_t pointIndex[10]);

  PndEmcApdHit(Int_t trackid, Int_t id, Float_t energy, Float_t time, Float_t X, Float_t Y, Float_t Z, std::vector<PndEmcApdPoint *> PointList);

  PndEmcApdHit(const PndEmcApdHit &);

  virtual ~PndEmcApdHit();

  virtual void Print(const Option_t *opt = "") const;

  virtual void SetEnergy(Double32_t energy) { fEnergy = energy; };
  virtual void SetTime(Double32_t time) { fTime = time; };

  virtual Double_t GetEnergy() const { return fEnergy; };
  virtual Double_t GetTime() const { return fTime; };
  Float_t GetTheta() const { return fX == 0.0 && fY == 0.0 && fZ == 0.0 ? 0.0 : TMath::ATan2(sqrt(fX * fX + fY * fY), fZ) * TMath::RadToDeg(); };
  Float_t GetPhi() const { return fX == 0.0 && fY == 0.0 ? 0.0 : TMath::ATan2(fY, fX) * TMath::RadToDeg(); };
  Short_t GetModule() const { return (fDetectorID / 100000000); };
  Short_t GetRow() const { return ((fDetectorID / 1000000) % 100); };
  Short_t GetCrystal() const { return (fDetectorID % 10000); };
  Short_t GetCopy() const { return ((fDetectorID / 10000) % 100); };

  Short_t GetXPad() const;
  Short_t GetYPad() const;

  Int_t GetNPoints() const { return fNPoints; };
  Int_t GetMCIndex(Short_t ind) const { return fPointIndex[ind]; };

  // Possible memory issue (FIXME).
  std::vector<PndEmcApdPoint *> &GetPointList() { return fPointList; }

  // TArrayI* GetPoindIndexes() { return fPointIndexes;};

 protected:
  Double32_t fTime;      // time
  Double32_t fEnergy;    // hit amplitude
  Int_t fPointIndex[10]; // index of ther first 10 particles entering

  // TArrayI *fPointIndexes;

  std::vector<PndEmcApdPoint *> fPointList; //  points contributed to hit
  Int_t fNPoints;                           // number of ApdPoints

  ClassDef(PndEmcApdHit, 1)
};

#endif // PNDEMCAPDHIT_H