HGCalTriggerCellThresholdCodec.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/fe_codecs/HGCalTriggerCellThresholdCodec.h"
#include <limits>

using namespace HGCalTriggerFE;

DEFINE_EDM_PLUGIN(HGCalTriggerFECodecFactory, 
        HGCalTriggerCellThresholdCodec,
        "HGCalTriggerCellThresholdCodec");

HGCalTriggerCellThresholdCodec::
HGCalTriggerCellThresholdCodec(const edm::ParameterSet& conf) : Codec(conf),
    codecImpl_(conf)
{
}


void
HGCalTriggerCellThresholdCodec::
setDataPayloadImpl(const HGCEEDigiCollection& ee,
        const HGCHEDigiCollection& fh,
        const HGCBHDigiCollection& bh) 
{
    data_.reset();
    std::vector<HGCDataFrame<DetId,HGCSample>> dataframes;
    std::vector<std::pair<DetId, uint32_t > > linearized_dataframes;
    // convert ee, fh and bh hit collections into the same object
    if(ee.size()>0)
    {
        for(const auto& eedata : ee)
        {
            dataframes.emplace_back(eedata.id());
            for(int i=0; i<eedata.size(); i++)
            {
                dataframes.back().setSample(i, eedata.sample(i));
            }
        }
    }
    else if(fh.size()>0)
    {
        for(const auto& fhdata : fh)
        {
            dataframes.emplace_back(fhdata.id());
            for(int i=0; i<fhdata.size(); i++)
            {
                dataframes.back().setSample(i, fhdata.sample(i));
            }
        }
    }
    else if(bh.size()>0)
    {
        for(const auto& bhdata : bh)
        {
            dataframes.emplace_back(bhdata.id());
            for(int i=0; i<bhdata.size(); i++)
            {
                dataframes.back().setSample(i, bhdata.sample(i));
            }
        }
    }
    // linearize input energy on 16 bits
    codecImpl_.linearize(dataframes, linearized_dataframes);
    // sum energy in trigger cells
    codecImpl_.triggerCellSums(*geometry_, linearized_dataframes, data_);
    // choose thresholds selected cells in the module
    codecImpl_.thresholdSelect(data_);
}

void
HGCalTriggerCellThresholdCodec::
setDataPayloadImpl(const l1t::HGCFETriggerDigi& digi)
{
    data_.reset();
    // decode input data with different parameters
    // (no selection, so NData=number of trigger cells in module)
    // FIXME:
    // Not very clean to define an alternative codec within this codec 
    // Also, the codec is built each time the method is called, which is not very efficient
    // This may need a restructuration of the FECodec
    edm::ParameterSet conf;
    conf.addParameter<std::string>("CodecName",     name());
    conf.addParameter<uint32_t>   ("CodecIndex",    getCodecType());
    conf.addParameter<uint32_t>   ("MaxCellsInModule", codecImpl_.nCellsInModule());
    conf.addParameter<uint32_t>   ("NData",         codecImpl_.nCellsInModule());
    // The data length should be the same for input and output, which is limiting
    conf.addParameter<uint32_t>   ("DataLength",    codecImpl_.dataLength());
    conf.addParameter<double>     ("linLSB",        codecImpl_.linLSB());
    conf.addParameter<uint32_t>   ("linnBits",      codecImpl_.linnBits());
    conf.addParameter<double>     ("adcsaturation", codecImpl_.adcsaturation());
    conf.addParameter<uint32_t>   ("adcnBits",      codecImpl_.adcnBits());
    conf.addParameter<double>     ("tdcsaturation", codecImpl_.tdcsaturation());
    conf.addParameter<uint32_t>   ("tdcnBits",      codecImpl_.tdcnBits());
    conf.addParameter<double>     ("tdcOnsetfC",    codecImpl_.tdcOnsetfC());
    conf.addParameter<double>     ("adcsaturationBH", codecImpl_.adcsaturationBH());
    conf.addParameter<uint32_t>   ("adcnBitsBH",      codecImpl_.adcnBitsBH());
    conf.addParameter<uint32_t>   ("triggerCellTruncationBits", codecImpl_.triggerCellTruncationBits());
    conf.addParameter<double>        ("TCThreshold_fC", codecImpl_.TCThreshold_fC());
    conf.addParameter<double>        ("TCThresholdBH_MIP", codecImpl_.TCThresholdBH_MIP());
    conf.addParameter<std::vector<double>>("ThicknessCorrections", codecImpl_.ThicknessCorrections());
    HGCalTriggerCellThresholdCodec codecInput(conf);
    codecInput.setGeometry(geometry_);
    digi.decode(codecInput,data_);
    // choose threshold selected cells in the module
    codecImpl_.thresholdSelect(data_);
}


std::vector<bool>
HGCalTriggerCellThresholdCodec::
encodeImpl(const HGCalTriggerCellThresholdCodec::data_type& data) const 
{
    return codecImpl_.encode(data, *geometry_);
}

HGCalTriggerCellThresholdCodec::data_type
HGCalTriggerCellThresholdCodec::
decodeImpl(const std::vector<bool>& data, const uint32_t module) const 
{
    return codecImpl_.decode(data, module, *geometry_);
}