HGCalTriggerCellThresholdCodecImpl Class Reference

#include <HGCalTriggerCellThresholdCodecImpl.h>

Public Types
typedef HGCalTriggerCellThresholdDataPayload	data_type

Public Member Functions
uint32_t	adcnBits () const
double	adcsaturation () const
size_t	dataLength () const
data_type	decode (const std::vector< bool > &, const uint32_t, const HGCalTriggerGeometryBase &) const
std::vector< bool >	encode (const data_type &, const HGCalTriggerGeometryBase &) const
	HGCalTriggerCellThresholdCodecImpl (const edm::ParameterSet &conf)
void	linearize (const std::vector< HGCDataFrame< HGCalDetId, HGCSample >> &, std::vector< std::pair< HGCalDetId, uint32_t > > &)
double	linLSB () const
size_t	nCellsInModule () const
int	TCThreshold_ADC () const
double	TCThreshold_fC () const
uint32_t	tdcnBits () const
double	tdcOnsetfC () const
double	tdcsaturation () const
void	thresholdSelect (data_type &)
uint32_t	triggerCellSaturationBits () const
void	triggerCellSums (const HGCalTriggerGeometryBase &, const std::vector< std::pair< HGCalDetId, uint32_t > > &, data_type &)
uint32_t	triggerCellTruncationBits () const

Private Attributes
double	adcLSB_
uint32_t	adcnBits_
double	adcsaturation_
size_t	dataLength_
double	linLSB_
size_t	nCellsInModule_
int	TCThreshold_ADC_
double	TCThreshold_fC_
double	tdcLSB_
uint32_t	tdcnBits_
double	tdcOnsetfC_
double	tdcsaturation_
uint32_t	triggerCellSaturationBits_
uint32_t	triggerCellTruncationBits_

Detailed Description

Definition at line 27 of file HGCalTriggerCellThresholdCodecImpl.h.

Member Typedef Documentation

typedef HGCalTriggerCellThresholdDataPayload HGCalTriggerCellThresholdCodecImpl::data_type

Definition at line 30 of file HGCalTriggerCellThresholdCodecImpl.h.

Constructor & Destructor Documentation

HGCalTriggerCellThresholdCodecImpl::HGCalTriggerCellThresholdCodecImpl

(

const edm::ParameterSet &

conf

)

Definition at line 6 of file HGCalTriggerCellThresholdCodecImpl.cc.

References adcLSB_, adcnBits_, adcsaturation_, dataLength_, encode(), createfilelist::int, linLSB_, funct::pow(), TCThreshold_ADC_, TCThreshold_fC_, tdcLSB_, tdcnBits_, tdcsaturation_, triggerCellSaturationBits_, and triggerCellTruncationBits_.

                                                                :
    dataLength_(conf.getParameter<uint32_t>("DataLength")),
    nCellsInModule_(conf.getParameter<uint32_t>("MaxCellsInModule")), 
    linLSB_(conf.getParameter<double>("linLSB")),
    adcsaturation_(conf.getParameter<double>("adcsaturation")),
    adcnBits_(conf.getParameter<uint32_t>("adcnBits")),
    tdcsaturation_(conf.getParameter<double>("tdcsaturation")), 
    tdcnBits_(conf.getParameter<uint32_t>("tdcnBits")), 
    tdcOnsetfC_(conf.getParameter<double>("tdcOnsetfC")),
    triggerCellTruncationBits_(conf.getParameter<uint32_t>("triggerCellTruncationBits")),
    TCThreshold_fC_(conf.getParameter<double>("TCThreshold_fC"))
{
  adcLSB_ =  adcsaturation_/pow(2.,adcnBits_);
  tdcLSB_ =  tdcsaturation_/pow(2.,tdcnBits_);
  triggerCellSaturationBits_ = triggerCellTruncationBits_ + dataLength_;
  TCThreshold_ADC_ = (int) (TCThreshold_fC_ / linLSB_);
}

Member Function Documentation

uint32_t HGCalTriggerCellThresholdCodecImpl::adcnBits ( ) const

inline

Definition at line 50 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return adcnBits_;}

double HGCalTriggerCellThresholdCodecImpl::adcsaturation ( ) const

inline

Definition at line 49 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return adcsaturation_;}

size_t HGCalTriggerCellThresholdCodecImpl::dataLength ( ) const

inline

