 |
Belle II Software
release-05-02-19
|
16 #include <mdst/dataobjects/ECLCluster.h>
17 #include <framework/logging/Logger.h>
25 B2ERROR(
"This cluster does not support the requested hypothesis,"
28 <<
LogVar(
" You requested the hypothesis bitmask: ", std::bitset<16>(
static_cast<unsigned short>(hypothesis)).to_string()));
29 return std::numeric_limits<double>::quiet_NaN();
37 B2ERROR(
"EHypothesisBit is not supported yet: " <<
static_cast<unsigned short>(hypothesis));
38 return std::numeric_limits<double>::quiet_NaN();
47 return TVector3(cluster_x, cluster_y, cluster_z);
52 TMatrixDSym covmatecl(3);
61 for (
int i = 0; i < 3; i++)
62 for (
int j = 0; j < i ; j++)
63 covmatecl(j, i) = covmatecl(i, j);
72 if (theta < 0.2164208)
return 0;
73 if (theta < 0.5480334)
return 1;
74 if (theta < 0.561996)
return 11;
75 if (theta < 2.2462387)
return 2;
76 if (theta < 2.2811453)
return 13;
77 if (theta < 2.7070057)
return 3;
85 #pragma GCC diagnostic push
86 #pragma GCC diagnostic ignored "-Wdeprecated-declarations"
88 #pragma GCC diagnostic pop
90 return 100000 * crid + 1000 * hypoid + showerid;
TVector3 getClusterPosition() const
Return TVector3 on cluster position (x,y,z)
EHypothesisBit
The hypothesis bits for this ECLCluster (Connected region (CR) is split using this hypothesis.
double getTheta() const
Return Corrected Theta of Shower (radian).
double getPhi() const
Return Corrected Phi of Shower (radian).
@ c_nPhotons
CR is split into n photons (N1)
int getHypothesisId() const
Return hypothesis id.
double getR() const
Return R.
int getUniqueId() const
Return (pseudo) unique Id based on CRId, ShowerId and HypothesisID.
bool hasHypothesis(EHypothesisBit bitmask) const
Return if specific hypothesis bit is set.
double getEnergy(const EHypothesisBit &hypothesis) const
Return Energy (GeV).
Abstract base class for different kinds of events.
Class to store variables with their name which were sent to the logging service.
Double32_t m_sqrtcovmat_22
Covariance matrix 22, sigma_theta, between 0 and 50 mrad.
@ c_neutralHadron
CR is reconstructed as a neutral hadron (N2)
Double32_t m_sqrtcovmat_00
Covariance entry 00 sigma_E (1% to 25% between 10 MeV and 8 GeV)
Double32_t m_covmat_21
Covariance matrix 21, not used yet.
Double32_t m_logEnergy
< Log.
Double32_t m_covmat_20
Covariance matrix 20, not used yet.
TMatrixDSym getCovarianceMatrix3x3() const
Return TMatrixDsym 3x3 covariance matrix for E, Phi and Theta.
int getDetectorRegion() const
Return detector region: 0: below acceptance, 1: FWD, 2: BRL, 3: BWD, 11: FWDGAP, 13: BWDGAP.
Double32_t m_covmat_10
Covariance matrix 10, not used yet.
int getConnectedRegionId() const
Return connected region id.
Double32_t m_sqrtcovmat_11
Covariance matrix 11, sigma_phi, between 0 and 50 mrad.
int getClusterId() const
Return cluster id.
Double32_t m_logEnergyRaw
Log.