mtd_digitizer Namespace Reference

Namespaces
	MTDHelpers

Classes
struct	MTDCellId
struct	MTDCellInfo
class	MTDDigitizer

Typedefs
typedef float	MTDSimData_t
typedef std::array< MTDSimData_t, nSamples >	MTDSimHitData
typedef std::unordered_map< MTDCellId, MTDCellInfo >	MTDSimHitDataAccumulator

Functions
void	loadSimHitAccumulator (MTDSimHitDataAccumulator &simData, const PMTDSimAccumulator &simAccumulator, const float minCharge, const float maxCharge)
void	saveSimHitAccumulator (PMTDSimAccumulator &simResult, const MTDSimHitDataAccumulator &simData, const float minCharge, const float maxCharge)

Variables
constexpr size_t	nSamples = 15
constexpr float	PREMIX_MAX_TOF = 26.0f

Typedef Documentation

typedef float mtd_digitizer::MTDSimData_t

Definition at line 13 of file MTDDigitizerTypes.h.

typedef std::array<MTDSimData_t, nSamples> mtd_digitizer::MTDSimHitData

Definition at line 15 of file MTDDigitizerTypes.h.

typedef std::unordered_map<MTDCellId, MTDCellInfo> mtd_digitizer::MTDSimHitDataAccumulator

Definition at line 38 of file MTDDigitizerTypes.h.

Function Documentation

void mtd_digitizer::loadSimHitAccumulator ( MTDSimHitDataAccumulator &  simData, const PMTDSimAccumulator &  simAccumulator, const float  minCharge, const float  maxCharge  )

inline

Definition at line 93 of file MTDDigitizer.h.

References newFWLiteAna::base, constexpr, dqm-mbProfile::log, PREMIX_MAX_TOF, PMTDSimAccumulator::Data::sampleOffset, logintpack::unpack16log(), and relativeConstraints::value.

Referenced by HGCDigitizer::accumulate(), and mtd_digitizer::MTDDigitizer< Traits >::accumulate().

                                                           {
    const float minPackChargeLog = minCharge > 0.f ? std::log(minCharge) : -2;
    const float maxPackChargeLog = std::log(maxCharge);
    constexpr uint16_t base = 1 << PMTDSimAccumulator::Data::sampleOffset;

    for (const auto& detIdIndexHitInfo : simAccumulator) {
      auto foo = simData.emplace(
          MTDCellId(detIdIndexHitInfo.detId(), detIdIndexHitInfo.row(), detIdIndexHitInfo.column()), MTDCellInfo());
      auto simIt = foo.first;
      auto& hit_info = simIt->second.hit_info;

      size_t iEn = detIdIndexHitInfo.energyIndex();
      size_t iSample = detIdIndexHitInfo.sampleIndex();

      float value;
      if (iEn == 1) {
        value = static_cast<float>(detIdIndexHitInfo.data()) / base * PREMIX_MAX_TOF;
      } else {
        value = logintpack::unpack16log(detIdIndexHitInfo.data(), minPackChargeLog, maxPackChargeLog, base);
      }

      if (iEn == 0) {
        hit_info[iEn][iSample] += value;
      } else if (hit_info[iEn][iSample] == 0) {
        // For iEn==1 the digitizers just set the TOF of the first SimHit
        hit_info[iEn][iSample] = value;
      }
    }
  }

void mtd_digitizer::saveSimHitAccumulator ( PMTDSimAccumulator &  simResult, const MTDSimHitDataAccumulator &  simData, const float  minCharge, const float  maxCharge  )

inline

Definition at line 57 of file MTDDigitizer.h.

References newFWLiteAna::base, constexpr, HTMLExport::elem(), PMTDSimAccumulator::emplace_back(), PMTDSimAccumulator::Data::energyMask, dqm-mbProfile::log, nSamples, logintpack::pack16log(), PREMIX_MAX_TOF, PMTDSimAccumulator::reserve(), PMTDSimAccumulator::Data::sampleMask, PMTDSimAccumulator::Data::sampleOffset, EgammaValidation_cff::samples, and relativeConstraints::value.

Referenced by HGCDigitizer::finalizeEvent(), and mtd_digitizer::MTDDigitizer< Traits >::finalizeEvent().

                                                           {
    constexpr auto nEnergies = std::tuple_size<decltype(MTDCellInfo().hit_info)>::value;
    static_assert(nEnergies <= PMTDSimAccumulator::Data::energyMask + 1,
                  "PMTDSimAccumulator bit pattern needs to updated");
    static_assert(nSamples <= PMTDSimAccumulator::Data::sampleMask + 1,
                  "PMTDSimAccumulator bit pattern needs to updated");

    const float minPackChargeLog = minCharge > 0.f ? std::log(minCharge) : -2;
    const float maxPackChargeLog = std::log(maxCharge);
    constexpr uint16_t base = 1 << PMTDSimAccumulator::Data::sampleOffset;

    simResult.reserve(simData.size());
    // mimicing the digitization
    for (const auto& elem : simData) {
      // store only non-zero
      for (size_t iEn = 0; iEn < nEnergies; ++iEn) {
        const auto& samples = elem.second.hit_info[iEn];
        for (size_t iSample = 0; iSample < nSamples; ++iSample) {
          if (samples[iSample] > minCharge) {
            unsigned short packed;
            if (iEn == 1) {
              // assuming linear range for tof of 0..26
              packed = samples[iSample] / PREMIX_MAX_TOF * base;
            } else {
              packed = logintpack::pack16log(samples[iSample], minPackChargeLog, maxPackChargeLog, base);
            }
            simResult.emplace_back(elem.first.detid_, elem.first.row_, elem.first.column_, iEn, iSample, packed);
          }
        }
      }
    }
  }

Variable Documentation

constexpr size_t mtd_digitizer::nSamples = 15

Definition at line 11 of file MTDDigitizerTypes.h.

Referenced by saveSimHitAccumulator().

constexpr float mtd_digitizer::PREMIX_MAX_TOF = 26.0f

Definition at line 23 of file MTDDigitizerTypes.h.

Referenced by loadSimHitAccumulator(), and saveSimHitAccumulator().