File indexing completed on 2024-04-06 12:05:20
0001 #include "DataFormats/TrackerRecHit2D/interface/SiPixelRecHitQuality.h"
0002
0003 #include "FWCore/MessageLogger/interface/MessageLogger.h"
0004
0005 SiPixelRecHitQuality::Packing::Packing() {
0006
0007
0008
0009 probX_width = 14;
0010 probY_width = 8;
0011 qBin_width = 3;
0012 edge_width = 1;
0013 bad_width = 1;
0014 twoROC_width = 1;
0015 hasFilledProb_width = 1;
0016 spare_width = 3;
0017
0018 if (probX_width + probY_width + qBin_width + edge_width + bad_width + twoROC_width + hasFilledProb_width +
0019 spare_width !=
0020 32) {
0021 throw cms::Exception("SiPixelRecHitQuality::Packing: ")
0022 << "\nERROR: The allocated bits for the quality word to not sum to 32."
0023 << "\n\n";
0024 }
0025
0026 probX_units = 1.0018;
0027 probY_units = 1.0461;
0028 probX_1_over_log_units = 1.0 / log(probX_units);
0029 probY_1_over_log_units = 1.0 / log(probY_units);
0030
0031
0032 probX_shift = 0;
0033 probY_shift = probX_shift + probX_width;
0034 qBin_shift = probY_shift + probY_width;
0035 edge_shift = qBin_shift + qBin_width;
0036 bad_shift = edge_shift + edge_width;
0037 twoROC_shift = bad_shift + bad_width;
0038 hasFilledProb_shift = twoROC_shift + twoROC_width;
0039
0040
0041
0042 QualWordType zero32 = 0;
0043
0044 probX_mask = ~(~zero32 << probX_width);
0045 probY_mask = ~(~zero32 << probY_width);
0046 qBin_mask = ~(~zero32 << qBin_width);
0047 edge_mask = ~(~zero32 << edge_width);
0048 bad_mask = ~(~zero32 << bad_width);
0049 twoROC_mask = ~(~zero32 << twoROC_width);
0050 hasFilledProb_mask = ~(~zero32 << hasFilledProb_width);
0051 }
0052
0053
0054 const SiPixelRecHitQuality::Packing SiPixelRecHitQuality::thePacking;
0055
0056 void SiPixelRecHitQuality::warningObsolete() {
0057 edm::LogWarning("ObsoleteVariable")
0058 << "Since 39x, probabilityX and probabilityY have been replaced by probabilityXY and probabilityQ";
0059 }
0060
0061 void SiPixelRecHitQuality::warningOutOfBoundQbin(int iValue, QualWordType const& iQualWord) {
0062 edm::LogWarning("OutOfBounds") << "Qbin outside the bounds of the quality word: Qbin = " << iValue
0063 << " and QualityWord = " << iQualWord << " --> Now defaulting to Qbin = 0.0";
0064 }
0065
0066 void SiPixelRecHitQuality::warningOutOfBoundProb(const char* iName, float iProb, QualWordType const& iQualWord) {
0067 edm::LogWarning("OutOfBounds") << "Prob" << iName << " outside the bounds of the quality word: Prob" << iName
0068 << " = 1+" << iProb - 1 << " and QualityWord = " << iQualWord
0069 << " --> Now defaulting to Prob" << iName << " = 0.0";
0070 }
0071
0072 void SiPixelRecHitQuality::warningOutOfBoundRaw(const char* iName, int iRaw, QualWordType const& iQualWord) {
0073 edm::LogWarning("OutOfBounds") << "Probability" << iName << " outside the bounds of the quality word: Raw = " << iRaw
0074 << " and QualityWord = " << iQualWord << " --> Now defaulting to Prob" << iName
0075 << " = 0.0";
0076 }