Definition at line 47 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return dataLength_;}

HGCalTriggerCellThresholdCodecImpl::data_type HGCalTriggerCellThresholdCodecImpl::decode	(	const std::vector< bool > &	data,
		const uint32_t	module,
		const HGCalTriggerGeometryBase &	geometry
	)		const

Definition at line 90 of file HGCalTriggerCellThresholdCodecImpl.cc.

References dataLength_, PV3DBase< T, PVType, FrameType >::eta(), Exception, HGCalTriggerGeometryBase::getOrderedTriggerCellsFromModule(), HGCalTriggerGeometryBase::getTriggerCellPosition(), mps_fire::i, diffTreeTool::index, linearize(), nCellsInModule_, p4, HGCalTriggerCellThresholdDataPayload::payload, PV3DBase< T, PVType, FrameType >::phi(), point, HGCalTriggerCellThresholdDataPayload::reset(), and mps_fire::result.

Referenced by HGCalTriggerCellThresholdCodec::decodeImpl(), and encode().

{
    data_type result;
    result.reset();
    // TODO: could eventually reserve result memory to the max size of trigger cells
   
    HGCalTriggerGeometryBase::geom_ordered_set trigger_cells_in_module = geometry.getOrderedTriggerCellsFromModule(module);
    size_t iselected = 0;
    size_t index = 0;
    for(const auto& triggercell : trigger_cells_in_module)
    {
        if(index>=nCellsInModule_)
        {
            throw cms::Exception("BadGeometry")
                << "Number of trigger cells in module too large for available data payload\n";
        }
        if(data[index])
        {
            uint32_t value = 0;
            for(size_t i=0;i<dataLength_;i++)
            {
                size_t ibit = nCellsInModule_+iselected*dataLength_+i; 
                if(data[ibit]) value |= (0x1<<i);
            }
            iselected++;
            // Build trigger cell
            if(value>0)
            {
                // Currently no hardware eta, phi and quality values
                result.payload.emplace_back(reco::LeafCandidate::LorentzVector(),
                        value, 0, 0, 0, triggercell); 
                GlobalPoint point = geometry.getTriggerCellPosition(triggercell);
                // 'value' is hardware, so p4 is meaningless, except for eta and phi
                math::PtEtaPhiMLorentzVector p4((double)value/cosh(point.eta()), point.eta(), point.phi(), 0.);
                result.payload.back().setP4(p4);
                result.payload.back().setPosition(point);
            }
        }
        index++;
    }
    return result;
}

std::vector< bool > HGCalTriggerCellThresholdCodecImpl::encode	(	const data_type &	data,
		const HGCalTriggerGeometryBase &	geometry
	)		const

Definition at line 28 of file HGCalTriggerCellThresholdCodecImpl.cc.

References dataLength_, decode(), Exception, HGCalTriggerGeometryBase::getModuleFromTriggerCell(), HGCalTriggerGeometryBase::getOrderedTriggerCellsFromModule(), mps_fire::i, diffTreeTool::index, nCellsInModule_, HGCalTriggerCellThresholdDataPayload::payload, mps_fire::result, findQualityFiles::size, triggerCellSaturationBits_, and triggerCellTruncationBits_.

Referenced by HGCalTriggerCellThresholdCodec::encodeImpl(), and HGCalTriggerCellThresholdCodecImpl().

