PndEmcHit.h

Go to the documentation of this file.

// -------------------------------------------------------------------------
// -----                     PndEmcHit header file                     -----
// -----               Created 14/08/06  by S.Spataro                  -----
// -------------------------------------------------------------------------

//#pragma once
#ifndef PNDEMCHIT_H
#define PNDEMCHIT_H

#include "FairHit.h"
#include "TMath.h"
#include <vector>
#include <map>
#include <algorithm>

class PndEmcPoint;

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

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

  PndEmcHit(Int_t trackid, Int_t id, Float_t energy, Float_t time, Float_t X, Float_t Y, Float_t Z, std::vector<Int_t> McTruth, FairMultiLinkedData enteringTrack,
            FairMultiLinkedData exitingTrack, std::set<Int_t> clusterList);

  PndEmcHit(const PndEmcHit &);

  virtual ~PndEmcHit();

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

  void SetEnergy(Double32_t energy) { fEnergy = energy; };
  void SetTime(Double32_t time) { fTimeStamp = time * 1E9; };
  void SetDepositedEnergyMap(const std::map<Int_t, Double_t> &depEnergy) { fDepositedEnergy = depEnergy; }
  void AddDepositedEnergy(Int_t mcId, Double_t depEnergy) { fDepositedEnergy[mcId] += depEnergy; }

  Double_t GetEnergy() const { return fEnergy; };
  Double_t GetTime() const { return fTimeStamp / 1E9; };
  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;
  FairMultiLinkedData GetTrackEntering() const { return fTrackEntering; }
  FairMultiLinkedData GetTrackExiting() const { return fTrackExiting; }

  const std::vector<Int_t> &GetMcList() { return fMcList; }
  std::set<Int_t> GetClusterList() { return fClusterList; }

  void AddClusterID(int clusterId) { fClusterList.insert(clusterId); }

  std::map<Int_t, Double_t> GetDepositedEnergyMap() const { return fDepositedEnergy; };

 protected:
  Double32_t fEnergy; // hit amplitude

  std::vector<Int_t> fMcList;                 //  Mc TrackIndex contributed to hit
  FairMultiLinkedData fTrackEntering;         // Links to tracks entering the crystal
  FairMultiLinkedData fTrackExiting;          // Links to tracks exiting the crystal
  std::set<Int_t> fClusterList;               // To which cluster belongs the hit based on MC information (could be more than one)
  std::map<Int_t, Double_t> fDepositedEnergy; // MC deposited energy per (primary) MC particle

  ClassDef(PndEmcHit, 4)
};

#endif // PNDEMCHIT_H