TTelescopeData.h

Go to the documentation of this file.

/**
 * @file    TTelescopeData.h
 * @brief
 * @author  Kodai Okawa <okawa@cns.s.u-tokyo.ac.jp>
 * @date    2024-01-17 17:11:50
 * @note    last modified: 2025-01-08 10:36:09
 * @details
 */
 
#ifndef _CRIB_TTELESCOPEDATA_H_
#define _CRIB_TTELESCOPEDATA_H_
 
#include "TArtTypes.h"
#include <TDataObject.h>
#include <TVector3.h>
#include <constant.h>
 
namespace art::crib {
class TTelescopeData;
} // namespace art::crib
 
class art::crib::TTelescopeData : public TDataObject {
  public:
    typedef enum { kID,
                   kTiming } ESortType;
    typedef enum { kASC,
                   kDESC } ESortOrder;
 
    TTelescopeData();
    ~TTelescopeData() override;
 
    TTelescopeData(const TTelescopeData &rhs);
    TTelescopeData &operator=(const TTelescopeData &rhs);
 
    TVector3 GetPosition() const { return fPos; }
    void SetPosition(TVector3 vec) { fPos = vec; }
    void SetPosition(Double_t x, Double_t y, Double_t z) { fPos.SetXYZ(x, y, z); }
 
    Int_t GetN() const { return fNE; }
    void SetN(Int_t arg) { fNE = arg; }
 
    Int_t GetTelID() const { return fTelID; }
    void SetTelID(Int_t arg) { fTelID = arg; }
    Int_t GetXID() const { return fXID; }
    void SetXID(Int_t arg) { fXID = arg; }
    Int_t GetYID() const { return fYID; }
    void SetYID(Int_t arg) { fYID = arg; }
 
    Double_t GetdE() const { return fdE; }
    void SetdE(Double_t arg) { fdE = arg; }
    Double_t GetdEX() const { return fdEX; }
    void SetdEX(Double_t arg) { fdEX = arg; }
    Double_t GetdEY() const { return fdEY; }
    void SetdEY(Double_t arg) { fdEY = arg; }
 
    Double_t GetE() const { return fE; }
    void SetE(Double_t arg) { fE = arg; }
    Double_t GetEtotal() const { return fEtotal; }
    void SetEtotal(Double_t arg) { fEtotal = arg; }
    Double_t GetTelXTiming() const { return fXTiming; }
    void SetTelXTiming(Double_t arg) { fXTiming = arg; }
    Double_t GetTelYTiming() const { return fYTiming; }
    void SetTelYTiming(Double_t arg) { fYTiming = arg; }
 
    Double_t GetTheta_L() const { return fTheta_L; }
    void SetTheta_L(Double_t arg) { fTheta_L = arg; }
 
    DoubleVec_t GetEnergyArray() const { return fEnergyArray; }
    Double_t GetEnergyArray(Int_t id) const { return fEnergyArray[id]; }
    void PushEnergyArray(Double_t arg) { fEnergyArray.emplace_back(arg); }
    DoubleVec_t GetTimingArray() const { return fTimingArray; }
    Double_t GetTimingArray(Int_t id) const { return fTimingArray[id]; }
    void PushTimingArray(Double_t arg) { fTimingArray.emplace_back(arg); }
 
    Double_t E() const { return fEtotal; } // get total energy
    Double_t E(Int_t id) const {
        if (fEnergyArray.size() == 0) {
            return kInvalidD;
        }
        if (id < 0 || id >= fNE) {
            return kInvalidD;
        }
        return fEnergyArray[id];
    } // get each layer energy
 
    Double_t T() const {
        if (fTimingArray.size() == 0) {
            return kInvalidD;
        }
        return fTimingArray[0];
    } // get timing
    Double_t T(Int_t id) const {
        if (fTimingArray.size() == 0) {
            return kInvalidD;
        }
        if (id < 0 || id >= fNE) {
            return kInvalidD;
        }
        return fTimingArray[id];
    } // get each layer timing
 
    Double_t A() const { return fTheta_L; } // get angle: A()-> Lab angle (deg)
    Double_t X() const { return fPos.X(); } // get hit X position (mm)
    Double_t Y() const { return fPos.Y(); } // get hit Y position (mm)
    Double_t Z() const { return fPos.Z(); } // get hit Z position (mm)
 
    void Copy(TObject &dest) const override;
    void Clear(Option_t *opt = "") override;
 
  protected:
    TVector3 fPos; // detected position (X, Y, Z)
 
    Int_t fTelID; // telescope ID
    Int_t fXID;   // X strip number
    Int_t fYID;   // Y strip number
    Int_t fNE;    // number of all SSDs
 
    Double_t fdE;      // energy at first layor
    Double_t fdEX;     // X side energy (=~ fdEY)
    Double_t fdEY;     // Y side energy (=~ fdEX)
    Double_t fE;       // added energy at thick SSDs
    Double_t fEtotal;  // all energy deposit in the telescope
    Double_t fXTiming; // timing information at the first layor (X side)
    Double_t fYTiming; // for case that X side have trouble (Y side)
 
    Double_t fTheta_L; // reaction angle in LAB system
 
    DoubleVec_t fEnergyArray; // energy array for each SSD
    DoubleVec_t fTimingArray; // timing array for each SSD
 
    ClassDefOverride(TTelescopeData, 3)
};
 
#endif // _TTELESCOPEDATA_H_