FullModuleSumAlgo.cc

Go to the documentation of this file.

#include "L1Trigger/L1THGCal/interface/HGCalTriggerBackendAlgorithmBase.h"
#include "L1Trigger/L1THGCal/interface/fe_codecs/HGCalTriggerCellBestChoiceCodec.h"
#include "L1Trigger/L1THGCal/interface/fe_codecs/HGCalTriggerCellThresholdCodec.h"
#include "DataFormats/ForwardDetId/interface/HGCalDetId.h"

#include "DataFormats/L1THGCal/interface/HGCalCluster.h"

using namespace HGCalTriggerBackend;

template<typename FECODEC, typename DATA>
class FullModuleSumAlgo : public Algorithm<FECODEC> 
{
    public:

        using Algorithm<FECODEC>::name;

    protected:
        using Algorithm<FECODEC>::codec_;

    public:    
        FullModuleSumAlgo(const edm::ParameterSet& conf):
        Algorithm<FECODEC>(conf),
        cluster_product_( new l1t::HGCalClusterBxCollection ){}

        virtual void setProduces(edm::EDProducer& prod) const override final 
        {
            prod.produces<l1t::HGCalClusterBxCollection>(name());
        }

        virtual void run(const l1t::HGCFETriggerDigiCollection& coll, const edm::EventSetup& es) override final
        {
            for( const auto& digi : coll ) 
            {
                DATA data;
                data.reset();
                const HGCalDetId& moduleId = digi.getDetId<HGCalDetId>();
                digi.decode(codec_, data);

                // Sum of trigger cells inside the module
                uint32_t moduleSum = 0;
                for(const auto& triggercell : data.payload)
                {
                    moduleSum += triggercell.hwPt();
                }
                // dummy cluster without position
                // moduleId filled in place of hardware eta
                l1t::HGCalCluster cluster( reco::LeafCandidate::LorentzVector(), moduleSum, 0, 0);
                cluster.setModule(moduleId.wafer());
                cluster.setLayer(moduleId.layer());
                cluster.setSubDet(moduleId.subdetId());
                cluster_product_->push_back(0,cluster);
            }
        }
        virtual void putInEvent(edm::Event& evt) override final 
        {
            evt.put(std::move(cluster_product_),name());
        }

        virtual void reset() override final 
        {
            cluster_product_.reset( new l1t::HGCalClusterBxCollection );
        }

    private:
        std::unique_ptr<l1t::HGCalClusterBxCollection> cluster_product_;

};

        typedef FullModuleSumAlgo<HGCalTriggerCellBestChoiceCodec, HGCalTriggerCellBestChoiceCodec::data_type> FullModuleSumAlgoBestChoice;
        typedef FullModuleSumAlgo<HGCalTriggerCellThresholdCodec, HGCalTriggerCellThresholdCodec::data_type> FullModuleSumAlgoThreshold;
 

DEFINE_EDM_PLUGIN(HGCalTriggerBackendAlgorithmFactory, 
                  FullModuleSumAlgoBestChoice,
                  "FullModuleSumAlgoBestChoice");

DEFINE_EDM_PLUGIN(HGCalTriggerBackendAlgorithmFactory,
                  FullModuleSumAlgoThreshold,
                  "FullModuleSumAlgoThreshold");