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 int	kInTimeBX = 9
constexpr int	kNumberOfBX = 15
constexpr size_t	nSamples = 15
constexpr float	PREMIX_MAX_TOF = 25.0f

Typedef Documentation

MTDSimData_t

typedef float mtd_digitizer::MTDSimData_t

Definition at line 13 of file MTDDigitizerTypes.h.

MTDSimHitData

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

Definition at line 15 of file MTDDigitizerTypes.h.

MTDSimHitDataAccumulator

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

Definition at line 39 of file MTDDigitizerTypes.h.

Function Documentation

loadSimHitAccumulator()

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

inline

Definition at line 92 of file MTDDigitizer.h.

References newFWLiteAna::base, PMTDSimAccumulator::Data::dataMask, PMTDSimAccumulator::Data::energyMask, Exception, dqm-mbProfile::log, PREMIX_MAX_TOF, PMTDSimAccumulator::Data::sampleMask, logintpack::unpack16log(), and relativeConstraints::value.

Referenced by 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 = PMTDSimAccumulator::Data::dataMask;

    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();

      if (iEn > PMTDSimAccumulator::Data::energyMask + 1 || iSample > PMTDSimAccumulator::Data::sampleMask)
        throw cms::Exception("MTDDigitixer::loadSimHitAccumulator")
            << "Index out of range: iEn = " << iEn << " iSample = " << iSample << std::endl;

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

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

saveSimHitAccumulator()

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

inline

Definition at line 56 of file MTDDigitizer.h.

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

Referenced by 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 be updated");
    static_assert(nSamples == PMTDSimAccumulator::Data::sampleMask,
                  "PMTDSimAccumulator bit pattern needs to be updated");

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

    simResult.reserve(simData.size());
    // mimicking 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 || iEn == 3) {
              // assuming linear range for tof of 0..25
              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

kInTimeBX

constexpr int mtd_digitizer::kInTimeBX = 9

Definition at line 42 of file MTDDigitizerTypes.h.

Referenced by BTLDeviceSim::getHitsResponse(), and BTLElectronicsSim::run().

kNumberOfBX

constexpr int mtd_digitizer::kNumberOfBX = 15

Definition at line 41 of file MTDDigitizerTypes.h.

Referenced by BTLDeviceSim::getHitsResponse().

nSamples

constexpr size_t mtd_digitizer::nSamples = 15

Definition at line 11 of file MTDDigitizerTypes.h.

Referenced by saveSimHitAccumulator().

PREMIX_MAX_TOF

constexpr float mtd_digitizer::PREMIX_MAX_TOF = 25.0f

Definition at line 24 of file MTDDigitizerTypes.h.

Referenced by loadSimHitAccumulator(), and saveSimHitAccumulator().