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 90 of file MTDDigitizer.h.

References runEdmFileComparison::base, constexpr, cmsBatch::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 56 of file MTDDigitizer.h.

References runEdmFileComparison::base, constexpr, HTMLExport::elem(), PMTDSimAccumulator::emplace_back(), PMTDSimAccumulator::Data::energyMask, cmsBatch::log, nSamples, logintpack::pack16log(), PREMIX_MAX_TOF, PMTDSimAccumulator::reserve(), PMTDSimAccumulator::Data::sampleMask, PMTDSimAccumulator::Data::sampleOffset, hcalTTPDigis_cfi::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().