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

                                                            {
     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;
       }
     }
   }

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

◆ saveSimHitAccumulator()

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

inline

Definition at line 57 of file MTDDigitizer.h.

                                                            {
     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);
           }
         }
       }
     }
   }