{
    // First nCellsInModule_ bits are encoding the map of selected trigger cells
    // Followed by size words of dataLength_ bits, corresponding to energy/transverse energy of
    // the selected trigger cells
  
    // Convert payload into a map for later search
    std::unordered_map<uint32_t, uint32_t> data_map; // (detid,energy)
    size_t size=0;
    for(const auto& triggercell : data.payload)
      {
        data_map.emplace(triggercell.detId(), triggercell.hwPt());
        if (triggercell.hwPt()>0) size++;
      }
    std::vector<bool> result(nCellsInModule_ + dataLength_*size, 0);
    // No data: return vector of 0
    if(data.payload.size()==0) return result;
    // All trigger cells are in the same module
    // Loop on trigger cell ids in module and check if energy in the cell
    size_t index = 0; // index in module
    size_t idata = 0; // counter for the number of non-zero energy values
    // Retrieve once the ordered list of trigger cells in this module
    uint32_t module = geometry.getModuleFromTriggerCell(data.payload.begin()->detId());
    HGCalTriggerGeometryBase::geom_ordered_set trigger_cells_in_module = geometry.getOrderedTriggerCellsFromModule(module);
    for(const auto& triggercell_id : trigger_cells_in_module)
      {
        // Find if this trigger cell has data
        const auto& data_itr = data_map.find(triggercell_id);
        // if not data, increase index and skip
        if(data_itr==data_map.end())
        {
            index++;
            continue;
        }
        // else fill result vector with data
        // (set the corresponding adress bit and fill energy if >0)
        if(index>=nCellsInModule_) 
        {
          throw cms::Exception("BadGeometry")
            << "Number of trigger cells in module too large for available data payload\n";
        }
        uint32_t value = data_itr->second; 
        if(value>0)
          {
            // Set map bit to 1
            result[index] =  1;
            // Saturate and truncate energy values
            if(value+1>(0x1u<<triggerCellSaturationBits_)) value = (0x1<<triggerCellSaturationBits_)-1;
            for(size_t i=0; i<dataLength_; i++)
              {
                // remove the lowest bits (=triggerCellTruncationBits_)
                result[nCellsInModule_ + idata*dataLength_ + i] = static_cast<bool>(value & (0x1<<(i+triggerCellTruncationBits_)));
              }
            idata++;
          }
        index++;
      }
    return result;
}

void HGCalTriggerCellThresholdCodecImpl::linearize	(	const std::vector< HGCDataFrame< HGCalDetId, HGCSample >> &	dataframes,
		std::vector< std::pair< HGCalDetId, uint32_t > > &	linearized_dataframes
	)

Definition at line 136 of file HGCalTriggerCellThresholdCodecImpl.cc.

References adcLSB_, CustomPhysics_cfi::amplitude, data, linLSB_, ALCARECOPromptCalibProdSiPixelAli0T_cff::mode, tdcLSB_, tdcOnsetfC_, and triggerCellSums().

Referenced by decode(), and HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{
    double amplitude; uint32_t amplitude_int;
   

    for(const auto& frame : dataframes) {//loop on DIGI
        if (frame[2].mode()) {//TOT mode
            amplitude =( floor(tdcOnsetfC_/adcLSB_) + 1.0 )* adcLSB_ + double(frame[2].data()) * tdcLSB_;
        }
        else {//ADC mode
            amplitude = double(frame[2].data()) * adcLSB_;
        }

        amplitude_int = uint32_t (floor(amplitude/linLSB_+0.5)); 
        if (amplitude_int>65535) amplitude_int = 65535;

        linearized_dataframes.push_back(std::make_pair (frame.id(), amplitude_int));
    }
}

double HGCalTriggerCellThresholdCodecImpl::linLSB ( ) const

inline

Definition at line 48 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return linLSB_;}

size_t HGCalTriggerCellThresholdCodecImpl::nCellsInModule ( ) const

inline

Definition at line 46 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return nCellsInModule_;}

int HGCalTriggerCellThresholdCodecImpl::TCThreshold_ADC ( ) const

inline

Definition at line 56 of file HGCalTriggerCellThresholdCodecImpl.h.

{return TCThreshold_ADC_;} 

double HGCalTriggerCellThresholdCodecImpl::TCThreshold_fC ( ) const

inline

Definition at line 57 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return TCThreshold_fC_;} 

uint32_t HGCalTriggerCellThresholdCodecImpl::tdcnBits ( ) const

inline

Definition at line 52 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return tdcnBits_;}

double HGCalTriggerCellThresholdCodecImpl::tdcOnsetfC ( ) const

inline

Definition at line 53 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return tdcOnsetfC_;}

double HGCalTriggerCellThresholdCodecImpl::tdcsaturation ( ) const

inline

Definition at line 51 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return tdcsaturation_;}

void HGCalTriggerCellThresholdCodecImpl::thresholdSelect

(

data_type &

data

)

Definition at line 190 of file HGCalTriggerCellThresholdCodecImpl.cc.

References mps_fire::i, HGCalTriggerCellThresholdDataPayload::payload, and TCThreshold_ADC_.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl(), and triggerCellSums().

{
  for (size_t i = 0; i<data.payload.size();i++){
    if (data.payload[i].hwPt() < TCThreshold_ADC_)  data.payload[i].setHwPt(0);
  }
  
}

uint32_t HGCalTriggerCellThresholdCodecImpl::triggerCellSaturationBits ( ) const

