PndEmcPoint.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

// ------------------------------------------------------------------------
// -----                     PndEmcPoint header file                  -----
// -----               Created 14/08/06  by S.Spataro                 -----
// ------------------------------------------------------------------------
//#pragma once
#ifndef PNDEMCPOINT_H
#define PNDEMCPOINT_H

#include "TObject.h"
#include "TVector3.h"
#include "TLorentzVector.h"
#include "FairMCPoint.h"
#include "PndDetectorList.h"

class PndEmcPoint : public FairMCPoint {

 public:
  PndEmcPoint();

  PndEmcPoint(Int_t trackID, Int_t detID, Int_t evtID, TVector3 pos, TVector3 mom, Double_t tof, Double_t length, Double_t eLoss, Short_t mod, Short_t row, Short_t crys,
              Short_t copy, Bool_t entering = kFALSE, Bool_t exiting = kFALSE, int clusteID = -1);

  PndEmcPoint(const PndEmcPoint &point);

  virtual ~PndEmcPoint();

  Double_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(); };
  Double_t GetPhi() const { return fX == 0.0 && fY == 0.0 ? 0.0 : TMath::ATan2(fY, fX) * TMath::RadToDeg(); };
  Short_t GetXPad() const;
  Short_t GetYPad() const;

  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); };
  Bool_t GetEntering() const { return fEntering; }
  Bool_t GetExiting() const { return fExiting; }
  Int_t GetClusterID() const { return fClusterID; }
  void SetModule(Short_t mod) { nModule = mod; };
  void SetRow(Short_t row) { nRow = row; };
  void SetCrystal(Short_t crys) { nCrystal = crys; };
  void SetTrackID(Int_t trackId)
  {
    FairMCPoint::SetTrackID(trackId);
    SetLink(FairLink("MCTrack", trackId));
  } // 14.09.10 Stefano FIX

  void SetClusterID(int clusterID) { fClusterID = clusterID; }

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

 protected:
  Short_t nModule;  // Module number
  Short_t nRow;     // Row number
  Short_t nCrystal; // Crystal number
  Short_t nCopy;    // Copy number
  Bool_t fEntering; //< Is particle entering into crystal
  Bool_t fExiting;  //< Is particle exiting the crystal
  Int_t fClusterID; //< MCParticle which created the cluster

  ClassDef(PndEmcPoint, 2)
};

#endif