ECLCrystalData.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
namespace Belle2 {
  class ECLCrystalData {
  private :
 
    //Crystal identifiers
    int m_cell_ID;
    int m_phi_ID;
    int m_theta_ID;
 
    //indexes (whatever those mean...)
    int m_phi_idx;
    int m_theta_idx;
 
    //Crystal properties
    float m_mass;
    float m_volume;
    float EvalMass();
 
    //Crystal Location
    float m_PosR;
    float m_PosP;
    float m_PosX;
    float m_PosY;
    float m_PosZ;
 
    double EvalPhi();
    double EvalR();
    double EvalZ();
    double EvalX();
    double EvalY();
 
    void Mapping(int cid);
 
  public :
 
 
    static const int Ring[69];
    static const int sumPrevious[24];
 
    const int nECLCrystalTot = 8736; 
    const int nECLCrystalBAR = 6624; 
    const int nECLCrystalECF = 1152; 
    const int nECLCrystalECB =  960; 
    static const float CrystalMassEndcap[132];
    static const float CrystalMassBarrel[46];
 
 
    static const double theta[69];
 
    static const double EndcapRadius[132];
 
    static const double BarrelZ[46];
 
    ECLCrystalData();
    explicit ECLCrystalData(int cid);
    ECLCrystalData(int tid, int pid);
    ~ECLCrystalData();
    void Eval();
    int GetCrystalIndex();
    bool IsEndCap();
    bool IsBarrel();
    int GetCellID();
    int GetCellID(int tid, int pid);
 
    double GetPhi() { return m_PosP; };
    inline double GetTheta() { return theta[m_theta_ID]; } ;
    double GetR() { return m_PosR; };
    double GetZ() { return m_PosZ; };
    double GetX() { return m_PosX; };
    double GetY() { return m_PosY; };
    double GetMass() { return m_mass; };
    inline int GetThetaID() { return m_theta_ID; } ;
    inline int GetPhiID() { return m_phi_ID; } ;
    inline int GetNperThetaID() { return 16 * Ring[m_theta_ID];};
 
 
  };
}