inline

Definition at line 55 of file HGCalTriggerCellThresholdCodecImpl.h.

{return triggerCellSaturationBits_;}

void HGCalTriggerCellThresholdCodecImpl::triggerCellSums	(	const HGCalTriggerGeometryBase &	geometry,
		const std::vector< std::pair< HGCalDetId, uint32_t > > &	linearized_dataframes,
		data_type &	data
	)

Definition at line 160 of file HGCalTriggerCellThresholdCodecImpl.cc.

References HGCalTriggerGeometryBase::getModuleFromTriggerCell(), HGCalTriggerGeometryBase::getOrderedTriggerCellsFromModule(), HGCalTriggerGeometryBase::getTriggerCellFromCell(), jetCorrFactors_cfi::payload, HGCalTriggerCellThresholdDataPayload::payload, thresholdSelect(), and relativeConstraints::value.

Referenced by linearize(), and HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{
    if(linearized_dataframes.size()==0) return;
    std::map<HGCalDetId, uint32_t> payload;
    // sum energies in trigger cells
    for(const auto& frame : linearized_dataframes)
    {
        HGCalDetId cellid(frame.first);
        // find trigger cell associated to cell
        uint32_t tcid = geometry.getTriggerCellFromCell(cellid);
        HGCalDetId triggercellid( tcid );
        payload.insert( std::make_pair(triggercellid, 0) ); // do nothing if key exists already
        uint32_t value = frame.second; 
        payload[triggercellid] += value; // 32 bits integer should be largely enough 

    }
    uint32_t module = geometry.getModuleFromTriggerCell(payload.begin()->first);
    HGCalTriggerGeometryBase::geom_ordered_set trigger_cells_in_module = geometry.getOrderedTriggerCellsFromModule(module);
    // fill data payload
    for(const auto& id_value : payload)
    {
        // Store only energy value and detid
        // No need position here
        data.payload.emplace_back(reco::LeafCandidate::LorentzVector(),
                        id_value.second, 0, 0, 0, id_value.first.rawId());
    }
}

uint32_t HGCalTriggerCellThresholdCodecImpl::triggerCellTruncationBits ( ) const

inline

Definition at line 54 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodec::setDataPayloadImpl().

{return triggerCellTruncationBits_;}

Member Data Documentation

double HGCalTriggerCellThresholdCodecImpl::adcLSB_

private

Definition at line 68 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl(), and linearize().

uint32_t HGCalTriggerCellThresholdCodecImpl::adcnBits_

private

Definition at line 64 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl().

double HGCalTriggerCellThresholdCodecImpl::adcsaturation_

private

Definition at line 63 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl().

size_t HGCalTriggerCellThresholdCodecImpl::dataLength_

private

Definition at line 60 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by decode(), encode(), and HGCalTriggerCellThresholdCodecImpl().

double HGCalTriggerCellThresholdCodecImpl::linLSB_

private

Definition at line 62 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl(), and linearize().

size_t HGCalTriggerCellThresholdCodecImpl::nCellsInModule_

private

Definition at line 61 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by decode(), and encode().

int HGCalTriggerCellThresholdCodecImpl::TCThreshold_ADC_

private

Definition at line 72 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl(), and thresholdSelect().

double HGCalTriggerCellThresholdCodecImpl::TCThreshold_fC_

private

Definition at line 73 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl().

double HGCalTriggerCellThresholdCodecImpl::tdcLSB_

private

Definition at line 69 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl(), and linearize().

uint32_t HGCalTriggerCellThresholdCodecImpl::tdcnBits_

private

Definition at line 66 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl().

double HGCalTriggerCellThresholdCodecImpl::tdcOnsetfC_

private

Definition at line 67 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by linearize().

double HGCalTriggerCellThresholdCodecImpl::tdcsaturation_

private

Definition at line 65 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by HGCalTriggerCellThresholdCodecImpl().

uint32_t HGCalTriggerCellThresholdCodecImpl::triggerCellSaturationBits_

private

Definition at line 71 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by encode(), and HGCalTriggerCellThresholdCodecImpl().

uint32_t HGCalTriggerCellThresholdCodecImpl::triggerCellTruncationBits_

private

Definition at line 70 of file HGCalTriggerCellThresholdCodecImpl.h.

Referenced by encode(), and HGCalTriggerCellThresholdCodecImpl().